举例如下：
假设所处的域为2^4的伽罗华域(本原多项式为x^4+x+1)
码型为(15,9),则n=15,k=9,t=3,d=7; 1 2 3 4 5 6 7 8 9
g(x) = (x+a1) * (x+a2) * (x+a3) * (x+a4) * (x+a5) * (x+a6);
g(x) = x^6+(a1+a2+a3+a4+a5+a6)*x^5+(a6*(a1+a2+a3+a4+a5)+a1*a2+a4*(a1+a2+a3)+a5*(a1+a2+a3+a4)+a3*(a1+a2))*x^4+(a6*(a1*a2+a4*(a1+a2+a3)+a5*(a1+a2+a3+a4)+a3*(a1+a2))+a4*(a1*a2+a3*(a1+a2))+a5*(a1*a2+a4*(a1+a2+a3)+a3*(a1+a2))+a1*a2*a3)*x^3+(a6*(a4*(a1*a2+a3*(a1+a2))+a5*(a1*a2+a4*(a1+a2+a3)+a3*(a1+a2))+a1*a2*a3)+a5*(a4*(a1*a2+a3*(a1+a2))+a1*a2*a3)+a1*a2*a3*a4)*x^2+(a6*(a5*(a4*(a1*a2+a3*(a1+a2))+a1*a2*a3)+a1*a2*a3*a4)+a1*a2*a3*a4*a5)*x+a1*a2*a3*a4*a5*a6
化简成GF(4)域的十进制元素后的码元发生多项式为：
g(x) = x^6+(a1+a2+a3+a4+a5+a6)*x^5+(a1*a2+a1*a3+a1*a4+a2*a3+a1*a5+a2*a4+a1*a6+a2*a5+a3*a4+a2*a6+a3*a5+a3*a6+a4*a5+a4*a6+a5*a6)*x^4+(a1*a2*a3+a1*a2*a4+a1*a2*a5+a1*a3*a4+a1*a2*a6+a1*a3*a5+a2*a3*a4+a1*a3*a6+a1*a4*a5+a2*a3*a5+a1*a4*a6+a2*a3*a6+a2*a4*a5+a1*a5*a6+a2*a4*a6+a3*a4*a5+a2*a5*a6+a3*a4*a6+a3*a5*a6+a4*a5*a6)*x^3+(a1*a2*a3*a4+a1*a2*a3*a5+a1*a2*a3*a6+a1*a2*a4*a5+a1*a2*a4*a6+a1*a3*a4*a5+a1*a2*a5*a6+a1*a3*a4*a6+a2*a3*a4*a5+a1*a3*a5*a6+a2*a3*a4*a6+a1*a4*a5*a6+a2*a3*a5*a6+a2*a4*a5*a6+a3*a4*a5*a6)*x^2+(a1*a2*a3*a4*a5+a1*a2*a3*a4*a6+a1*a2*a3*a5*a6+a1*a2*a4*a5*a6+a1*a3*a4*a5*a6+a2*a3*a4*a5*a6)*x+a1*a2*a3*a4*a5*a6
系数分别为
x^6系数 1
x^5系数 a1+a2+a3+a4+a5+a6
x^4系数 a1*a2+a1*a3+a1*a4+a2*a3+a1*a5+a2*a4+a1*a6+a2*a5+a3*a4+a2*a6+a3*a5+a3*a6+a4*a5+a4*a6+a5*a6
x^3系数 a1*a2*a3+a1*a2*a4+a1*a2*a5+a1*a3*a4+a1*a2*a6+a1*a3*a5+a2*a3*a4+a1*a3*a6+a1*a4*a5+a2*a3*a5+a1*a4*a6+a2*a3*a6+a2*a4*a5+a1*a5*a6+a2*a4*a6+a3*a4*a5+a2*a5*a6+a3*a4*a6+a3*a5*a6+a4*a5*a6
x^2系数 a1*a2*a3*a4+a1*a2*a3*a5+a1*a2*a3*a6+a1*a2*a4*a5+a1*a2*a4*a6+a1*a3*a4*a5+a1*a2*a5*a6+a1*a3*a4*a6+a2*a3*a4*a5+a1*a3*a5*a6+a2*a3*a4*a6+a1*a4*a5*a6+a2*a3*a5*a6+a2*a4*a5*a6+a3*a4*a5*a6
x^1系数 a1*a2*a3*a4*a5+a1*a2*a3*a4*a6+a1*a2*a3*a5*a6+a1*a2*a4*a5*a6+a1*a3*a4*a5*a6+a2*a3*a4*a5*a6
x^0系数 a1*a2*a3*a4*a5*a6
--------------------------------------------
这个是以前直观的展开后的多项式情况，但是发现这个变量设置有问题。
当初为了省掉一个”^“，结果使得结果不能化简。
下边为设置变量如下，将a1 a2等转换成a^1 a^2,则计算结果如下：
x^6系数 1
x^5系数 a^6 + a^5 + a^4 + a^3 + a^2 + a
x^4系数 a^11 + a^10 + 2*a^9 + 2*a^8 + 3*a^7 + 2*a^6 + 2*a^5 + a^4 + a^3
x^3系数 a^15 + a^14 + 2*a^13 + 3*a^12 + 3*a^11 + 3*a^10 + 3*a^9 + 2*a^8 + a^7 + a^6
x^2系数 a^18 + a^17 + 2*a^16 + 2*a^15 + 3*a^14 + 2*a^13 + 2*a^12 + a^11 + a^10
x^1系数 a^20 + a^19 + a^18 + a^17 + a^16 + a^15
x^0系数 a^21
进一步化简可得：
x^6系数 1
x^5系数 7
x^4系数 9
x^3系数 3
x^2系数 12
x^1系数 10
x^0系数 12
然后，发现×……&×……×&%78%*&%*&)(&*)(&*%&^&*^*(()........
可以直接用matlab计算，rsgenpoly(15,9)，狂汗。。。。
========================================
本来以为告一段落，但是发现结果和一篇论文里的结果有点出路，
本来应该怀疑论文可能打错了，但是还是对比了一下：
rsgenpoly(15,11)结果为 1 13 12 8 7，即 1 a^13 a^6 a^4 a^10
论文中的计算如下：
g(x)=(x+1)(x+2)(x+4)(x+8)
=(x^2+3x+2)(x+4)(x+8)
=(x^3+7x^2+14x+8)(x+8)
=x^4 + 15x^3 + 3x^2 + x + 12
对系数用指数形式也可以表达为：
g(x)=a^0*x + a^12*x + a^4*x + a^0*x + a^6
则结果为 1 15 3 1 12，即 a^0 a^12 a^4 a^0 a^6
居然不一样，通过一步一步展开发现
x^3的系数应该为 1 + a^1 + a^2 + a^3
按照这里的元素计算与论文的结果吻合，但是matlab的结果像是
a^1 + a^2 + a^3 + a^4 这样算出来的
这个还需要进一步确认一下。
|
本文地址：http://114er.blogspot.com/2010/11/rs.html 原创文章如转载，请注明链接： 转自Welcome Funny Guys |
0 评论: