基于MATLAB——汉明码设计与实现

1、实验报告书-汉明码设计与实现汉明码编译码器系统班级： 姓名： 学号：1. 实验原理描述1.1汉明码编码原理一般来说，若汉明码长为n，信息位数为k，则监督位数r=n-k。若希望用r个监督位构造出r个监督关系式来指示一位错码的n种可能位置，则要求或下面以（7，4）汉明码为例说明原理： 设汉明码（n,k）中k=4，为了纠正一位错码，由式（1）可知，要求监督位数r3。若取r=3,则n=k+r=7。我们用来表示这7个码元，用的值表示3个监督关系式中的校正子，则的值与错误码元位置的对应关系可以规定如表1所列。表1 校正子和错码位置的关系错码位置错码位置001101010110 100111011000无错

2、码则由表1可得监督关系式：S1=a6a5a4a2 (2)S2=a6a5a3a1(3) S3=a6a4a3a0(4)在发送端编码时，信息位的值决定于输入信号，因此它们是随机的。监督位、应根据信息位的取值按监督关系来确定，即监督位应使式（2）式（4）中、的值为0（表示编成的码组中应无错码） a6a5a4a2=0a6a5a3a1=0a6a4a3a0=0(5)式（5）经过移项运算，接触监督位a2=a6a5a4a1=a6a5a3a0=a6a4a3(6)式（5）其等价形式为：1 1 1 0 1 0 01 1 0 1 0 1 01 0 1 1 0 0 1a6a5a4a3a2a1a0=000(7)式（6）还可

3、以简记为HAT=0T或AHT=0 (8)其中H=1 1 1 0 1 0 01 1 0 1 0 1 01 0 1 1 0 0 1A=a6a5a4a3a2a1a0P=1 1 1 01 1 0 11 0 1 1Ir=1 0 00 1 00 0 10=0 0 0所以有H=PIr式（6）等价于a2a1aa=a6 a5 a4 a31 1 11 1 01 0 10 1 1=a6 a5 a4a3Q(10)其中Q为P的转置，即式（10）表示,信息位给定后，用信息位的行矩阵乘矩阵Q就产生出监督位。我们将Q的左边加上一个k×k阶单位方阵，就构成一个矩阵GG=IkQ=1 0 0 0 1 1 10 1 0 0

4、 1 1 00 0 1 0 1 0 10 0 0 1 0 1 1(12)G称为生成矩阵，因为由它可以产生整个码组，即有a6a5a4a3a2a1a0=a6a5a4a3G(13) 或者A=a6a5 a4 a3G（14）式(13)即汉明码的编码原理1.2汉明码纠错原理当数字信号编码成汉明码形式（本文中即A）后在信道中传输，由于信道中噪声的干扰，可能由于干扰引入差错，使得接收端收到错码，因此在接收端进行汉明码纠错，以提高通信系统的抗干扰能力及可靠性。 一般来说接收码组与A不一定相同。若设接收码组为一n列的行矩阵B，即B=a6a5a4a3a2a1a0(15)则发送码组和接收码组之差为B-A=E(16)E

5、就是传输中产生的错码行矩阵E=e6e5e4e3e2e1e0(17)若ei=0，表示接收码元无错误，若ei=1，则表示该接收码元有错。式（16）可改写成B=A+E (18)若E=0，即接收码组无错,则，将它代人式（8），该是仍成立，即有BHT=0(19)当接收码组有错时，E0，将B带入式（8）后，该式不一定成立。在未超过检错能力时，式（19）不成立。假设此时式（19）的右端为S,即BHT=S(20)将B=A+E代入式（20），可得S=A+EHT+EHT由式（8）可知，所以S=EHT(21)此处S与前面的有着一一对应关系，则S能代表错码位置。因此，纠错原理即，接收端收到码组后按式（20）计算出S,

6、再根据表1判断错码情况，进行差错纠正。二.实验仪器1.通信原理综合实验系统一台2.电脑-MATlab一台三.实验目的1.熟悉掌握汉明码的原理与实现2.观察了解汉明距离的作用3.通过已经知道的汉明码监督方程为传输的编码进行纠错也验证四.实验内容1.汉明码编码规则汉明码是1950年由美国贝尔实验室提出来的，是第一个设计用来纠正一位误码的线性分组码，汉明码及其变型已广泛应用于数字通信和数据存储系统中作为差错控制码。汉明码的原始设计思想来自于前面讨论的奇偶监督码。通过一个例子来说明如何具体构造这些监督关系式。 设分组码 (n,k)中k = 4。为了纠正一位错码，按汉明不等式可得r3，若取r =3，则n

7、 = kr =7。我们用 a6 a5 a4 a3 a2 a1 a0 a2 a6 a5 a4 a1 a6 a5 a3 a0 a6 a4 a32.编码的检验与纠错用S1，S2，S3表示三个监督关系式式中的校正子，监督方程为： S1 =a6 a5 a4 a2 S2 =a6 a5 a3 a1 S3 =a6 a4 a3 a0则S1，S2，S3的值与错码位置的对应关系可以规定如下表：五.实验要求1.编写7位汉明码的程序，输出汉明码function f=hammingencod(a) G=1 0 0 0 1 1 1;0 1 0 0 1 1 0;0 0 1 0 1 0 1;0 0 0 1 0 1 1; a=i

8、nput('输入信息元序列:'); c=mod(a*G,2); disp('编码后序列为:'); disp(c);x=.01:.01:4;m,n=size(a'*ones(1,100);y=reshape(a'*ones(1,100)',1,m*n);plot(x,y)axis(0 4 0 1.5);set(gca,'XTick',0:1:4);set(gca,'YTick',0:0.5:1.5);title('hanmingencode')xlabel('value')yl

9、abel('value')end输入信息元序列:1101编码后序列为: 1 0 0 0 1 1 1 0 1 0 0 1 1 0 0 0 1 0 1 0 1 0 0 0 1 0 1 12.任意输入一个7位编码使用程序判断编码是否正确，如果错误，指出错位并纠正。function g=hammingdecod(B) H=1 1 1 0 1 0 0 ;1 1 0 1 0 1 0;1 0 1 1 0 0 1; B=input('输入接收序列B='); S=mod(B*H',2); %计算B的伴随式if S=0 disp('接收到的码字无错误。');

10、 E=dec2bin(0,7); endfori=1:1:7if S=H(:,i)' E=dec2bin(2(7-i),7); %计算R的错误图样fprintf('错误出现在第%1.0f位n',i); break; endend a=mod(B-E,2); %计算原发送码序列disp('原发送码字为:'); disp(a) x=.01:.01:7;m,n=size(a'*ones(1,100);y=reshape(a'*ones(1,100)',1,m*n);m,n=size(B'*ones(1,100);z=reshape(B'*ones(1,100)',1,m*n);plot(x,y)holdon;plot(x,z,'-r')axis(0 7 0 1.5);set(gca,'XTick',0:1:7);set(gca,'YTick',0:0.5:2.5);set(gca,'ZTick',0:0.5:2.5);title('hanmingdecode')xlabel('value')ylabel(