n维欧拉公式

欧拉公式在n维的推广

​ 我们知道 ez=er(cosθ+sinθ)e^{z} = e^r(cos \theta + sin\theta) 为最熟悉的欧拉公式形式,而复数的表示可以用二元的矩阵取代,因此我们可以将复数形式的欧拉公式改写为矩阵形式,因此很自然的一个想法是,这个公式是否能推广到 nn 维欧式空间

欧拉公式的意义

​ 以一个线性微分方程组为例:

{ddtx1(t)=a11x1(t)+a12x2(t)++a1nxn(t)ddtx2(t)=a21x1(t)+a22x2(t)++a2nxn(t)ddtxn(t)=an1x1(t)+an2x2(t)++annxn(t)\left\{\begin{matrix} \frac{d}{dt}x_{1}(t) = a_{11}x_{1}(t) + a_{12}x_{2}(t) + {\cdots} + a_{1n}x_{n}(t) \\ \frac{d}{dt}x_{2}(t) = a_{21}x_{1}(t) + a_{22}x_{2}(t) + {\cdots} + a_{2n}x_{n}(t) \\ {\cdots} \\ \frac{d}{dt}x_{n}(t) = a_{n1}x_{1}(t) + a_{n2}x_{2}(t) + {\cdots} + a_{nn}x_{n}(t) \end{matrix} \right.

可以把它等价地写为 ddtx(t)=Ax(t)\frac{d}{dt}\boldsymbol x(t) = A\boldsymbol x(t),它的解就像普通微分方程一样:x(t)=eAtx(0)\boldsymbol x(t) = e^{\boldsymbol At}\boldsymbol x(0),它的几何意义是一个 nn 维向量在空间内的运动情况,我们把 x\boldsymbol x 视为运动轨迹,则 A\boldsymbol A 则视为对空间内的位置到对应的速度的变换:

A的含义

该方程的解 eAte^{\boldsymbol At} 视为让一个初始值沿着这个场运动 tt 时间(它仍然是一个矩阵,还是对空间的一个变换):

e^T

​ 如果能写出沿某个平面的沿着某个东西旋转的速度场,那么我们是否能利用这个公式来拓展我们的欧拉公式,让他描绘出旋转的过程?

​ 令 R=[0,11, 0]R = \left[ \begin{array}{c} ​ 0, -1\\ ​ 1,\ 0\\ \end{array} \right],则 R\boldsymbol R 代表在这个空间逆时针旋转 90°90\degree 的速度场,eRte^{\boldsymbol Rt} 就是旋转矩阵:

旋转

那么按理来说 eRte^{\boldsymbol Rt} 具有周期性,我们用数学的方式论证一下:

​ 由于 RR 的特殊性:

R4=[0110]4=[1001]2=I\boldsymbol{R}^4=\left[ \begin{matrix} 0& -1\\ 1& 0\\ \end{matrix} \right] ^4=\left[ \begin{matrix} -1& 0\\ 0& -1\\ \end{matrix} \right] ^2=\boldsymbol{I}

我们可以按照矩阵指数函数的定义直接展开:

eRt=n=0Rnn!tn=eRt=n=0Rnn!tn=n0(mod4)1n![tn00tn]+n1(mod4)1n![0tntn0]+n2(mod4)1n![tn00tn]+n3(mod4)1n![0tntn0]=[1t22+t44!t+t33!t55!+tt33!+t55!1t22+t44]=[costsintsintcost]e^{\boldsymbol{Rt}}=\sum_{n=0}^{\infty}{\frac{\boldsymbol{R}^n}{n!}t^n} =e^{\boldsymbol{Rt}}=\sum_{n=0}^{\infty}{\frac{\boldsymbol{R}^n}{n!}t^n} \\ =\sum_{n\equiv 0\left( mod4 \right)}{\frac{1}{n!}\left[ \begin{matrix} t^n& 0\\ 0& t^n\\ \end{matrix} \right]}+\sum_{n\equiv 1\left( mod4 \right)}{\frac{1}{n!}\left[ \begin{matrix} 0& -t^n\\ t^n& 0\\ \end{matrix} \right]}+\sum_{n\equiv 2\left( mod4 \right)}{\frac{1}{n!}\left[ \begin{matrix} -t^n& 0\\ 0& -t^n\\ \end{matrix} \right]}+\sum_{n\equiv 3\left( mod4 \right)}{\frac{1}{n!}\left[ \begin{matrix} 0& t^n\\ -t^n& 0\\ \end{matrix} \right]} \\ =\left[ \begin{matrix} 1-\frac{t^2}{2}+\frac{t^4}{4!}-\cdots& -t+\frac{t^3}{3!}-\frac{t^5}{5!}+\cdots\\ t-\frac{t^3}{3!}+\frac{t^5}{5!}-\cdots& 1-\frac{t^2}{2}+\frac{t^4}{4}-\cdots\\ \end{matrix} \right] =\left[ \begin{matrix} \cos t& -\sin t\\ \sin t& \cos t\\ \end{matrix} \right]

​ 这不就是线性代数里面熟悉的旋转矩阵吗(其实这也是二维的欧拉公式),有了这个理解,我们可以继续理解我们的高维欧拉公式

高维欧拉公式

​ 很遗憾的是,高维的欧拉公式只能规定在某一个平面里面旋转 nn 维向量,由于复平面的特殊性,我们在复平面内往往意识不到这一点

​ 我们以三维空间为例,设 u v\boldsymbol{\vec{u}\ \vec{v}} 为相互正交的单位向量(这样可以先不考虑它的伸缩性),考虑 w=u×v    r=w×r\boldsymbol{\vec{w}}=\boldsymbol{\vec{u} \times \vec{v}\ \ \ \ \vec{r}'=\vec{w} \times \vec{r}}(理解它的几何含义),这样就能把 r\boldsymbol {\vec{r}}u v\boldsymbol{\vec{u}\ \vec{v}} 平面内顺时针旋转 π/2\pi/2我们发现三重叉积就是我们想要的旋转,下面我们用矩阵的语言描述它:

considering:a×b=[0azayaz0axayax0]   a=[axayaz]Then:(u×v)×b=[0uyvxuxvyuzvxuxvzuxvyuyvx0uzvyuyvzuxvzuzvxuyvzuzvy0]bconsidering: \boldsymbol{\vec{a}}\times \boldsymbol{\vec{b}}=\left[ \begin{matrix} 0& -a_z& a_y\\ a_z& 0& -a_x\\ -a_y& a_x& 0\\ \end{matrix} \right] \ \ \ \boldsymbol{\vec{a}}=\left[ \begin{array}{c} a_x\\ a_y\\ a_z\\ \end{array} \right] \\ Then:\left( \overrightarrow{\mathbf{u}} {\times} \overrightarrow{\mathbf{v}} \right) {\times} \overrightarrow{\mathbf{b}} = \begin{bmatrix} 0 & u_{y}v_{x} {-} u_{x}v_{y} & u_{z}v_{x} {-} u_{x}v_{z} \\ u_{x}v_{y} {-} u_{y}v_{x} & 0 & u_{z}v_{y} {-} u_{y}v_{z} \\ u_{x}v_{z} {-} u_{z}v_{x} & u_{y}v_{z} {-} u_{z}v_{y} & 0 \end{bmatrix}\overrightarrow{\mathbf{b}}

注意到:

[0uyvxuxvyuzvxuxvzuxvyuyvx0uzvyuyvzuxvzuzvxuyvzuzvy0]=[uxvxuyvxuzvxuxvyuxvyuzvyuxvzuyvzuzvz][uxvxuxvyuxvzuyvxuyvyuyvzuzvxuzvyuzvz]=[vxvyvz][uxuyuz][uxuyuz][vxvyvz]=vuTuvT\left[ \begin{matrix} 0 & u_{y}v_{x} {-} u_{x}v_{y} & u_{z}v_{x} {-} u_{x}v_{z} \\ u_{x}v_{y} {-} u_{y}v_{x} & 0 & u_{z}v_{y} {-} u_{y}v_{z} \\ u_{x}v_{z} {-} u_{z}v_{x} & u_{y}v_{z} {-} u_{z}v_{y} & 0 \end{matrix} \right] \\= \left[ \begin{matrix} u_xv_x & u_yv_x & u_zv_x \\ u_xv_y & u_xv_y & u_zv_y \\ u_xv_z & u_yv_z & u_zv_z \end{matrix} \right] - \left[ \begin{matrix} u_xv_x & u_xv_y & u_xv_z \\ u_yv_x & u_yv_y & u_yv_z \\ u_zv_x & u_zv_y & u_zv_z \end{matrix} \right] \\ =\left\lbrack \begin{array}{l} v_{x} \\ v_{y} \\ v_{z} \end{array} \right\rbrack\left\lbrack \begin{array}{lll} u_{x} & u_{y} & u_{z} \end{array} \right\rbrack {-} \left\lbrack \begin{array}{l} u_{x} \\ u_{y} \\ u_{z} \end{array} \right\rbrack\left\lbrack \begin{array}{lll} v_{x} & v_{y} & v_{z} \end{array} \right\rbrack \\ = \boldsymbol {vu^T -uv^T}

​ 其实到这一步我们可以发现这就是外积(outer product)或更一般的楔积(wedge product),这里不展开讲。我们把这一步运算记为:

F(u,v)=vuTuvTF(\boldsymbol {u,v})=\boldsymbol {vu^T -uv^T}

虽然上面的公式是用叉乘来导出的,而不是所有的 nn 维欧式空间都可以定义叉乘,但是这个公式实际上是可以推广到任意 nn 维空间的,只是在三维空间内它恰好可以被翻译为叉乘

欧拉公式:

eθF(u,v)re^{\theta F(\boldsymbol {u, v})}\boldsymbol{\vec{r}}

代表着将 r\vec{r}u v\boldsymbol{\vec{u}\ \vec{v}} 平面内顺时针旋转 θπ/2\theta\pi/2 度(取弧度制)

eθF(u,v)=I+F(u,v)sinθ+F2(u,v)(1cosθ)e^{\theta F(\boldsymbol{u, v})}=I+F(\boldsymbol{u, v})\sin\theta + F^2(\boldsymbol {u, v})(1-\cos \theta)

可以自行对比一下二维复数情形,取 u,v\boldsymbol{u, v}[1,0]T,[0,1]T[1,0]^T, [0,1]^T

F(u,v)=[01][10][10][01]=[0010][0100]=[0110]F(\boldsymbol{u, v}) =\left[ \begin{matrix} 0 \\ 1 \\ \end{matrix}\right]\left[ \begin{matrix} 1 & 0 \end{matrix}\right] - \left[ \begin{matrix} 1 \\ 0 \\ \end{matrix}\right]\left[ \begin{matrix} 0 & 1 \end{matrix}\right] = \left[ \begin{matrix} 0 & 0 \\ 1 & 0\\ \end{matrix}\right] - \left[ \begin{matrix} 0 & 1 \\ 0 & 0 \end{matrix}\right] = \left[ \begin{matrix} 0 & -1 \\ 1 & 0 \end{matrix}\right]

​ 这其实就是前面举的例子!我们再继续把二维欧拉公式写下去:

eθF([10],[01])=I+[0110]sinθ+[1001](1cosθ)=[0110](sinθ+cosθ)eiθ=cosθ+isinθe^{\theta F\left(\left[\begin{matrix} 1 \\ 0 \\ \end{matrix}\right],\left[ \begin{matrix} 0 \\ 1 \\ \end{matrix} \right] \right)} = I + \left[ \begin{matrix} 0 & -1 \\ 1 & 0 \end{matrix}\right] \sin \theta + \left[ \begin{matrix} -1 & 0 \\ 0 & -1 \end{matrix}\right] ( 1-\cos \theta) \\ = \left[ \begin{matrix} 0 & -1 \\ 1 & 0 \end{matrix}\right] (\sin \theta + \cos \theta) \\ e^{i\theta} = \cos \theta + i\sin\theta

球坐标