1 有限域基础知识
1.1 有限域(Galois域)的构造
令
注:任意一个有限域,其元素的个数一定为
例1(有限域
GF(p) ) 令p 为一个素数,集合
GF(p)=Zp={0,1,2,…,p−1}.
在GF(p) 上定义加法⊕ 和乘法⊙ 分别为模p 加法和模p 乘法,即任意的a,b∈GF(p) ,
a⊕b=(a+b)modp, a⊙b=(a⋅b)modp
则<GF(p),⊕,⊙> 为一个有p 个元素的有限域,其中零元素为0 ,单位元为1 .令
a 为GF(p) 中的一个非零元素. 由于gcd(a,p)=1 ,因此,存在整数b,c ,使得ab+pc=1 . 由此得到a 的逆元为a−1=bmodp .域
GF(p) 称为一个素域(prime field).
例注1: 给定
例2(有限域
GF(pn) ) 从GF(p) 出发,对任意正整数n ,n≥2 ,我们可以构造元素元素个数为pn 的有限域GF(pn) 如下:
令g(x) 为一个GF(p) 上次数为n 的不可约多项式,集合
GF(pn)=GF(p)[x]/⟨g(x)⟩={a0+a1x+a2x2+⋯+an−1xn−1 | ai∈GF(p),0≤i≤n−1}
在GF(pn) 上定义加法⊕ 和乘法⊙ 分别为模g(x) 加法和模g(x) 乘法,即任意的a(x),b(x)∈GF(pn) ,
a(x)⊕b(x)=a(x)+b(x), a(x)⊙b(x)=(a(x)⋅b(x))modg(x)
则<GF(pn),⊕,⊙> 为一个有pn 个元素,特征为p 的有限域,其中零元素为GF(p) 中的0 ,单位元为GF(p) 中的1 .令
a(x) 为GF(pn) 中的一个非零元素. 由于gcd(a(x),g(x))=1 ,因此,存在GF(p) 上的多项式b(x),c(x) ,使得a(x)b(x)+g(x)c(x)=1 . 由此得到a(x) 的逆元为a−1(x)=b(x)modg(x) .域
GF(pn) 称为GF(p) 的(n 次)扩域(extension field),而GF(p) 称为GF(pn) 的子域(subfield).
例注2.1: 给定
例注2.2:设
其中
1.2 有限域的性质
令
若
并且
定义:设
GF(q) 是一个含有q 个元素的有限域,GF(p) 是GF(q) 的一个含有p 个元素的子域(p 不一定为素数),α∈GF(q) . 则GF(p) 上以α 为根,首项系数为1 ,并且次数最低的多项式称为α 在GF(p) 上的极小多项式(minimal polynomial ofα overGF(p) ).特别地,若
α∈GF(q) 为GF(q) 的一个本原元,则α 在GF(p) 上的极小多项式称为GF(p) 上的一个本原多项式(primitive polynomial forGF(q) overGF(p) ).
定义注1:对任意的
定义注2:设
中的元素具有相同的极小多项式. 设
定理:设
GF(q) 是一个含有q 个元素的有限域,GF(p) 是GF(q) 的一个含有p 个元素的子域. 设α∈GF(q) ,r 为满足αpr=α 的最小正整数. 则α 在GF(p) 上的极小多项式g(x) 是一个r 次不可约多项式,并且
B(α)={α,αp,αp2,…,αpr−1}
中的元素为g(x) 在GF(q) 上的所有不同的根,即
g(x)=(x−α)(x−αp)(x−αp2)⋯(x−αpr−1).
注:
在模
推论:设
更进一步,
注:设
例3 考虑二元域
GF(2) 上的不可约多项式p(α)=α3+α+1 ,构造有限域
GF(23)=GF(2)[α]/⟨p(α)⟩={0,1,α,α+1,α2,α2+1,α2+α,α2+α+1}.
容易验证,α,α2,α3,α4,α5,α6 都是GF(23) 的本原元.GF(2) 上的首项系数为1 的3 次本原多项式有两个,分别为
(i)α,α2,α4 在GF(2) 上的极小多项式
g(x)=(x+α)(x+α2)(x+α4)=x3+x+1
(ii)α3,α5,α6 在GF(2) 上的极小多项式
g(x)=x3+x2+1
有限域
定理:设
GF(q) 是一个含有q 个元素的有限域,GF(p) 是GF(q) 的一个含有p 个元素的子域,g(x) 是GF(p) 上的一个不可约多项式. 则g(x) 为GF(p) 上的本原多项式当且仅当g(x) 在GF(q) 上的根都是GF(q) 的本原元.
下面例子说明不可约多项式不一定是本原多项式.
例4 考虑二元域
GF(2) 上的不可约多项式p(x)=x4+x3+x2+x+1 ,构造有限域
GF(24)=GF(2)[x]/⟨p(x)⟩={a+bx+cx2+dx3 | a,b,c,d∈GF(2)}.
显然,x∈GF(24) . 由于x5=1 ,即x 的阶为5 ,因此,x 不是GF(24) 的本原元. 于是,p(x) 不是GF(2) 上的本原多项式. 另外,可以验证x+1 是GF(24) 的本原元.
2 Matlab 中的有限域计算函数
Matlab 中自带的有限域的计算是在
由 “1.1 有限域的构造” 的 “例2” 可知,我们只需先找到一个
然而,这样得到的有限域
Matlab 中
GF(2m) 的元素: 在 Matlab 中GF(2m):=GF(2)[D]/⟨p(D)⟩ ,其中p(D) 为一个GF(2) 上的m 次本原多项式.
GF(2m)={am−1Dm−1+am−2Dm−2+⋯+a1D+a0, | ai∈GF(2),0≤i≤m−1}
因此,每个GF(2m) 中的元素本质上是一个次数小于m 的多项式,每个元素和多项式之间有“1-1”对应关系. 例如,取m=3 和本原多项式p(D)=D3+D+1 ,则我们得到有限域GF(23) ,其中的元素和多项式之间的对应关系如下:
二进制 | ||
---|---|---|
0 | 000 | |
1 | 001 | |
2 | 010 | |
3 | 011 | |
4 | 100 | |
5 | 101 | |
6 | 110 | |
7 | 111 |
GF(2) 上的多项式由系数组成的二进制所对应的(十进制)数字来表示. 例如,多项式p(D)=D3+D+1 的系数组成的二进制为1011 ,因此,多项式p(D) 表示为数字11 .
2.1 定义有限域数组
在 Matlab 中,函数 gf 用来定义一个有限域数组,函数申明如下:
X_GF = GF(X,M,PRIM_POLY)
函数创建有限域
例如,生成有限域
>> GF8 = gf(0:7,3,13)
GF8 = GF(2^3) array. Primitive polynomial = D^3+D^2+1 (13 decimal)
Array elements =
0 1 2 3 4 5 6 7
如果不指定本原多项式,则 Matlab 将使用默认本原多项式. 例如
>> gf(0:7,3)
ans = GF(2^3) array. Primitive polynomial = D^3+D+1 (11 decimal)
Array elements =
0 1 2 3 4 5 6 7
在这里例子中,Matlab 使用了
如果不指定次数 M 和本原多项式 PRIM_POLY,则生成二元域
>> gf(0:1)
ans = GF(2) array.
Array elements =
0 1
生成的有限域中的数组可以参与运算(+、、.、.^、\等). 注意:参与运算的操作数必须来自同一个有限域,用于生成有限域的本原多项式也必须相同!
一个典型的例子是计算有限域的乘法表如下:
>> GF8 = gf(0:7,3)
GF8 = GF(2^3) array. Primitive polynomial = D^3+D+1 (11 decimal)
Array elements =
0 1 2 3 4 5 6 7
>> GF8'*GF8
ans = GF(2^3) array. Primitive polynomial = D^3+D+1 (11 decimal)
Array elements =
0 0 0 0 0 0 0 0
0 1 2 3 4 5 6 7
0 2 4 6 3 1 7 5
0 3 6 5 7 4 1 2
0 4 3 7 6 2 5 1
0 5 1 4 2 7 3 6
0 6 7 1 5 3 2 4
0 7 5 2 1 6 4 3
>> GF8 = gf(0:7,3,13)
GF8 = GF(2^3) array. Primitive polynomial = D^3+D^2+1 (13 decimal)
Array elements =
0 1 2 3 4 5 6 7
>> GF8'*GF8
Warning: Lookup tables not defined for this order 2^3 and
primitive polynomial 13. Arithmetic still works
correctly but multiplication, exponentiation, and
inversion of elements is faster with lookup tables.
Use gftable to create and save the lookup tables.
> In gf.gettables at 35
In gf.mtimes at 20
ans = GF(2^3) array. Primitive polynomial = D^3+D^2+1 (13 decimal)
Array elements =
0 0 0 0 0 0 0 0
0 1 2 3 4 5 6 7
0 2 4 6 5 7 1 3
0 3 6 5 1 2 7 4
0 4 5 1 7 3 2 6
0 5 7 2 3 6 4 1
0 6 1 7 2 4 3 5
0 7 3 4 6 1 5 2
在这里我们用两个不同的本原多项式构造有限域
注1:当我们计算
GF(2)[D]/⟨D3+D2+1⟩ 的乘法表时,Matlab 给产生一个警告 “Warning: Lookup tables not defined for this order 2^3 and primitive polynomial 13.” 从警告中我们可以看出,Matlab 中有限域的乘法是通过查表来完成的,这样可以显著地提高计算的速度. 我们可以通过命令 gftable 来创建并保存查找表格.
注2:用本原多项式D3+D+1 和D3+D2+1 生成两个不同的元素个数为8 的有限域,然而这两个有限域是同构的. 一般地,我们有如下有限域同构定理:定理: 任意两个元素个数相同的有限域一定同构.
与本原元多项式相关的函数
primpoly
函数 primpoly 用于计算
PR = PRIMPOLY(M, OPT, 'nodisplay')
其中 M 为本原多项式的次数,其取值为
OPT = 'min' 给出一个权值最小的本原多项式
OPT = 'max' 给出一个权值最大的本原多项式
OPT = 'all' 给出所有的本原多项式
OPT = L 给出所有权值为L的本原多项式
字符串 ‘nodisplay’ 用于关闭默认的本原多项式显示方式.
例如,输出
>> primpoly(3,'all')
Primitive polynomial(s) =
D^3+D^1+1
D^3+D^2+1
ans =
11
13
>> primpoly(3,'all','nodisplay')
ans =
11
13
isprimitive
函数 isprimitive 用来检查
CK = ISPRIMITIVE(A)
其中 A 为一个表示多项式的数字,并且表示的多项式的次数不能超过
例如,检查多项式
>> isprimitive(13)
ans =
1
>> isprimitive(15)
ans =
0
其它函数
见 Matlab 帮助.