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