%產生比特信號
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（‘