viterbi 卷积译码代码MATLAB语言

/* * 硬判决viterbi译码9约束长度 ************************************************************************************************ * MODULE NAME: test9.cpp ************************************************************************************************* * DESCRIPTION :BER测试程序 * * FUNCTIONS :信息位编码，加白噪声，译码，测BER * ************************************************************************************************* */ #include "decode.h" #include "pn_seq.h" #include "channel.h" #include <iomanip> 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++; //加信道噪声 awgn_channel(b1_polar , channel_trans1 ,pi,l+m); awgn_channel(b2_polar , channel_trans2 ,pi,l+m); //接收端量化判决 大于等于0->1 小于0->0 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++; }; if (error)error_frame++; //统计误帧 } cout<<endl<<"总帧数:"<<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; }