程序.rar_MFSK 误码率_MFSK解调_MFSK调制解调_mfsk 调制解调_mfsk误码率

% MFSK信号仿真 M=2,4,8可以修改 M=2; %调制阶数 freq_sep=16; %频率间隔（Hz） nsamp=32; %每个符号的采样点数目 SymbolRate=8; %符号速率 Fs=nsamp*SymbolRate; %采样速率 注意：M*freq_sep<=Fs multi=4; %载波频率为符号速率的倍数 fc=multi*SymbolRate; %设置载波频率 N_symbol=200; %每帧包含的符号数目 frame=1000; %帧数 EbNo=0:2:16; %设置信噪比 bit_err=zeros(1,length(EbNo)); FSKMod=comm.FSKModulator('ModulationOrder',M,'BitInput',1,'SymbolMapping','Binary','FrequencySeparation',freq_sep,'SamplesPerSymbol',nsamp,'SymbolRate',SymbolRate); FSKDemod=comm.FSKDemodulator('ModulationOrder',M,'BitOutput',1,'SymbolMapping','Binary','FrequencySeparation',freq_sep,'SamplesPerSymbol',nsamp,'SymbolRate',SymbolRate); data=randi([0 1],log2(M)*N_symbol,1); %产生发送比特数据 MFSK_Mod=step(FSKMod,data); %产生MFSK基带数据符号 %加入载波，得到MFSK基带信号 I_mod=real(MFSK_Mod); %分成I路和Q路两部分 Q_mod=imag(MFSK_Mod); %乘以载波 I_temp=interp(I_mod,multi); %插值 Q_temp=interp(Q_mod,multi); fs2=fc*nsamp; %采样速率 t=1/fs2:1/fs2:length(I_temp)*1/fs2; %采样点 MFSK_I=I_temp.'.*cos(2*pi*fc*t); MFSK_Q=Q_temp.'.*sin(2*pi*fc*t); %得到最终的MFSK信号 MFSK_signal=MFSK_I-MFSK_Q; %画出MFSK信号的波形 figure(1) plot(MFSK_signal); hold on; title('MFSK调制信号波形'); axis([0 length(MFSK_signal) min(MFSK_signal)-0.5 max(MFSK_signal)+0.5]) %画出MFSK信号的频谱图 h=dsp.SpectrumAnalyzer('SampleRate',Fs); step(h,MFSK_Mod); %MFSK信号的解调 data_demod=step(FSKDemod,MFSK_Mod); %画出MFSK调制信号 figure(3) stairs(1:length(data),data.'); title('MFSK基带信号'); axis([0 length(data) min(data)-0.5 max(data)+0.5]) %画出MFSK解调信号 figure(4) stairs(1:length(data_demod),data_demod.'); title('MFSK解调信号'); axis([0 length(data_demod) min(data_demod)-0.5 max(data_demod)+0.5]) %加噪声信号 EbNo0=10; en=EbNo0-10*log10(nsamp)+10*log10(log2(M)); hChan=comm.AWGNChannel('NoiseMethod','Signal to noise ratio (Eb/No)','SignalPower',1,'SamplesPerSymbol',1); hChan.EbNo=en; MFSK_signal_awgn=step(hChan,MFSK_signal); %画出加噪声之后的信号波形 figure(5) plot(MFSK_signal_awgn); title('MFSK信号加入噪声信号波形');