fs?=?100;
dt?=?1/fs;
f1?=?10;
f2?=?30;??%A兩個信號的頻率
a?=?round?（rand（110））%隨機信號
g1?=?a;
g2?=?~a;?%信號翻轉，和g1反向
g11?=?（ones（1fs））‘*?g1?;%抽樣
g11（:）;
g1a?=?g11（:）‘?;
g21?=?（ones（1fs））‘*g2;
g2a?=?g21（:）‘;
t?=?0?:?dt:10-dt;
t1?=?length（t）;
fsk1?=?g1a.*cos（2*pi?*f1.*t）;
fsk2?=?g2a.*cos（2*pi?*f2.*t）;
fsk?=?fsk1+fsk2;????%??產生信號
n0??=?0.01*randn（1t1）;
sn?=?fsk+n0;
subplot（311）;
plot（tn0）;%噪聲波形
title?（‘噪聲波形‘）
ylabel（‘幅度‘）

subplot（312）;
plot（tfsk）;
title?（‘產生波形‘）
ylabel（‘幅度‘）
subplot（313）;
plot（tsn）;
title?（‘將要通過濾波器的波形‘）
xlabel（‘t‘）
ylabel（‘幅度‘）

figure（2）?%?fsk解調------------------------------------圖二
b1?=?fir1（10?[8/50?12/50]）;
b2?=?fir1（10?[28/50?32/50]）;%設置帶通參數
h1?=?filter（b1?1?sn）;
h2?=?filter（b2?1?sn）;??%?經過帶通濾波器后的信號

subplot（211）;
plot（th1）;
title?（‘經過帶通濾波器f1后的波形‘）

ylabel（‘幅度‘）
subplot（212）;
plot（th2）;
ti