矩阵论期末考试（二）

1⃣️求微分方程的解（一般结合矩阵函数考）

1.齐次微分方程（x^’(t)=A·x(t)）的解
记住：
满足x(0)=x₀的特解为： x(t)=e^tAx₀

2.非齐次微分方程（x^’(t)=A·x(t)+b(t)）的解
满足x(0)=x₀的特解为:
：x(t)=e^tA·[x₀+${\int }_0^t$e^-uAb(u)d_u]

例：

步骤：
1.微分方程的等式右边部分=A·(x1,x2,x3)，求出坐标矩阵A（a₁₁ a₂₁ … a_n1)^T,（a₁₂ a₂₂ … a_n2)^T…（a_1n a_2n … a_nn)^T
2.求出矩阵A最小特征多项式m(λ)（即，使得（A-λ₁E)^m…(A-λ_nE)^s=0阵）
3.f(λ)=m(λ)·g(λ)+r(λ)=e^λt(通解)，而 r(λ) = a+bλ+cλ²+…+gλ^n-1,根据最小多项式m(λ)的次数和=确定r(λ)的项数（比如m(λ）=（λ-1)(λ-2)中 m(λ）的次数和为2，因此对应的r(λ)为两项，也就是 r(λ)=a+bλ ）,求出f(λ₁)…f(λ_n)（如果解出来的特征值为重根，则求出f(λ)和f^’(λ)），并使用待定系数法解出r(λ)的系数 a,b,…,g
4.求出f(A)=e^At
5.求出特解

2⃣️证明题

1.证明T是线性变换，只需要证明:∀a，b∈子空间V₁，∀l₁,l₂∈R，有T(l₁a+l₂b) = l₁(Ta)+l₂(Tb)

2.证明T是正交变换，只需要证明：∀a∈欧氏空间V，(Ta,Ta)=(a,a)

3.证明T是对称变换，只需要证明：∀x,y∈欧氏空间V，(Tx,y)=(x,Ty)

3.和下面这个例子类似的证明：
其中题目已知条件：T₁T₂=T₂T₁
要证明：R(T₂)是T₁的不变子空间，即证明：∀x∈R(T₂)，有T₁x∈R(T₂)

要证明：N(T₂)是T₁的不变子空间，即证明：∀x∈N(T₂)，有T₁(x)=0∈N(T₂)
4.关于谱半径的证明
(1)证明：∀A_n*n，有谱半径ρ(A)<||A||_M

（2）证明：∀A_n*n，ε>0，∃范数||·||_M，使得||A||_M<ρ(A)+ε

3⃣️使用盖尔圆隔离特征值

步骤：
1. |z-a_ii|<=R_i，其中a_ii为主对角线元素，R_i是指第i 行除主对角线元素之外的所有元素的绝对值之和。
2.画出盖尔圆
3.构造隔离矩阵D = diag（a₁,a₂,…,a_n），构造的原则是，对原矩阵主对角线元素放大缩小。比如，第二个盖尔圆和第三个盖尔圆连通，第一个盖尔圆孤立，则应该将第一个盖尔圆变大，从而使得其他两个盖尔圆相对变小而不连通。将第一个盖尔圆变大的具体做法是将隔离矩阵D的主对角线的第一个元素a₁设为>1的数（具体为几，可自定义，能分开盖尔圆即可），而其他位置的主对角线元素设置为1。
4.新矩阵B = DAD^-1，并对B矩阵画出盖尔圆，发现盖尔圆不相交，特征值已隔离

上题答案：

4⃣️QR分解（Givens，Householder）

目的：将矩阵A分解为QR（Q为正交矩阵，R为可逆上三角矩阵），A = QR

1.householder 有限个H矩阵之积S，使得R=SA为可逆上三角矩阵，则有A=S^-1R=S^TR=QR
步骤：
(1)

1. 取出矩阵A的第一列β⁽⁰⁾
2. 求k = β⁽⁰⁾-|β⁽⁰⁾|e₁
3. 求u = k/||k||₂
4. 求H₀ = E - 2uu^T
5. 求H₀A
6. 取矩阵H₀A（nn）的第一行为R的第一行，剩下右下角的（n-1)(n-1）矩阵为A⁽¹⁾

(2)

8. 取出矩阵A⁽¹⁾的第一列β⁽¹⁾
   重复（1）的2-5步
9. 矩阵H₁A⁽¹⁾的上三角部分，和第（1）步的6小步中取出的第一行拼成一个上三角矩阵R
10. S = diag（1，H⁽¹⁾）H₀
11. Q = S^-1=S^T
12. H = QR

上题答案：

2.Givens
步骤:
（1）
1.取出矩阵A的第一列β⁽⁰⁾，β⁽⁰⁾=（v1,v2,v3,…,vn)^T
2.如果矩阵是4阶的，先计算T₁₄，T_ij 矩阵的特征为（注意看第i行和第j行，第i，j列，且i<j ，i=1，j从n开始取）,且c=vi/sqrt(vi²+vj²)，s=vj/sqrt(vi²+vj²)，且满足c²+s²=1,如果不满足则求T₁₃，直到有满足的c,s出现为止，然后求出T_ij
β⁽⁰⁾ (验证=(|β⁽⁰⁾|,0,0,…,0)^T)，令T₀=T_ij

3.求出T₀A
4.取矩阵T₀A（nn）的前（n-j+1）行为R的前n-j+1行，剩下右下角的（j-1)(j-1）矩阵为A⁽¹⁾

(2）
1.取出矩阵A⁽¹⁾的第一列β⁽¹⁾
重复（1）中的2步
2.计算T₁A⁽¹⁾(如果是2阶，则将其上三角部分和（1）步第4小步的R拼接起来)
3.而令T=diag(E,T₁)T₀
4.Q=T_-1=T^T
5.A=QR

例题：