2FSK调制解调MATLAB 程序

Fc=10;   %载频      
Fs=40;   %系统采样频率      
Fd=1;    %码速率     
N=Fs/Fd;        
df=10;      
numSymb=25;%进行仿真的信息代码个数     
M=2;       %进制数      
SNRpBit=60;%信噪比     
SNR=SNRpBit/log2（M）;%60   
seed=[12345 54321];     
numPlot=15;    
x=randsrc（numSymb1[0:M-1]）;%产生25个二进制随机码     
figure（1）     
stem（[0:numPlot-1]x（1:numPlot）‘bx‘）;%显示15个码元，杆图，从x的前十五个随机数中选取     
title（‘二进制随机序列 杆图‘）      
xlabel（‘T‘）;     
ylabel（‘幅度‘）; 

%调制
 y=dmod（xFcFdFs‘fsk‘Mdf）;%数字带通调制 
 %x是一个随机序列，如果M=2则x=[1 1 0 0]之类的，Fd为码元速率，Fc为载频，Fs为采样频率，fsk表示为频移键控调制。     
 numModPlot=numPlot*Fs; %15*40     
 t=[0:numModPlot-1]./Fs;%数组除法（仿真时间）     
 figure（2）      
 plot（ty（1:length（t））‘b-‘）;     
 axis（[min（t） max（t） -1.5 1.5]）;      
 title（‘调制后的信号‘）    
 xlabel（‘Time‘）;     
 yl