Givens矩阵与Givens变换

Givens矩阵与Givens变换

1.引入

Givens矩阵与Givens变换
由图可得[x2y2]=[cos(θ)sin(θ)sin(θ)cos(θ)][x1y1] \left[ \begin{matrix} x_2 \\ y_2 \end{matrix} \right] = \left[ \begin{matrix} cos(\theta) & sin(\theta)\\ -sin(\theta) & cos(\theta) \end{matrix} \right] \left[ \begin{matrix} x_1 \\ y_1 \end{matrix} \right] 形如[cos(θ)sin(θ)sin(θ)cos(θ)]\left[ \begin{matrix} cos(\theta) & sin(\theta)\\ -sin(\theta) & cos(\theta) \end{matrix}\right]类型的矩阵被称为Givens矩阵。

2.Givens矩阵的定义

Givens矩阵与Givens变换
c2+s2=1c^2+s^2=1,称这样的矩阵为Givens矩阵(初等旋转矩阵),也记做Tij=Tij(c,s)T_{ij}=T_{ij}(c,s)
说明:

  • 因为c2+s2=1c^2+s^2=1,所以存在θ\theta,使得c=cos(θ),s=sin(θ)c=cos(\theta),s=sin(\theta)
  • 二阶情况下,对应的是一个平面直角坐标系中的一个旋转变换。

性质:

  • [Tij(c,s)]1=[Tij(c,s)]T=[Tij(c,s)],det[Tij(c,s)]=1[T_{ij}(c,s)]^{-1}=[T_{ij}(c,s)]^T=[T_{ij}(c,-s)],det[T_{ij}(c,s)]=1
  • x=[ξ1ξ2ξn]Tx=\left[ \begin{matrix}\xi_1 & \xi_2\cdots \xi_n\end{matrix}\right]^T,即x是一个列向量,有y=Tijx=[η1η2ηn]Ty=T_{ij}x=\left[ \begin{matrix}\eta_1 & \eta_2\cdots \eta_n\end{matrix}\right]^T,则有ηk={ξkki,jcξi+sξjk=isξi+cξjk=j\eta_k= \begin{cases} \xi_k & k\neq i,j \\ c\xi_i+s\xi_j & k=i \\ -s\xi_i+c\xi_j & k=j \end{cases}总可以选c=ξiξi2+ξj2,s=ξjξi2+ξj2c=\frac{\xi_i}{\sqrt{\xi_i^2+\xi_j^2}},s=\frac{\xi_j}{\sqrt{\xi_i^2+\xi_j^2}},使得y=[η1η2ξi2+ξj20ηn]Ty=\left[ \begin{matrix}\eta_1 & \eta_2&\sqrt{\xi_i^2+\xi_j^2}&\cdots & 0&\eta_n\end{matrix}\right]^T,其中ξi2+ξj2\sqrt{\xi_i^2+\xi_j^2}为第i项,0为第j项

3.两个定理

  • 定理一
    x=[ξ1ξ2ξn]Tx=\left[ \begin{matrix}\xi_1 & \xi_2\cdots \xi_n\end{matrix}\right]^T,存在有限个Givens矩阵的乘积TT,使得Tx=xe1Tx=|x|e_1,其中e1=[1,0...0]e_1=[1,0...0]
    证明:很容易,对x做T12,T13...T1nT_{12},T_{13}...T_{1n},即可得证

  • 定理二
    对任何非列零向量x,及任何单位向量z(z=1|z|=1),存在有限个Givens矩阵使得Tx=xzTx=|x|z
    证明:对左右两面分别进行上面定理的变换,则可得到T(1)x=xe1T^{(1)}x=|x|e_1T(2)xz=xe1T^{(2)}|x|z=|x|e_1 令两边相等即可得到结果