卷积码编码和维特比译码.doc

卷积码编码维特比 SUN 一、实验目的 掌握卷积码编码和维特比 二、实验原理 1.卷积码是由连续输入的信息序列得到连续输出的已编码序列。其编码器将k个信息码元编为n个码元时，这n个码元不仅与当前段的k个信息有关，而且与前面的（m-1）段信 息有关（m为编码的约束长度）。 2.一般地，最小距离d表明了卷积码在连续m段以内的距离特性，该码可以在m个连续码流内纠正(d-1)/2个错误。卷积码的纠错能力不仅与约束长度有关，还与采用的译码方式有关。 3. 维特比译码算法基本原理是将接收到的信号序列和所有可能的发送信号序列比较，选择其中汉明距离最小的序列认为是当前发送序列。 卷积码的Viterbi 译码是根据接收码字序列寻找编码时通过网格图最佳路径的过程,找到最佳路径即完成了译码过程,并可以纠正接收码字中的错误比特。 4.所谓“最佳”, 是指最大后验条件概率:P( C/ R) = max [ P ( Cj/ R) ] , 一般来说, 信道模型并不使用后验条件概率,因此利用Beyes 公式、根据信道特性出结论:max[ P ( Cj/ R) ]与max[ P ( R/ Cj) ]等价。考虑到在系统实现中往往采用对数形式的运算,以求降低运算量,并且为求运算值为整数加入了修正因子a1 、a2 。 令M ( R/ Cj) = log[ P ( R/ Cj) ] =Σa1 (log[ P( Rm/ Cmj ) ] + a2) 。 其中, M 是组成序列的码字的个数。因此寻找最佳路径, 就变成寻找最大M( R/ Cj) , M( R/ Cj) 称为Cj 的分支路径量度,含义为发送Cj 而接收码元为R的似然度5.卷积码的viterbi译码是根据接收码字序列寻找编码时通过网格图最佳路径的过程，找到最佳路径即完成了译码过程并可以纠正接收码字中的错误比特。 三、实验代码 #includestdio.h #include Conio.h #define N 7 #include math.h #include stdlib.h #includetime.h #define randomize() 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;jnbytes;j++) { input=*data; data++; *symbols = input^a1^a2^a3^a6; //c1(171) symbols++; *symbols = input^a2^a3^a5^a6; //c2(133) symbols++; a2=a1; a1=input; } return 0; } int trandistance(int m, int state1, int state2) /*符号m与从state1到state2时输出符号的汉明距离,如果state1无法到state2则输出度量值为100*/ { int c; int sym,sym1,sym2; sym1=((state21)1)^(state21)^(state11); sym2=((state21)1)^(state11); sym=(sym11) | sym2; if ( ((st