四元数

                     

贡献者: JierPeter; addis; 叶月2_

预备知识 复数,矢量叉乘,域

  

未完成:补充四元数与 S0(4),SO(3) 关系。新增覆叠映射的词条

1. 四元数的历史

   复数可以定义为具有特定乘法有序实数对,也可以看成是在实数域上添加一个新的元素 i 所得到的扩域。

   复数应用极为广泛,因为它具有非常良好的代数性质和分析性质,同时可以用来表示二维矢量空间,极大地方便了许多几何运算。John Derbyshire 在他的科普作品Unknown Quantity: Real and Imaginary History of Algebra中写道,“如果从一维的实数到二维的复数的变更会给我们带来这样巨大的力量和见识,那么为什么要止步呢?难道就没有正等待着被发现的其他种类的数,比如说超复数,它们本身就是三维的吗?还有,难道那些数就不能让我们进一步认识数学吗?”

   19 世纪初的爱尔兰数学家 William Hamilton 就曾试图寻找一种具有特定乘法的三元数组,期望能和三维空间的坐标一一对应。然而他努力多年却未能找到可以满足结合律、分配律等的三元数乘法。事实上,这样的三元数乘法是不存在的,但 Hamilton 当时不知道。

   1843 年 10 月 16 日,Hamilton 在和妻子一同前往爱尔兰皇家学会议会的路上,突然灵光一现,意识到了问题的本质。他 “很不理智” 地拿出刀来,在布鲁姆桥的一块石头上刻下了以下公式:

(1)i2=j2=k2=ijk=1 .
这一公式包含了一个四元数的良好的乘法定义。

2. 四元数的定义和基本运算

   现代四元数的定义是把每一个四元数看成一个实数和一个三维实矢量的组合:(a0,a)。在这个语境下,ijk 被看成是一个标准正交基中的矢量,为了避免将来的混淆,把这三个基矢量重新表示成 1,i^,j^k^ 后,四元数的运算可以定义如下:

定义 1 四元数的加法和乘法

   记四元数集合为 H={(a0,a)|a0R,aR3}(a0,a),(b0,b)H,加法定义为

(2)(a0,a)+(b0,b)=(a0+b0,  a+b) .
乘法定义为(乘号可以省略或者用其他符号)
(3)(a0,a)×(b0,b)=(a0b0ab,  a0b+b0a+a×b) .

   注意,虽然四元数的加法显然满足交换律,但由于一般来说矢量的叉积不可交换,因此四元数乘法一般不具有交换律。只有 a 平行于 b(包括其中之一为 0)的时候,才会有 a×b=0=b×a。四元数的乘法满足分配律习题 1 )。

   四元数 A=(a0,a) 的实数部分 a0 被称为标部,矢量部分 a 称为矢部。为了方便,有时候也把这个四元数记为 a0+a,类比复数的表示法 (a,b)=a+bi。当标部为 0 时,可以直接将四元数记为矢量 a,同理,当矢部为零时也可以记为标量 a0。这样,a0+a 就可以看作两个四元数相加。和复数同理,我们认为矢部为 0 的四元数就是实数。

例 1 

   两个标量相乘也可以看作是四元数相乘,易证

(4)ab=(ab,0)=(a,0)(b,0)=(b,0)(a,0) .
类似地,矢量的数乘(子节 2 λv 也可以认为是四元数相乘,易证
(5)λv=(0,λv)=(λ,0)(0,v)=(0,v)(λ,0) ,
标量与四元数相乘也类似
(6)λ(a0,a)=(λ,0)(a0,a)=(a0,a)(λ,0)=(λa0,λa) ,
而两个矢量(的四元数)相乘却是以前没有见过的,要区别于点乘和叉乘
(7)ab=(0,a)(0,b)=(ab, a×b) .
另外若令 a 的模长为 a,单位矢量 a^=a/a,使用四元数乘法有
(8)an=(0,a)n={(1)n/2an(n=0,2,4,) (1)(n1)/2ana^(n=1,3,5,) .
这里 an 的记号可能会产生歧义(区分于例如 a2=aa),但一般不难通过语境判断。也就是说在四元数乘法下,矢量的偶数次方是标量,奇数次方是矢量。

3. 四元数的性质

引理 1 与 Hamilton 定义的相容性

   给定四元数 A=(a,a)B=(x,x),如果将 ax 分别表示为 bi+cj+dkyi+zj+wk,其中 i2=j2=k2=ijk=1,那么

(9)(a+bi+cj+dk)×(x+yi+zj+wk)=(axbyczdw)+(ay+bx+cwdz)i+(az+cx+dybw)j+(aw+dx+bycz)k. 

   证明

   首先计算出 ax=by+cz+dwa×x=(cwdz)i+(dybw)j+(bzcy)k,然后可得

(10)A×B=(a,a)×(xx)=(axax,ax+xa+a×x)=(ax(by+cz+dw))+(ayi+azj+awk)+(xbi+xcj+xdk)+a×x=(cwdz)i+(dybw)j+(bzcy)k=(axbyczdw)+(ay+xb+cwdz)i+(az+xc+dybw)j+(aw+xd+bzcy)k. 

   证毕

定理 1 四元数结合性

   给定四元数 A=a0+a1i+a2j+a3kB=b0+b1i+b2j+b3kC=c0+c1i+c2j+c3k,则有

(11)A(B+C)=AB+AC, 
以及
(12)A(BC)=(AB)C. 

   证明

(13)A(B+C)=(a0+a1i+a2j+a3k)×((b0+c0)+(b1+c1)i+(b2+c2)j+(b3+c3)k)=(a0(b0+c0)a1(b1+c1)a2(b2+c2)a3(b3+c3))+(a0(b1+c1)+a1(b0+c0)+a2(b3+c3)a3(b2+c2))i+(a0(b2+c2)+a2(b0+c0)+a3(b1+c1)a1(b3+c3))j+(a0(b3+c3)+a3(b0+c0)+a1(b2+c2)a2(b1+c1))k=(a0b0a1b1a2b2a3b3)+(a0b1+a1b0+a2b3a3b2)i+(a0b2+a2b0+a3b1a1b3)j+(a0b3+a3b0+a1b2a2b1)k+(a0c0a1c1a2c2a3c3)+(a0c1+a1c0+a2c3a3c2)i+(a0c2+a2c0+a3c1a1c3)j+(a0c3+a3c0+a1c2a2c1)k=AB+AC. 
由此得证分配性。

   由于每个四元数都可以写成 1,i,j,k 四个元素分别乘以一个实数后相加,且四元数乘法具有分配性,故要证明结合性,只需要证明这四个元素的乘法总满足结合性即可;又因为 1 是乘法单位元,与任何其它两个四元数做乘法的时候必定满足结合性,故我们只需要证明 i,j,k 之间有结合性即可。最后,因为 ij=kjk=iki=j 的 “3-轮换对称性1,我们只需要证明以下 9 个结合性:

(14){i(jj)=i×(1)=ki=(ij)j;j(ij)=jk=kj=(ji)j;j(ji)=j×(k)=i=(jj)i;i(ij)=ik=j=(ii)j;i(ji)=i×(k)=ki=(ij)i;j(ii)=j=ki=(ji)i;i(jk)=ii=1=kk=(ij)k;i(kj)=i×(i)=1=jj=(ik)j;i(ii)=i=(ii)i. 
即可得证四元数乘法的结合性2

   证毕

   如果我们把四元数的矢部限制在一个方向上,那么得到的运算结构实际上就是复数。你可以通过比较四元数乘法规则和复数乘法规则来看出这一点——注意在此限制下任何矢部的向量积都是 0

   除了把四元数看成一个标量加一个三维矢量以外,我们也可以把四元数看成一个四维矢量,就像复数既可以看成一个实数配上一个一维矢量,也可以看成一个二维矢量一样。这种情况下,四元数实际上就是实数域上的一个四维代数3

   两个复数相乘,就是把它们的模长相乘,辐角相加(式 13 ),类似地两个四元数相乘的模长也等于各自的模长相乘(证明留做习题)

(15)|AB|=|A||B| ,
与复数类似,任何四元数都可以表示为模长乘以一个模长等于 1 的四元数:A=|A|q。其中 q 类似于复数的 eiθeiθ 可以代表一个角度为 θ 的平面旋转,而在 “四元数与旋转矩阵” 中会看到,任意模长为 1 的四元数 q 都可以表示一个唯一的三维旋转。

定义 2 共轭、模方和模

  

  • 令四元数为 A=a0+a,记 A~=a0a,称为 A共轭(conjugate)
  • A=a02+|a|2,称为 A模方
  • |A|=A,称为 A

习题 1 

  

  • 证明四元数满足乘法分配律
  • 证明四元数满足乘法结合律
  • 证明 ||A||=||A~||=A~A=AA~(类比复数)
  • 证明 |AB|=|A||B|,即 AB=AB(类比复数)
  • 若将单位矢量 e^ 看作四元数(标部为零),证明 e^2=1, e^3=e^

   和复数类似,定义非零四元数 A 的逆 A1 满足

(16)AA1=A1A=1(A0) .
易得每个非零四元数都有唯一的逆
(17)A1=A~|A|2(A0) .
代入式 16 即可证明。证明逆元的唯一性:假设 B1,B2 都是 A 的逆元,那么 AB1=AB2,等式两边左乘 B1,得 B1=B2。证毕。

   另外易证

(18)(AB)1=B1A1 .

   由于四元数的乘法不可交换,我们不定义除法,因为除法可能有两种不同的定义 AB1B1A。我们直接使用这两种符号即可。

   根据以上讨论不难发现,所有的非零四元数及其乘法运算构成一个定义 3 )。

4. 四元数的指数表示

   记 e 是自然对数的底数,那么对于任意实数 x,有 ex=n=0xn/n!。如果把实数 x 替换为其它能够进行乘方的对象,比如说矩阵、复数等,那么还可以把指数的概念拓展到更多的对象上,构造出复数的指数、矩阵的指数等。

定义 3 四元数的指数表示

   令四元数为 A=a0+a,则使用四元数的乘法(式 3 ),定义(类比复数的指数函数定义

(19)eA=n=0Ann! .

   可以证明(留做习题)

(20)eA=ea0ea=ea0n=0ann!=ea0(cosa+a^sina) .
其中 a^ 是和 a 同方向的单位矢量,aa 作为四元数的模或者作为三维实矢量的模。

习题 2 与复数的类比

   比较复数和四元数的指数表示,说明为什么可以将复数看成矢部限制在一个方向上的四元数。

   由于四元数的加法满足交换律,而乘法一般不满足,因此一般情况下四元数并不满足 eA+B=eAeB。因此,四元数的指数表示法在计算中其实用得不多。式 19 中的第一个等号成立是因为 a0a=aa0 满足交换律。

5. 四元数的拓展

   我们以上定义的四元数,可以看成四个基矢量 1i^j^k^ 的线性组合的集合,而这些线性组合的系数取自实数域 R。如果系数还可以取自复数域 C,定义出来的数被称为双四元数,或复四元数。一个双四元数可以看成一个复数和一个复数域上的三维矢量的组合。

   对于一个双四元数 A=a0+a1i^+a2j^+a3k^,共轭的概念就被细分了:

定义 4 双四元数的共轭

  • 四元共轭:A~=a0a1i^a2j^a3k^
  • 反共轭:AC=a0+a1i^+a2j^+a3k^
  • 转置:AT=a0+a1i^a2j^+a3k^
  • 复共轭:A=a0a1i^+a2j^a3k^
  • 厄密(hermitian)共轭:A=a0a1i^a2j^a3k^

   双四元数的模方是用四元共轭定义的:||A||=AA~=A~A=n=03an2。但是和四元数不同的是,双四元数的模方并不总是非负实数,而是有可能取任何复数。

习题 3 双四元数共轭模方的关系

   证明:||A||=||A~||=||AT||=||AC||=||A||=||A||

习题 4 双四元数共轭的乘法规律

   设 AB 是两个双四元数,证明:

(21){ACBC=(AB)CAB=(AB)A~B~=BA~ATBT=(BA)TAB=(BA) .

   和普通四元数不同的是,双四元数中有可能存在非零但模为 0 的数,比如下面例 2 所举例的 uv

   类似地,让四元数的系数取值为四元数,所得到的数被称为二级四元数。一般地,让四元数的系数取值为 n级四元数,则能得到 n+1级四元数。但是由于四元数不再是一个域,也没有交换律,用其做系数的效果远远不如复数系数,应用范围相对有限,在此不作详细介绍。

6. 四元数的矩阵表示

四元数与 SU(2)

   如果将四元数的四个基 1,i^,j^k^ 分别对应到以下四个矩阵

(22)1=(1001) ,i^=(i00i) ,j^=(0110) ,k^=(0ii0) .
并设该映射是线性的,那么一个任意的四元数或双四元数 A=a0+a 可以用矩阵表示为
(23)f(A)=(a0+ia1a2+ia3a2+ia3a0ia1) .

   这是数乘特殊酉矩阵的一般形式,满足 f(A)=f(A)1,且行列式为四元数的模方。显然,若 A 是单位四元数,则 detf(A)=1,对应 SU(2) 群中唯一一个矩阵。也就是说单位四元数与 SU(2) 群之间有一一对应。

   利用四元数的性质可以验证,单位四元数构成群。式 22 的映射是保四元数乘法的,即 f(i^j^)=f(i^)f(j^),其他同理。因此在群意义上,单位四元数与 SU(2) 同构:QSU(2)

例 2 矩阵表示的例子

  

(24)1+(10i)=1+1i^+ik^=(1+i111i) ,
(25)u=12(1ii^)=(1000) ,
(26)v=12(j^+ik^)=(0010) .

   在矩阵表示中,A 的转置共轭就是矩阵 A 的转置,复共轭就是 A 的各矩阵元取复共轭,厄密共轭就是 A 的厄密共轭矩阵(矩阵转置后再取各矩阵元的复共轭),而四元数的模方就是其矩阵表示的行列式:||A||=detA。因此,模方不为零的(双)四元数对应的矩阵是可逆的,其逆矩阵所表示的就是该四元数的逆。

单位四元数与 SO(3)

   设单位四元数 q=(cosθ2,n^sinθ2),v 为任意向量,则 qvq1 表示 vn^ 旋转 θ4。易证,qq 都代表同一个旋转,因此这是一个二对一的覆叠映射。当我们讨论单位四元数同态于 SO(3) 时,qvq1 所代表的旋转过程实际上是把四元数的矢量部分 ImH={a1i^+a2j^+a3k^} 同构为一般意义上的矢量空间 R3。 又因为 (q1q2)v(q1q2)1=q1(q2vq21)q11,因此单位四元数与 SO(3) 群同态。

推论 1 

(27)SU(2)/Z2SO(3) .

   证明:

   设所有单位四元数构成的集合为 Qf:QSU(2) 代表式 17 的同构映射。对于任意 qQg:QSO(3) 表示 q 对应的三维旋转,由前文讨论可知,这是一个同态映射。

   若 qvq1=v,即 q 对应 SO(3) 中的单位元 I,则 θ=0,2π。代入单位四元数的表示可得:q=1,1。因此 kerg={1,1}

   又因为 QSU(2),易证 f(kerg) 也是 SU(2) 的正规子群。所以 Q/{1,1}SU(2)/Z2SO(3),得证。

单位四元数与 SO(4)

   因为四元数乘法是 R 线性的,因此对于任意 a,bHab 可视作对 a 或者对 b 的线性变换。设 a=(a0,a)=(a0,a1,a2,a3),b=(b0,b)=(b0,b1,b2,b3),这两种线性变换可以分别对应

(28)ab=Lab=(a0a1a2a3a1a0a3a2a2a3a0a1a3a2a1a0)(b0b1b2b3)=Rba=(b0b1b2b3b1b0b3b2b2b3b0b1b3b2b1b0)(a0a1a2a3) .
显然,无论 La 还是 Rb,不同列和不同行都是相互正交的,而且行列式和列矢量平方都对应为四元数的模方。所以映射是 SO(4) 的群元当且仅当对应的四元数是单位四元数。

   利用 Lab=LaLb,Rab=RbRa,我们可以诱导 Q×QSO(4) 的同态。

   为满足群直积运算 (q1,q2)(q3,q4)=(q1q2,q3q4) 和运算保同态,旋转 u 对应的四元数乘法为

(29)Lq1Rq21u=q1uq21 .

定理 2 

   存在 2 对 1 的覆叠同态 ϕ:SU(2)×SU(2)SO(4),核元素为 Z25

   由于 SU(2)Q,且式 29 满足 ϕ:Q×QSO(4),我们只需证明核元素只有 2 个即可。

   设 (p,q)ker(Q×Q),则对于任意 uH 都有 puq1=u,对应单位矩阵。代入特例 u=1 得:p=q。结合前面的讨论易证,pup1 代表的是固定 u 的实轴部分,只对 ImH 部分进行旋转的线性变换。因此,若 pup1=u,则 p=1,1,所以 ker(Q×Q)={(1,1),(1,1)}。证毕。


1. ^ 即把 i,j,k 放在等边三角形的三个顶点上,通过逆时针旋转三角形 060120 来改变这几个元素的名字,其乘法结构不变。比如,如果 i,j,k 依次按逆时针排列,则三角形逆时针旋转 60 后,i 变为 kk 变为 j,而 j 变为 i,此时 ij=k 就变成了 ki=j,变换结果依然是成立的。
2. ^ 比如说,要证明 k(kj)=(kk)j,只需要利用轮换对称性,把 i 变为 kj 变为 ik 变为 j,即可将问题化为证明 j(ji)=(jj)i 即可。
3. ^ 从代数角度来说,Hamilton 爵士当年苦苦思索多年的问题就是 “是否存在实数域上的三维代数”,今天我们已经知道,实数域上只存在一维、二维、四维和八维的 “具有积性的可除代数”,分别对应实数、复数、四元数和八元数,其中四元数已经不具有交换性,而八元数连结合性都丢了。这里的 “积性” 是指 |a||b|=|ab|,而 “可除性” 就是 “乘法单位元存在性”。
4. ^ 证明见四元数与旋转矩阵一节。
5. ^ 可以证明,本节涉及的同态关系都是李群同态。


致读者: 小时百科一直以来坚持所有内容免费无广告,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此,我们请求广大读者热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 20 元,我们一周就能脱离亏损, 并在接下来的一年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。

                     

友情链接: 超理论坛 | ©小时科技 保留一切权利