当前位置:HOME - 线性代数入门答案 - 6.2 正定矩阵

6.2 正定矩阵

6.2.1-3
\\ 略。\\
6.2.4 \quadAA对称正定,BB是实矩阵。\\ 1. 证明,对任意整数,AkA^k也正定。\\ 2. 若存在正整数rr,使得ArB=BArA^rB=BA^r,证明AB=BA.AB=BA.\\
1. 证明:\\ 显然AA有谱分解A=QΛQTA=Q\varLambda Q^{\mathrm{T}},则Ak=QΛkQTA^k=Q\varLambda^k Q^{\mathrm{T}}。而Λ>0    Λk>0\varLambda \gt 0 \implies \varLambda^k \gt 0,所以AkA^k也正定。\\ 2. 证明:\\ 因为AA对称正定,所以存在正交矩阵QQ,使QAQT=Λ=[λ1IλnI]QAQ^{\mathrm{T}}=\varLambda=\begin{bmatrix} \lambda_1 I&&\\ &\ddots&\\ &&\lambda_n I \end{bmatrix},其中λ1,λ2,,λn\lambda_1,\lambda_2,\cdots,\lambda_n均为正数且互不相等。于是 ArB=BAr    QTΛrQB=BQTΛrQ    ΛrQBQT=QBQTΛrA^rB=BA^r\implies Q^{\mathrm{T}}\varLambda^rQB=BQ^{\mathrm{T}}\varLambda^rQ \implies \varLambda^r QBQ^{\mathrm{T}}=QBQ^{\mathrm{T}}\varLambda^r 这说明QBQTQBQ^{\mathrm{T}}Λr\varLambda^r可交换。所以QBQT=[B11B22Bnn]QBQ^{\mathrm{T}}=\begin{bmatrix} B_{11}&&&\\ &B_{22}&&\\ &&\ddots&\\ &&&B_{nn} \end{bmatrix},显然其与Λ\varLambda也可交换。于是 AB=QTΛQB=QTΛQBQTQ=QTQBQTΛQ=BQTΛQ=BA.AB=Q^{\mathrm{T}}\varLambda QB=Q^{\mathrm{T}}\varLambda QBQ^{\mathrm{T}}Q=Q^{\mathrm{T}}QBQ^{\mathrm{T}}\varLambda Q=BQ^{\mathrm{T}}\varLambda Q=BA.
6.2.5 \quadnn阶实对称矩阵AA,证明,当实数tt充分大时,tIn+AtI_n+A正定。
\\ 证明:\\ xT(tIn+A)x=txTx+xTAx=tx+xTAx\bm{x}^{\mathrm{T}}(tI_n+A)\bm{x}=t\bm{x}^{\mathrm{T}}\bm{x}+\bm{x}^{\mathrm{T}}A\bm{x}=t\Vert \bm{x} \Vert + \bm{x}^{\mathrm{T}}A\bm{x}显然当tt充分大时,tx+xTAx>0t\Vert \bm{x} \Vert + \bm{x}^{\mathrm{T}}A\bm{x} \gt 0对任意x0\bm{x}\neq \bm{0}都成立,所以tIn+AtI_n+A正定。\\
6.2.6 \quadnn阶实对称矩阵AA,证明,存在正实数cc,使得对任意xRn\bm{x}\in\mathbb{R}^n,都有xTAxcxTx.\vert\bm{x}^{\mathrm{T}}A\bm{x}\vert\leq c \bm{x}^{\mathrm{T}}\bm{x}.\\
证明:\\ xTAxcxTx    xTAxxTxc\vert\bm{x}^{\mathrm{T}}A\bm{x}\vert\leq c \bm{x}^{\mathrm{T}}\bm{x} \implies \cfrac{\vert\bm{x}^{\mathrm{T}}A\bm{x}\vert}{\bm{x}^{\mathrm{T}}\bm{x}} \leq c。显然只需要cλmaxc \geq {\vert\lambda\vert}_{\max}即可。\\
6.2.7 \quad (Hadamard不等式) 给定nn阶对称正定矩阵AA,求证:
\\
1.对任意y\bm{y}det([AyyT0])0;\det(\begin{bmatrix} A&\bm{y}\\ \bm{y}^{\mathrm{T}}&0 \end{bmatrix})\leq 0;\\
证明:\\ 显然AA可逆,于是AyyT0=A0yTyTA1y=(yTA1y)det(A)\begin{vmatrix} A&\bm{y}\\ \bm{y}^{\mathrm{T}}&0\\ \end{vmatrix} = \begin{vmatrix} A&\bm{0}\\ \bm{y}^{\mathrm{T}}&-\bm{y}^{\mathrm{T}}A^{-1}\bm{y} \end{vmatrix} = -(\bm{y}^{\mathrm{T}}A^{-1}\bm{y}) \det(A),而A1A^{-1}显然也正定,即对任意y\bm{y}yTA1y0\bm{y}^{\mathrm{T}}A^{-1}\bm{y} \geq 0,等号成立当且仅当y=0\bm{y}=\bm{0}。而det(A)>0\det(A) \gt 0,所以对任意y\bm{y}det([AyyT0])0\det(\begin{bmatrix} A&\bm{y}\\ \bm{y}^{\mathrm{T}}&0 \end{bmatrix})\leq 0 \\
2.记A=[aij]A=\begin{bmatrix} a_{ij} \end{bmatrix},则det(A)anndet(An1)\det(A)\leq a_{nn}\det(A_{n-1}),其中An1A_{n-1}AAn1n-1阶顺序主子阵;\\
证明:\\ 显然An1A_{n-1}可逆,于是 det(A)=An1aaTann=A0aTannaTAn11a=(annaTAn11a)det(An1)\det(A)=\begin{vmatrix} A_{n-1}&\bm{a}\\ \bm{a}^{\mathrm{T}}&a_{nn}\\ \end{vmatrix} = \begin{vmatrix} A&\bm{0}\\ \bm{a}^{\mathrm{T}}&a_{nn}-\bm{a}^{\mathrm{T}}A_{n-1}^{-1}\bm{a} \end{vmatrix} = (a_{nn}-\bm{a}^{\mathrm{T}}A_{n-1}^{-1}\bm{a}) \det(A_{n-1}) 显然aTAn11a0\bm{a}^{\mathrm{T}}A_{n-1}^{-1}\bm{a} \geq 0对任意a\bm{a}成立,所以det(A)anndet(An1)\det(A)\leq a_{nn}\det(A_{n-1})\\
3.det(A)a11a22ann.\det(A) \leq a_{11}a_{22}\cdots a_{nn}.\\
证明:\\ 由上面结论立得。\\
利用上述结论证明:如果实矩阵T=[tij]T=\begin{bmatrix} t_{ij} \end{bmatrix}可逆,那么det(T)2i=1n(t1i2+t2i2++tni2).\det(T)^2\leq\overset{n}{\underset{i=1}{\prod}}(t^2_{1i}+t^2_{2i}+\cdots+t^2_{ni}).\\
证明:\\ det(T)2=det(TTT)\det(T)^2=\det(T^{\mathrm{T}}T),而矩阵TTTT^{\mathrm{T}}T显然正定,且其第ii行第ii列的对角元为t1i2+t2i2++tni2t_{1i}^2+t_{2i}^2+\cdots+t_{ni}^2。由上述结论得det(T)2i=1n(t1i2+t2i2++tni2).\det(T)^2\leq\overset{n}{\underset{i=1}{\prod}}(t^2_{1i}+t^2_{2i}+\cdots+t^2_{ni}).\\
6.2.8 \quad 证明A=[1i+j]n×nA=\begin{bmatrix} \cfrac{1}{i+j} \end{bmatrix}_{n \times n}正定。
\\ 证明:\\
法一:
\\ 注意到1i+j=01ti+j1dt\cfrac{1}{i+j}=\int_0^1t^{i+j-1}dtx=[x1x2xn]\bm{x}=\begin{bmatrix} x_1\\ x_2\\ \vdots\\ x_n \end{bmatrix},则 xTAx=1in1jnxixji+j=01t(1in1jnxixjti+j2)=01t(x1+x2t++xntn1)2dt\bm{x}^{\mathrm{T}}A\bm{x}=\sum_{\substack{1 \leq i \leq n\\1 \leq j \leq n}}\cfrac{x_ix_j}{i+j}=\int_0^1 t(\sum_{\substack{1 \leq i \leq n\\1 \leq j \leq n}}x_ix_jt^{i+j-2})=\int_0^1 t(x_1+x_2t+\cdots+x_nt^{n-1})^2 dt 显然对任意x0\bm{x}\neq\bm{0}xTAx>0\bm{x}^{\mathrm{T}}A\bm{x} \gt 0恒成立     \implies AA正定。\\
法二:
\\An=[1i+j]n×nA_n=\begin{bmatrix} \cfrac{1}{i+j} \end{bmatrix}_{n \times n},于是 det(An)=11+111+211+n12+112+212+n1n+11n+21n+n=n1(1+1)(n+1)n1(1+2)(n+2)n1(1+n)(n+n)n2(2+1)(n+1)n2(2+2)(n+2)n2(2+n)(n+n)1n+11n+21n+n=1in1n+i1i<n(ni)11+111+211+n12+112+212+n111=1in1n+i1i<n(ni)n1(1+1)(1+n)n2(1+2)(1+n)11+nn1(2+1)(2+n)n2(2+2)(2+n)12+n001=1i<n(ni)1in(n+i)n1(1+1)(1+n)n2(1+2)(1+n)n(n1)(1+n1)(1+n)n1(2+1)(2+n)n2(2+2)(2+n)n(n1)(2+n1)(2+n)n1(n1+1)(n1+n)n2(n1+2)(n1+n)n(n1)(n1+n1)(n1+n)=1i<n(ni)21in(n+i)1i<n(i+n)11+111+211+n112+112+212+n11n1+11n1+21n1+n1=1i<n(ni)2(n+n)1i<n(n+i)1i<n(i+n)An1 \begin{align*} \det(A_n)&=\begin{vmatrix} \cfrac{1}{1+1}&\cfrac{1}{1+2}&\cdots&\cfrac{1}{1+n}\\ \cfrac{1}{2+1}&\cfrac{1}{2+2}&\cdots&\cfrac{1}{2+n}\\ \vdots&\vdots&&\vdots\\ \cfrac{1}{n+1}&\cfrac{1}{n+2}&\cdots&\cfrac{1}{n+n}\\ \end{vmatrix}\\ &= \begin{vmatrix} \cfrac{n-1}{(1+1)(n+1)}&\cfrac{n-1}{(1+2)(n+2)}&\cdots&\cfrac{n-1}{(1+n)(n+n)}\\ \cfrac{n-2}{(2+1)(n+1)}&\cfrac{n-2}{(2+2)(n+2)}&\cdots&\cfrac{n-2}{(2+n)(n+n)}\\ \vdots&\vdots&&\vdots\\ \cfrac{1}{n+1}&\cfrac{1}{n+2}&\cdots&\cfrac{1}{n+n}\\ \end{vmatrix}\\ &= \prod_{1 \leq i \leq n}\cfrac{1}{n+i} \prod_{1 \leq i \lt n} (n-i) \begin{vmatrix} \cfrac{1}{1+1}&\cfrac{1}{1+2}&\cdots&\cfrac{1}{1+n}\\ \cfrac{1}{2+1}&\cfrac{1}{2+2}&\cdots&\cfrac{1}{2+n}\\ \vdots&\vdots&&\vdots\\ 1&1&\cdots&1\\ \end{vmatrix}\\ &= \prod_{1 \leq i \leq n}\cfrac{1}{n+i} \prod_{1 \leq i \lt n} (n-i) \begin{vmatrix} \cfrac{n-1}{(1+1)(1+n)}&\cfrac{n-2}{(1+2)(1+n)}&\cdots&\cfrac{1}{1+n}\\ \cfrac{n-1}{(2+1)(2+n)}&\cfrac{n-2}{(2+2)(2+n)}&\cdots&\cfrac{1}{2+n}\\ \vdots&\vdots&&\vdots\\ 0&0&\cdots&1\\ \end{vmatrix}\\ &= \cfrac{\displaystyle\prod_{1 \leq i \lt n} (n-i)}{\displaystyle\prod_{1 \leq i \leq n}(n+i)} \begin{vmatrix} \cfrac{n-1}{(1+1)(1+n)}&\cfrac{n-2}{(1+2)(1+n)}&\cdots&\cfrac{n-(n-1)}{(1+n-1)(1+n)}\\ \cfrac{n-1}{(2+1)(2+n)}&\cfrac{n-2}{(2+2)(2+n)}&\cdots&\cfrac{n-(n-1)}{(2+n-1)(2+n)}\\ \vdots&\vdots&&\vdots\\ \cfrac{n-1}{(n-1+1)(n-1+n)}&\cfrac{n-2}{(n-1+2)(n-1+n)}&\cdots&\cfrac{n-(n-1)}{(n-1+n-1)(n-1+n)}\\ \end{vmatrix}\\ &= \cfrac{\displaystyle\prod_{1 \leq i \lt n} (n-i)^2}{\displaystyle\prod_{1 \leq i \leq n}(n+i)\displaystyle\prod_{1 \leq i \lt n}(i+n)} \begin{vmatrix} \cfrac{1}{1+1}&\cfrac{1}{1+2}&\cdots&\cfrac{1}{1+n-1}\\ \cfrac{1}{2+1}&\cfrac{1}{2+2}&\cdots&\cfrac{1}{2+n-1}\\ \vdots&\vdots&&\vdots\\ \cfrac{1}{n-1+1}&\cfrac{1}{n-1+2}&\cdots&\cfrac{1}{n-1+n-1}\\ \end{vmatrix}\\ &=\cfrac{\displaystyle\prod_{1 \leq i \lt n} (n-i)^2}{(n+n) \displaystyle\prod_{1 \leq i \lt n}(n+i)\displaystyle\prod_{1 \leq i \lt n}(i+n)} A_{n-1} \end{align*} A1=12A_1=\frac{1}{2},于是det(An)=1i<jn(ij)21in1jn1(i+j)\det(A_n)=\displaystyle\prod_{1\leq i \lt j \leq n}(i-j)^2 \displaystyle\prod_{\substack{1\leq i \leq n\\1\leq j \leq n}}\cfrac{1}{(i+j)},显然任意顺序主子式都为正,所以AA正定。\\
6.2.9 \quad 证明Hilbert矩阵Hn=[1i+j1]n×nH_n=\begin{bmatrix} \cfrac{1}{i+j-1} \end{bmatrix}_{n \times n}正定。
\\ 证明:\\ 与6.2.8类似,也有两种方法。\\
法一:
\\ 注意到1i+j1=01ti+j2dt\cfrac{1}{i+j-1}=\displaystyle\int_0^1t^{i+j-2}dt, 设x=[x1x2xn]\bm{x}=\begin{bmatrix} x_1\\ x_2\\ \vdots\\ x_n \end{bmatrix},则 xTAx=1in1jnxixji+j1=01(1in1jnxixjti+j2)=01(x1+x2t++xntn1)2dt\bm{x}^{\mathrm{T}}A\bm{x}=\sum_{\substack{1 \leq i \leq n\\1 \leq j \leq n}}\cfrac{x_ix_j}{i+j-1}=\int_0^1 (\sum_{\substack{1 \leq i \leq n\\1 \leq j \leq n}}x_ix_jt^{i+j-2})=\int_0^1 (x_1+x_2t+\cdots+x_nt^{n-1})^2 dt 显然对任意x0\bm{x}\neq\bm{0}xTAx>0\bm{x}^{\mathrm{T}}A\bm{x} \gt 0恒成立     \implies AA正定。\\
法二:
\\ 法二与6.2.8中法二几乎完全一致,故从略。\\
6.2.10 \quadA,BA,B是实对称矩阵,如果ABA-B正定,则记作ABA\succ B。求证:\\ 1.AB,BCA\succ B,B\succ C可以推出 AC.A\succ C.\\ 2.ABA\succ BBAB\succ A不可能同时成立。\\ 3.对任意实对称矩阵AA,都存在实数k1,k2k_1,k_2使得k1InAk2In.k_1I_n\succ A\succ k_2I_n.\\
1. 证明:\\ ABA-B正定    \implies 对任意x0,xT(AB)x\bm{x}\neq\bm{0}, \bm{x}^{\mathrm{T}}(A-B)\bm{x}\\ BCB-C正定    \implies 对任意x0,xT(BC)x\bm{x}\neq\bm{0}, \bm{x}^{\mathrm{T}}(B-C)\bm{x}\\ 两式相加立得xT(AC)x\bm{x}^{\mathrm{T}}(A-C)\bm{x},所以ACA-C正定,即AC.A \succ C.\\ 2. 证明:\\ 假设同时成立,则{xT(AB)x>0xT(BA)x>0\begin{cases} \bm{x}^{\mathrm{T}}(A-B)\bm{x} \gt 0\\ \bm{x}^{\mathrm{T}}(B-A)\bm{x} \gt 0 \end{cases},两式相加得xT(AB+BA)x=0>0\bm{x}^{\mathrm{T}}(A-B+B-A)\bm{x}=0 \gt 0,矛盾。所以ABA\succ BBAB\succ A不可能同时成立。\\ 3. 证明:\\ 先证存在实数k1k_1使得k1InAk_1I_n\succ A,即(k1InA)(k_1I_n-A)正定。\\ 对任意非零向量x\bm{x}xT(k1InA)x=k1xTxxTAxk1xTxλmaxxTx \bm{x}^{\mathrm{T}}(k_1I_n-A)\bm{x}=k_1\bm{x}^{\mathrm{T}}\bm{x}-\bm{x}^{\mathrm{T}}A\bm{x}\geq k_1\bm{x}^{\mathrm{T}}\bm{x}-\lambda_{\max}\bm{x}^{\mathrm{T}}\bm{x}。显然当k1k_1足够大时,上式大于0对任意非零向量x\bm{x}均成立。于是存在实数k1k_1使得k1InAk_1I_n\succ A\\ 同理可得存在足够小的实数k2k_2使得Ak2In A\succ k_2I_n\\
6.2.11 \quad 举例说明,实对称矩阵AA的所有顺序主子式都非负,但AA并不半正定.
\\ 解:\\ [0001]\begin{bmatrix} 0&0\\ 0&-1 \end{bmatrix}所有顺序主子式都为0,但不半正定。\\
6.2.12 \quad 证明命题 6.2.4. 对实对称矩阵AA,以下叙述等价:\\ 1. AA半正定。\\2. AA的特征值都是非负数。\\3. 存在列满秩矩阵TT,使得A=TTTA=TT^{\mathrm{T}}\\4. AA存在LDLTLDL^{\mathrm{T}}分解,且DD的对角元都是非负数。\\
证明:\\ 采用轮转证法。\\ 1    21 \implies 2 \quad 任取AA的一个特征向量x\bm{x},则xTAx=λxTx0\bm{x}^{\mathrm{T}}A\bm{x}=\lambda\bm{x}^{\mathrm{T}}\bm{x}\geq 0,于是λ0.\lambda\geq 0.\\ 2    32 \implies 3 \quadA=QΛQTA=Q\varLambda Q^{\mathrm{T}},其中Λ=diag(λ1,,λr,0,,0)\varLambda=\operatorname{diag}(\lambda_1,\cdots,\lambda_r,0,\cdots,0)。因为Λ\varLambda对角元非负,所以设Λ=D2\varLambda=D^2,则A=(QD)(QD)T=[λ1q1λrqr00][λ1q1TλrqrT0T0T]A=(QD)(QD)^{\mathrm{T}}=\begin{bmatrix} \sqrt{\lambda_1}\bm{q}_1&\cdots&\sqrt{\lambda_r}\bm{q}_r&\bm{0}&\cdots&\bm{0} \end{bmatrix} \begin{bmatrix} \sqrt{\lambda_1}\bm{q}_1^{\mathrm{T}}\\ \vdots\\ \sqrt{\lambda_r}\bm{q}_r^{\mathrm{T}}\\ \bm{0}^{\mathrm{T}} \\ \vdots\\ \bm{0}^{\mathrm{T}} \end{bmatrix},设T=[λ1q1λrqr]T=\begin{bmatrix} \sqrt{\lambda_1}\bm{q}_1&\cdots&\sqrt{\lambda_r}\bm{q}_r \end{bmatrix},显然TT列满秩且A=TTT.A=TT^{\mathrm{T}}.\\ 3    43 \implies 4 \quadTTT^{\mathrm{T}}的QR分解为QL~TQ\widetilde{L}^{\mathrm{T}},其中QQ是列正交矩阵,L~\widetilde{L}是下三角矩阵。于是 A=TTT=L~QTQL~T=L~L~TA=TT^{\mathrm{T}}=\widetilde{L}Q^{\mathrm{T}}Q\widetilde{L}^{\mathrm{T}}=\widetilde{L}\widetilde{L}^{\mathrm{T}},令L~=LD~\widetilde{L}=L\widetilde{D},其中L~\widetilde{L}是单位下三角阵,D~\widetilde{D}是对角阵则A=L~L~T=LD~2LTA=\widetilde{L}\widetilde{L}^{\mathrm{T}}=L\widetilde{D}^2L^{\mathrm{T}},令D=D~2D=\widetilde{D}^2立得。 \\ 4    14 \implies 1 \quad xTAx=(LTx)TD(LTx)0.\bm{x}^{\mathrm{T}}A\bm{x}=(L^{\mathrm{T}}\bm{x})^{\mathrm{T}}D(L^{\mathrm{T}}\bm{x})\geq 0.\\
6.2.13 \quad 证明,实对称矩阵半正定,当且仅当它的所有主子式都非负,其中主子式的定义见练习 4.3.22.
\\ 证明:\\
法一:
\\ 先证:实对称矩阵AA半正定    \implies AA所有主子式都非负。\\ 显然,AA的某个kk阶主子阵可表示为PTAP,P=[ei1ei2eik]P^{\mathrm{T}}AP, P=\begin{bmatrix} \bm{e}_{i_1}&\bm{e}_{i_2}&\cdots&\bm{e}_{i_k} \end{bmatrix},其中i1<i2<<iki_1 \lt i_2 \lt \cdots \lt i_kei\bm{e}_i是标准基。于是对任意x0\bm{x}\neq\bm{0},都有xTPTAPx=(Px)TA(Px)0    \bm{x}^{\mathrm{T}}P^{\mathrm{T}}AP\bm{x}=(P\bm{x})^{\mathrm{T}}A(P\bm{x})\geq 0 \implies所有主子阵都半正定    A\implies A所有主子式都非负。\\ 再证:实对称矩阵AA所有主子式都非负    \implies AA半正定。\\ 用数学归纳法,按照阶数归纳。\\ 阶数n=1n=1时显然成立。现假设阶数小于nn时命题成立,考虑nn阶实对称矩阵AA。分两种情况讨论:\\ (a)若rank(A)<n\operatorname{rank}(A) \lt n\\x1,,xnr\bm{x}_1,\cdots,\bm{x}_{n-r}N(A)\mathcal{N}(A)的一组基。逐一考虑标准基e1,,en\bm{e}_1,\cdots,\bm{e}_n,若ei\bm{e}_ix1,,xnr\bm{x}_1,\cdots,\bm{x}_{n-r}线性无关,则将其加入基中。这样可以得到Rn\mathbb{R}^{n}的一组基(ei1,,eir,x1,,xnr)(\bm{e}_{i_1},\cdots,\bm{e}_{i_r},\bm{x}_1,\cdots,\bm{x}_{n-r}),其中i1<i2<<iki_1 \lt i_2 \lt \cdots \lt i_k。设X=[ei1eirx1xnr]X=\begin{bmatrix} \bm{e}_{i_1}&\cdots&\bm{e}_{i_r}&\bm{x}_1&\cdots&\bm{x}_{n-r} \end{bmatrix},则XTAX=[Cr×rOOO].X^{\mathrm{T}}AX=\begin{bmatrix} C_{r \times r}&O\\ O&O \end{bmatrix}.\\ 因为CCAA的一个主子阵,且AA所有主子式都非负,所以CC所有主子式也都非负。根据归纳假设,CC半正定,于是XTAXX^{\mathrm{T}}AX也半正定,而合同变换不改变正定性,所以AA也半正定。\\ 这说明若阶数小于nn时成立,矩阵阶数为nn时也成立,于是n1n \geq 1时命题均成立。\\ (b)若rank(A)=n\operatorname{rank}(A)=n\\ 显然此时det(A)>0\det(A) \gt 0,这说明AA只可能有0个或偶数个负数特征值。\\ 根据归纳假设,AAn1n-1阶顺序主子阵An1A_{n-1}半正定。假设AA有偶数个负数特征值,任意取其中两个负特征值λi,λj\lambda_i,\lambda_j,对应的特征向量是x,y\bm{x},\bm{y},注意xTy=yTx=0.\bm{x}^{\mathrm{T}}\bm{y}=\bm{y}^{\mathrm{T}}\bm{x}=0.\\u=αx+βy\bm{u}=\alpha\bm{x}+\beta\bm{y},存在一组α,β\alpha,\beta使得u\bm{u}的最后一个元素是为0,但u0\bm{u}\neq\bm{0}。即u=[u10],u10\bm{u}=\begin{bmatrix} \bm{u}_1\\ 0 \end{bmatrix},\bm{u}_1\neq\bm{0}。于是 uTAu=u1TAn1u1=(αx+βy)TA(αx+βy)=α2λixTx+αβλjxTy+αβλiyTx+β2λjyTy=α2λixTx+β2λjyTy<0 \begin{align*} \bm{u}^{\mathrm{T}}A\bm{u}&=\bm{u}_1^{\mathrm{T}}A_{n-1}\bm{u}_1 \\ &=(\alpha\bm{x}+\beta\bm{y})^{\mathrm{T}}A(\alpha\bm{x}+\beta\bm{y})\\ &= \alpha^2\lambda_i\bm{x}^{\mathrm{T}}\bm{x}+\alpha\beta\lambda_j\bm{x}^{\mathrm{T}}\bm{y}+\alpha\beta\lambda_i\bm{y}^{\mathrm{T}}\bm{x}+\beta^2\lambda_j\bm{y}^{\mathrm{T}}\bm{y}\\ &= \alpha^2\lambda_i\bm{x}^{\mathrm{T}}\bm{x}+\beta^2\lambda_j\bm{y}^{\mathrm{T}}\bm{y}\\ & \lt 0 \end{align*} 即存在u1\bm{u}_1使得u1TAn1u1<0\bm{u}_1^{\mathrm{T}}A_{n-1}\bm{u}_1 \lt 0,这与归纳假设相矛盾。所以AA没有负数特征值    \implies AA半正定。\\ 这说明若阶数小于nn时成立,矩阵阶数为nn时也成立,于是n1n \geq 1时命题均成立。\\
法二:
\\ 先证:实对称矩阵AA半正定    \implies AA所有主子式都非负。\\ 考虑矩阵A+ϵIA+\epsilon I,其中ϵ\epsilon是任意小的正实数。于是对任意x0\bm{x}\neq\bm{0},有xT(A+ϵI)x=xTAx+ϵxTx>0    A+ϵI\bm{x}^{\mathrm{T}}(A+\epsilon I)\bm{x}=\bm{x}^{\mathrm{T}}A\bm{x}+\epsilon\bm{x}^{\mathrm{T}}\bm{x} \gt 0 \implies A+\epsilon I正定    A+ϵI\implies A+\epsilon I的任意kk阶主子式det(Ak+ϵI)>0\det(A_k+\epsilon I) \gt 0。而det(Ak+ϵI)\det(A_k+\epsilon I)是关于ϵ\epsilon的连续函数,所以det(Ak)0\det(A_k)\geq 0,即AA所有主子式都非负。\\ 再证:实对称矩阵AA所有主子式都非负    \implies AA半正定。\\AkA_kAA的任一kk阶主子阵,det(Ak)0\det(A_k)\geq 0\\ 考虑矩阵A+ϵIA+\epsilon I,其中ϵ\epsilon是任意小的正实数。其任一kk阶主子式为det(Ak+ϵIk)\det(A_k+\epsilon I_k)。设λ\lambdaAkA_k的特征值,则det(AkλIk)=0    det(Ak+ϵIk(λ+ϵ)Ik)=0.\det(A_k-\lambda I_k)=0 \implies \det(A_k+\epsilon I_k - (\lambda+\epsilon)I_k)=0.这说明若λ\lambdaAkA_k的特征值,则λ+ϵ\lambda+\epsilonAk+ϵIkA_k+\epsilon I_k的特征值。因为ϵ\epsilon是任意小的正实数,所以若λ0\lambda\neq 0λ+ϵ\lambda+\epsilon不改变符号,若λ=0\lambda=0λ+ϵ>0\lambda+\epsilon \gt 0,于是det(Ak+ϵIk)>0    A+ϵI\det(A_k+\epsilon I_k) \gt 0 \implies A + \epsilon I任意顺序主子阵为正    A+ϵI\implies A+\epsilon I正定。\\ 所以对任意x0\bm{x}\neq\bm{0},有xT(A+ϵI)x=xTAx+ϵxTx>0\bm{x}^{\mathrm{T}}(A+\epsilon I)\bm{x}=\bm{x}^{\mathrm{T}}A\bm{x}+\epsilon\bm{x}^{\mathrm{T}}\bm{x} \gt 0,而xTAx+ϵxTx\bm{x}^{\mathrm{T}}A\bm{x}+\epsilon\bm{x}^{\mathrm{T}}\bm{x}是关于ϵ\epsilon的连续函数,于是xTAx0    A\bm{x}^{\mathrm{T}}A\bm{x}\geq 0 \implies A半正定。\\
6.2.14 \quadAA是实对称矩阵,证明:\\ 1.AA半正定,当且仅当存在实对称矩阵BB,使得A=B2.A=B^2.\\ 2.若AA半正定,则存在唯一的半正定实对称矩阵BB,使得A=B2.A=B^2.\\
1.证明:\\ AA是实对称矩阵且半正定    A=QTΛQ\iff A=Q^{\mathrm{T}}\varLambda Q,其中Λ\varLambda是对角元均为非负数的对角阵    \iff存在对角阵DD,使得A=QTD2Q=QTDQQTDQ=B2.A=Q^{\mathrm{T}} D^2 Q=Q^{\mathrm{T}}DQQ^{\mathrm{T}}DQ=B^2.\\ 2.证明:\\ A=QTΛQA=Q^{\mathrm{T}}\varLambda Q,其中Λ\varLambda是对角元均为非负数的对角阵,Λ=diag(λ1,λ2,,λn)\varLambda=\operatorname{diag}(\lambda_1,\lambda2,\cdots,\lambda_n),显然存在唯一的分解Λ=D2\varLambda=D^2使得DD也是对角元均为非负数的对角阵,即D=diag(λ1,λ2,,λn)D=\operatorname{diag}(\sqrt{\lambda_1},\sqrt{\lambda2},\cdots,\sqrt{\lambda_n})。此时A=QTD2Q=QTDQQTDQ=B2A=Q^{\mathrm{T}} D^2 Q=Q^{\mathrm{T}}DQQ^{\mathrm{T}}DQ=B^2,且B=QTDQB=Q^{\mathrm{T}}DQ是半正定矩阵。\\
6.2.15 \quad 证明,若实对称矩阵对角占优,且对角元素全为正数,则该矩阵正定。
\\ 证明:\\
法一:
\\AA对角占优,则det(A)0\det(A)\neq 0.\\ 考虑矩阵A+tInA+tI_n,其中tt是充分大的正实数。显然A+tInA+tI_n也对角占优,且xT(A+tIn)x=xTAx+txTx>0\bm{x}^{\mathrm{T}}(A+tI_n)\bm{x}=\bm{x}^{\mathrm{T}}A\bm{x}+t\bm{x}^{\mathrm{T}}\bm{x} \gt 0,于是A+tInA+tI_n正定。\\AkA_kAAkk阶顺序主子阵,则对任意kk,都有det(Ak+tIk)>0\det(A_k+tI_k) \gt 0。又因为det(Ak+tIk)\det(A_k+tI_k)是关于tt的连续函数,且对任意t0t\geq 0det(Ak+tIk)0\det(A_k+tI_k)\neq 0,而tt充分大时,det(Ak+tIk)>0\det(A_k+tI_k) \gt 0。所以det(Ak)>0\det(A_k) \gt 0对任意kk均成立    A\implies A正定。\\
法二:
\\x=[x1x2xn]\bm{x}=\begin{bmatrix} x_1\\ x_2\\ \vdots\\ x_n \end{bmatrix},则 xTAx=1in1jnaijxixj=1i<jn2aijxixj+1inaiixi2(由对称性)>1i<jn2aijxixj+1inxi2ijaij(由对角占优)=1i<jn2aijxixj+1i<jnaij(xi2+xj2)(由对称性)1i<jn2aijxixj+1i<jn2aijxixj(由Cauchy不等式)0 \begin{align*} \bm{x}^{\mathrm{T}}A\bm{x}&=\sum_{\substack{1\leq i \leq n\\ 1\leq j \leq n}}a_{ij}x_ix_j\\ &=\sum_{1\leq i \lt j \leq n}2a_{ij}x_ix_j+\sum_{1\leq i \leq n}a_{ii}x_i^2 \quad \text{(由对称性)}\\ & \gt \sum_{1\leq i \lt j \leq n}2a_{ij}x_ix_j+\sum_{1\leq i \leq n}x_i^2\sum_{i \neq j}\vert a_{ij} \vert \quad \text{(由对角占优)}\\ &=\sum_{1\leq i \lt j \leq n}2a_{ij}x_ix_j+\sum_{1\leq i \lt j \leq n}\vert a_{ij} \vert (x_i^2+x_j^2) \quad \text{(由对称性)}\\ &\geq \sum_{1\leq i \lt j \leq n}2a_{ij}x_ix_j+\sum_{1\leq i \lt j \leq n}2\vert a_{ij} \vert x_ix_j \quad \text{(由Cauchy不等式)}\\ &\geq 0 \end{align*} xTAx>0\bm{x}^{\mathrm{T}}A\bm{x} \gt 0对任意x0\bm{x}\neq\bm{0}恒成立    A\implies A正定。\\
6.2.16 \quad 证明或举出反例:
\\
1.如果AA对称正定,则A1A^{-1}正定。\\
证明:\\ AA对称正定则其特征值全是正数,而A1A^{-1}的特征值是AA特征值的倒数,显然也全为正数,所以A1A^{-1}正定。\\
2.如果A,BA,B对称正定,则A+BA+B正定。\\
显然正确。\\
3.如果A,BA,B对称半正定,则A+BA+B半正定。\\
显然正确。\\
4.如果A,BA,B对称不定,则A+BA+B不定。\\
不正确。\\A=[1111]A=\begin{bmatrix} -1&1\\1&1\\ \end{bmatrix},其特征值是±2\pm\sqrt{2},不定。 B=[1112]B=\begin{bmatrix} 1&-1\\ -1&-2 \end{bmatrix},其特征值是1±132\cfrac{-1\pm\sqrt{13}}{2},不定。而A+B=[0001]A+B=\begin{bmatrix} 0&0\\ 0&-1 \end{bmatrix},其特征值是001-1,半负定。\\
5.如果AA列满秩,BB对称正定,则ATBAA^{\mathrm{T}}BA正定。\\
证明:\\ AA列满秩,所以当x0\bm{x}\neq\bm{0}时,Ax0A\bm{x}\neq\bm{0}。而BB对称正定,即对任意x0,xTBx>0\bm{x}\neq\bm{0}, \bm{x}^{\mathrm{T}}B\bm{x} \gt 0。所以xTATBAx=(Ax)TB(Ax)>0    ATBA\bm{x}^{\mathrm{T}}A^{\mathrm{T}}BA\bm{x}=(A\bm{x})^{\mathrm{T}}B(A\bm{x}) \gt 0 \implies A^{\mathrm{T}}BA正定。\\
6.如果S=ATAS=A^{\mathrm{T}}AAA有简化 QR 分解A=QRA=QR,则S=RTRS=R^{\mathrm{T}}RSS的 Cholesky 分解。\\
证明:\\ S=ATA=(QR)T(QR)=RTQTQR=RTRS=A^{\mathrm{T}}A=(QR)^{\mathrm{T}}(QR)=R^{\mathrm{T}}Q^{\mathrm{T}}QR=R^{\mathrm{T}}R,其中RR是上三角阵。\\
7.如果A,BA,B对称正定,则ABAB的特征值都是正数。\\
ABAB的特征值是λ\lambda,对应的特征向量是x\bm{x},则ABx=λxAB\bm{x}=\lambda\bm{x}。两边左乘xTBT\bm{x}^{\mathrm{T}}B^{\mathrm{T}}xTBTABx=λxTBTx    (Bx)TA(Bx)=λxTBTx\bm{x}^{\mathrm{T}}B^{\mathrm{T}}AB\bm{x}=\lambda\bm{x}^{\mathrm{T}}B^{\mathrm{T}}\bm{x}\implies (B\bm{x})^{\mathrm{T}}A(B\bm{x})=\lambda\bm{x}^{\mathrm{T}}B^{\mathrm{T}}\bm{x},而A,BA,B都正定,(Bx)TA(Bx)>0,xTBTx>0(B\bm{x})^{\mathrm{T}}A(B\bm{x}) \gt 0,\bm{x}^{\mathrm{T}}B^{\mathrm{T}}\bm{x} \gt 0,所以λ>0.\lambda \gt 0.\\
6.2.17 \quad 证明命题 6.2.7. 方阵的合同关系是等价关系。
\\ 证明:\\ 反身性:A=InTAInA=I_n^{\mathrm{T}}AI_n,因此AAAA合同。\\ 对称性:若A,BA,B合同,则存在可逆矩阵XX,使得XTAX=BX^{\mathrm{T}}AX=B,显然同时也有XTBX1=AX^{-\mathrm{T}}BX^{-1}=A。所以BBAA也合同。\\ 传递性:若A,BA,B合同,B,CB,C合同,则存在可逆矩阵X,YX,Y,使得XTAX=B,YTBY=CX^{\mathrm{T}}AX=B,Y^{\mathrm{T}}BY=C。即(XY)TA(XY)=C(XY)^{\mathrm{T}}A(XY)=C,而XYXY是可逆矩阵,因此A,CA,C也合同。\\ 所以方阵的合同关系是等价关系。\\
6.2.18 \quad 考虑实对称矩阵S=[ABBTC].S=\begin{bmatrix}A&B\\B^\mathrm{T}&C\end{bmatrix}.证明SS正定当且仅当AA及其 Schur 补都正定。
\\ 证明:\\ 注意到[ABBTC]=[IBTA1I][ACBTA1B][IA1BI]\begin{bmatrix} A&B\\ B^{\mathrm{T}}&C \end{bmatrix} = \begin{bmatrix} I&\\ B^{\mathrm{T}}A^{-1}&I \end{bmatrix} \begin{bmatrix} A&\\ &C-B^{\mathrm{T}}A^{-1}B \end{bmatrix} \begin{bmatrix} I&A^{-1}B\\ &I \end{bmatrix} 于是SS正定    S\iff SLDLTLDL^{\mathrm{T}}分解    \iff存在下三角阵L1,L2L_1,L_2使得A=L1DL1T,CBTA1B=L2DL2T    A=L_1DL_1^{\mathrm{T}},C-B^{\mathrm{T}}A^{-1}B=L_2DL_2^{\mathrm{T}} \iff AA及其Schur补都正定。\\
6.2.19 \quadA,BA,Bnn阶实对称矩阵,AA正定。证明,存在可逆矩阵TT,使得TTATT^\mathrm{T}ATTTBTT^\mathrm{T}BT同时是对角矩阵。
\\ 证明:\\ AA是正定的实对称矩阵,所以存在正交矩阵QAQ_A使得QATAQA=ΛAQ_A^{\mathrm{T}}AQ_A=\varLambda_A。设ΛA=D2\varLambda_A=D^2,则(QAD1)TA(QAD1)=In.(Q_AD^{-1})^{\mathrm{T}}A(Q_AD^{-1})=I_n.\\ 注意到(QAD1)TB(QAD1)(Q_AD^{-1})^{\mathrm{T}}B(Q_AD^{-1})也是实对称矩阵,存在谱分解QT(QAD1)TB(QAD1)Q=ΛQ^{\mathrm{T}}(Q_AD^{-1})^{\mathrm{T}}B(Q_AD^{-1})Q=\varLambda,而QT(QAD1)TA(QAD1)Q=QTQ=InQ^{\mathrm{T}}(Q_AD^{-1})^{\mathrm{T}}A(Q_AD^{-1})Q=Q^{\mathrm{T}}Q=I_n。所以存在可逆矩阵T=QAD1QT=Q_AD^{-1}Q使得TTATT^\mathrm{T}ATTTBTT^\mathrm{T}BT同时是对角矩阵。\\
6.2.20 \quadA,BA,Bnn阶实对称矩阵,A,BA,B半正定。证明,存在可逆矩阵TT,使得TTATT^\mathrm{T}ATTTBTT^\mathrm{T}BT同时是对角矩阵。
\\ 证明:\\ A,BA,B半正定    \implies任意x0,xTAx0,xTBx0\bm{x}\neq\bm{0},\bm{x}^{\mathrm{T}}A\bm{x}\geq 0, \bm{x}^{\mathrm{T}}B\bm{x}\geq 0。分两种情况讨论:\\ (a)若N(A)N(B)={0}.\mathcal{N}(A)\cap\mathcal{N}(B)=\set{\bm{0}}.\\ 对任意x0,xT(A+B)x=xTAx+xTBx>0    A+B\bm{x}\neq\bm{0},\bm{x}^{\mathrm{T}}(A+B)\bm{x}=\bm{x}^{\mathrm{T}}A\bm{x}+\bm{x}^{\mathrm{T}}B\bm{x} \gt 0 \implies A+B正定。显然A+BA+B也是对称矩阵,由6.2.19,存在可逆矩阵TT,使得TT(A+B)T=Λ1T^\mathrm{T}(A+B)T=\varLambda_1TTAT=Λ2T^\mathrm{T}AT=\varLambda_2同时是对角矩阵。而TTBT=TT(A+B)TTTAT=Λ1Λ2T^\mathrm{T}BT=T^\mathrm{T}(A+B)T-T^\mathrm{T}AT=\varLambda_1-\varLambda_2也是对角矩阵。所以存在可逆矩阵TT,使得TTATT^\mathrm{T}ATTTBTT^\mathrm{T}BT同时是对角矩阵。\\ (b)若N(A)N(B){0}.\mathcal{N}(A)\cap\mathcal{N}(B)\neq\set{\bm{0}}.\\N(A)N(B)\mathcal{N}(A)\cap\mathcal{N}(B)的一组正交基(q1,,qs)(\bm{q}_1,\cdots,\bm{q}_s),并扩张为正交矩阵QQ,则QTAQ=[OA1],QTBQ=[OB1].Q^{\mathrm{T}}AQ=\begin{bmatrix} O&\\ &A_1 \end{bmatrix}, Q^{\mathrm{T}}BQ=\begin{bmatrix} O&\\ &B_1 \end{bmatrix}.\\ 显然此时N(A1)N(B1)={0}\mathcal{N}(A_1)\cap\mathcal{N}(B_1)=\set{\bm{0}},问题即转化为(a)中情况。