test91.cpp

截短长度72viterbi约束长度9测试

/*
*                          viterbi译码9约束长度截短测试
************************************************************************************************
* MODULE NAME:  test91.cpp  
*************************************************************************************************
* DESCRIPTION     :BER测试程序 截短长度72
*
* FUNCTIONS       :信息位编码，加白噪声，译码，测BER 
*                                                   
*************************************************************************************************
*/

#include "decode.h"
#include "pn_seq.h"
#include "channel.h"
void main()
{   int i,j;
    int  error_count=0;    //译码后错误信息比特
    int  trans_error=0;    //译码前传输错误码元
	int  error_frame=0;    //错误帧
    int num;               //记录结束个数
    bool error=false;     
    unsigned a[192];  //信源信息帧
    double channel_trans1[192]; //信道经噪声污染的已编码帧
    double channel_trans2[192];
    int b1_polar[192];   //已编码帧加极性
	int b2_polar[192];

	l=184;	m=8; //信息位长度184 归零位长度8

AWGN p,*pi;
cout<<"信噪比:";
cin>>p.snr;
cout<<"跳出次数:";
cin>>num;

p.rate=0.5;   //设置白噪声
p.reset=1;
p.ix=1;
p.iy=1;
p.iz=1;
p.I_sigma=2*(p.rate)*pow(10.0,(p.snr)/10);
p.sigma=sqrt(1/p.I_sigma);
pi=&p;
	PN_seq(a,l,true); //伪随机序列生成信息位
   
	for (i=l;i<l+m;i++) //加归零位
	 a[i]=0;
    
   conv_encode_m9(a,l+m); //编码
    
	j=0;
	error_count=0;
    trans_error=0;
    error_frame=0;
 
 init_trans();//建立状态网格图


  for ( i=0;i<l+m;i++)   //编码加极性
	{ 
		if (b1[i]==1) b1_polar[i]=1;
		              else b1_polar[i]=-1;
        if (b2[i]==1) b2_polar[i]=1;
		              else b2_polar[i]=-1;
}

//while	((error_count<num)||(j<1000))
 while (error_count<num)
 {  j++;
  

  

  //加信道噪声
awgn_channel(b1_polar , channel_trans1 ,pi,l+m);
awgn_channel(b2_polar , channel_trans2 ,pi,l+m);

  //接收端量化判决
 for ( i=0;i<m+l;i++)
 {    
	   if (channel_trans1[i]>0)  receiv1[i]=1;
	                    else receiv1[i]=0;
       if (channel_trans2[i]>0) receiv2[i]=1;
	                    else receiv2[i]=0;
   }

  
 search();    //译码
 
 for (i=0;i<l+m;i++) //统计接收错误
 {    
	 if (b1[i]!=receiv1[i]) {trans_error++;}
	 if (b2[i]!=receiv2[i]) {trans_error++;}
  }

   error=false;
   for (i=0;i<l+m;i++)  //统计译码错误
   if (a[i]!=result[i])
   {      error=true;
	      error_count++;
	      cout<<error_count; 
	  };    
  if (error)error_frame++; //统计误帧

} 
cout<<"总帧数:"<<j<<endl;	
cout<<"错误比特:"<<error_count<<"BER:"<<(double)(error_count)/(192*j)<<endl;
cout<<"传输错误码元:"<<trans_error<<"误码元率:"<<(double)(trans_error)/(192*j*2)<<endl; 
cout<<"传输错误帧:"<<error_frame<<"误帧率:"<<(double)error_frame/j<<endl; 


}