dpsk8.m

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

clear all
echo off
SNRindB1=0:2:20;
SNRindB2=0:0.1:20;
for i=1:length(SNRindB1)
     smld_err_prb(i)=8dpsk_2(SNRindB1(i));  %仿真误码率
end;
for i=1:length(SNRindB2),
     SNR=exp(SNRindB2(i)*log(10)/10);       %信噪比
     theo_err_prb(i)=exp((-2)*SNR*(sin(pi/16).^2));   %理论符号误码率
end;
%随后为绘图曲线
semilogy(SNRindB1,smld_err_prb,'*');
hold
semilogy(SNRindB2,theo_err_prb);
grid
xlabel('Eb/No in dB')
ylabel('error probability')
title('8DPSK通信系统的蒙特卡洛仿真（多径信道）')