echo off
clear all
close all
SNRindB1=-4:0.5:6;%定義信噪比的測試點
SNRindB2=-4:0.2:6;
for i=1:length(SNRindB1),
    [class1_err_rate(i),sig_err_rate(i)]=gbwhole(SNRindB1(i));%改善后的誤碼率，調(diào)用函數(shù),卷積加交織，高斯加突發(fā)干擾
    sig_err_rate0(i)=gbconv(SNRindB1(i));%改善后的誤碼率，調(diào)用函數(shù)，只有卷積，高斯加突發(fā)干擾
end;
for i=1:length(SNRindB2),
    snr=exp(SNRindB2(i)*log(10)/10);
    theo_err_rate(i)=Qfunct(sqrt(4*snr));%理論上的原始誤碼率，調(diào)用函數(shù)
end;
semilogy(SNRindB1,sig_err_rate,':r*',SNRindB1,class1_err_rate,':*');
hold
semilogy(SNRindB1,sig_err_rate0,':g+');
semilogy(SNRindB2,theo_err_rate);