clear
n=8;%隨機碼元個數
N=100;%模擬一個碼元的點數
K=4;%figure1畫四個小圖
a=randint（1n）;%碼元生成
bita=[];%定義空數組，存放基帶信號
for?i=1:length（a）
????if?a（i）==0
????????bit1=zeros（1N）;
????else?
????????bit1=ones（1N）;
????end
?bita=[bitabit1];%基帶信號
end
figure（1）;
subplot（411）;
plot（bita‘LineWidth‘1.5）;
axis（[0N*length（a）-22]）;
title（‘基帶信號‘）;
grid?on;
bitRate=1e3;%每一個碼元中采樣點的間隔寬度0.001s
fc=1e3;%載頻1KHZ
t=linspace（01/bitRateN）;
tz=[];
????c1=sin（2*pi*t*fc）;%載波
????c2=sin（2*pi*t*fc*2）;%載波
????for?i=1:length（a）
???????if?a（i）==1
??????????tz=[tzc1];
???????else
?????????tz=[tzc2];
???????end
????end
subplot（412）;
plot（tz‘LineWidth‘1.5）;
title（‘2FSK已調信號‘）;
grid?on;
signal=awgn（tz.2‘measured‘）;%加噪
subplot（413）;
plot（signal‘LineWidth‘1.5）;
grid?on;
title（‘信號+噪聲‘）;
Fs=5e3;%采樣頻率
[b1a1]=ellip（40.140[999.91000.1]*2/Fs）;%設計IIR帶通濾波器，階數為4，通帶紋波0.1，阻帶衰減40DB
[b2a2]=ellip（40.140[1999.92000.1]*2/Fs）;
sa=filter（b1a1signal）;%信號通過該濾波器
sb=filter（b2a2signal）;
figure（2）;
sia=[];
sib=[];
????for?i=1:n
????????sia=[sia