mqam_simulation.m

matlab通信仿真

% Performance of 16QAM system from the Monte Carlo Simulation
clear;
% clf;
echo off
SNRindB1=0:2:15;
SNRindB2=0:0.1:15;
M=16;
k=log2(M);
for i=1:length(SNRindB1)
    simu_err_prb(i)=cm_sm41(SNRindB1(i));     % simulated error rate
end

for i=1:length(SNRindB2)
    SNR=exp(SNRindB2(i)*log(10)/10);        % signal to noise ratio
    
    % theoretical symbol error rate
    theo_err_prb(i)=3*Qfunct(sqrt(3*k*SNR/(M-1)));
end

% plotting commands
semilogy(SNRindB1,simu_err_prb,'*');
hold on;
semilogy(SNRindB2,theo_err_prb,'linewidth',2);
xlabel('Eb/n0 (dB)');

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