dpsk2.m

通原实验课的作业

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