mfsk调制解调matlab源程序

      Fc=10;   %载频
      Fs=40;   %系统采样频率
      Fd=1;    %码速率
      N=Fs/Fd;
      df=10;
      numSymb=25;%进行仿真的信息代码个数
      M=8;       %进制数
      SNRpBit=60;%信噪比
      SNR=SNRpBit/log2(M);
      seed=[12345 54321];
      numPlot=15;
      %产生25个二进制随机码
      x=randsrc(numSymb,1,[0:M-1]);%产生25个二进制随机码
      figure(1)
      stem([0:numPlot-1],x(1:numPlot),'bx');
      title('二进制随机序列')
      xlabel('Time');
      ylabel('Amplitude');
      %调制
      y=dmod(x,Fc,Fd,Fs,'fsk',M,df);
      numModPlot=numPlot*Fs;
      t=[0:numModPlot-1]./Fs;
      figure(2)
      plot(t,y(1:length(t)),'b-');
      axis([min(t) max(t) -1.5 1.5]);
      title('调制后的信号')
      xlabel('Time');