《卷积码编码和维特比译码(共9页).doc》由会员分享，可在线阅读，更多相关《卷积码编码和维特比译码(共9页).doc（9页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上卷积码编码维特比译码实验设计报告SUN一、实验目的掌握卷积码编码和维特比译码的基本原理，利用了卷积码的特性, 运用网格图和回溯以得到译码输出。二、实验原理1.卷积码是由连续输入的信息序列得到连续输出的已编码序列。其编码器将k个信息码元编为n个码元时，这n个码元不仅与当前段的k个信息有关，而且与前面的（m-1）段信 息有关（m为编码的约束长度）。2.一般地，最小距离d表明了卷积码在连续m段以内的距离特性，该码可以在m个连续码流内纠正(d-1)/2个错误。卷积码的纠错能力不仅与约束长度有关，还与采用的译码方式有关。3. 维特比译码算法基本原理是将接收到的信号序列和所有可能
2、的发送信号序列比较，选择其中汉明距离最小的序列认为是当前发送序列。 卷积码的Viterbi 译码是根据接收码字序列寻找编码时通过网格图最佳路径的过程,找到最佳路径即完成了译码过程,并可以纠正接收码字中的错误比特。4.所谓“最佳”, 是指最大后验条件概率:P( C/ R) = max P ( Cj/ R) , 一般来说, 信道模型并不使用后验条件概率,因此利用Beyes 公式、根据信道特性出结论:max P ( Cj/ R) 与max P ( R/ Cj) 等价。考虑到在系统实现中往往采用对数形式的运算,以求降低运算量,并且为求运算值为整数加入了修正因子a1 、a2 。 令M ( R/ Cj)
3、= log P ( R/ Cj) =a1 (log P( Rm/ Cmj ) + a2) 。 其中, M 是组成序列的码字的个数。因此寻找最佳路径, 就变成寻找最大M( R/ Cj) , M( R/ Cj) 称为Cj 的分支路径量度,含义为发送Cj 而接收码元为R的似然度。5.卷积码的viterbi译码是根据接收码字序列寻找编码时通过网格图最佳路径的过程，找到最佳路径即完成了译码过程并可以纠正接收码字中的错误比特。三、实验代码#include #include Conio.h #define N 7 #include math.h #include #include #define rando
4、mize() srand(unsigned)time(NULL) encode( unsigned int *symbols, /*编码输出*/ unsigned int *data, /*编码输入*/ unsigned int nbytes, /*nbytes=n/16,n为实际输入码字的数目*/ unsigned int startstate /*定义初始化状态*/ ) /卷积码编码/ unsigned int j; unsigned int input,a1=0,a2=0,a3=0,a4=0,a5=0,a6=0; for(j=0;j1)&1)(state2&1)(state1&1); s
5、ym2=(state21)&1)(state1&1); sym=(sym11)=(state2&1) c=(m&1)(sym&1)+(m 1)&1)(sym 1)&1); else c=10000; return(c); int traninput(int a,int b) /*状态从a到b时输入卷积码的符号*/ int c; c=(b&2)1); return(c); int tranoutput(int a,int b) /*状态从a到b时卷积码输出的符号*/ int c,s1,s2; s1=(a&1)(a&2)1)(b&2)1); s2=(a&1)(b&2)1); c=(s11)|s2;
6、 return(c); /维特比译码/void viterbi( int initialstate, /*定义解码器初始状态*/ int *viterbiinput, /*解码器输入码字序列*/ int *viterbioutput /*解码器输出码字序列*/ ) struct sta /*定义网格图中每一点为一个结构体,其元素包括*/ int met; /*转移到此状态累计的度量值*/ int value; /*输入符号 */ struct sta *last; /*及指向前一个状态的指针*/ ; struct sta state4N; struct sta *g,*head; int i,
7、j,p,q,t,r,u,l; for(i=0;i4;i+) /* 初始化每个状态的度量值*/ for(j=0;jN;j+) stateij.met=0;for(l=0;l4;l+) statel0.met=trandistance(*viterbiinput,initialstate,l); statel0.value=traninput(initialstate,l); statel0.last=NULL; viterbiinput+; /*扩展第一步幸存路径*/ for(t=1;tN;t+) for(p=0;p4;p+) statept.met=state0t-1.met+trandist
8、ance(*viterbiinput,0,p); statept.value=traninput(0,p); statept.last=&state0t-1; for(q=0;q4;q+) if(stateqt-1.met+trandistance(*viterbiinput,q,p)0;u-) /*向前递归的找出最大似然路径 */ *(viterbioutput+(u-1)=g-value; g=g-last; /* for(u=0;u8;u+) *(viterbioutput+u)=stateu2.met; */ /*此行程序可用于检测第n列的 度量值*/ void decode(unsi
9、gned int *input, int *output,int n) int viterbiinput100; int j; for(j=0;jn+2;j+) viterbiinputj=(inputj*21)|inputj*2+1; viterbi(0,viterbiinput,output); void main() unsigned int encodeinput100,wrong10=0,0,0,0,0,0,0,0,0,0,encodeoutput100; int n=5,i,m,j=0,decodeinput100,decodeoutput100; randomize(); for
10、(i=0; in; i+) encodeinputi=rand()%2; encodeinputn= encodeinputn+1=0; encode(encodeoutput,encodeinput,n+2,0); printf(the input of encoder is :n); /信息源输入的信息码（随机产生） for(i=0;in; i+) printf(%2d,encodeinputi); printf(n); printf(the output of encoder is :n); /编码之后产生的卷积码 for(i=0;i(n+2)*2;i+) printf(%2d,enco
11、deoutputi); if(i%20=19) printf(n); printf(n); printf(please input the number of the wrong bitn); /信道传输收到干扰而产生的错误码 scanf(%d,&m); printf(please input the positions of the wrong bit(0-9)n); for(i=0;im;i+) scanf(%d,&wrongm); if(encodeoutputwrongm=0) encodeoutputwrongm=1; else encodeoutputwrongm=0; print
12、f(the input of decoder is :n); for(i=0;i(n+2)*2;i+) printf(%2d,encodeoutputi); if(i%20=19) printf(n); printf(n); decode(encodeoutput,decodeoutput,n+2); printf(the output of decoder is :n); for(i=0;in;i+) printf(%2d,decodeoutputi); printf(n); for(i=0;in;i+) if(encodeinputi!=decodeoutputi) j+; printf(the number of incorrect bit is:%dn,j); 四、实验总结（1）了解实验原理，分析实验所占数组变量很重要，也是相对考虑较多的；（2）对于读写文件，通过本实验更加熟悉；（3）记录实验程序最佳路径是本实验的难点；专心-专注-专业