2fsk调制解调MATLAB仿真代码

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