2PSK的调制与解调

clear clc;
fs=9e5;               %采样频率
fm=21e3;              %基带频率
n=2*（6*fs/fm）;
final=（1/fs）*（n-1）;
fc=2e5;               % 载波频率
t=0:1/fs:（final）;
Fn=fs/2;              %奈奎斯特频率
%产生方波
wct=2*pi*fm*t; 
A=1;
p=0;
x = A * cos（wct + p）;
am=1;
x（x>0）=am;
x（x<0）=-1;
figure（1）
subplot（321）;
plot（tx）;
axis（[0 2e-4 -2 2]）;
title（‘基带信号‘）;
grid on
Fc=sin（2*pi*fc*t）;   %载波
ask=x.*Fc;           %载波调制
subplot（322）;
plot（task）;
axis（[0 200e-6 -2 2]）;
title（‘2PSK信号‘）;
grid on;
%解调
%1 带通滤波
fBW=40e3;
f=[0:3e3:4e5];
w=2*pi*f/fs;
z=exp（w*j）;
BW=2*pi*fBW/fs;
a=.8547;              
p=（j^2*a^2）;
gain=.135;
Hz=gain*（z+1）.*（z-1）./（z.^2-（p））;
Hz（Hz==0）=10^（8）;   
%滤波器参数
a=[1 0 0.7305];        
b=[0.135 0 -0.135];    
faskn=fi