clear
n=8;%隨機(jī)碼元個數(shù)
N=100;%模擬一個碼元的點數(shù)
K=4;%figure1畫四個小圖
a=randint（1n）;%碼元生成
bita=[];%定義空數(shù)組，存放基帶信號
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=-c1;%載波
????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已調(diào)信號‘）;
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）;%設(shè)計IIR帶通濾波器，階數(shù)