dpsk8.m

来自「根据蒙特卡罗仿真原理」· M 代码 · 共 30 行

%********************************************************
% dpsk8.m
% 8dPSK通信系统蒙特卡罗仿真              
% 在不同信噪比Eb/N0下发送10000符号的蒙特卡罗仿真，
% 其中Eb=Es/2是比特能量，图中给出了误码率和相应的理论误码率
% Programmed by Liu Dejian
% Dtat:2008-3-8
%********************************************************
clear
clc
echo on
SNRindB1=0:2:12;
SNRindB2=0:0.1:12;
for i=1:length(SNRindB1),
     smld_err_prb(i)=cm_sm8(SNRindB1(i));  % 仿真符号误码率
end;
for i=1:length(SNRindB2),
     SNR=exp(SNRindB2(i)*log(10)/10);       % 信噪比
     theo_err_prb(i)=2*Qfunct(sqrt(SNR));   % 理论符号误码率
    theo_err_prb2(i)=2*Qfunct(sqrt(6*SNR)*sin(pi/(8*sqrt(2))));
    % 此理论符号误码率是根据曹志刚编著《现代通信原理》书P292（10-169）
    % 式得到理论符号误码率计算公式
end;
% 随后为绘图曲线
semilogy(SNRindB1,smld_err_prb,'*');
hold
semilogy(SNRindB2,theo_err_prb);
 semilogy(SNRindB2,theo_err_prb2,'r');
title('8dPSK通信系统的蒙特卡洛仿真')
legend('仿真符号误码率','理论符号误码率','理论符号误码率(曹)')