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);