mtkl.m

来自「通原实验课的作业」· M 代码 · 共 21 行

close all
clear all
SNRindB1=0:1:8;
SNRindB2=0:0.1:8;
for i=1:length(SNRindB1),
	%误码率仿真
	smld_err_prb(i)=smldPe(SNRindB1(i));
end
for i=1:length(SNRindB2),
	SNR=exp(SNRindB2(i)*log(10)/10);
	%理论误码率
	theo_err_prb(i)=Qfunct(sqrt(2*SNR));
end
%随后为绘图命令
semilogy(SNRindB1,smld_err_prb,'*');
hold
semilogy(SNRindB2,theo_err_prb);	
grid
xlabel('Eb/No in dB')
ylabel('error probability')
title('2PSK通信系统的蒙特卡洛仿真(理想信道)')