课件：Malab四.pptVIP

4.4.2 循环码的仿真模块 4.4.3 循环码的仿真模块参数 循环码的解码仿真模块 循环码的解码仿真模块参数 cyclpoly 功能：生成循环码的生成多项式 语法：p=cyclpoly(N,K); p=cyclpoly(N,K,fd_flag); fd_flag=’min’时，寻找最小阶数循环生成多项式； fd_flag=’max’时，寻找最大阶数循环生成多项式； fd_flag=’all’时，对所有给定阶数的生成多项式进行搜索； fd_flag=L时，搜索有L项的生成多项式； 如果没有找到满足条件的多项式，p为空输出。 4.4.4 循环码的仿真函数 c1 = cyclpoly(15,4,all) c2 = cyclpoly(15,4,max) 4.4.4 循环码的仿真编码函数 code=encode(msg，N，K，‘cyclic’，p) [h,g]=cyclgen(n,p) 用生成多项式生成循环码的生成矩阵g和校验矩阵h code = mod(msg*g,2); clear all n = 3; k = 2; N=10000; msg = randint(N,k); pol=cyclpoly(n,k); [h,g]=cyclgen(n,pol); code1 = encode(msg,n,k,cyclic/binary); code2 = mod(msg*g,2); noisy=randerr(N,n,[0 1;0.7 0.3]); noisycode1 = mod(code1 + noisy, 2); noisycode2 = mod(code2 + noisy,2); newmsg1 = decode(noisycode1,n,k,cyclic); newmsg2 = decode(noisycode2,n,k,cyclic); [number,ratio1] = biterr(newmsg1,msg); [number,ratio2] = biterr(newmsg2,msg); disp([The bit error rate1 is ,num2str(ratio1)]) disp([The bit error rate2 is ,num2str(ratio2)]) 作业 3 仿真（7，3）循环编码/译码过程。（用siumulink模块做）。 4 仿真（7，4）循环编码/译码过程，仿真中使用最小阶数循环生成多项式。（用matlab函数做）。 THANK YOU SUCCESS * * 可编辑 4.3 线性分组码 4.4.1 线性分组码原理 反映线性分组码纠错能力的重要参数是最小（汉明）码距，即任何两个不同码字间的最小汉明距离。在线性码中，最小距离等于该码的最小码重。 4.3.1 线性分组码解码原理 从长为N的码字中恢复出长为K的信元需要经过一个矩阵运算－－伴随式计算。 假定接收到的信号是码字v和信道中产生的错误序列e的叠加，即r＝v+e，则译码过程包括下面四个步骤： 计算r的伴随式的值；s=rHT 利用伴随式的值确定错误的位置； 将接收到的矢量还原为码字矢量 v＝r-e； 从处理过的码字中恢复出原始的信息码元。 4.4.2 线性分组码解码原理 伴随式计算公式为 s=rHT 由于r＝v+e并且v＝uG，得到 S=rHT=(v+e) HT。 =(uG+e) HT =uG HT+e HT =e HT 4.4.2 线性分组码解码原理 当生成矩阵 ，校正子与真值表的关系 序号 错误码位 真值表E 校正子S e6 e5 e2 e3 e2 e1 e0 s2 s1 s0 0 无 0 0 0 0 0 0 0 0 0 0 1 b4 0 0 1 0 0 0 0 0 0 1 2 b5 0 1 0 0 0 0 0 0 1 0 3 b2 0 0 0 0 1 0 0 0 1 1 4 b6 1 0 0 0 0 0 0 1 0 0 5 b0 0 0 0 0 0 0 1 1 0 1 6 b3 0 0 0 1 0 0 0 1 1 0 7 b1 0 0 0 0 0 1 0 1 1 1 4.3.1 线性分组码解码原理 通信工具箱为对称线性分组码提供了一个函数gen2par，该函数可以由生成矩阵G计算出校验矩阵H