function FSKFc=10;   %載頻Fs=40;   %系統采樣頻率Fd=1;    %碼速率N=Fs/Fd;df=10;numSymb=25;%進行仿真的信息代碼個數M=2;       %進制數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');ylabel('Amplitude');%在已調信號中加入高斯白噪聲randn('state',seed(2));y=awgn(y,SNR-10*log10(0.5)-10*log10(N),'measured',[],'dB');%在已調信號中加入高斯白噪聲figure(3)plot(t,y(1:length(t)),'b-');%畫出經過信道的實際信號axis([min(t) max(t) -1.5 1.5]);title('加入高斯白噪聲后的已調信號')xlabel('Time');ylabel('Amplitude');%相干解調figure(4)z1=ddemod(y,Fc,Fd,Fs,'fsk/eye',M,df);title('相干解調后的信號的眼圖')%帶輸出波形的相干M元頻移鍵控解調figure(5)stem([0:numPlot-1],x(1:numPlot),'bx');hold on;stem([0:numPlot-1],z1(1:numPlot),'ro');hold off;axis([0 numPlot -0.5 1.5]);title('相干解調后的信號原序列比較')legend('原輸入二進制隨機序列','相干解調后的信號')xlabel('Time');ylabel('Amplitude');%非相干解調figure(6)z2=ddemod(y,Fc,Fd,Fs,'fsk/eye/noncoh',M,df);title('非相干解調后的信號的眼圖')%帶輸出波形的非相干M元頻移鍵控解調figure(7)stem([0:numPlot-1],x(1:numPlot),'bx');hold on;stem([0:numPlot-1],z2(1:numPlot),'ro');hold off;axis([0 numPlot -0.5 1.5]);title('非相干解調后的信號')legend('原輸入二進制隨機序列','非相干解調后的信號')xlabel('Time');ylabel('Amplitude');%誤碼率統計[errorSym ratioSym]=symerr(x,z1);figure(8)simbasebandex([0:1:5]);title('相干解調后誤碼率統計')[errorSym ratioSym]=symerr(x,z2);      figure(9)simbasebandex([0:1:5]);title('非相干解調后誤碼率統計')%濾除高斯白噪聲Delay=3;R=0.5;PropD=0;   %滯后3s[yf,tf]=rcosine(Fd,Fs,'fir',R,Delay);     %升余弦函數[yo2,to2]=rcosflt(y,Fd,Fs,'filter',yf);%加入高斯白噪聲后的已調信號和經過升余弦濾波器后的已調信號t=[0:numModPlot-1]./Fs;figure(10)plot(t,y(1:length(t)),'r-');hold on;plot(to2,yo2,'b-');hold off;axis([0 30 -1.5 1.5]);xlabel('Time');ylabel('Amplitude');legend('加入高斯白噪聲后的已調信號','經過升余弦濾波器后的已調信號')title('升余弦濾波前后波形比較')eyediagram(yo2,N);%眼圖title('加入高斯白噪聲后的已調信號的眼圖')