5.9_课后习题

5.9 课后习题

练习5.1设 ARn×nA\in \mathbb{R}^{n\times n} ,若 (λ1,x1)(\lambda_1,x_1)(λ2,x2)(\lambda_{2},x_{2})AA 的特征值对,且 λ1λ2\lambda_1\neq \lambda_2 ,证明: x1x_{1}x2x_{2} 线性无关.进一步,若 AA 对称,则 x1x_{1}x2x_{2} 正交

练习5.2设 x,yRnx,y\in \mathbb{R}^n ,试证明: det(I+xy)=1+yx.\operatorname *{det}(I + xy^{\top}) = 1 + y^{\top}x. (注:在复数域也成立)

练习5.3设 A=D+uuTA = D + uu^{\mathsf{T}} ,其中 D=diag(d1,d2,,dn)D = \mathrm{diag}(d_1,d_2,\ldots ,d_n) 满足 d1d2dn,u=[u1,u2,,un]TRn.d_{1}\geq d_{2}\geq \dots \geq d_{n},u = [u_{1},u_{2},\ldots ,u_{n}]^{\mathsf{T}}\in \mathbb{R}^{n}.

(1) 证明: did_{i}AA 的特征值的充要条件是 di=di+1d_{i} = d_{i + 1}di=di1d_{i} = d_{i - 1}ui=0u_{i} = 0 ;
(2) 若 ui=0u_{i} = 0 , 则 eie_{i} 是与 did_{i} 对应的特征向量;
(3) 若 di1>di=di+1>di+2d_{i-1} > d_i = d_{i+1} > d_{i+2}ui0u_i \neq 0 , 证明: 对应于 λ=di\lambda = d_i 的特征向量 xxxix_ixi+1x_{i+1} 外, 其余分量全部为 0, 且 xiui+xi+1ui+1=0x_i u_i + x_{i+1} u_{i+1} = 0 .

思考: 如果 di1>di=di+1=di+2>di+3d_{i-1} > d_i = d_{i+1} = d_{i+2} > d_{i+3} , 则结论如何?

练习 5.45.4^{*}αR,D=diag(d1,d2,,dn)\alpha \in \mathbb{R}, D = \mathrm{diag}(d_1, d_2, \ldots, d_n) . 若 u,u^Rnu, \hat{u} \in \mathbb{R}^n 非负, 矩阵 D+αuuD + \alpha uu^{\top}D+u^u^D + \hat{u}\hat{u}^{\top} 具有相同的特征值, 记为 λ1,λ2,,λn\lambda_1, \lambda_2, \ldots, \lambda_n , 且满足交错性质 λ1>d1>λ2>d2>>λn>dn\lambda_1 > d_1 > \lambda_2 > d_2 > \dots > \lambda_n > d_n , 则它们具有相同的特征向量.

思考: 如果没有交错性质, 则结论如何?

练习 5.55.5^{*}qRnq\in \mathbb{R}^n 满足 q2=1\| q\| _2 = 1 .对任意与 qq 正交的向量 dRnd\in \mathbb{R}^n ,试证明:

(q+d)qI2=q+d2.\| (q + d) q ^ {\intercal} - I \| _ {2} = \| q + d \| _ {2}.

练习5.6设 SCn×nS\in \mathbb{C}^{n\times n} 是skew-Hermite矩阵,即 S=SS^{*} = -S 证明:

(1) SS 的非零特征值是纯虚数;
(2) I+SI + S 非奇异;
(3) 矩阵 (I+S)1(IS)(I + S)^{-1}(I - S) 是酉矩阵. (该矩阵称为 SS 的 Cayley 变换)

练习5.7 设 BRm×n,mnB \in \mathbb{R}^{m \times n}, m \geq nB2<1\| B \|_2 < 1 . 若 A=[IBBI]A = \left[ \begin{array}{cc} I & B \\ B^\top & I \end{array} \right] , 证明:

κ2(A)=1+B21B2.\kappa_ {2} (A) = \frac {1 + \| B \| _ {2}}{1 - \| B \| _ {2}}.

练习 5.85.8^{*}BB 是二对角矩阵

B=[a1b1bn1an].B = \left[ \begin{array}{c c c c} a _ {1} & b _ {1} & & \\ & \ddots & \ddots & \\ & & \ddots & b _ {n - 1} \\ & & & a _ {n} \end{array} \right].

证明: 存在正交矩阵 Q1Q_{1}Q2Q_{2} , 使得 Q1TBQ2Q_{1}^{\mathrm{T}} B Q_{2} 仍然是二对角矩阵且所有元素都非负.

练习5.9设 x,yRnx,y\in \mathbb{R}^n ,若 yTxy^{\mathrm{T}}x 只有零特征值,证明: xyxy^{\top} 也只有零特征值

X,YRn×2X, Y \in \mathbb{R}^{n \times 2} , 若 YTXY^{\mathsf{T}}X 只有零特征值, 则 XYTXY^{\mathsf{T}} 是否也只有零特征值?

练习 5.105.10^{*} (极分解) 设 ACn×nA \in \mathbb{C}^{n \times n} . 证明:

(1)存在酉矩阵 UU 和唯一的Hermite半正定矩阵 PP ,使得 A=PUA = PU

(2) 进一步, 若 AA 非奇异, 则 UU 也唯一.

练习 5.115.11^{*}ACn×nA \in \mathbb{C}^{n \times n} . 证明: AA 可对角化当且仅当存在 Hermite 正定矩阵 PP 使得 P1APP^{-1}AP 是正规矩阵.

(提示: 利用极分解, 但不是对 AA 进行极分解)

以下为可选题

练习5.12 设 λR\lambda \in \mathbb{R} 是对称矩阵 ARn×nA \in \mathbb{R}^{n \times n} 的一个特征值, 对应的特征向量为 xRnx \in \mathbb{R}^n . 若 x~Rn\tilde{x} \in \mathbb{R}^nxx 的一个 O(ε)\mathcal{O}(\varepsilon) 近似, 即 x~=x+O(ε)\tilde{x} = x + \mathcal{O}(\varepsilon) , 证明:

x~TAx~x~Tx~=λ+O(ε2),\frac {\tilde {x} ^ {\mathsf {T}} A \tilde {x}}{\tilde {x} ^ {\mathsf {T}} \tilde {x}} = \lambda + \mathcal {O} (\varepsilon^ {2}),

x~\tilde{x} 对应的Rayleigh商是 λ\lambdaO(ε2)\mathcal{O}(\varepsilon^2) 逼近

练习5.13 设 A,ERn×nA, E \in \mathbb{R}^{n \times n} 都是对称矩阵,它们的特征值分别为 λ1λ2λn\lambda_1 \geq \lambda_2 \geq \dots \geq \lambda_nθ1θ2θn\theta_1 \geq \theta_2 \geq \dots \geq \theta_n . 设 A+EA + E 的特征值为 λ^1λ^2λ^n\hat{\lambda}_1 \geq \hat{\lambda}_2 \geq \dots \geq \hat{\lambda}_n ,试证明

λi+θnλ^iλi+θ1,i=1,2,,n.\lambda_ {i} + \theta_ {n} \leq \hat {\lambda} _ {i} \leq \lambda_ {i} + \theta_ {1}, \quad i = 1, 2, \dots , n.

并由此可知, 若 EE 对称正定, 则 λ^iλi\hat{\lambda}_i \geq \lambda_i .

练习5.14设 ARn×nA\in \mathbb{R}^{n\times n} 是对称矩阵, An1A_{n - 1}AAn1n - 1 阶顺序主子矩阵,它们的特征值分别为 λ1λ2λn\lambda_1\geq \lambda_2\geq \dots \geq \lambda_nθ1θ2θn1\theta_{1}\geq \theta_{2}\geq \dots \geq \theta_{n - 1} .试证明

λ1θ1λ2θ2λn1θn1λn.\lambda_ {1} \geq \theta_ {1} \geq \lambda_ {2} \geq \theta_ {2} \geq \dots \geq \lambda_ {n - 1} \geq \theta_ {n - 1} \geq \lambda_ {n}.

更一般地, 记 λi(B)\lambda_{i}(B) 为对称矩阵 BB 的第 ii 个特征值 (按降序排列), 设 AkA_{k}AA 任意一个 kk 阶主子矩阵 (1kn1)(1 \leq k \leq n - 1) , 则有

λi(A)λi(Ak)λnk+i(A),i=1,2,,k.\lambda_ {i} (A) \geq \lambda_ {i} \left(A _ {k}\right) \geq \lambda_ {n - k + i} (A), \quad i = 1, 2, \dots , k.

练习5.15 设 A=D+uuA = D + uu^{\top} , 其中 D=diag(d1,d2,,dn)D = \mathrm{diag}(d_1, d_2, \ldots, d_n) 满足 d1>d2>>dnd_1 > d_2 > \dots > d_n . 由此可知, AAn1n - 1 个特征值满足 λi+1(di+1,di),i=1,2,,n1\lambda_{i + 1} \in (d_{i + 1}, d_i), i = 1, 2, \ldots, n - 1 . 试给出特征值 λ1\lambda_1 的取值范围.

练习 5.165.16^{*} 证明以下结论:

(1) 设 xRnx \in \mathbb{R}^n 是一个正向量, 即 xi>0x_i > 0 . 证明: 由 xx 定义的 Cauchy 矩阵 AA

aij=1xi+xja _ {i j} = \frac {1}{x _ {i} + x _ {j}}

是对称半正定的. 进一步, 若 xix_{i} 互不相等, 则 AA 对称正定. (参见 [153])

(2) 证明: Hilbert 矩阵是对称正定的.

以下为实践题

练习5.17 编写程序, 实现对称矩阵 ARn×nA \in \mathbb{R}^{n \times n} 的三对角化

练习5.18 编写程序, 实现矩阵 ARm×nA \in \mathbb{R}^{m \times n} ( mnm \geq n ) 的二对角化

练习5.19 编写程序, 实现计算对称三对角矩阵特征值的带Wilkinson位移的QR算法.

练习5.20 编写程序, 实现计算二对角矩阵奇异值的 dqds 算法