echo off
clear all
close all
snr_in_DB=-4:1:6;
SNR_in_DB=-4:0.1:6;

for i=1:length(snr_in_DB)
   [py(i),pb(i)]=AwgnPeG(snr_in_DB(i));%py為沒有經(jīng)過信道編碼的誤碼率，pb為經(jīng)過信道編碼的誤碼率
end;
for i=1:length(SNR_in_DB),
   SNR=10^(SNR_in_DB(i)/10);
   p(i)=Qfunct(sqrt(2*SNR));%計(jì)算理想誤碼率p
end; 
%畫出誤比特率圖
semilogy(snr_in_DB,py,'b*');
hold on
semilogy(snr_in_DB,pb,'r+');
semilogy(SNR_in_DB,p);
py
pb