705 字
4 分钟

线性代数基础

2025-05-28
无标签
文章目录

Vector#

Vector Representation#

a=(a1,a2,,a3)T=[a1a2an]\vec{a} = (a_{1}, a_{2},\dots,a_{3})^T=\begin{bmatrix} a_{1}\\ a_{2} \\ \vdots \\ a_{n} \end{bmatrix}a=ax2+ay2+az2||\vec{a}||=\sqrt{ a_{x}^2+a_{y}^2+a_{z}^2 }

Dot product#

ab=a1b1+a2b2++anbn\vec{a}\cdot\vec{b}=a_{1}b_{1}+a_{2}b_{2}+\cdots +a_{n}b_{n}
  • In Eucidean Space,
    ab=a bcosθ\vec{a}\cdot\vec{b}=||\vec{a}||\space||\vec{b}||\cos \theta image.png

Cross product#

image.png
image.png

a×b=absin(θ)n\vec{a} \times \vec{b} = ||\vec{a}|| ||\vec{b}|| \sin (\theta) \vec{n}

注意, 叉积满足反交换律:

a×b=b×a\vec{a} \times \vec{b} = -\vec{b} \times \vec{a}

叉乘不满足结合律。

How to find the rotation between vector#

image.png

Dodrigues’s Rotation Formula#

image.png
Dodrigues’ Rotation Fomula:

b=a+(sinθ)u×a+(1cosθ)u×(u×a)\vec{b}=\vec{a}+(\sin \theta)\vec{u} \times \vec{a}+(1-\cos \theta)\vec{u}\times (\vec{u}\times \vec{a})

注意, 旋转方向的正负由左右手轴决定。针对左右手轴, 公式都是一致的。

同时这个公式也可以写成矩阵的形式(见 Matrix Form of Cross Product)。
image.png

Orthogonal Coordinate#

对于满足下面要求的轴满足右手系:

ex=ey=ez=1exey=eyez=ezex=0ex×ey=ezey×ez=exez×ex=ey\begin{split} ||e_{x}||=||e_{y}||=||e_{z}||=1 \\ e_{x}\cdot e_{y}=e_{y}\cdot e_{z}=e_{z}\cdot e_{x}=0 \\ e_{x}\times e_{y}=e_{z} \\ e_{y}\times e_{z}=e_{x} \\ e_{z}\times e_{x}=e_{y} \end{split}

这样可以得出以下的结论:

ab=axbx+ayby+azbza×b=(aybzazby,azbxaxbz,axbyaybz)\begin{align} \vec{a}\cdot \vec{b}= a_{x}b_{x}+a_{y}b_{y}+a_{z}b_{z} \\ \vec{a} \times \vec{b} = (a_{y}b_{z}-a_{z}b_{y}, a_{z}b_{x}-a_{x}b_{z},a_{x}b_{y}-a_{y}b_{z}) \end{align}

Matrix#

Matrix Operation#

image.png

Matrix Form of Dot Produnct#

image.png

Matrix Form of Cross Product#

叉乘可以表达为一个反对称阵的形式:

image.png
叉乘也可以用行列式进行表示:
image.png

叉乘满足下面的一些性质:
image.png

Orthogonal Matrix#

如果矩阵的列向量互相正交,且长度为1,那么该矩阵称之为标准正交矩阵,不要求矩阵满秩。
如果满秩,即 QQ 是方阵,称之为正交矩阵(Orthogonal Matrix)。标准正交矩阵有很多好的性质:

  • QTQ=IQ^TQ=I
  • QxQx 不改变 xx 的长度。Qx2=(Qx)TQx=xTQTQx=xTx=x2|Qx|^2=(Qx)^TQx=x^TQ^TQx=x^Tx=|x|^2
  • QQ 不改变向量点积。QxQy=(Qx)TQy=xTQTQy=xTyQx\cdot Qy=(Qx)^TQy=x^TQ^TQy=x^Ty
  • 如果 QQ 是方阵, QTQ=QQT=I    QT=Q1Q^TQ=QQ^T=I\implies Q^T=Q^{-1}

Determinant of a Matrix#

image.png
行列式的性质:

detI=1detAB=detAdetBdetAT=detAif A is invertible,detA1=(detA)1if U is orthogonal,detU=±1\begin{align} \det I = 1 \\ \det AB = \det A*\det B \\ \det A^{T}= \det A \\ if \ A \ is \ invertible, \det A^{-1} = (\det A)^{-1} \\ if \ U \ is \ orthogonal, \det U = \pm 1 \end{align}

Eigenvalues and Eigenvetors#

对于矩阵 AA, 存在一个非零向量 xx 满足:

Ax=λxAx=\lambda x

则称 λ\lambdaAA 的特征值, xxAA 的特征向量。

特别地, 对于 3×33\times {3} 的正交矩阵 UU, 它至少有一个特征值 λ=detU=±1\lambda=\det U=\pm {1}

Rigid Transformation#

Scaling#

image.png

Rotation#

向量 aa 经过旋转得到 aa', 对于 a=Raa'=Ra, 则称 RRa    aa\implies a' 的旋转矩阵。

旋转矩阵是一个正交矩阵, 满足一些性质:

R1=RTRTR=RRT=IdetR=+1Rx=x\begin{align} R^{-1}=R^{T} \\ R^{T}R = RR^{T} = I \\ \det R = +1 \\ ||Rx||=||x|| \end{align}

Combination of Rotations#

image.png

Rotation around Coordinate Axes#

image.png

Rotation Axis and Angle#

旋转矩阵始终有一个特征值 +1+1

也就是存在一个特征向量 Ru=uRu=u

这样, 旋转矩阵 RR 可以看做是 aa 绕轴 uu 旋转 θ\theta 角度。

image.png

这样反对称矩阵可以和叉乘建立关系(见线性代数基础)
image.png

这样 uu'uu 是平行的, 对 uu' 进行标准化后就能得到旋转轴 uu
image.png

同时根据线性代数基础 可以计算得出旋转角 θ\theta
image.png

View of Rotation#

image.png
image.png
对于原点相同的正交基, 可以看做 objectobject 的正交基是由 globalglobal 的正交基旋转 RR 得到的。
也就是说旋转矩阵给出了从局部坐标系到全局坐标系的转换公式。

image.png

加上位移后就是上图所示的坐标变换公式。