完整word版基于MATLAB的7 4汉明码编译码设计与仿真结果分析

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.

3、1 H矩阵根据(7, 4)汉明码可知一般有363534320(1)3635333103634333001361 351340331320310301361 350341330321310301360 35134133032031130简写成“ +”0000a ”现在将上面它改写为上式中已经将上式可以表示成如下矩阵形式:3上式还可以简记为H AT式中a6a51110100110101010110011110100H 11010101011001a4a3a2a1ao00(模 20A = a6 a5 a4 a3 a2 a1 a0 0 = 000右上标“T”表示将矩阵转置。例如，HT是H的转置，即HT

4、的第一行为 H的第一列，HT的第二行为H的第二列等等。将H称为监督矩阵。只要监督矩阵H给定，编码时监督位和信息位的关系就完全确定了。H矩阵的性质：1) H的行数就是监督关系式的数目，它等于监督位的数目 r。H的每行H的第一行1110100中“1”的位置表示相应码元之间存在的监督关系。例如, 表示监督位a2是由a6 a5 a4之和决定的。H矩阵可以分成两部分，例如1110100H11010101011001k阶矩阵，Ir为rPIr(8)式中，P为r式的H矩阵称为典型阵。2)由代数理论可知，H矩阵的各行应该是线性无关的，否则将得不到r个线性无关的监督关系式，从而也得不到r个独立的监督位。若一矩阵能

5、写成典型阵形式PI r，则其各行一定是线性无关的。因为容易验证I r的各行是线 性无关的，故P I r的各行也是线性无关的。3.1.2 G矩阵：上面汉明码例子中的监督位公式为r阶单位方阵。我们将具有P I r形a2a6a5 a4a1a6a5a3a。a63 a4a3(9)5也可以改写成矩阵形式:a2a1a。111011011011a6a5(10)或者写成a2 a1a0a6a5a4a3111110101011a4a3a6a5a4a3 Q(11)式中,Q就产生出Q为一个k r阶矩阵，它为P的转置，即Q= PT上式表示，在信息位给定后，用信息位的行矩阵乘矩阵监督位。我们将Q的左边加上1个k k阶单位方

6、阵，就构成1个矩阵G1000 1110100 1100010 1010001 011G称为生成矩阵，因为由它可以产生整个码组，即有36 35 34 3332313036353433 GG I kQ(12)(13)A26858483 G或者(14)因此，如果找到了码的生成矩阵 G则编码的方法就完全确定了。具有 IkQ形式的生成矩阵称为典型生成矩阵。由典型生成矩阵得出的码组A中，信息位的位置不变，监督位附加于其后。这种形式的码称为系统码。G矩阵的性质：1) G矩阵的各行是线性无关的。因为由上式可以看出，任一码组A都是G的各行的线性组合。G共有k行，若它们线性无关，则可以组合出2k种不同的 码组A,

7、它恰是有k位信息位的全部码组。若G的各行有线性相关的，则不可能 由G生成2k种不同的码组了。2) 实际上，G的各行本身就是一个码组。因此，如果已有k个线性无关的码组，则可以用其作为生成矩阵 G,并由它生成其余码组。3.1.3校正子S当接收码组有错时，E 0，将B当作A代入公式(A H T = 0)后，该 式不一定成立。在错码较多，已超过这种编码的检错能力时，B变为另一许用码(15)组，则该式仍能成立。这样的错码是不可检测的。在未超过检错能力时，上式不 成立，即其右端不等于0。假设这时该式的右端为S,即B H T = S将B = A + E代入上式，可得(16)S = ( A + E) H T

8、= A H T + E H T由于A HT = 0,所以(17)S = E H T式中S称为校正子。它能用来指示错码的位置。S和错码E之间有确定的线性变换关系。若 S和E之间一一对应，则S 将能代表错码的位置。3.2 (7, 4)汉明码的构造原理aO,它和信息位an-1a1 一起由于使用了一位监督位在偶数监督码中，构成一个代数式：在接收端解码时,an 1 an 2实际上就是在计算a。(18)San 1an 2a。(19)若S = 0，就认为无错码；若系式，S称为校正子。由于校正子两种信息，而不能指出错码的位置。若监督位增加一位，即变成两位，个校正子的可能值有4中组合：00,若用其中1种组合表示

9、无错，则其余S = 1S只有两种取值，故它只能代表有错和无错这，就认为有错码。现将上式称为监督关则能增加一个类似的监督关系式。由于两 01，10，11,故能表示4种不同的信息。 3种组合就有可能用来指示一个错码的 3#种不同位置。同理，r个监督关系式能指示1位错码的(2r - 1)个可能位置。一般来说，若码长为n,信息位数为k,则监督位数r = n k。如果希望用r个监督位构造出r个监督关系式来指示1位错码的n种可能位置，则要求(20)2r 1 n 或 2 k r 1然后根据(7, 4)汉明码的编码原理，画出程序设计的流程图：然后根据流程图进行编写程序。4、程序设计clear allclose

10、 allN=30000;K=4*N; %言息长度Q=8; %调制形式对（7,4）汉明码进行编码输出编码后的（7,4）汉明码对编码后的（7,4） 汉明码进行调制加入高斯白噪声对加入噪声后的（7,4） 汉明码进行译码I输出译码码后的（7,4）汉明码画出误码率/数与 信噪比的仿真图图1编码流程图7n um=zeros(20,1);%误比特数9ber=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

11、=zeros(7*N,1);%错误图样I=zeros(7*N,1);% 解调纠错后的矩阵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

12、* n;s=B(M:L,1)*H：S(n ,:)=mod(s,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

13、1 0E(M: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);end11plot( num.b:s); %画岀误比特数xlabel(信噪比);ylabel(误比特数);figurepl ot(ber.k-*);%画岀误比特率xlabel(信噪比);ylabel(误比特率);5、设计结果及分析5.1通过对以上程序的编译运

14、行仿真可以得到以下两个仿真图。4x 10数特比误008 10 12信噪比口 E Q 由14161820图2汉明码译码误比特数图率特比误0.150.10.058101214 W W 20信噪比图35.2编码后信号A与经过校验子汉明码译码误比特率图S校验后译码得到的信号截图Cffuand TindovD 1T.* to HM1AE7 Wst cL thi* Vi占一or :rgJ. G-tting t,Th图4编码后信号A的截图Cauand Vindav_to MATLAS? fatcli this Ti dwo, see IkniQs, qr Tsad Gattiag Staj-E -图5译码得到的信号截图5.3仿真结果分析（1）观察图2和图3我们可以看到误比特数和误比特率随信噪比的增加而 减小，并最后趋于0.由图4和图5,再联系前面校正子s与各码元之间的对应关 系，我们可以发现校正子s的重要作用。并且我们还可以推算出当程序中 N取值 很小时画出来的图形将是锯齿减函数，也就是说随着N的增大观察到的图形越光 滑。（2）与其它检错纠错码相比较（7,4）汉明码是能够纠正1位错码且编码 效率较高的一种线性分组码6主要仪器与设备安装了 matlab R2009a程序