MATLAB ask psk fsk dpsk调制解调

%产生比特信号
t=0:0.01:7.99;
a=randint（18）;
figure（1）
m=a（ceil（t+0.01））;
plot（tm）
title（‘产生随机八比特二进制比特序列‘）;
axis（[08-1.51.5]）;
%*************调制************%
fc=4000;          %载波频率
fc1=8000;
fs=80000;         %采样频率
ts=0:1/fs:（800-1）/fs;
ts1=0:1/fs:（100-1）/fs;
tzxh1=cos（2*pi*fc*ts）;
tzxh2=cos（2*pi*fc1*ts）;%2FSK调制情况下另一载波
%**************2ASK调制************%
ask=m.*tzxh1;
figure（2）
plot（task）
title（‘2ASK调制波形‘）;
axis（[08-1.51.5]）;
%**************2FSK调制************%
fsk=cos（2*pi*（fc+fc*m）.*ts）;
figure（3）
plot（tfsk）
title（‘2FSK调制波形‘）;
axis（[08-1.51.5]）;
%**************2PSK调制************%
psk=cos（2*pi*fc.*ts+pi*m）;
figure（4）
plot（tpsk）
title（‘2PSK调制波形‘）;
axis（[08-1.51.5]）;
%**************2DPSK调制************%
k=1;
for i=0:7
    if m（1i*100+50）>0        
   dpsk（i*100+1:i*100+100）=cos（2*pi*fc.*ts1-pi/2*k）;
   k=-k;     
    else      
  dpsk（i*100+1:i*100+100）=cos（2*pi*fc.*ts1+pi/2*k）; 
    end
end 
figure（5）
plot（tdpsk）
title（‘2DPSK调制波形‘）;
axis（[08-1.51.5]）;
%**************叠加噪声************%
e_2ask=awgn（ask10）;
e_2fsk=awgn（fsk10）;
e_2psk=awgn（psk10）;
e_2dpsk=awgn（dpsk10）;
figure（6）
subplot（221）
plot（te_2ask）
title（‘2ASK调制信号叠加噪声波形‘）;
axis（[08-1.51.5]）;
subplot（222）
plot（te_2fsk）
title（‘2FSK调制信号叠加噪声波形‘）;
axis（[08-1.51.5]）;
subplot（223）
plot（te_2psk）
title（‘2PSK调制信号叠加噪声波形‘）;
axis（[08-1.51.5]）;
subplot（224）
plot（te_2dpsk）
title（‘2DPSK调制信号叠加噪声波形‘）;
axis（[08-1.51.5]）;

%**************2ASK相干解调************%
[b11a11]=ellip（50.560[20006000]*2/80000）;%带通椭圆滤波器设计
e_ask1=filter（b11a11e_2ask）;%相乘器
e_ask2=2*e_ask1.*tzxh1*2;%低通滤波器
[b12a12]=ellip（50.5601000*2/80000）;       %低通滤波器设计
ask_xgjt=filter（b12a12e_ask2）;       %相干解调后，抽样判决前的结果
figure（7）
subplot（211）
plot（task_xgjt）
title（‘2ASK调制信号相干解调后通过低通滤波器，抽样判决前的信号‘）;
axis（[08-1.51.5]）;
%**************2ASK非相干解调************%  
e_ask1_f=abs（e_ask1）;%包络检波
ask_fxgjt=filter（b12a12e_ask1_f）;       %抽样判决前的结果
subplot（212）
plot（task_xgjt）
title（‘2ASK调制信号非相干解调后通过低通滤波器，抽样判决前的信号‘）;
axis（[08-1.51.5]）;
%**************2FSK相干解调************%
[b21a21]=ellip（50.560[20006000]*2/80000）;%带通椭圆滤波器1设计
e_fsk11=filter（b21a21e_2fsk）; 
[b22a22]=ellip（50.560[600010000]*2/80000）;%带通椭圆滤波器2设计
e_fsk12=filter（b22a22e_2fsk）;   
fsk1_m1=e_fsk11.*tzxh1*2;    %相干解调
fsk1_m2=e_fsk12.*tzxh2*2;
fsk1_xgjt=filter（b12a12fsk1_m1）;           %相干解调后，抽样判决前的结果  
fsk2_xgjt=filter（b12a12fsk1_m2）;
figure（8）
subplot（211）;
plot（tfsk1_xgjt）;
title（‘2FSK调制信号相干解调后通过低通滤波器，抽样判决前的信号（发0时）‘）;
subplot（212）;
plot（tfsk2_xgjt）;
title（‘2FSK调制信号相干解调后通过低通滤波器，抽样判决前的信号（发1时）‘）;
%**************2FSK非相干解调************%
e_fsk1_f=abs（fsk1_m1）; %包络检波
e_fsk2_f=abs（fsk1_m2）;
fsk1_fxgjt=filter（b12a12e_fsk1_f）;           %相干解调后，抽样判决前的结果  
fsk2_fxgjt=filter（b12a12e_fsk2_f）;
figure（9）
subplot（211）;
plot（tfsk1_fxgjt）;
title（‘2FSK调制信号包络检波后通过低通滤波器，抽样判决前的信号（发0时）‘）;
subplot（212）;
plot（tfsk2_fxgjt）;
title（‘