qpsk_simulation.m

matlab通信仿真

clear;
clf;
echo off
SNRindB1=0:2:10;
SNRindB2=0:0.1:10;
for i=1:length(SNRindB1)
    [pb,ps]=cm_sm32(SNRindB1(i));    % simulated bit and symbol error rates
    bit_err_prb(i)=pb;
    symbol_err_prb(i)=ps;
    %echo off;
end

echo off;
for i=1:length(SNRindB2)
    SNR=exp(SNRindB2(i)*log(10)/10);     % Signal to Noise ratio
    theo_err_prb(i)=Qfunct(sqrt(2*SNR));    %theoretical bit error rate
   % echo off;
end;

echo off;
semilogy(SNRindB1,bit_err_prb,'*','linewidth',2);
hold on;
semilogy(SNRindB1,symbol_err_prb,'o','linewidth',2);
hold on;
semilogy(SNRindB2,theo_err_prb,'linewidth',2);
xlabel('Eb/n0 (dB)');

legend('仿真的误比特率曲线','仿真的误信元率曲线','误比特率曲线的理论值');