pskmotocalo.m

qpsk通信系统的MonteCalo仿真程序

echo on
SNRindB1=0:2:10;   %定义信噪比序列，共有六个值
SNRindB2=0:0.1:10; %扫描用的信噪比序列，采用dB作单位
for i=1:length(SNRindB1)              %计算误比特率
    [pb,ps]=pskmoto(SNRindB1(i));
    smld_bit_err_prb(i)=pb;
    smld_symbol_err_prb(i)=ps;
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_bit_err_prb,'*');%以对数形式作y坐标绘图
hold                                    %将上一曲线保持住
semilogy(SNRindB1,smld_symbol_err_prb,'o');%作出实际的信噪比——误比特率点
semilogy(SNRindB2,theo_err_prb);           %作出理论的信噪比——误比特率点