基于MATLAB的(7_4)汉明码编译码设计与仿真结果分析(DOC)

1、探2009级通信工程专业%通信原理课程设计通信原理课程设计报告课题名称基于MATLAB 的（7 , 4 ）汉明码编译码设计与仿真结果分析通信与电子工程学院通信工程指导教师2011年12月23日一、设计任务及要求:设计任务:利用MATLA编程，实现汉明码编译码设计。理解（7,4 ）汉明码的构造原 理，掌握（7, 4）汉明码的编码和译码的原理和设计步骤。并对其性能进行分析。要 求:通过MATLABS程，设计出（7,4 ）汉明码的编码程序，编码后加入噪声,然后译码，画出信噪比与误比特数和信噪比与误比特率的仿真图，然后对其结果 进行分析指导教师签名：2011年12月23日二、指导教师评语:指导教师签名

2、：年 月 日三、成绩验收盖章 年 月 日基于MATLAB£(乙4)汉明码编译码设计的与仿真结果分析1设计目的(1) 熟悉掌握汉明码的重要公式和基本概念。(2) 利用matlabs程，实现汉明码编译码设计。(3) 理解(7,4 )汉明码的构造原理，掌握(7, 4)汉明码的编码和译码的原理 和设计步骤。(4) 对其仿真结果进行分析。2设计要求(1) 通过MATLAB程，设计出(7,4 )汉明码的编码程序。(2) 编码后加入噪声，然后译码，画出信噪比与误比特数和信噪比与误比特率 的仿真图。(3) 然后对其结果进行分析。3设计步骤3.1 线性分组码的一般原理线性分组码的构造3.1.1 H矩阵

3、根据(7, 4)汉明码可知一般有(1)a © a © a © a 0Ia © a © a © a1 006 © a4 © as © ao = 0现在将上面它改写为1 06 +1 05 +1 04 +0 £3 +1 02 +0 q +0 a。=01 as +1 a5 +0 a4 +1 q +0 32+1 6+0 ao =0 »1 as +0 as +1 a4 +1 £3 +0 a0 a<1 a0 =0” 上式中已经将“ ©”简写成“ +”。上式可以表示成如下矩阵

4、形式：4上式还可以简记为H - AT11101001 1101010 1011001.站asa4asa2a1|_a0-00式中fl110110111011:100-010-001=Ipl r (8)(1110100111010101011001 .A = a6 a5 a4 a3 a2 a1 a00 = 000右上标“T”表示将矩阵转置。例如，HT是H的转置，即HT的第一行为 H的第一列，HT的第二行为H的第二列等等。将H称为监督矩阵。只要监督矩阵H给定，编码时监督位和信息位的关系就完全确定了。H矩阵的性质：1) H的行数就是监督关系式的数目，它等于监督位的数目 r。H的每行 中“1”的位置表示

5、相应码元之间存在的监督关系。例如，H的第一行1110100表示监督位a2是由a6 a5 a4之和决定的。H矩阵可以分成两部分，例如为r咒r阶单位方阵。我们将具有P I r形式中，P为r X k阶矩阵，Ir式的H矩阵称为典型阵。2)由代数理论可知，H矩阵的各行应该是线性无关的，否则将得不到r个线性无关的监督关系式，从而也得不到r个独立的监督位。若一矩阵能写成典型阵形式PI r，则其各行一定是线性无关的。因为容易验证I r的各行是线 性无关的，故P I r的各行也是线性无关的。3.1.2 G矩阵：上面汉明码例子中的监督位公式为$2 = as © a5 © 34a1 as 田 a

6、s 田 asIb。=a63 a4© as(9)5也可以改写成矩阵形式：82-81 1_80 .111011011011仟a5a4香3.(10)或者写成乩3!3。】=863534331111110=136858483 Q101L011(11)式中,Q为一个k X r阶矩阵，它为P的转置，即Q= PT 上式表示，在信息位给定后，用信息位的行矩阵乘矩阵Q就产生出监督位。我们将Q的左边加上1个k X k阶单位方阵，就构成1个矩阵G1000111 0100110 0010101 0001011.G称为生成矩阵，因为由它可以产生整个码组，即有a6a5a4a3a2a1a0 = a6a5a4a3 G

7、(12)(13)或者A= 86853433 G因此，如果找到了码的生成矩阵 G则编码的方法就完全确定了。具有 IkQ形式的生成矩阵称为典型生成矩阵。由典型生成矩阵得出的码组A中，信息位的位置不变，监督位附加于其后。这种形式的码称为系统码。G矩阵的性质：1) G矩阵的各行是线性无关的。因为由上式可以看出，任一码组A都是G的各行的线性组合。G共有k行，若它们线性无关，则可以组合出2k种不同的 码组A,它恰是有k位信息位的全部码组。若G的各行有线性相关的，则不可能 由G生成2k种不同的码组了。2) 实际上，G的各行本身就是一个码组。因此，如果已有 k个线性无关 的码组，则可以用其作为生成矩阵 G,并

8、由它生成其余码组。3.1.3校正子S当接收码组有错时，E H 0，将B当作A代入公式(A ' H T = 0)后，该 式不一定成立。在错码较多，已超过这种编码的检错能力时，B变为另一许用码(14)(15)组，则该式仍能成立。这样的错码是不可检测的。在未超过检错能力时，上式不 成立，即其右端不等于0。假设这时该式的右端为S,即B ” H T = S将B = A + E代入上式，可得(16)S = ( A + E) H T = A ” H T + E ” H T由于A - HT = 0,所以(17)S = E -H T式中S称为校正子。它能用来指示错码的位置。S和错码E之间有确定的线性变换

9、关系。若 S和E之间一一对应，则S 将能代表错码的位置。由于使用了一位监督位aO,它和信息位an-1a1 一起3.2 (7, 4)汉明码的构造原理在偶数监督码中，构成一个代数式：在接收端解码时,anj© an工©© a0 =0实际上就是在计算(18)7(19)则能增加一个类似的监督关系式。由于两 01，10，11,故能表示4种不同的信息。 3种组合就有可能用来指示一个错码的 3S = a© an _2 ©© a0若S = 0，就认为无错码；若S = 1 ,就认为有错码。现将上式称为监督关 系式，S称为校正子。由于校正子S只有两种取值，

10、故它只能代表有错和无错这 两种信息，而不能指出错码的位置。若监督位增加一位，即变成两位， 个校正子的可能值有4中组合：00, 若用其中1种组合表示无错，则其余 种不同位置。同理，r个监督关系式能指示1位错码的(2r - 1)个可能位置。一般来说，若码长为n,信息位数为k,则监督位数r = n k。如果希望用r 个监督位构造出r个监督关系式来指示1位错码的n种可能位置，则要求(20)2rT 沙或 2 1然后根据(7, 4)汉明码的编码原理，画出程序设计的流程图：然后根据流程图进行编写程序。4、程序设计clear allclose allN=30000;K=4*N; %言息长度Q=8; %调制形式

11、对（7,4）汉明码进行编码输出编码后的（7,4）汉明码对编码后的（7,4） 汉明码进行调制加入高斯白噪声对加入噪声后的（7,4）汉明码进行译码输出译码码后的（7,4）汉明码画出误码率/数与 信噪比的仿真图图1编码流程图9n um=zeros(20,1);%误比特数10ber=zeros(20,1);%误比特率for SNR=1:1:20 %生成矩阵G=0 0 0 1 0 1 1 ； %监督矩阵H=1 1 1 0 1 0 01 0 1 1 0 0 1;A=zeros(7*N,1);编码后信号S=zeros(N,3);%校验子E=zeros(7*N,1);%错误图样I=zeros(7*N,1);%

12、 解调纠错后的矩阵NEW=zeros(K,1);%收到的信号a=ra ndsrc(K,1,0,1);%信源%编码for k=1:Ncode=a(4*k-3) a(4*k-2) a(4*k-1) a(4*k)*G;for l=0:6A(7*k-I)=mod(code(7-I),2);end end %8PSK调制Xa=pskmod(A,Q);%信道传输加噪Xb=awg n(Xa,SNR);%解调B=p skdemod(Xb,Q);%编码for n=1:NM=7* n-6;Z=4* n-3;T=7* n-3;V=4* n;L=7* n;s=B(M:L,1)'*H：S(n ,:)=mod(s

13、,2);ifS( n,:)=0 0 0E(M:L,1)=0 0 0 0 0 0 0'elseifS( n,:)=0 0 1E(M:L,1)=0 0 0 0 0 0 1'elseifS( n,:)=0 1 0E(M:L,1)=0 0 0 0 0 1 0'elseifS( n,:)=0 1 1E(M:L,1)=0 0 0 0 1 0 0'elseif S(n ,:)=1 0 0E(M:L,1)=0 0 0 1 0 0 0'elseif S(n ,:)=1 0 1E(M:L,1)=0 0 1 0 0 0 0'elseifS( n,:)=1 1 0E(M

14、:L,1)=0 1 0 0 0 0 0'elseifS( n,:)=1 1 1E(M:L,1)=1 0 0 0 0 0 0'endP=B(M:L,1)+E(M:L,1);I(M:L,1)=mod( P,2);NEW(Z:V,1)=I(M:T,1);end n um(SNR,:),ber(SNR,:)=biterr(NEW,a);end13plot( num.'b:s'); %画岀误比特数xlabel('信噪比);ylabel('误比特数');figurepl ot(ber.'k-*');%画岀误比特率xlabel('

15、;信噪比);ylabel('误比特率);5、设计结果及分析5.1通过对以上程序的编译运行仿真可以得到以下两个仿真图。x 10数 特 比 误6Lit%bE(ht5- 4f口3E0211-LJn0rrrr.E厂匕5=口1口0口由024648101214161820信噪比图2汉明码译码误比特数图Cogjijd：普317率特比误0.10.058101214161820信噪比汉明码译码误比特率图5.2编码后信号A与经过校验子S校验后译码得到的信号截图Cojuzmd TindoTO to HAILAE' latch this VIdg s唱e 卩帥og or rewi Gwttitt宣 E

16、tart#图4编码后信号A的截图toi IMAILA5? If at ch this ¥i dmo. see EwmoN or read Getting StM；，图5译码得到的信号截图5.3仿真结果分析（1）观察图2和图3我们可以看到误比特数和误比特率随信噪比的增加而 减小，并最后趋于0.由图4和图5,再联系前面校正子s与各码元之间的对应关 系，我们可以发现校正子s的重要作用。并且我们还可以推算出当程序中 N取值 很小时画出来的图形将是锯齿减函数，也就是说随着N的增大观察到的图形越光 滑。（2）与其它检错纠错码相比较（7,4）汉明码是能够纠正1位错码且编码 效率较高的一种线性分组码6主要仪器与设