FSK的调制和解调matlab实现

clear all;
clc;

%2FSK调制
Fs=20;   %插入的个数
f1=2;    %对应1
f2=6;    %对应0
ts=1/Fs;
fd=1;
code=randint（120）;
[xm]=pulse（codeFsfd）;%插入，得到矩形波
t=ts:ts:length（code）;
x1=（x-1）;
x1=-x1;
y=x.*cos（2*pi*f1*t）+x1.*cos（2*pi*f2*t）;
X=fft（x）;
Y=fft（y）;

n=length（x）;
f=[0:Fs/n:Fs-Fs/n]-Fs/2;%表示f是从[-106]，间隔为4

figure（1）
subplot（221）;
plot（x）;
axis（[0 length（x） -1 2]）;
title（‘基带信号‘）;
grid;
subplot（222）;
plot（fabs（fftshift（X）））;
title（‘基带信号频谱‘）;
grid;

subplot（223）;
plot（tytx‘-r‘）;
title（‘2FSK信号‘）;
grid;
subplot（224）;
plot（fabs（fftshift（Y）））;
title（‘2FSK信号频谱‘）;
grid;

%加入噪声
figure（2）
snr=20;
y=awgn（ysnr）;
plot（y）;
title（‘加入噪声后的信号‘）;
xlabel（‘t‘）;
ylabel（‘Amplitude‘）;

%相干解调
fp1=[13];
fs1=[0.53.5];
fp2=[57];
fs2=[4.57.5];
wp1=fp1/（Fs/2）;
ws1=f