examinawgn.m

来自「依据GSM 05.03建议」· M 代码 · 共 17 行

echo off
clear all
close all
SNRindB1=-4:1:6;%定义信噪比的测试点
SNRindB2=-4:0.1:6;
for i=1:length(SNRindB1),
    cnv_err_rate(i)=cnv_awgn_pe(SNRindB1(i));%改善后的误码率，调用函数
    err_rate(i)=smldpe(SNRindB1(i));%原始误码率，调用函数
end;
for i=1:length(SNRindB2),
    snr=exp(SNRindB2(i)*log(10)/10);
    theo_err_rate(i)=Qfunct(sqrt(4*snr));%理论上的原始误码率，调用函数
end;
semilogy(SNRindB1,err_rate,'*');
hold
semilogy(SNRindB1,cnv_err_rate,'+');
semilogy(SNRindB2,theo_err_rate);