四元数法
[cos(b), sin(b)*sin(a), -sin(b)*cos(a);
sin(c)*sin(b), cos(c)*cos(a)-sin(c)*cos(b)*sin(a), cos(c)*sin(a)+sin(c)*cos(b)*cos(a);
cos(c)*sin(b), -sin(c)*cos(a)-cos(c)*cos(b)*sin(a), -sin(c)*sin(a)+cos(c)*cos(b)*cos(a)];
一阶算法：
（2-31）
二阶算法：
（2-32）
三阶算法：
（2-33）
四阶算法：
（2-34）
求解四元数微分方程除了以上的增量算法，还可以利用龙格—库塔法或阿达姆斯方法等数值积分算法求解微分方程。但迭代一定时间后，四元数会不满足 ，为了让其满足此式，必须进行归一化处理。
初始值的确定：
正负号的确定
可取
其中的 由初始的姿态矩阵求得
四元数是哈密顿于1843年建立的数学概念，但只有在近四十年中才在刚体运动学中得到实际应用。四元数是由1个实数单位1和3个虚数单位 ， ， 组成的包含4个实元的超复数。若将 ， ， 视为基矢量，则可以把四元数分为标量和矢量两部分组成。其形式为：
（2-22）
且 ，其中 为标量， 为矢量。
引入四元数后，则方向余弦矩阵可由四元数表示为
（2-23）
由（2-12），我们可得
（2-24）
利用毕卡逼近法求解可得： （2-25）
令：
（2-26）
（2-25）可简写ຫໍສະໝຸດ ：（2-27）将 展开可得：
（2-28）
由于：
（2-29）
将式（2-29）代入式（2-28）整理可得：
（2-30）
在实际解算中，把 和 展为级数形式并取有限项，得四元数的各阶近似算法。
四元数法希腊哲学
四元数法希腊哲学1.为什么要用四元数可能四元数的由来大家都看过很多遍。
很久以前，一位老者坐在大桥边上，看着过往船只，突然灵光一闪，在桥边石碑上洋洋洒洒刻上几行大字，四元数诞生了！故事大家都爱听，那么为什么我们需要四元数？一种说法是解决向量乘法，我们知道向量之间乘法有内积和外积，但这两个运算均不完美，即不满足群的条件（当然四元数诞生的时候也还没有内积外积的说法）。
那向量之间是否存在这样一个非常完美的乘法，于是三维空间无法解决的问题就映射到四维空间。
这便是四元数诞生的契机。
那么问题又来了，既然四元数只是为了解决矩阵乘法，那为什么我们现在要用四元数进行旋转，甚至替代了欧拉角、轴角等形式？首先，四元数并不是生来为了解决三维旋转，而是它的性质非常有利于表达旋转信息（后面会详述），所以了解四元数的性质要先于了解四元数在旋转中的应用。
至于四元数替代欧拉角等形式，就需要牵扯到一些别的知识点，我先罗列一下四元数相比其他形式的优点：✧解决万向节死锁（Gimbal Lock）问题✧仅需存储4个浮点数，相比矩阵更加轻量✧四元数无论是求逆、串联等操作，相比矩阵更加高效所以综合考虑，现在主流游戏或动画引擎都会以缩放向量+旋转四元数+平移向量的形式进行存储角色的运动数据。
2.四元数是怎么想出来的平庸的教程会直接提出四元数的定义、运算规则等等，然后读者不知所云。
相反，较为系统的教程一般会从复数（Complex Number）进行引导，逐步提出四元数的定义，这样会让读者更容易理解，在脑中也更好形成画面。
那复数与四元数之间的关系、以及如何从复数这样一个概念扩展到四元数是我们需要理清的一个思路。
2.先说几个概念。
空间中的子空间：一般而言，空间（维度>2）都存在更低维的子空间，比如二维空间中一维子空间，也就是直线；三维空间中的一维子空间和二维子空间，也就是直线和面。
当超过三维的概念我们就很难去想象是什么样子，但四维空间一定会存在三维子空间或二维子空间。
四元数详解
四元数详解四元数是一种数学概念，它在多个领域都有广泛的应用。
在计算机图形学中，四元数用于表示旋转变换。
下面我将以人类的视角来介绍四元数的定义、性质和应用。
四元数是一种扩展了复数的数学结构。
它由一个实部和三个虚部组成，可以写成q = a + bi + cj + dk的形式，其中a、b、c、d分别是实数，i、j、k是虚数单位。
与复数一样，四元数也有加法和乘法运算。
我们来看四元数的定义。
四元数的实部a对应于实数部分，而虚部bi + cj + dk对应于虚数部分。
四元数的加法定义很简单，就是将实部和虚部分别相加。
而乘法则稍微复杂一些，需要使用四元数的乘法规则：i² = j² = k² = ijk = -1。
通过这个规则，我们可以计算出两个四元数的乘积。
接下来，我们来探讨一下四元数的性质。
首先，四元数的加法满足交换律和结合律。
然而，四元数的乘法不满足交换律，即ab ≠ ba。
此外，四元数的乘法满足结合律，但不满足分配律。
这些性质使得四元数的运算有一些独特的特点。
四元数在计算机图形学中有广泛的应用。
由于四元数可以用于表示旋转变换，因此在三维游戏和动画中经常被用到。
与传统的欧拉角相比，四元数具有很多优点，例如不存在万向锁问题和旋转插值更加平滑。
因此，使用四元数可以提高计算机图形学的效率和质量。
除了计算机图形学，四元数还在其他领域有着重要的应用。
例如，在航空航天领域，四元数可以用于表示飞行器的姿态和旋转控制。
在物理学中，四元数可以用于描述粒子的自旋。
此外，四元数还可以用于解决某些数学问题，例如解四次方程和计算曲线的弯曲度。
四元数是一种重要的数学概念，具有广泛的应用。
它在计算机图形学、航空航天和物理学等领域都发挥着重要作用。
通过深入理解四元数的定义、性质和应用，我们能够更好地应用它们解决实际问题，推动科学技术的发展。
四元数计算位移公式
四元数计算位移公式四元数是一种用来表示旋转的数学工具，它可以进行旋转的合成和插值运算。
但是，四元数也可以用来表示位移，尽管它不如矢量或矩阵表示法直观。
在本文中，我们将介绍如何使用四元数进行位移计算。
首先，让我们回顾一下四元数的基本定义。
四元数可以表示为q=(w,x,y,z)，其中w是实数部分，而x、y和z是虚数部分。
四元数的基本属性是它们可以进行加法和乘法运算。
在位移计算中，我们首先需要将位移向量转换为四元数形式。
假设我们有一个三维位移向量p = (px, py, pz)，我们可以将其转换为四元数形式q' = (0, px, py, pz)。
注意，w部分为0，以防止位移影响旋转。
接下来，假设我们有一个代表当前旋转的四元数r = (rw, rx, ry, rz)。
我们希望将位移向量与当前旋转相结合，并计算得到一个新的位移向量p' = (px', py', pz')。
要计算p'，我们需要使用四元数的乘法运算。
首先，我们需要将p'转换为四元数形式q = (0, px', py', pz')。
然后，我们可以计算p' = rq，即将当前旋转应用于位移向量。
这可以通过以下公式来实现：p' = (rw, rx, ry, rz) * (0, px, py, pz) * (rw, -rx, -ry, -rz)这个乘法过程可以通过以下公式进行展开：p' = ((rw * 0 - rx * px - ry * py - rz * pz) + (rw * px + rx * 0 + ry * pz - rz * py) + (rw * py - rx * pz + ry * 0 + rz * px) + (rw * pz + rx * py - ry * px + rz * 0)) * (rw, rx, ry, rz)根据上面的展开，我们可以计算出p'的各个分量px'、py'和pz'。
四元数
我们举个最简单的例子：把点 P(1, 0, 1)绕旋转轴 u = (0, 1, 0)旋转 90°，求旋转后的顶点坐标。首先将 P 扩充到四元数， 即 p = (P, 0)。而 q = (u*sin45°, cos45°)。求 p′=qpq−1 的值。建议大家一定要在纸上计算一边，这样才能加深印象， 连笔都懒得动的人还是不要往下看了。最后的结果 p` = ((1, 0, -1), 0)，即旋转后的顶点位置是(1, 0, -1)。
在 Unity 里，tranform 组件有一个变量名为 rotation，它的类型就是四元数。很多初学者会直接取 rotation 的 x、y、 z，认为它们分别对应了 Transform 面板里 R 的各个分量。当然很快我们就会发现这是完全不对的。实际上，四元数的 x、y、z 和 R 的那三个值从直观上来讲没什么关系，当然会存在一个表达式可以转换，在后面会讲。
那么，四元数又是什么呢？简单来说，四元数本质上是一种高阶复数（听不懂了吧。。。），是一个四维空间，相对于 复数的二维空间。我们高中的时候应该都学过复数，一个复数由实部和虚部组成，即 x = a + bi，i 是虚数单位，如果你 还记得的话应该知道 i^2 = -1。而四元数其实和我们学到的这种是类似的，不同的是，它的虚部包含了三个虚数单位， i、j、k，即一个四元数可以表示为 x = a + bi + cj + dk。那么，它和旋转为什么会有关系呢？
[csharp] view plain copy print?
1.
transform.rotation = Quaternion.AngleAxis(degrees, transform.right) * transform.rotation;
四元数简介
可表示为 d + ai + bj + ck，a、b、c、d是 实数。 • 四元数的其它表示方法： 1、q=[w,v] 其中v=（x,y,z）是矢量，w是标量。 2、q=( cosθ/2 , (x ,y,z )sinθ/2) 3、q=[d,a,b,c]T
• • • • •
纠正陀螺仪积分得出的姿态。 长期融合用到的传感器是加速度计和电子罗盘 假设正确的姿态四元数为Q,那么可以利用四元数旋转将参考坐标系和 体坐标系下的向量互相转换,将Eh和Eg转换到体坐标系下 BEh=Q×Eh×Q* BEg=Q×Eg×Q * 理论上Bg=BEg,Bh=BEh 即 { Q×Eh×Q* - Bh = 0} , { Q×Eg×Q* - Bg = 0} 这两个方程成立,联立这两个方程就可以解得姿态四元数Q。但是，由 于各种误差的存在,这个方程组只能找到最优解,找最优解的问题有许 多方法可以采用,如[梯度下降][高斯牛顿]
• q1 q2=
w2 x 2 y2 z2
四元数在姿态运算上的应用
• 一个单位化的四元数可以描述一个三维旋转的过程。例如：
p0以原点为旋转中心，旋转的轴是(α, β, γ) （ α^2 + β^2 + γ^2 = 1）， 转θ角的旋转，用四元数表示就是： p1= qp0q-1
四元数
齐朋冲
四元数的引入
• 复数的加、乘运算可以表示平面向量的合成、伸
缩和旋转变换，这些知识我们已经在中学课程中 学过了。 那么，很自然的问题就是，在三维，或更高维空 间中是否也有复数的类似物？其中的元素还可以 像复数域那样做加、减、乘、除运算，并满足通 常复数的那些运算律，包括加法和乘法的交换律 与结合律、乘法对加法的分配律等待？更进一步， 我们是否可以期望用这样的数来表示三维或更高 维空间中的伸缩和旋转，就像用复数表示平面向 量的伸缩旋转那样方便？
四元数小总结
四元数⼩总结四元数记法：⼀个四元数包含⼀个标量分量和⼀个3D向量分量。
记标量为w，记向量为v或分开的x,y,z。
如下：[w,v][w,(x,y,z)]四元数与复数：四元数扩展了复数系统，它使⽤三个虚部i,j,k。
它们的关系如下：i2=j2=k2=-1ij=k,ji=-kjk=i,kj=-iki=j,ik=-j⼀个四元数[w,(x,y,z)]定义了复数 w+xi+yj+zk。
四元数和轴-⾓对：四元数能被解释为⾓位移的轴-⾓对⽅式。
其公式为下：设向量n为旋转轴，θ为绕轴旋转的量。
q=[cos(θ/2) sin(θ/2)n]=[cos(θ/2) (sin(θ/2)n x sin(θ/2)n y sin(θ/2)n z)]负四元数：-q=[-w (-x -y -z)]=[-w -v]q和-q代表的实际⾓位移是相同的，很奇怪吧！如果我们将θ加上360度的倍数，不会改变q代表的⾓位移，但它使q的四个分量变负了。
因此，3D中的任意⾓位移都有两种不同的四元数表⽰⽅式，它们互相为负。
单位四元数：⼏何上存在2个单位四元数：[1,0]和[-1,0]。
它们的意义是：当旋转⾓为360度的整数倍时，⽅位并没有改变，并且旋转轴也是⽆关紧要的。
数学上只有⼀个单位四元数：[1,0]。
任意四元数q乘以单位四元数[1,0]仍为q。
四元数的模：公式如下：||q||=||[w (x y z)]||=sqrt(w2+x2+y2+z2) =||[w v]||=sqrt(w2+||v||2)⼏何意义：||q||=sqrt(cos(θ/2)2+sin(θ/2)2||n||2)若n为单位向量，则：||q||=1四元数共轭：q*=[w -v]=[w (-x -y -z)]四元数的逆：q-1=q*/||q||但我们只使⽤单位四元数，故q-1=q*⼏何解释：使向量v反向，则旋转⽅向也反向了。
因此q绕轴旋转θ⾓，⽽q*沿相反的⽅向旋转相同的⾓度。
四元数乘法（叉乘）：[w1 v1][w2 v2]=[w1w2-v1v2 w1v2+w2v1+v2×v1]四元数叉乘满⾜结合律但不满⾜交换律：(ab)c=a(bc)ab!=ba四元数乘积的模等于模的乘积：||q1q2||=||q1|| ||q2||四元数乘积的逆等于各个四元数的逆以相反的顺序相乘：(ab)-1=b-1a-1如何⽤四元数将3D点绕轴旋转：让我们“扩展”⼀个标准3D点(x,y,z)到四元数空间，通过定义四元数p=[0, (x,y,z)]即可。
旋转矩阵和四元素法
旋转矩阵和四元素法引言：旋转矩阵和四元素法是计算机图形学中常用的两种方法，用于描述和计算三维空间中的旋转变换。
本文将详细介绍旋转矩阵和四元素法的原理、应用以及优缺点。
一、旋转矩阵的原理和应用：1. 原理：旋转矩阵是一个3x3的矩阵，用于描述三维空间中的旋转变换。
旋转矩阵的每一列代表了旋转后的坐标轴在原始坐标系中的表示。
通过将一个向量与旋转矩阵相乘，可以实现对该向量的旋转变换。
2. 应用：旋转矩阵广泛应用于计算机图形学、计算机视觉等领域。
在三维建模中，使用旋转矩阵可以实现物体的旋转、变形和姿态控制。
在游戏开发中，旋转矩阵常用于计算相机的朝向和角度变化。
此外，旋转矩阵还可以用于计算两个坐标系之间的转换。
二、四元素法的原理和应用：1. 原理：四元素法，又称为四元数法，是一种用四个实数表示旋转的方法。
四元数由实部和虚部组成，虚部是一个三维向量。
通过将旋转变换表示为一个旋转轴和旋转角度，可以通过四元数的乘法来实现旋转变换。
2. 应用：四元素法在计算机图形学中被广泛应用于旋转插值和动画的计算。
通过插值计算两个旋转变换之间的中间状态，可以实现平滑的动画过渡效果。
此外，四元素法还可以用于防止万向锁现象的发生，提高旋转变换的稳定性和精确性。
三、旋转矩阵和四元素法的优缺点比较：1. 旋转矩阵的优点：（1）计算简单直观，易于理解和实现；（2）可以直接应用于三维坐标系的变换；（3）可以通过矩阵的乘法来实现多个旋转变换的复合。
2. 旋转矩阵的缺点：（1）存在数值误差累积的问题，当进行多次旋转变换时，可能导致结果不准确；（2）矩阵的运算比较耗时，特别是在计算资源有限的设备上；3. 四元素法的优点：（1）无数值误差累积问题，旋转变换精确度高；（2）计算速度较快，适用于实时计算和动画插值；（3）可以方便地进行旋转插值和平滑动画的计算。
4. 四元素法的缺点：（1）计算过程相对复杂，需要使用四元数的乘法和插值计算；（2）不直观，难以理解和调试；（3）在某些特定情况下，可能出现奇异性和计算不稳定性。
四元数 四旋翼飞行器姿态解算
四元数四旋翼飞行器姿态解算四元数是用来描述空间旋转的数学工具，在飞行器姿态解算中具有重要的应用。
四旋翼飞行器是一种采用四个电动马达驱动的多旋翼飞行器，通过调节四个马达的转速实现飞行器的姿态控制。
在四旋翼飞行器的飞行过程中，需要实时获取飞行器的姿态信息，以便进行飞行控制。
四元数作为一种有效的姿态描述方法，被广泛应用于四旋翼飞行器的姿态解算中。
四元数是一种具有四个元素的数学结构，通常表示为q = w + xi + yj + zk，其中w、x、y、z分别表示四元数的实部和三个虚部。
四元数可以表示为一个旋转矩阵，通过四元数乘法运算可以实现空间旋转的复合。
在四旋翼飞行器的姿态解算中，通常使用四元数来描述飞行器的姿态状态。
四旋翼飞行器的姿态解算涉及到四元数的插值、积分和旋转等计算。
在飞行器的姿态控制过程中，需要将传感器获取的姿态信息进行融合处理，得到飞行器的姿态状态。
四元数插值可以实现飞行器姿态信息的平滑过渡，提高飞行的稳定性和平顺性。
四元数积分可以实现对飞行器姿态状态的更新，保持飞行器的正确姿态。
四元数旋转可以实现飞行器的姿态控制，使飞行器按照指定的姿态进行飞行。
在四旋翼飞行器的姿态解算中，需要考虑传感器误差、系统延迟和控制精度等因素。
传感器误差会影响到飞行器的姿态感知精度，需要通过滤波算法和校准方法来降低误差影响。
系统延迟会导致飞行器姿态状态的延迟更新，需要通过合理的控制策略来补偿延迟效应。
控制精度是指飞行器姿态控制的准确性，需要通过优化控制算法来提高飞行器的稳定性和精确性。
总的来说，四元数是一种有效的姿态描述方法，被广泛应用于四旋翼飞行器的姿态解算中。
通过四元数插值、积分和旋转等计算，可以实现对飞行器姿态状态的准确解算和控制。
在实际的飞行应用中，需要综合考虑传感器误差、系统延迟和控制精度等因素，全面提高飞行器的姿态解算精度和控制性能。
四旋翼飞行器的姿态解算是飞行控制领域的重要课题，将带来对未来飞行器飞行性能的提升和发展。
用四元数法的捷联惯性导航姿态解算程序
;no dirg ;)')s/m( /yticoleV'(lebaly ;)')sm01( / emiT'(lebalx ;)'noitaluclaC yticoleV'(eltit ;)'Z :etanidrooC noitagivaN','Y :etanidrooC noitagivaN','X :etanidrooC noitagivaN'(dnegel )yticolev,mt(tolp erugif ;no dirg ;)')eerged( /legnA'(lebaly ;)')sm01( / emiT'(lebalx ;)'noitaluclaC erutseG'(eltit ;)'legnA waY','legnA lloR','legnA hctiP'(dnegel ;)'.-b',):,3(ELGNA,mt,'.g',):,2(ELGNA,mt,'-r',):,1(ELGNA,mt(tolp ;ip/081*ELGNA=ELGNA erugif dne
值比度速加/压电的向轴个三计度速加% ?=za_orez% ?=ya_orez% ?=xa_orez% 压电点零的向轴个三计度速加%
*********************************************************% % 值度速加的应相为换转据数将% 3*4201/cca_atad=cca_atad ;]z_cca;y_cca;x_cca[=cca_atad 压电拟模为换转字数度速加% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ;)'txt.zetarelecca'(daol=z_cca ;)'txt.yetarelecca'(daol= y_cca ;)'txt.xetarelecca'(daol= x_cca ;4201/3=etar_cca% **************据数的计度速加入读****************% ;)'txt.zoryg'(daol=z_oryg ;)'txt.yoryg'(daol=y_oryg ;)'txt.xoryg'(daol=x_oryg 据数的仪螺陀入读**********% 据数入读****************************% ;']g-,0,0[=G 2^s/m8.9 位单%;24108.9=g 量矢度速加的生产力重% ;lla raelc ;lla esolc
四元数,矢量运算规则,场论基础,并矢,算符,场量的Taylor展开,正交曲线坐标系,Delta函数
四元数quaternions复数对四则运算,代数运算,极限自封.四元数是复数的扩展.四元数有四个单元:k ,j ,i ,1.四元数定义dk cj bi a +++=α,其中R d ,c ,b ,a ∈ 另一四元数R 'd ,'c ,'b ,'a ,k 'd j 'c i 'b 'a ∈+++=β,则四元数加减法定义对应分量相加减;四元数乘法定义为)k 'd j 'c i 'b 'a )(dk cj bi a (++++++=αβk )'cb 'bc 'da 'ad (j )'bd 'db 'ca 'ac (i )'dc 'cd 'ba 'ab ()'dd 'cc 'bb 'aa (-+++-+++-+++---=四元数的单元间的运算规则: j ik ki ,i kj jk ,k ji ij ,1k j i 222=-==-==-=-===四元数加法适合结合律,交换律;,即)()(βγαγαβ=而一般βααβ≠.(βααβα=⇒∈R ) 对实数有效的运算规则对复数总有效,但对复数有效的运算规则对四元数不总有效,(如上述的乘法的交换律)！！！ 四元数的共轭: dk cj bi a :---=α,若dk cj bi a +++=α 性质:αβαβ=四元数的迹: R a 2:)(S ∈=+=ααα性质: )(S )(S )(S βαβα+=+四元数的模:R d c b a :)(N 2222∈+++==ααα性质: )(N )(N )(N βααβ⋅=,0)(N 0=⇔=αα证明: 0,oder ,00==⇒=βααβ00)(N 00)(N 00,und ,0=⇒⎭⎬⎫≠⇒≠=⇒=⇒≠=βααβαααβαααβ,即00,und ,0=⇒≠=βααβ,同理00,,0=⇒≠=αβαβund证明:若α是方程0)(N )(S x 2=+-αα的根,则α也是其根.因为,α是方程0)(N )(S x 2=+-αα的根0)()(2=++-⇒ααααα⇒=++-⇒0)()(2αααααα也是其根)四元数域内二次方程一般不止两个根,如最简单的方程1x 2-=就最少存在k ,j ,i ±±±6个根,实际上1x 2-=有无穷多个根,因为使1r q p 222=++成立的实数r ,q ,p 有无穷多个,而1)r q p ()rk qj pi (2222-=++-=++Halmiton 四元素体;第一个非交换体,1843 年 W.R.Hamilton 为建立三维复数空间,把复数x+iy 作为有序偶的实数,并定义规则,使i 在有明确意义: 4阶实方阵集H内方阵型如⎪⎪⎪⎪⎭⎫ ⎝⎛------a b c d b a d c c d a b d c b a ,令⎪⎪⎪⎪⎭⎫ ⎝⎛--=⎪⎪⎪⎪⎭⎫ ⎝⎛--=⎪⎪⎪⎪⎭⎫ ⎝⎛--=⎪⎪⎪⎪⎭⎫ ⎝⎛=0110011010011001011001101111K ;J ,E ,I ,则集H 内任意方阵可唯一表为dK cJ bE aI +++,即}R d ,c ,b ,a |dK cJ bE aI {H ∈+++=,H 对矩阵减法封闭;且I K J E -===222,;J KE ,E JK ,K EJ ===J EK ,E KJ ,K JE -=-=-=,矩阵乘法在H 内封闭,故H 对矩阵加,乘法构成环;H 的元素个数>1;I 是H 的单位元,又因I )d c b a ()dK cJ bE aI )(dK cJ bE aI (2222+++=---+++,且当0≠+++dK cJ bE aI 时,d ,c ,b ,a 不全为零,故02222>+++d c b a ,所以H 中非零元在H 内存在逆元,综上所述H 是非交换体,常称H 为四元数环,称H 内的元为四元数Quaterion : t+xi+yj+zk,其中t为数量部分/纯量部分,xi+yj+zk 为向量部分.四元数系构成了以实数域为系数域的有限维可除代数,是向量代数和向量分析基础. 矢量运算规则两矢的内积:)b ,a cos(|b ||a |b a ∧=⋅ R V ,V →两矢的外积: )b ,a sin(|b ||a ||b a |∧=⨯, b ,a )b a ( ⊥⨯ V )V ,V (→ 物理意义: b ,a 两矢内积是功; b ,a 两矢外积的模是以b ,a两矢的为边平行四边形的面积. 故内积可交换,外积可反交换外积和内积的关系:)b ,a (sin |b ||a |))b ,a (cos 1(|b ||a |)b a (|b ||a ||b a |2222222222∧∧=-=⋅-=⨯ 即)b ,a sin(|b ||a ||b a |∧=⨯推论 22b a b b a a )b a ()b a (;b a )b a ()b a (-=⋅-⋅=+⋅-⨯=+⨯-四元数和两重积间的联系:两四元数k a j a i a 321++=α,k b j b i b 321++=β;两矢量)a ,a ,a (a 321=,)b ,b ,b (b 321= 间关系βα↔↔b ,a 两矢内积和四元数间的关系:两量积)Re()Re()(21)(21b a αββαβαβααββα-==+=+=⋅ ,即两矢内积b ,a 对应于四元数βα的实部.两矢外积和四元数间的关系:矢量内积)Im()Re()(21b a αβαβαβαβαβ=-=-=⨯ ,即两矢外积b a ⨯对应于四元数αβ的非实部.两矢内积,外积和四元数间的关系:αβ=⋅-⨯b a b a三矢内积)c b (a c )b a (:]c ,b ,a [ ⨯⋅≡⋅⨯=,R V ,V ,V → 物理意义: c ,b ,a 三矢的内积是以c ,b ,a三矢为边的平行六面体的体积性质:b )c a (a )b c (c )a b (b )a c (a )c b (c )b a (⋅⨯-=⋅⨯-=⋅⨯-=⋅⨯=⋅⨯=⋅⨯推论:0]q c ,p b ,r a []p c ,r b ,q a []r c ,q b ,p a [=⨯⨯⨯+⨯⨯⨯+⨯⨯⨯2]c ,b ,a []a c ,c b ,b a [ =⨯⨯⨯ 三矢外积c )b a ()c a (b )c b (a⋅-⋅=⨯⨯V )V ,V ,V (→c)b a (b )c a (c c c )b a b a b a (b b b )c a c a c a (c )b a b a b a (b )c a c a c a (c )b a b a b a (b )c a c a c a (c )b a b a b a (b )c a c a c a ()c b c b (a )c b c b (a )c b c b (a )c b c b (a )c b c b (a )c b c b (a c b c b c b c b c b c b a a a )c b (a 321332211321332211333221133322112332211233221113322111332211233223113112211233233113312212122131132332321⋅-⋅=⎪⎪⎪⎭⎫ ⎝⎛++-⎪⎪⎪⎭⎫ ⎝⎛++=⎪⎪⎪⎭⎫⎝⎛++-++++-++++-++=⎪⎪⎪⎭⎫ ⎝⎛---------=⎪⎪⎪⎭⎫ ⎝⎛---⨯⎪⎪⎪⎭⎫ ⎝⎛=⨯⨯推论0)b a (c )a c (b )c b (a=⨯⨯+⨯⨯+⨯⨯四矢内积:)c b )(d a ()d b )(c a (db cb d a ca )d c ()b a (⋅⋅-⋅⋅=⋅⋅⋅⋅=⨯⋅⨯R )V ,V ,V ,V (→)c b )(d a ()d b )(c a (a )d )c b (c )d b ((a ))d c (b ()d c ()b a (三矢外积三矢内积 ⋅⋅-⋅⋅=⋅⋅-⋅=⋅⨯⨯=⨯⋅⨯四矢外积:a ]b ,d ,c [b ]a ,d ,c [d ]c ,b ,a [c ]d ,b ,a [)d c ()b a (-=-=⨯⨯⨯ V )V ,V ,V ,V (→ a ]b ,d ,c [b ]a ,d ,c [a )b )d c ((b )a )d c (()d c ()b a (;d ]c ,b ,a [c ]d ,b ,a [d )c )b a ((c )d )b a (()d c ()b a (三矢外积三矢外积 -=⋅⨯-⋅⨯=⨯⨯⨯-=⋅⨯-⋅⨯=⨯⨯⨯推论c ]c ,b ,a []d ,b ,a [b ]c ,b ,a []c ,d ,a [a ]c ,b ,a []c ,b ,d [d 0]c ,b ,a [ ++=→≠ )d a )(c b ()c a )(d b ()}d c (b {a ⨯⋅-⨯⋅=⨯⨯⨯流线 等X 面/线 通量 环流量 散度 旋度 方向导数 梯度为形象描述矢量场)z ,y ,x (f 定义)z ,y ,x (f 的流线f.为形象描述标量场)z ,y ,x (ϕ定义)z ,y ,x (ϕ的等X 面/线.S d 为开/闭有向曲面S 上一面元,矢量f 在面元S d 上的元通量S d f d f⋅=Φ,面积分得矢量场)z ,y ,x (f 在曲面S 上的通量(标)⎰⋅=ΦSf S d fl d 为开/闭有向曲线l 上一面元,矢量f 在线元l d 上的元环量l d f d f⋅=Θ,线积分得矢量场)z ,y ,x (f 在曲线l 上的环量(标)⎰⋅=Θlfl d f 矢量场)z ,y ,x (f的散度(标):描述有源场源/汇强度. 正/负/零散度对应于源/汇/无源无汇闭合曲面S 包围体积V ∆,0V →∆时f 在S上的通量与V ∆比的极限称为矢量场)z ,y ,x (f 的散度V /S d f lim f div S 0V ∆⋅=⎰→∆矢量场)z ,y ,x (f的旋度(矢):描述有旋场旋涡强度和旋涡法矢方向. 旋度的法向分量的模的大小顺比于涡旋场旋涡程度.闭合曲线l 包围有向曲面S ∆,0S |S |→∆=∆ 时f 在l 上的环量与S ∆的比的极限称为矢量场)z ,y ,x (f 的旋度f rot 沿S∆法向的分量S /l d f lim )f rot (l 0S n ∆⋅=⎰→∆ 等效于0S →∆时S )f rot (l d f ∆⋅=⋅⎰标量场)z ,y ,x (ϕ的梯度(矢):描述标量场各点空间变化率及方向.某场点的梯度的方向是标量场变化最快的方向,其模是标量场单位长度的变化率.场)z ,y ,x (ϕ沿l d 向改变ϕd ,称dl d ϕ为ϕ沿l d 向的方向导数,dl d ϕ等于ϕ的梯度的l d 向分量l d )grad (d dld )grad (l ⋅=⇔=ϕϕϕϕ积分变换公式Gauss 定律: ⎰⎰∂⋅=⋅∇VV dS f n dV f(f 的散度对体积V 体积分 ←转换→ f 对V 的包面的闭面积分)Stokes 定律: ⎰⎰∂⋅=⋅⨯∇SS l d f dS n )f ((f 对有向曲线S ∂的闭线积分 ←转换→ f 的旋度对以S∂为边的有向曲面S 的面积分)Green 恒等式:⎰⎰⎰⎰⎰∂∂∂∂∂=⋅∇=⋅∇=∇⋅∇=∇+∇⋅∇V V V V V 2dS ndS n )(S d )(dV )(dV )(ψφψφψφψφψφψφ (n ∂∂:外法线方向导数)Green 定理:⎰⎰⎰⎰⎰∂∂∂∂∂-∂∂=⋅∇-∇=⋅∇-∇=∇-∇⋅∇=∇-∇VV V V V 22dS )n n (dS n )ˆˆ(S d )ˆˆ(dV )ˆˆ(ˆdV )ˆˆ(φψψφφψψφφψψφφψψφφψψφ ⎰⎰⨯∇=⨯∂V V dV A A S d (⎰⎰⎰⎰⎰⎰⨯∇⋅=⨯⋅⇔⨯∇⋅=⨯⋅∇=⨯⋅=⨯⋅∂∂∂VV V V V V dV )A (C )A S d (C dV )A (C dV )C A ()C A (S d )A S d (C) ⎰⎰∇=∂VVdV S d ψψ (⎰⎰⎰⎰⎰⎰∇⋅=⋅⇔∇⋅=⋅∇=⋅=⋅∂∂∂VVVVVVdV C )S d (C dV C dV )C (S d C )S d (C ψψψψψψ)⎰⎰∂=∇⨯SSl d S d ψψ(⎰⎰⎰⎰⎰⎰∂∂⋅=⋅∇⨯⇔⋅=⋅⨯∇=⋅⨯∇=⋅∇⨯SSSSSSl d C C )S d (l d C S d )C (S d )C (C )S d (ψψψψψψ)并矢 及其运算标量是零阶张量,矢量是一阶张量,二阶三维张量借助于直角坐标转动矩阵定义:矢量i T 在坐标架转动满足变换关系i mi m T R T =,坐标转动矩阵miR 即二阶张量.二阶张量ij T 满足变换关系ij nj mi mn T R R T =.由两矢B A ,并列放置且之间无运算则构成并矢B A,含9个分量,记为j i B A ,由于i A 和j B 分别满足:i mi m A R A =,j nj m B R B =,故并矢B A满足j i nj mi B A R R B A = ,故并矢是二阶张量的一种形式,显然三个矢量的并矢具有三阶张量的变换关系. ()⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=≡><≡⊗≡z z z z z z y y y y y y x x x x x x z y x z y x v v g f g f g f g f g f g f g f g f g f g g g f f f g f g f g f g f || , 单位并矢(单位二阶张量)ij ij kk jj ii r δ=I =++=∇=I,性质:X X X =⋅I =I ⋅;(X 为矢量或算符); 2:∇=∇∇ I ; ϕϕ∇=⋅∇)(I ; )(:T Spur T T I ii ==; g f g f I ⋅=:;并矢-矢量点乘区分左右:右点乘p g f p g f )(⋅=⋅;左点乘)(f p g f p g⋅=⋅,这样,三矢外积可用并矢表示)()(p g g p f p g f -⋅=⨯⨯两二阶张量B A ,间的双点乘:ji ij B A B A = :(或))(()(:)(q f p g q p g f⋅⋅=)双点乘得到的标量是两矩阵积的迹.并矢的微分运算要注意是对那个张量进行的,一般需加括号.g f f g g f)()()(∇⋅+⋅∇=⋅∇g f g f g f)()()(∇⨯-⨯∇=⨯∇f r f r r f ⋅+⋅∇=⋅∇2)(22f r f r f +⋅∇=⋅∇)()(f r r f r r f r r f++⋅∇=⋅∇)()( ⎰⎰⎰∂∂⋅I =I ⋅=I ⋅∇VV V S d dS n dV⎰⎰∂I ⨯=I ⨯∇VV dS n dV⎰⎰∂=∇V V dS f n dV f⎰⎰∂⋅=⋅∇VVS d g f dV g f)()(根据以上矢量运算定理,可把Gauss 定理⎰⎰∂=∇VV n S d dV 和Stokes 定理⎰⎰∂=∇⨯SS l d n dS 的运算推广到对标量,矢量,张量的各种运算∇算符具有:矢量性和算符性.∇对矢量左/右/点/叉乘不可交换,矢量运算规则也适于∇,但需调整∇在结果中的位置,使等式左右量同型. f )g (g )f ()g f ( ⋅⨯∇-⋅⨯∇=⨯⋅∇ (第一式点和叉换位,取正;第二式交换第一式中的两矢量次序,取负)f )f ()f (2 ∇-⋅∇∇=⨯∇⨯∇ (按c )b a ()c a (b )c b (a⋅-⋅=⨯⨯写结果,再调整次序,使右端得矢量)ψϕϕψϕψ∇+∇=∇ )( f f )()f ( ⋅∇+⋅∇=⋅∇ϕϕϕf f )()f (⨯∇+⨯∇=⨯∇ϕϕϕ 2f 21f )f ()f f (f )f (f )f ( ∇-⋅∇=⋅∇-⋅∇=⨯⨯∇]g )f ()f (g []f )g ()g (f [)g f (∇⋅+⋅∇-∇⋅+⋅∇=⨯⨯∇ )f (g )g (f f )g (g )f ()g f (⨯∇⨯+⨯∇⨯+∇⋅+∇⋅=⋅∇3r =⋅∇ ; I r r =⊗∇≡∇; r e r r r ==∇;0e r r =⨯∇=⨯∇; r 2e r =⋅∇ 2r 3re r r r 1 -=-=∇ r e dr df r ˆdr df )r (f ˆ =∇=∇ ⎩⎨⎧=∞≠=-=-∇=-∇=∇)0(,)0(,0)(4ˆˆ1ˆ232r r r r e r r r r πδ Coulomb 定理的微分式:Green 函数|'r r |141)'r ,r (G ),r (4r e ˆ02r -==∇πεπδ标量场的梯度场无旋0)(≡∇⨯∇ϕ无旋场必可表为一标量场的梯度ϕ∇=⇒=⨯∇f 0f矢量场的旋度场无源0)f (≡⨯∇⋅∇ 涡旋场必可表为一矢量场的旋度A f 0f⨯∇=⇒=⋅∇a,0E ,k 为常矢a r )a ()r a ( =∇⋅=⋅∇ r r a r a ⋅=∇⋅ r a]e )e a (a [r 1e )a (r r r ⊥=⋅-=∇⋅ 533r r )r a (3r a r r )a ( ⋅-=∇⋅ 0a )r (r )a ()r a (=⋅⨯∇-⋅⨯∇=⨯⋅∇a 2r )a ()r (a )a (r a )r ()r a ()r a ()r a (r a =∇⋅-⋅∇+⋅∇-∇⋅=⨯⨯∇+⨯⨯∇=⨯⨯∇ 5333333r r )r a (3r a a )r r ()a (r r r r )a ()r r (a )r r a ( ⋅-=∇⋅+⨯∇⨯+∇⋅+⨯∇⨯=⋅∇ 0a )rr (r r )a ()r r a (333=⋅⨯∇-⋅⨯∇=⨯⋅∇ 35333333r a r )r a (3r r )a (a )r r ()a (r r r r )a ()r r (a )r r a ( -⋅=∇⋅-=∇⋅+⋅∇-∇⋅-⋅∇=⨯⨯∇33rr a )r r (a r 1a ⋅-=-⋅=∇⋅ r k i 0e E E ⋅=:E k i E ⋅=⋅∇;E k i )k i (e E )]r k i (e [E e E )e E (E r k i 0r k i 0r k i 0r k i 0⨯=⨯-=⋅∇⨯-=⨯∇-=⨯∇=⨯∇⋅⋅⋅⋅其中1:=⋅=r r r r e e I e e 故r r e e I r r I r r :'':''22==.Taylor 展开:...)(...'...''!)(...)(''!)('|'|,...,,+∂∂∂∂-++∂∂∂+∂∂-=-∑∑∑kj i k j i nk j i ji j i j i ii i rx x x x x x n r x x x x r x x r r r 111211112 其中k j i ,...,,取1,2,3; i x 代表直角坐标系的三个分量,注意:1 上式是对'r 展开; 2 对'r 的展开和对'r的展开相差一个负号. 曲线正交坐标系(Krummlinigen Koordinaten)三维空间里确定一点P 的位置需3个坐标321u ,u ,u .若P 点坐标在直角坐标系中表为)u ,u ,u (z ),u ,u ,u (y y ),u ,u ,u (x x 321321321===,则)z ,y ,x (u u ),z ,y ,x (u u ),z ,y ,x (u u 332211===,两坐标系等价.=i u 常数)3,2,1i (=的曲面是坐标面,他们的单位法向矢量为)3,2,1i (,e i =,其指向为iu 增加的方向.当过P 点的三坐标曲面两两垂直时,三坐标面的三交线也两两垂直,称此类坐标系为正交曲线坐标系.正交条件)j i (,0)u z )(u z ()u y)(u y ()u x )(u x (h ji j i j i 2ij≠=∂∂∂∂+∂∂∂∂+∂∂∂∂=.由i 31i i i 31i i i 31i i du u z dz ,du u ydy ,du u x dx ∑∑∑===∂∂=∂∂=∂∂=得2222)dz ()dy ()dx ()ds (++=232322222121j,i j i 2ij )du (h )du (h )du (h du du h ++==∑,其中2i 2i 2i 2ii 2i )u z ()u y()u x (h h ∂∂+∂∂+∂∂==,称)3,2,1i (,h i =为Lame 系数或度量因子.Delta 函数定义)a (f )x (f )a x (=-δ⎰性质1 偶函数)x ()x (δ=-δ2 采样性)a (f )x (f )a x (=-δ⎰3 函数下的面积⎩⎨⎧∉∈=-δ⎰])b ,c [a (,])b ,c [a (,dx )a x (bc 01; ⎩⎨⎧=∞≠=-δ)a x (,)a x (,)a x (04 缩放 )x (|a |)ax (δ=δ1证明)x (|a |)ax ()a (,|a |dz )z (|a |dz )z ()a (,a dz)z (dx )ax (|a ||a |a a a a δ=δ⇒⎪⎩⎪⎨⎧<δ=-δ>δ=δ⎰⎰⎰⎰∆∆-∆∆-∆∆-∆∆-100 5 若)(x f 为连续函数,且∆为包含a 电的任意长度区间,则a )x (g |)]x ('g /)x (f [dx ]a )x (g [)x (f =∆=-δ⎰证明dy a y g g dx dy x g dx dx x g dy a x g y a y g x a x g y )](['1)('1)(')()()(11+=→⎪⎭⎪⎬⎫=→=→-=+=→-=-- 若)x (f 为单值连续函数,且有N 个过零点N ,...,,i ,x i 321=,则a x g x g x f a g g a g f dy a y g g y a y g f dx a x g x f =--∆--∆==++=-⎰⎰)(1111|)('1)()](['1)]([)](['1)()]([])([)(δδ 6 复合函数∑-δ=δii i )x x (|x /)x (df |)]x (f [1 证明: )x (f y =单值连续,则)x (f 在每个过零点的邻域内可逆;且)x (g 为任意品优(gutartig)函数,则dy )]y (f ['f 'x dy )x ('f 'x dy )]'y (f ['x )y (f x )x (f y 11111---=→=→=→=→=⎰∑∑∑∑⎰⎰+∞∞---∆+∆---+∞∞--δ===δ=δii i ii i iix x dx )x x (|)x ('f |)x (g |)x ('f |)x (g )](f ['f )](f [g dy )]y (f ['f ]y [)]y (f [g dx )]x (f [)x (g i i1101011111 被积函数须相等,再由)x (g 的任意性,得∑-δ=δii i )x x (|x /)x (df |)]x (f [1 6 导数x)a x (:)a x ('∂-δ∂=-δ则)a ('f dx )a x (')x (f -=-δ⎰∆证明)a ('f dx )x ('f )a x ()x (df )a x (|)]a x ()x (f [)a x (d )x (f dx )a x (')x (f =-δ-=-δ--δ=-δ=-δ⎰⎰⎰⎰∆∆=+∆-∆∆∆ 07)a (f )(dx )a x ()x (f )n (n )n (1-=-δ⎰∆证明)a (f )(dx )x (f )a x ()(dx )x (f )a x ()(...dx )x (''f )a x ()x ('df )a x (|)]a x ()x ('f [)a x (d )x ('f dx )x ('f )a x ()x (df )a x (|)]a x ()x (f [)a x (d )x (f dx )a x ()x (f )n (n )n (n)m ()m n (m)n ()n ()n ()n ()n ()n ()n ()n ()n (1112202211011-=-δ-=-δ-==-δ=-δ+-δ-=-δ-=-δ-=-δ--δ=-δ=-δ⎰⎰⎰⎰⎰⎰⎰⎰⎰∆∆-∆-∆-=+∆-∆-∆-∆-∆-=+∆-∆-∆-∆⎰∆=-=-δa)x (g |]})x ('g )x (f [dx d )x ('g {dx ]a )x (g [')x (f 1 三维δ函数⎩⎨⎧=≠=-δ-δ-δ=-δ)r r (,)r r (,)z z ()y y ()x x ()r r (00000010曲线系下的三维δ函数3213020103020100h h h )u u ()u u ()u u (|)u ,x (J |)u u ()u u ()u u ()r r (i i -δ-δ-δ=-δ-δ-δ=-δ ,(其中)u ,x (J i i 为Jacobi 行列式)柱坐标系下)z z ()()()r r (00001-δϕ-ϕδρ-ρδρ=-δ球坐标系下)()()r r (sin r )r r (000201ϕ-ϕδθ-θδ-δθ=-δ注意:n 维δ函数的量纲为n m -,即n -米δ函数的逼近钟形曲线: 2201xa a lim )x (a +π=δ→ Gauss 曲线;)x n exp(n lim )x (n 220π-=δ→sinc 函数: )kx (c sin k lim x kx sin lim )x (k k π=π=δ∞→∞→1sinc 函数平方: )kx (c sin k lim kx kx sin lim )x (k k 2221π=π=δ∞→∞→ 复指函数:⎰⎰+∞∞-+∞∞-π=±π=δdk )kx cos(dk )ikx exp()x (121盒子函数:∑+∞-∞==δn )L /inx exp(L )x (21*********************.cn。