当前位置:HOME - 线性代数入门答案 - 3.2 正交矩阵和QR分解

3.2 正交矩阵和QR分解

3.2.1-2
\\ 略。\\
阅读3.2.3 (Hadamard矩阵) \quad 给定nn阶矩阵AA,如果AA的元素都是111-1,且ATA=nInA^{\mathrm{T}}A=nI_{n},则称AA是一个nn 阶 Hadamard 矩阵。\\ 显然 Hadamard 是所有元素绝对值相同的正交矩阵的倍数. 可以证明 Hadamard 矩阵的阶只能是1,21,24k,k=1,2,.4k,k=1,2,\cdots.\\ 然而是否存在4k4k阶 Hadamard 矩阵,还是一个悬而未决的问题,称为 Hadamard 猜想。Hadamard 矩阵在信号处理中有应用。
\\
3.2.4
\\ 1.列举所有的一,二阶 Hadamard 矩阵。\\ 2.说明不存在三阶 Hadamard 矩阵。\\ 3.找出一个4阶Hadamard矩阵。\\ 4.证明如果AA是 Hadamard 矩阵,则[AAAA]\begin{bmatrix}A&A\\A&-A\end{bmatrix}也是 Hadamard 矩阵、以此说明存在2n2^n阶Hadamard 矩阵。\\ 1.解:\\ 一阶 Hadamard 矩阵:[1],[1].\begin{bmatrix} 1 \end{bmatrix},\begin{bmatrix} -1 \end{bmatrix}.\\ 二阶 Hadamard 矩阵:[1111],[1111],[1111],[1111],[1111],[1111],[1111],[1111].\begin{bmatrix} 1&1\\ 1&-1 \end{bmatrix},\begin{bmatrix} -1&1\\ 1&1 \end{bmatrix}, \begin{bmatrix} 1&1\\ -1&1 \end{bmatrix}, \begin{bmatrix} 1&-1\\ 1&1 \end{bmatrix}, \begin{bmatrix} -1&-1\\ -1&1 \end{bmatrix}, \begin{bmatrix} 1&-1\\ -1&-1 \end{bmatrix}, \begin{bmatrix} -1&-1\\ 1&-1 \end{bmatrix}, \begin{bmatrix} -1&1\\ -1&-1 \end{bmatrix}.\\ 2.证明:\\ 设三阶 Hadamard 矩阵A=[a1a2a3],a1=[x1x2x3],a2=[y1y2y3]A=\begin{bmatrix} \bm{a}_1&\bm{a}_2&\bm{a}_3 \end{bmatrix},\bm{a}_1=\begin{bmatrix} x_1\\ x_2\\ x_3 \end{bmatrix},\bm{a}_2=\begin{bmatrix} y_1\\ y_2\\ y_3 \end{bmatrix}。显然a1,a2\bm{a}_1,\bm{a}_2互相正交    a1Ta2=x1y1+x2y2+x3y3=0\implies \bm{a}_1^{\mathrm{T}}\bm{a}_2=x_1y_1+x_2y_2+x_3y_3=0。因为xiyix_iy_i只能取±1\pm 1,所以方程无解,即不存在三阶 Hadamard 矩阵。\\ 3.解:\\ 其中一个四阶 Hadamard 矩阵为[1111111111111111].\begin{bmatrix} 1&1&-1&1\\ -1&-1&-1&1\\ 1&-1&1&1\\ -1&1&1&1\\ \end{bmatrix}.\\ 4.证明:\\ [ATATATAT][AAAA]=[2nIn2nIn]=2nI2n\begin{bmatrix} A^{\mathrm{T}}&A^{\mathrm{T}}\\ A^{\mathrm{T}}&-A^{\mathrm{T}} \end{bmatrix} \begin{bmatrix} A&A\\ A&-A \end{bmatrix} = \begin{bmatrix} 2nI_n&\\ &2nI_n \end{bmatrix} =2nI_{2n}
3.2.5 \quad 证明,分块上三角矩阵[caT0Q]\begin{bmatrix}c&\bm{a}^{\mathrm{T}}\\ \bm{0}&Q\end{bmatrix} 是正交矩阵时,必有 c=±1,a=0,Qc=\pm 1,\bm{a}=\bm{0},Q 是正交矩阵。
\\ 利用QTQ=InQ^{\mathrm{T}}Q=I_n即可。略。\\
3.2.6 \quad 证明,上三角矩阵是正交矩阵时,必是对角矩阵,且对角元素是±1\pm 1
\\ 证明:\\A=[a11a12a1na22a2nann]A=\begin{bmatrix} a_{11}&a_{12}&\cdots&a_{1n}\\ &a_{22}&\cdots&a_{2n}\\ &&\ddots&\vdots\\ &&&a_{nn} \end{bmatrix}。考察列向量[a1100]\begin{bmatrix} a_{11}\\ 0\\ \vdots\\ 0 \end{bmatrix},由于列向量长度为1且互相正交,于是a11=±1,a12==a1n=0a_{11}=\pm 1,a_{12}=\cdots=a_{1n}=0。再考察 [0a220]\begin{bmatrix} 0\\ a_{22}\\ \vdots\\ 0 \end{bmatrix},可得a22=±1,a23==a2n=0a_{22}=\pm 1, a_{23}=\cdots=a_{2n}=0。由此类推可得AA是对角阵,且对角元为±1\pm 1\\
3.2.7 \quad对标准基e1,e2,,en\bm{e}_1,\bm{e}_2,\cdots,\bm{e}_n,显然i=1neieiT=In\displaystyle\sum_{i=1}^n\bm{e}_i\bm{e}_i^{\mathrm{T}}=I_n。对任意标准正交基q1,q2,,qn\bm{q}_1,\bm{q}_2,\cdots,\bm{q}_n,求证i=1nqiqiT=In.\displaystyle\sum_{i=1}^n\bm{q}_i\bm{q}_i^{\mathrm{T}}=I_n.\\
证明:\\Q=[q1qn]Q=\begin{bmatrix} \bm{q}_1&\cdots&\bm{q}_n \end{bmatrix},则i=1nqiqiT=[q1qn][q1TqnT]=QQT=In.\displaystyle\sum_{i=1}^n\bm{q}_i\bm{q}_i^{\mathrm{T}}=\begin{bmatrix} \bm{q}_1&\cdots&\bm{q}_n \end{bmatrix} \begin{bmatrix} \bm{q}_1^{\mathrm{T}}\\ \vdots\\ \bm{q}_n^{\mathrm{T}} \end{bmatrix}=QQ^{\mathrm{T}}=I_n.\\
3.2.8 \quada1,a2,,an\bm{a}_1,\bm{a}_2,\cdots,\bm{a}_nb1,b2,,bn\bm{b}_1,\bm{b}_2,\cdots,\bm{b}_nRn\mathbb{R}^n的两组标准正交基,证明存在正交矩阵QQ,使得Qai=bi,i=1,2,,n.Q\bm{a}_i=\bm{b}_i,i=1,2,\cdots,n.\\
证明:\\A=[a1an],B=[b1bn]A=\begin{bmatrix} \bm{a}_1&\cdots&\bm{a}_n \end{bmatrix},B=\begin{bmatrix} \bm{b}_1&\cdots&\bm{b}_n \end{bmatrix}A,BA,B显然都是n×nn\times n的正交矩阵。设存在QQ使Qai=biQ\bm{a}_i=\bm{b}_i,则QA=B    Q=BATQA=B \implies Q=BA^{\mathrm{T}}。显然此时QQ是正交矩阵。\\
3.2.9\quad回顾命题3.2.4,设a1,a2,,an\bm{a}_1,\bm{a}_2,\cdots,\bm{a}_nRn\mathbb{R}^n的一组基,考虑下列放松的条件。
\\ 1.AA在一组基上保距,即如果对i=1,2,,n,Aai=aii=1,2,\cdots,n,\Vert A\bm{a}_i \Vert = \Vert \bm{a}_i \Vert。那么AA是否一定是正交矩阵?\\ 2.AA在一组基上保内积,即如果对i,j=1,2,,n,Aaii,j=1,2,\cdots,n, A\bm{a}_iAajA\bm{a}_j的内积等于ai\bm{a}_iaj\bm{a}_j的内积。那么AA是否一定是正交矩阵?\\ 1.一定是正交矩阵。\\ 证明:\\ Aai=ai    aiTATAai=aiTai    aiT(ATAIn)ai=0\Vert A\bm{a}_i \Vert=\Vert \bm{a}_i \Vert \implies \bm{a}_i^{\mathrm{T}}A^{\mathrm{T}}A\bm{a}_i=\bm{a}_i^{\mathrm{T}}\bm{a}_i \implies \bm{a}_i^{\mathrm{T}}(A^{\mathrm{T}}A-I_n)\bm{a}_i=0。上式对所有基都成立,于是(ATAIn)ai=0    N(ATAIn)=Rn    ATA=In    A(A^{\mathrm{T}}A-I_n)\bm{a}_i=\bm{0} \implies \mathcal{N}(A^{\mathrm{T}}A-I_n)=\mathbb{R}^n \implies A^{\mathrm{T}}A=I_n \implies A是正交矩阵。\\ 2.一定是正交矩阵。\\ 证明:\\ aiTATAaj=aiTaj    aiT(ATAIn)aj\bm{a}_i^{\mathrm{T}}A^{\mathrm{T}}A\bm{a}_j=\bm{a}_i^{\mathrm{T}}\bm{a}_j \implies \bm{a}_i^{\mathrm{T}}(A^{\mathrm{T}}A-I_n)\bm{a}_j。上式对所有基都成立,于是(ATAIn)aj=0    N(ATAIn)=Rn    ATA=In    A(A^{\mathrm{T}}A-I_n)\bm{a}_j=\bm{0} \implies \mathcal{N}(A^{\mathrm{T}}A-I_n)=\mathbb{R}^n \implies A^{\mathrm{T}}A=I_n \implies A是正交矩阵。\\
3.2.10 \quadHv:=In2vvTH_v:=I_n-2vv^\mathrm{T}Rn\mathbb{R}^n上反射变换的表示矩阵,QQnn阶正交矩阵,证明 Q1HvQQ^{-1}H_vQ 也是某个反射变换的表示矩阵。
\\ 证明: Q1HvQ=QTQ2QTvvTQ=In2(QTv)(QTv)T.Q^{-1}H_vQ=Q^{\mathrm{T}}Q-2Q^{\mathrm{T}}\bm{v}\bm{v}^{\mathrm{T}}Q=I_n-2(Q^{\mathrm{T}}\bm{v})(Q^{\mathrm{T}}\bm{v})^{\mathrm{T}}.
3.2.11 \quad 证明命题 3.2.6. 给定Rn\mathbb{R}^n中向量x,y\bm{x},\bm{y},满足x=y\Vert \bm{x} \Vert = \Vert \bm{y} \Vert,则存在反射Hv\bm{H_v},其中v=yxyx\bm{v}=\cfrac{\bm{y}-\bm{x}}{\Vert \bm{y}-\bm{x} \Vert},使得Hv(x)=y.\bm{H_v}(\bm{x})=\bm{y}.\\
证明:\\ 显然反射变换的表示矩阵为H=In2vvT=In2(yx)(yx)Tyx2H=I_n-2\bm{v}\bm{v}^{\mathrm{T}}=I_n-\cfrac{2(\bm{y}-\bm{x})(\bm{y}-\bm{x})^{\mathrm{T}}}{\Vert \bm{y}-\bm{x} \Vert^2},且yTy=xTx\bm{y}^{\mathrm{T}}\bm{y}=\bm{x}^{\mathrm{T}}\bm{x}。于是Hv(x)=Hx=x2(yx)(yx)TxyTy+xTx2xTy=x2(yx)(yTxxTx)2xTx2xTy=x(yx)(xTyxTx)xTxxTy=x(xy)=y\begin{align*} \bm{H}_v(\bm{x})=H\bm{x}&=\bm{x}-\cfrac{2(\bm{y}-\bm{x})(\bm{y}-\bm{x})^{\mathrm{T}}\bm{x}}{\bm{y}^{\mathrm{T}}\bm{y}+\bm{x}^{\mathrm{T}}\bm{x}-2\bm{x}^{\mathrm{T}}\bm{y}}\\ &=\bm{x}-\cfrac{2(\bm{y}-\bm{x})(\bm{y}^{\mathrm{T}}\bm{x}-\bm{x}^{\mathrm{T}}\bm{x})}{2\bm{x}^{\mathrm{T}}\bm{x}-2\bm{x}^{\mathrm{T}}\bm{y}}\\ &=\bm{x}-\cfrac{(\bm{y}-\bm{x})(\bm{x}^{\mathrm{T}}\bm{y}-\bm{x}^{\mathrm{T}}\bm{x})}{\bm{x}^{\mathrm{T}}\bm{x}-\bm{x}^{\mathrm{T}}\bm{y}}\\ &=\bm{x}-(\bm{x}-\bm{y})\\ &=\bm{y} \end{align*}
3.2.12
\\ 略。\\
3.2.13 \quad 设向量组v1,v2,,vk\bm{v}_1,\bm{v}_2,\cdots,\bm{v}_k线性无关,首先令q1\bm{q}_1为与v1\bm{v}_1平行的单位向量,然后令q2\bm{q}_2为二维子空间 span(v1,v2)\operatorname{span}(\bm{v}_1,\bm{v}_2)中垂直于直线 span(v1)\operatorname{span}(\bm{v}_1)的单位向量,再令q3\bm{q}_3为三维子空间span(v1,v2,v3)\operatorname{span}(\bm{v}_1,\bm{v}_2,\bm{v}_3)中垂直于平面 span(v1,v2)\operatorname{span}(\bm{v}_1,\bm{v}_2)的单位向量,以此类推。这样得到的q1,q2,,qk\bm{q}_1,\bm{q}_2,\cdots,\bm{q}_k与 Gram-Schmidt 正交化得到的结果是否一致?如果有区别的话,区别在哪里?从QR分解的角度如何解释?
\\ 解:\\ 不一定一致。可能与GS正交化得到的向量方向相反。从QR分解的角度看,GS正交化得到的RR对角元全是正数,但通过这种方法得到的RR对角元不一定全为正。\\
3.2.14
\\ 略。\\
3.2.15 \quad 证明若矩阵列满秩,则其简化QR分解唯一。
\\ 证明:\\ 先证明一个引理,若AA是上三角矩阵且AA可逆,则A1A^{-1}也是上三角矩阵。\\ 显然,倍乘变换和将第ii行的kk倍加到第jj(i>j)(i \gt j)的倍加变换对应的矩阵是上三角阵。而AA经过有限步骤上述行变换即可得到InI_n,即E1E2EsA=InE_1E_2 \cdots E_sA=I_n,所以A1=E1E2EsA^{-1}=E_1E_2 \cdots E_s是上三角矩阵的乘积,也是上三角矩阵。\\ 引理证毕。\\ 设列满秩矩阵有两个QR分解A=Q1R1=Q2R2A=Q_1R_1=Q_2R_2,于是Q1R1=Q2R2    Q2TQ1=R2R11Q_1R_1=Q_2R_2 \implies Q_2^{\mathrm{T}}Q_1=R_2R_1^{-1}。注意R2,R11R_2,R_1^{-1}均为上三角矩阵且对角元全为正,所以R2R11R_2R_1^{-1}也是对角元全为正的上三角矩阵。而Q2TQ1=R2R11Q_2^{\mathrm{T}}Q_1=R_2R_1^{-1}是正交矩阵,由3.2.6中结论立得Q2TQ1=In    Q1=Q2    R1=R2Q_2^{\mathrm{T}}Q_1=I_n \implies Q_1=Q_2 \implies R_1=R_2。于是列满秩矩阵的QR分解唯一。\\
3.2.16 \quad 证明任意nn阶正交矩阵可以表示成不多于nn个反射的乘积。
\\ 证明:\\ 用数学归纳法,按照阶数归纳。\\ 阶数n=1n=1时显然成立。现假设阶数n=k1n=k-1时成立,考虑kk阶正交矩阵QQ\\Q=[q1q2qk]Q=\begin{bmatrix} \bm{q}_1&\bm{q}_2&\cdots&\bm{q}_k \end{bmatrix}。由于反射矩阵保距,所以存在反射矩阵H1H_1使得 H1Q=[e1H1q2H1qk]=[1aT0Q1]H_1Q=\begin{bmatrix} \bm{e}_1&H_1\bm{q}_2&\cdots&H_1\bm{q}_k \end{bmatrix} = \begin{bmatrix} 1&\bm{a}^{\mathrm{T}}\\ \bm{0}&Q_1 \end{bmatrix} 显然H1QH_1Q也是正交矩阵,由练习3.2.5中结论可得a=0\bm{a}=\bm{0}Q1Q_1也是正交矩阵,即H1Q=[1Q1]H_1Q= \begin{bmatrix} 1&\\ &Q_1 \end{bmatrix}。由归纳假设,Q1=H2~H3~Hk~Q_1=\widetilde{H_2}\widetilde{H_3}\cdots\widetilde{H_k},于是 H1Q=[1Q1]=[1H2~][1H3~][1Hk~].H_1Q=\begin{bmatrix} 1&\\ &Q_1 \end{bmatrix} = \begin{bmatrix} 1&\\ &\widetilde{H_2} \end{bmatrix} \begin{bmatrix} 1&\\ &\widetilde{H_3} \end{bmatrix}\cdots \begin{bmatrix} 1&\\ &\widetilde{H_k} \end{bmatrix}. Hi~=Ik12viviT,viRk1\widetilde{H_i}=I_{k-1}-2\bm{v}_i\bm{v}_i^{\mathrm{T}}, \bm{v}_i\in\mathbb{R}^{k-1},则[1Hi~]=In2[0vi][0vi]T\begin{bmatrix} 1&\\ &\widetilde{H_i} \end{bmatrix}=I_n-2\begin{bmatrix} 0\\ \bm{v}_i \end{bmatrix} \begin{bmatrix} 0\\ \bm{v}_i \end{bmatrix}^{\mathrm{T}}也是反射矩阵。令Hi=[1Hi~],i=2,3,,kH_i=\begin{bmatrix} 1&\\ &\widetilde{H_i} \end{bmatrix},i=2,3,\cdots,k,注意H1H1=IkH_1H_1=I_k,有H1Q=H2H3Hk    Q=H1H2Hk.H_1Q=H_2H_3\cdots H_k \implies Q=H_1H_2\cdots H_k.\\ 这说明若n=k1n=k-1时命题成立,则n=kn=k时命题也成立。于是n1n\geq 1时,任意nn阶正交矩阵可以表示成不多于nn个反射的乘积。\\
3.2.17 \quada1,a2,,as\bm{a}_1,\bm{a}_2,\cdots,\bm{a}_sb1,b2,,bs\bm{b}_1,\bm{b}_2,\cdots,\bm{b}_sRn\mathbb{R}^n的两个向量组,证明,存在正交矩阵QQ,使得Qai=bi,i=1,2,,sQ\bm{a}_i=\bm{b}_i,i=1,2,\cdots,s,当且仅当aiTaj=biTbj,i,j=1,2,,s.\bm{a}_i^\mathrm{T}\bm{a}_j=\bm{b}_i^\mathrm{T}\bm{b}_j,i,j=1,2,\cdots,s.\\
证明:\\ 先证存在正交矩阵,使得Qai=bi    aiTaj=biTbj(i,j=1,2,,s).Q\bm{a}_i=\bm{b}_i \implies \bm{a}_i^{\mathrm{T}}\bm{a}_j=\bm{b}_i^{\mathrm{T}}\bm{b}_j (i,j=1,2,\cdots,s).\\ 由于正交矩阵保内积的特性,命题显然成立。\\ 再证若aiTaj=biTbj(i,j=1,2,,s)    \bm{a}_i^{\mathrm{T}}\bm{a}_j=\bm{b}_i^{\mathrm{T}}\bm{b}_j (i,j=1,2,\cdots,s) \implies存在正交矩阵QQ使得Qai=bi,i=1,2,,s.Q\bm{a}_i=\bm{b}_i,i=1,2,\cdots,s.\\A=[a1,a2,,as],B=[b1,b2,,bs]A=\begin{bmatrix} \bm{a}_1,\bm{a}_2,\cdots,\bm{a}_s \end{bmatrix}, B=\begin{bmatrix} \bm{b}_1,\bm{b}_2,\cdots,\bm{b}_s \end{bmatrix}。取a1,a2,,as\bm{a}_1,\bm{a}_2,\cdots,\bm{a}_s中的一组极大线性无关向量组ai1,ai2,,aik,ks\bm{a}_{i_1},\bm{a}_{i_2},\cdots,\bm{a}_{i_k}, k\leq s。设A1=[ai1,ai2,,aik],B1=[bi1,bi2,,bik]A_1=\begin{bmatrix} \bm{a}_{i_1},\bm{a}_{i_2},\cdots,\bm{a}_{i_k} \end{bmatrix}, B_1=\begin{bmatrix} \bm{b}_{i_1},\bm{b}_{i_2},\cdots,\bm{b}_{i_k} \end{bmatrix}。显然A1A_1列满秩且a1,a2,,asR(A1)\bm{a}_1,\bm{a}_2,\cdots,\bm{a}_s \in \mathcal{R}(A_1)\\ 接下来我们证明B1B_1也列满秩且b1,b2,,bsR(B1)\bm{b}_1,\bm{b}_2,\cdots,\bm{b}_s \in \mathcal{R}(B_1)也成立。\\x,yRk\bm{x},\bm{y} \in \mathbb{R}^k满足A1x=0,B1y=0A_1\bm{x}=\bm{0},B_1\bm{y}=\bm{0}。显然x=0\bm{x}=\bm{0},由aiTaj=biTbj\bm{a}_i^{\mathrm{T}}\bm{a}_j=\bm{b}_i^{\mathrm{T}}\bm{b}_jaiTA1=biTB1\bm{a}_i^{\mathrm{T}}A_1=\bm{b}_i^{\mathrm{T}}B_1,于是0=aiTA1x=biTB1y=aiTA1y    A1y0=\bm{a}_i^{\mathrm{T}}A_1\bm{x}=\bm{b}_i^{\mathrm{T}}B_1\bm{y}=\bm{a}_i^{\mathrm{T}}A_1\bm{y} \implies A_1\bm{y}与所有ai\bm{a}_i都垂直    y=0    B1\implies \bm{y}=\bm{0} \implies B_1列满秩。\\aiTaj=biTbj(i,j=1,2,,s)\bm{a}_i^{\mathrm{T}}\bm{a}_j=\bm{b}_i^{\mathrm{T}}\bm{b}_j (i,j=1,2,\cdots,s)进一步可得ATA=BTBA^{\mathrm{T}}A=B^{\mathrm{T}}B,所以ATA=BTB    xRs,xTATAx=xTBTBx    Ax=Bx    N(A)=N(B)    rank(A)=rank(B) A^{\mathrm{T}}A=B^{\mathrm{T}}B \implies \forall \bm{x}\in\mathbb{R}^s, \bm{x}^{\mathrm{T}}A^{\mathrm{T}}A\bm{x}=\bm{x}^{\mathrm{T}}B^{\mathrm{T}}B\bm{x} \implies \Vert A\bm{x} \Vert=\Vert B\bm{x} \Vert \implies \mathcal{N}(A)=\mathcal{N}(B) \implies \operatorname{rank}(A)=\operatorname{rank}(B) 。这说明bi1,bi2,,bik\bm{b}_{i_1},\bm{b}_{i_2},\cdots,\bm{b}_{i_k}也是b1,b2,,bs\bm{b}_1,\bm{b}_2,\cdots,\bm{b}_s的一组极大线性无关向量组,即b1,b2,,bsR(B1)\bm{b}_1,\bm{b}_2,\cdots,\bm{b}_s \in \mathcal{R}(B_1)\\ 又设aj=A1v,bj=B1w\bm{a}_j=A_1\bm{v},\bm{b}_j=B_1\bm{w},则aiTaj=biTbj=aiTA1v=biTB1w=aiTA1w    aiTA1(vw),i=1,2,,s.    A1(vw)=0    v=w\bm{a}_i^{\mathrm{T}}\bm{a}_j=\bm{b}_i^{\mathrm{T}}\bm{b}_j=\bm{a}_i^{\mathrm{T}}A_1\bm{v}=\bm{b}_i^{\mathrm{T}}B_1\bm{w}=\bm{a}_i^{\mathrm{T}}A_1\bm{w} \implies \bm{a}_i^{\mathrm{T}}A_1(\bm{v}-\bm{w}),i=1,2,\cdots,s. \implies A_1(\bm{v}-\bm{w})=\bm{0} \implies \bm{v}=\bm{w}。于是ai=A1v    bi=B1v.\bm{a}_i=A_1\bm{v} \iff \bm{b}_i=B_1\bm{v}.\\A1A_1的简化QR分解为A1=QARAA_1=Q_AR_A,则QA=A1RA1Q_A=A_1R^{-1}_A。利用A1TA1=B1TB1A_1^{\mathrm{T}}A_1=B_1^{\mathrm{T}}B_1QATQA=RATA1TA1RA1=RATB1TB1RA1=IkQ_A^{\mathrm{T}}Q_A=R_A^{-\mathrm{T}}A_1^{\mathrm{T}}A_1R_A^{-1}=R_A^{-\mathrm{T}}B_1^{\mathrm{T}}B_1R_A^{-1}=I_k。所以QB=B1RA1Q_B=B_1R_A^{-1}也是正交矩阵    B1=QBRA=QBQATQARA=QBQATA1\implies B_1=Q_BR_A=Q_BQ_A^{\mathrm{T}}Q_AR_A=Q_BQ_A^{\mathrm{T}}A_1。显然QBQATQ_BQ_A^{\mathrm{T}}是正交矩阵,令Q=QBQATQ=Q_BQ_A^{\mathrm{T}},结合上面讨论就能得到bi=B1vi=QA1vi=Qai. \bm{b}_i=B_1\bm{v}_i=QA_1\bm{v}_i=Q\bm{a}_i. 证毕。\\
3.2.18 (保角变换) \quad 设可逆矩阵AA对应的线性变换保保持向量之间的夹角不变。
\\ 1.对AA进行QR分解,证明RR也保持向量之间的夹角不变。\\ 2.证明RR为对角矩阵。\\ 3.证明R=kInR=kI_n,这里kk为常数。由此得到,AA必是某个正交矩阵的倍数。\\ 1.证明:\\ a,bab=Aa,AbAaAb=QRa,QRbQRaQRb=Ra,RbRaRb.\cfrac{\langle \bm{a},\bm{b} \rangle}{\Vert \bm{a}\Vert \Vert \bm{b}\Vert}=\cfrac{\langle A\bm{a},A\bm{b} \rangle}{\Vert A\bm{a}\Vert \Vert A\bm{b}\Vert}=\cfrac{\langle QR\bm{a},QR\bm{b} \rangle}{\Vert QR\bm{a}\Vert \Vert QR\bm{b}\Vert}=\cfrac{\langle R\bm{a},R\bm{b} \rangle}{\Vert R\bm{a}\Vert \Vert R\bm{b}\Vert}. 2.证明:\\R=[a11a12a1na22a2nann]R=\begin{bmatrix} a_{11}&a_{12}&\cdots&a_{1n}\\ &a_{22}&\cdots&a_{2n}\\ &&\ddots&\vdots\\ &&&a_{nn} \end{bmatrix}e1,e2,,en\bm{e}_1,\bm{e}_2,\cdots,\bm{e}_n是标准基。\\ 由保角性,iji\neq j时有ei,ej=Rei,Rej=0    aij=0(ij).\langle \bm{e}_i, \bm{e}_j \rangle = \langle R\bm{e}_i, R\bm{e}_j \rangle=0 \implies a_{ij}=0 (i\neq j).于是RR是对角矩阵。\\ 3.证明:\\R=[a1a2an]R=\begin{bmatrix} a_{1}&&&\\ &a_{2}&&\\ &&\ddots&\\ &&&a_{n} \end{bmatrix}a1,a2,,an>0a_1,a_2,\cdots,a_n \gt 0e1,e2,,en\bm{e}_1,\bm{e}_2,\cdots,\bm{e}_n是标准基。\\ 由保角性,iji\neq j时有eiej,ei+ej=R(eiej),R(ei+ej)=ai2aj2=0.\langle \bm{e}_i-\bm{e}_j, \bm{e}_i+\bm{e}_j \rangle = \langle R(\bm{e}_i-\bm{e}_j), R(\bm{e}_i+\bm{e}_j) \rangle=a_i^2-a_j^2=0. 于是a1=a2==an    R=kIn.a_1=a_2=\cdots=a_n \implies R=kI_n.\\
3.2.19 \quada1,a2,,am\bm{a}_1,\bm{a}_2,\cdots,\bm{a}_mRn\mathbb{R}^n中的mm个向量,定义矩阵 G(a1,a2,,am):=[a1Ta1a1Ta2a1Tama2Ta1a2Ta2a2TamamTa1amTa2amTam],G(\bm{a}_1,\bm{a}_2,\cdots,\bm{a}_m):= \begin{bmatrix} \bm{a}_1^\mathrm{T}\bm{a}_1&\bm{a}_1^\mathrm{T}\bm{a}_2&\cdots&\bm{a}_1^\mathrm{T}\bm{a}_m\\ \bm{a}_2^\mathrm{T}\bm{a}_1&\bm{a}_2^\mathrm{T}\bm{a}_2&\cdots&\bm{a}_2^\mathrm{T}\bm{a}_m\\ \vdots&\vdots&&\vdots\\ \bm{a}_m^\mathrm{T}\bm{a}_1&\bm{a}_m^\mathrm{T}\bm{a}_2&\cdots&\bm{a}_m^\mathrm{T}\bm{a}_m \end{bmatrix}, 称为a1,a2,,am\bm{a}_1,\bm{a}_2,\cdots,\bm{a}_m的Gram矩阵。证明:
\\ 1.a1,a2,,am\bm{a}_1,\bm{a}_2,\cdots,\bm{a}_m是正交单位向量组当且仅当G(a1,a2,,am)=Im.G(\bm{a}_1,\bm{a}_2,\cdots,\bm{a}_m)=I_m.\\ 2.Gram矩阵G=G(a1,a2,,am)G=G(\bm{a}_1,\bm{a}_2,\cdots,\bm{a}_m)mm阶对称矩阵,且对任意xRm\bm{x}\in\mathbb{R}^{m},都有xTGx0.\bm{x}^{\mathrm{T}}G\bm{x}\geqslant 0.\\ 3.a1,a2,,am\bm{a}_1,\bm{a}_2,\cdots,\bm{a}_m线性无关当且仅当G=G(a1,a2,,am)G=G(\bm{a}_1,\bm{a}_2,\cdots,\bm{a}_m)可逆,也等价于对任意非零xRm\bm{x}\in\mathbb{R}^{m},都有xTGx>0.\bm{x}^{\mathrm{T}}G\bm{x} \gt 0.\\ 1.\\ 显然,略。\\ 2.证明:\\ 由内积的对称性得GG是对角矩阵。设A=[a1am]A=\begin{bmatrix} \bm{a}_1&\cdots&\bm{a}_m \end{bmatrix},则G=ATAG=A^{\mathrm{T}}A,于是xTGx=xTATAx=(Ax)T(Ax)0.\bm{x}^{\mathrm{T}}G\bm{x}=\bm{x}^{\mathrm{T}}A^{\mathrm{T}}A\bm{x}=(A\bm{x})^{\mathrm{T}}(A\bm{x})\geq 0.\\ 3.证明:\\ 考察方程Gx=[a1TAxamTAx]=0G\bm{x}=\begin{bmatrix} \bm{a}_1^{\mathrm{T}}A\bm{x}\\ \vdots\\ \bm{a}_m^{\mathrm{T}}A\bm{x}\\ \end{bmatrix}=\bm{0}。由于aiR(A)\bm{a}_i \in \mathcal{R}(A),所以aiTAx=0    Ax=0\bm{a}_i^{\mathrm{T}}A\bm{x}=\bm{0}\iff A\bm{x}=\bm{0}。即Gx=0G\bm{x}=\bm{0}只有零解当且仅当Ax=0A\bm{x}=\bm{0}只有零解。于是a1,a2,,am\bm{a}_1,\bm{a}_2,\cdots,\bm{a}_m线性无关    Ax=0\iff A\bm{x}=\bm{0}只有零解    Gx=0\iff G\bm{x}=\bm{0}只有零解    G\iff G可逆。\\ 显然a1,,am\bm{a}_1,\cdots,\bm{a}_m线性无关    Ax=0\iff A\bm{x}=\bm{0}只有零解     x0,xTGx=(Ax)T(Ax)=Ax2>0.\iff \forall \bm{x}\neq\bm{0}, \bm{x}^{\mathrm{T}}G\bm{x}=(A\bm{x})^{\mathrm{T}}(A\bm{x})=\Vert A\bm{x}\Vert^2 \gt 0.