fsk.zip_2FSK_2FSK 调制_2FSK 调制解调_2FSK调制解调

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(numSymb,1,[0:M-1]);%产生25个二进制随机码 figure(1) stem([0:numPlot-1],x(1:numPlot),'bx');%显示15个码元，杆图，从x的前十五个随机数中选取 title('二进制随机序列 杆图') xlabel('T'); ylabel('幅度'); %调制 y=dmod(x,Fc,Fd,Fs,'fsk',M,df);%数字带通调制 %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(t,y(1:length(t)),'b-'); axis([min(t) max(t) -1.5 1.5]); title('调制后的信号') xlabel('Time'); ylabel('Amplitude'); %在已调信号中加入高斯白噪声 randn('state',seed(2)); %生成-2到+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('相干解调后的信号的眼图') figure(5) stem([0:numPlot-1],x(1:numPlot),'rx'); hold on; stem([0:numPlot-1],z1(1:numPlot),'go'); hold off; axis([0 numPlot -0.5 1.5]); title('相干解调后的信号原序列比较') legend('原输入二进制随机序列','相干解调后的信号') xlabel('Time'); ylabel('Amplitude');