%蒙特卡羅仿真（AWGN信道，硬判決，2PSK）
%
Eb=2;     %每比特的能量,能量越大，卷積譯碼的誤碼率越小
k=1;
G=[1 0 1;1 1 1];
input=[1 0 1 0 1 1 1 0];
dsource=cnv_encd(G,k,input);
N=length(dsource);
n=size(G,1);    %得到碼字長度
Ec=k*Eb/n;      %得到碼元能量

% 通過AWGN信道
for i=1:N
%matched filter output
%soft-decision
   %產生的高斯噪聲均值為0，方差為1
   if(dsource(i)==0)
       r(i)=-Ec+gngauss(0,1);
   else 
       r(i)=Ec+gngauss(0,1);
   end;
   %detector follows
%hard decision
   if(r(i)>0)
       decis(i)=1;
   else
       decis(i)=0;
   end;
end;
%硬判決維特比譯碼
[decoder_output,survivor_state,cumulated_metric]=viterbi(G,k,decis);