课程设计任务书2011—2012学年第一学期专业：通信工程学号：080110501 姓名：李琼课程设计名称：信息论与编码课程设计设计题目：线性分组码编码的分析与实现完成期限：自2011 年12 月19 日至2011 年12 月25 日共 1 周一．设计目的1、深刻理解信道编码的基本思想与目的；2、理解线性分组码的基本原理与编码过程；3、提高综合运用所学理论知识独立分析和解决问题的能力；4、使用MATLAB或其他语言进行编程。
二．设计内容给定消息组M及生成矩阵G，编程求解其线性分组码码字。
三．设计要求编写的函数要有通用性。
四．设计条件计算机、MATLAB或其他语言环境五．参考资料[1]曹雪虹，张宗橙.信息论与编码.北京：清华大学出版社，2007.[2]王慧琴.数字图像处理.北京：北京邮电大学出版社，2007.指导教师（签字）：教研室主任（签字）：批准日期：年月日该系统是（6，3）线性分组码的编码的实现，它可以对输入的三位的信息码进行线性分组码编码。
当接收到的六位码字中有一位发生错误时，可以纠正这一位错码；当接收到的码字有两位发生错误时，只能纠正一位错误，但同时能检测出另一位错误不能纠正。
只有特定位有两位错误时，才能纠正两位错误。
这样就译出正确的信息码组，整个过程是用MATLAB语言实现的。
关键词：编码；MA TLAB；纠错1课程描述 02 设计原理 (1)2.1 线性分组码的编码 (1)2.1.1 生成矩阵 (1)2.1.2 校验矩阵 (3)2.2 伴随式与译码 (4)2.2.1 码的距离及纠检错能力 (4)2.2.2 伴随式与译码 (4)3 设计过程 (5)3.1 编码过程 (5)3.2 仿真程序 (7)3.4 结果分析 (11)总结 (13)致谢 (14)参考文献 (15)1课程描述线性分组码具有编译码简单，封闭性好等特点，采用差错控制编码技术是提高数字通信可靠性的有效方法，是目前较为流行的差错控制编码技术。
对线性分组码的讨论都在有限域GF(2)上进行，域中元素为{0,1},域中元素计算为模二加法和模二乘法。
分组码是一组固定长度的码组，可表示为（n , k ），通常它用于前向纠错。
在分组码中，监督位被加到信息位之后，形成新的码。
在编码时，k 个信息位被编为n 位码组长度，而n-k 个监督位的作用就是实现检错与纠错。
对于长度为n 的二进制线性分组码，它有种2n 可能的码组，从2n 种码组中，可以选择M=2k 个码组（k<n ）组成一种码。
这样，一个k 比特信息的线性分组码可以映射到一个长度为n 码组上，该码组是从M=2k 个码组构成的码集中选出来的，这样剩下的码组就可以对这个分组码进行检错或纠错。
要设计一个（6，3）线性分组码的编译码程序，最基本的是要具备对输入的信息码进行编码，让它具有抗干扰的能力。
同时，还要让它具有对接收到的整个码组中提取信息码组的功能。
但是，在实际的通信系统中，由于信道传输特性不理想以及加性噪声的影响，接收到的信息中不可避免地会发生错误，影响通信系统的传输可靠性，因而，本设计还要让该程序具有纠正错误的能力，当接收到的码组中有一位码，发生错误时可以检测到这一位错码，并且可以纠正这一位错码，并且让系统从纠正后的码组中提取正确的信息码组。
针对给定的矩阵Q=完成如下的工作：1 完成对任意信息序列的编码2 根据生成矩阵，形成监督矩阵；3 根据得到的监督矩阵，得到伴随式，并根据它进行译码；4 验证工作的正确性。
2 设计原理2.1 线性分组码的编码2.1.1 生成矩阵线性分组码（n，k）中许用码字（组）为2k个。
定义线性分组码的加法为模二加法，乘法为二进制乘法。
即1+1=0、1+0=1、0+1=1、0+0=0；1×1=1、1×0=0、0×0=0、0×1=0。
且码字与码字的运算在各个相应比特位上符合上述二进制加法运算规则。
线性分组码具有如下性质（n，k）的性质：1、封闭性。
任意两个码组的和还是许用的码组。
2、码的最小距离等于非零码的最小码重。
对于码组长度为n、信息码元为k位、监督码元为r＝n－k位的分组码，常记作（n，k）码，如果满足2r－1≥n，则有可能构造出纠正一位或一位以上错误的线性码。
下面我们通过（7，3）分组码的例子来说明如何具体构造这种线性码。
设分组码（n，k）中，k = 3，为能纠正一位误码，要求r≥3。
现取r＝4，则n＝k＋r ＝7。
该例子中，信息组为(c6c5c4),码字为(c6c5c4c3c2c1c0).当已知信息组时，按以下规则得到四个校验元，即c3=c6+c4c2=c6+c5+c4 （2-1）c1=c6+c5c0=c5+c4这组方程称为校验方程。
（7，3）线性分组码有23（8）个许用码字或合法码字，另有27-23个禁用码字。
发送方发送的是许用码字，若接收方收到的是禁用码字，则说明传输中发生了错误。
为了深化对线性分组码的理论分析，可将其与线性空间联系起来。
由于每个码字都是一个二进制的n重，及二进制n维线性空间Vn中的一个矢量，因此码字又称为码矢。
线性分组码的一个重要参数是码率r=k/n,它说明在一个码字中信息位所占的比重，r 越大，说明信息位所占比重越大，码的传输信息的有效性越高。
由于(n,k)线性分组，线性分组码的2k 个码字组成了n 维线性空间Vn 的一个K 维子空间。
因此这2k 个码字完全可由k 个线性无关的矢量所组成。
设此k 个矢量为c 1,c 2,…,c k ,有生成矩阵形式为（2-2）(n,k)码字中的任一码字c i ,均可由这组基底的线性组合生成，即 c i =m i ·G=[m n-1 m n-2 … m n-k ]·G式中，mi =[m n-1 m n-2 …m n-k ]是k 个信息元组成的信息组。
表2-1 （7，3）线性分组码对于表2-1给出的（7，3）线性分组码，可将写成矩阵形式[c 6 c 5 c 4 c 3 c 2 c 1 c 0]=[c 6 c 5 c 4]·101110011100100111001故（7，3）码的生成矩阵为G= 101110011100100111001可以看到，从（7，3）码的8个码字中，挑选出k=3个线性无关的码字（1001110）（0100111），（00111101）作为码的一组基底，用c=m ·G 计算得码字。
一个系统码的生成矩阵G ，其左边k 行k 列应是一个k 阶单位方阵I k ，因此生成矩阵G 表示为G=[I k P] （2-3） 式中，P 是一个k ×(n-k)阶矩阵。
2.1.2 校验矩阵表2-1所示的（7，3）线性分组码的四个校验元由式（2-1）所示的线性方程组决定的。
把（2-1）移相，有c6+c4+c3=0c6+c5+c4+c2=0 c6+c1+c5=0 （2-4） c5+c4+c0=0上式的矩阵形式为1000110110001100101110001101 · 0123456c c c c c c c =000这里的四行七列矩阵称为（7，3）码的一致校验矩阵，用H 表示，即1000110110001100101110001101H= （2-5）由H 矩阵得到(n,k)线性分组码的每一码字c i,(i=1,2,…,2k ),都必须满足由H 矩阵各行所确定的线性方程组，即 c i ·H T =0.（7，3）码的生成矩阵G 中每一行及其线性组合都是（n,k ）码的码字，所以有G ·H T =0。
由G 和H 构成的行生成的空间互为零空间，即G 和H 彼此正交。
H=[P T I r ]其右边r 行r 列组成一个单位方阵。
2.2 伴随式与译码2.2.1 码的距离及纠检错能力 1．码的距离两个码字之间，对应位取之不同的个数，称为汉明距离，用d 表示。
一个吗的最小距离d min 定义为d min =min{d(c i ,c j ),i ≠j,c i ,c j ∈(n,k)},两个码字之间的距离表示了它们之间差别的大小。
距离越大，两个码字的差别越大，则传送时从一个码字错成另一码字的可能性越小。
码的最小距离愈大，其抗干扰能力愈强。
2. 线性码的纠检错能力对于任一个（n,k ）线性分组码，若要在码字内（1） 检测出e 个错误，则要求码的最小距离d ≥e+1;(2) 纠正t 个错误，则要求码的最小距离d ≥2t+1;(3)纠正t 个错误同时检测e(≥t)个错误，则要求 d ≥t+e+1； 2.2.2 伴随式与译码假设接收端收到的码字为B ，那么它和原来发送端发送的码字A 之间就有可能存在着误差。
即在码组A={a 6 a 5 a 4 a 3 a 2 a 1 a 0 }中的任意一位就有可能出错。
这样我们在接收端接收到一个码组是就有可能判断错发送端原来应该要表达的意思。
为了描述数据在传输信道中出现错误的情况，引入了错误图样E ，在错误图样中，0代表对应位没有传错，1代表传输错误。
实际上错误图样E就是收序列与发送序列的差。
所以在译码中用接收到的码字B模尔加错误图样E就可以得到发送端的正确码字A。
因此译码的过程就是要找到错误图样E。
定义：校正子SS = B * H T= ( A + E ) * H T= A * H T+ E * H T= E * H T因为A是编得的正确码字。
根据前面所叙述，它和监督矩阵的转置相乘为0。
显然，S仅与错误图样有关，它们之间是一一对应的关系。
找到了校正子S，也就可以找到E。
而与发送的码字无关。
若E=0，则S=0；因此根据S是否为0可进行码字的检错。
如果接收码字B中只有一位码元发生错误，又设错误在第i位。
即E i-1=1,其他的E i均为0。
在后面的译码程序中，建立了一个校正子S与错误图样E对应的表。
也就是收到一个B序列，就可以通过计算得到一个校正子，而每一个校正子都对应着一个错误图样E，再通过B模尔加上E，就可以得到正确的码字A。
因为在不同的错误序列B中，同一位码元错误时对应的E是一样的，所以可以利用0000000这个正确的码字让它每位依次错误，来求得它的八个校正子。
而这时的矩阵B就是错误图样E。
这样就算得了8个校正子S。
而这时的错误序列B，就是错误图样E，所以有：E与S都已经得到，这时就可以建立一个表来将它们一一对应起来。
3 设计过程3.1 编码过程监督矩阵H与生成矩阵G的关系：由H与G的分块表示的矩阵形式H =[P In-kG=[IkP=Q T则有G ·H T =0 （3-4） 或 H ·G T =0 （3-5） 已知给出的（6，3）码的Q 矩阵（3-6）则可以根据G=[I k Q]求出生成矩阵（3-7）由 P=Q T 和 H =[P I n-k ] 即 MATLAB 函数为：C=rem(I*G ,2); （3-8）其中C 为编码后的结果，I 为信息矩阵，G 为生成矩阵。
