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

%產生加性高斯噪聲
for i=1:N
awgn(i)=gngauss(0,0.5);
end;

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