~信道编码~无差错传输条件下的线性分组码反推分析程序.docVIP

“线性分组码反推分析”实验报告 一、实验题目 线性分组码的反推分析 二、实验目的 1. 练习应用线性分组码理论分析和解决问题 2. 初步体验线性分组码的反推分析 三、实验环境 Matlab仿真实验环境。 四、实验原理 推断：n=3？（推断码长为3,4,5…） 从编码素材第1位起，将码元每3个分一组；构造【3列】【足够多的行】 矩阵，求秩；绘制函数图像（横坐标为所取的码字起点，纵坐标为求得的秩）； 从编码素材第2位（从不同位置开始用于推断码字起点），3个分一组；构造【3列】【足够多的行】 矩阵，求秩；绘图； 取不同的起始位，不断重复上述操作，得到n=3对应的S-R曲线 类似的得到n=4,n=5,……,n=35对应的S-R曲线 并让它们在同一幅图中显示出来 根据判断原则： 当推断到码长或码长的倍数时，所得曲线为规则折线，其余为直线。 从第一条规则折线上读取最小秩的值，即为k值；该秩对应的横坐标即为码字起点。 可以找出规则折线对应的码长n，信息位数k（即最小的秩），码字起点d。 从码字起点开始，将码元每n个分为一组，取k组线性无关的码元，作为生成矩阵的k行，即得到生成矩阵g，对其标准化为标准生成矩阵G，在对G部分转置即得出标准校验矩阵H 五、实验步骤 1.绘制S-R曲线 根据所给的编码素材code ,假设码长n=3~35之间的一个值，起始位置d=1~100之间的一个数。 建立两重循环，求出n,d为定值时对应的码矩阵x，并将x的秩存入y(d,n-2)得到一个存贮了n=3:35,d=1:100的x的秩矩阵，再用指令plot（y）画出矩阵的图像即S-R图 2.求n,k,d 从y矩阵第一列开始，连续取i列，用rank()函数对其求秩，由于推断的码长不等于实际码长时，该列的S-R图像为平行于横轴的直线，体现在矩阵中即该列的每一行为相同的常数。因而有当i+2=n-1时，rank(y(1:100,1:i)) =1,直到i+2=n时， rank(y(1:100,1:i))~=1，跳出循环，此时n=i+2,k即第i列的最小值，在第i列里找到第一个值等于k的元素的行号，即起始位置d 3，求标准生成矩阵G，标准校验矩阵H 1) 从code中以d为起点，每n个一组，取 k组线性无关的码字，作为g的行向量，即得到非标准的生成矩阵g 2) 对矩阵g标准化，即可得到标准生成矩阵G. （ 这里将原函数库里的函数rref()稍作修改，在其中运算中加上mod(…,2)建立一个新的函数gfrref()可对二元域上的任意矩阵进行标准化的函数） 3) G部分转置得到标准校验矩阵H. (这里使用函数gen2 par()) 六、程序代码 clear load(D:\实验@大作业\编码素材\6\code.mat); for n=3:35 %假设码字长度为n for d=1:100 %假设码字起点为d for l=1:100 %从 code中以d为起点，每n个码元为一组取100组，得到码矩阵x x(l,1:n)=code(n*(l-1)+d:n*l+d-1); end y(d,n-2)=gfrank(x); %y为存储矩阵x的秩的矩阵，其中y(i，j)对应n=j+2，d=i是的矩阵x的秩 end end plot(y) %求码字长度n，及信息位个数k for i=2:48 if(rank(y(1:100,1:i))~=1) break end end n=i+2 k=min(y(1:100,n-2)) %求码字起点d for d=1:100 if(y(d,i)==k) break end end %求非标准生成矩阵g p=0; q=1; while q~=k+1; l=p+q; g(q,1:n)=code((l-1)*n+d:l*n+d-1); q=q+1; if(gfrank(g)==q-2) p=p+1; q=q-1; end end G=gfrref(g) %求标准生成矩阵G H=gen2par(G) %G部分转置得到标准校验矩阵H 七、实验结果 运行在matlab中运行程序，得出实验结