awgn_hard_2fsk.m

该MATLAB程序评估了viterbi算法在衰落信道下的卷积译码的性能

%蒙特卡罗仿真（AWGN信道，硬判决，2FSK）
%
Eb=1;     %每比特的能量,能量越大，卷积译码的误码率越小
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，方差为0.5     
   %高斯噪声方差越小，信噪比越大，则卷积译码的误码率越小
   if(dsource(i)==0)
       r0(i)=sqrt(Ec)+gngauss(0,0.5);
       r1(i)=gngauss(0,0.5);
   else 
       r0(i)=gngauss(0,0.5);
       r1(i)=sqrt(Ec)+gngauss(0,0.5);
   end;
   %detector follows
%hard decision
   if(r0(i)>r1(i))
       decis(i)=0;
   else
       decis(i)=1;
   end;
end;
%硬判决维特比译码
[decoder_output,survivor_state,cumulated_metric]=viterbi(G,k,decis);