文章目录
- 【一(18)】填空题
- 【二(10)】范数证明
- 【三(15)】矩阵函数
- 1 计算 e A t e^{At} eAt
- 2 求微分方程的解
- 【四(10)】QR分解
- 【五(10)】Gerschgorim隔离特征值
- 【六(15)】 A + A^+ A+计算及求解线性方程组
- 1 计算满秩分解
- 2 计算 A + A^+ A+
- 3 判断线性方程组解是否存在
- 【七(15)】线性空间与线性变换
- 求线性变换T在一组基下的矩阵
- 生成子空间的基与维数
- 线性变换矩阵的化简(相似对角化)
- 证明线性变换是对称变换
- 不变子空间证明
【一(18)】填空题
-
计算矩阵的Jordan标准形
- 能分块先分块( 但反求对应的广义特征向量不可分块);
- 分块对角阵,可以按块分别求Jordan标准形【11、13、16、17】;
- 2,3阶用特征向量法(P9),4阶及以上用Smith标准形法;
- 矩阵求特征多项式,2,3阶套公式,4阶以上利用初等变换(包括行、列变换,不改变行列式的值)化简计算;
- 计算相似变换矩阵P(P15-16)【14】。
-
计算向量/矩阵范数
-
向量范数
∥ x ∥ 1 = ∑ k = 1 n ∣ ξ k ∣ ∥ x ∥ 2 = ∑ k = 1 n ∣ ξ k ∣ 2 = x H x ∥ x ∥ ∞ = max k ∣ ξ k ∣ ∥ x ∥ p = ( ∑ k = 1 n ∣ ξ k ∣ p ) 1 p ∥ x ∥ A = x H A x , A 为 H e r m i t e 正定矩阵 . \|x\|_{1}=\sum_{k=1}^{n}|\xi_{k}|\\ \|x\|_{2}=\sqrt{\sum_{k=1}^{n}|\xi_{k}|^{2}}=\sqrt{x^{H}x}\\ \|x\|_{\infty}=\max_{k}|\xi_{k}|\\ \left.\|x\|_p=\left(\begin{array}{c}\sum_{k=1}^n|\xi_k|^p\\\end{array}\right.\right)^{\frac{1}{p}}\\ \|x\|_{A}=\sqrt{x^{H}Ax},\,A为Hermite正定矩阵. ∥x∥1=k=1∑n∣ξk∣∥x∥2=k=1∑n∣ξk∣2=xHx∥x∥∞=kmax∣ξk∣∥x∥p=(∑k=1n∣ξk∣p)p1∥x∥A=xHAx,A为Hermite正定矩阵. -
矩阵范数
∥ A ∥ m 1 = ∑ i = 1 m ∑ j = 1 n ∣ a i j ∣ ∥ A ∥ F = ∑ i = 1 n ∑ j = 1 n ∣ a i j ∣ 2 = t r ( A H A ) ∥ A ∥ m ∞ = n max i , j ∣ a i j ∣ ∥ A ∥ 1 = max j ∑ i = 1 n ∣ a i j ∣ ∥ A ∥ 2 = A H A 的最大特征值 ∥ A ∥ ∞ = max i ∑ j = 1 n ∣ a i j ∣ \|A\|_{m_{1}}=\sum_{i=1}^m\sum_{j=1}^n|a_{ij}|\\ \|A\|_F=\sqrt{\sum_{i=1}^n\sum_{j=1}^n\left|a_{ij}\right|^2}=\sqrt{tr(A^HA)}\\ \|A\|_{m_{\infty}}=n\max_{i,j}|a_{ij}|\\ \|A\|_{1}=\max_j\sum_{i=1}^n|a_{ij}|\\ \|A\|_{2}=\sqrt{A^{H}A\text{的最大特征值}}\\ \|A\|_\infty=\max_i\sum_{j=1}^n|a_{ij}| ∥A∥m1=i=1∑mj=1∑n∣aij∣∥A∥F=i=1∑nj=1∑n∣aij∣2=tr(AHA)∥A∥m∞=ni,jmax∣aij∣∥A∥1=jmaxi=1∑n∣aij∣∥A∥2=AHA的最大特征值∥A∥∞=imaxj=1∑n∣aij∣若A是正规矩阵(例如实对称矩阵), ∥ A ∥ 2 \|A\|_{2} ∥A∥2为 A A A的最大特征值。
1 + 2 + 3 + . . . + n = n ( n + 1 ) 2 1 2 + 2 2 + 3 2 + . . . + n 2 = n ( n + 1 ) ( 2 n + 1 ) 6 1+2+3+...+n=\frac{n(n+1)}2 \\1^2+2^2+3^2+...+n^2=\frac{n(n+1)(2n+1)}6 1+2+3+...+n=2n(n+1)12+22+32+...+n2=6n(n+1)(2n+1)
看清到底是向量范数还是矩阵范数!!!
-
-
计算矩阵的M-P逆
-
( A ⊗ B ) + = A + ⊗ B + (A\otimes{B})^+=A^+\otimes{B^+} (A⊗B)+=A+⊗B+。考点!!!;
-
A ( 1 ) A^{(1)} A(1)的计算【21】。
S A P = [ I r K O O ] , A ∈ C r m × n ( r > 0 ) , S ∈ C m m × m S 为变换矩阵, P 为 n 阶 ( 列 ) 置换矩阵。 SAP=\left[\begin{array}{ccc} I_r & K \\ O & O \end{array}\right],\,A\in{C_r^{m\times{n}}}(r>0),\,S\in{C_m^{m\times{m}}}\\S为变换矩阵,P为n阶(列)置换矩阵。 SAP=[IrOKO],A∈Crm×n(r>0),S∈Cmm×mS为变换矩阵,P为n阶(列)置换矩阵。
A ( 1 ) A^{(1)} A(1)为:
X = P [ I r O O L ] S , X ∈ C n × m , L ∈ C ( n − r ) × ( m − r ) X=P\left[\begin{array}{ccc} I_r & O \\ O & L \end{array}\right]S,\,X\in{C^{n\times{m}}},\,L\in{C^{(n-r)\times{(m-r)}}} X=P[IrOOL]S,X∈Cn×m,L∈C(n−r)×(m−r)
当上式中 L = O L=O L=O时, X X X为 A ( 1 , 2 ) A^{(1,2)} A(1,2)。
-
-
线性空间
-
矩阵的直积【11、14、15、16、20、21】
- 不会单独考直积的计算和性质,一般结合范数、广义逆出题;
- 正交投影矩阵 A A A, A 2 = A , A H = A A^2=A,A^H=A A2=A,AH=A, λ ( A ) = 0 o r 1 \lambda(A)=0\,or\,1 λ(A)=0or1;
- 投影矩阵 A A A, A 2 = A A^2=A A2=A;
- 酉矩阵A, A H A = I , A H = A − 1 A^HA=I,A^H=A^{-1} AHA=I,AH=A−1;
-
矩阵函数
-
已知 f ( A t ) f(At) f(At)求 A A A。【11、12】
已知 e A t e^{At} eAt求 A A A。
首先计算 d e A t d t = A e A t \frac{\text{d}e^{At}}{\text{d}t}=Ae^{At} dtdeAt=AeAt(对矩阵中每个函数分别求导),令 t = 0 t=0 t=0,即可求得 A A A。( sin A t \sin{At} sinAt同思路)
-
公式计算【15、17】
-
【二(10)】范数证明
-
Holder不等式
∑ k = 1 n ∣ ξ k ∣ ∣ η k ∣ ≤ ( ∑ k = 1 n ∣ ξ k ∣ p ) 1 p ( ∑ k = 1 n ∣ η k ∣ q ) 1 q \sum_{k=1}^n|\xi_k||\eta_k|\leq(\sum_{k=1}^n|\xi_k|^p)^{\frac{1}{p}}(\sum_{k=1}^n|\eta_k|^q)^{\frac{1}{q}} k=1∑n∣ξk∣∣ηk∣≤(k=1∑n∣ξk∣p)p1(k=1∑n∣ηk∣q)q1
上式中, p > 1 , q > 1 , 且 1 p + 1 q = 1. p>1,q>1,且\frac{1}{p}+\frac{1}{q}=1. p>1,q>1,且p1+q1=1.当 p = q = 2 p=q=2 p=q=2 时,即得到 Cauchy-Schwarz 不等式。
∑ k = 1 n ∣ ξ k ∣ ∣ η k ∣ ≤ ( ∑ k = 1 n ∣ ξ k ∣ 2 ) ( ∑ k = 1 n ∣ η k ∣ 2 ) \sum_{k=1}^n|\xi_k||\eta_k|\leq\sqrt{(\sum_{k=1}^n|\xi_k|^2)(\sum_{k=1}^n|\eta_k|^2)} k=1∑n∣ξk∣∣ηk∣≤(k=1∑n∣ξk∣2)(k=1∑n∣ηk∣2) -
C = A B , c i j = ∑ k = 1 n a i k b k j C=AB,\,c_{ij}=\sum_{k=1}^na_{ik}b_{kj} C=AB,cij=∑k=1naikbkj
-
对称正定矩阵 A A A满足:存在可逆矩阵 P P P使得 A = P H P A=P^HP A=PHP
-
几种常见的不等式
-
( ∑ k = 1 n ∣ ξ k ∣ 2 ) ( ∑ k = 1 n ∣ η k ∣ 2 ) ≤ ( ∑ k = 1 n ∣ ξ k ∣ ) ( ∑ k = 1 n ∣ η k ∣ ) \sqrt{(\sum_{k=1}^n|\xi_k|^2)(\sum_{k=1}^n|\eta_k|^2)}\leq(\sum_{k=1}^n|\xi_k|)(\sum_{k=1}^n|\eta_k|) (∑k=1n∣ξk∣2)(∑k=1n∣ηk∣2)≤(∑k=1n∣ξk∣)(∑k=1n∣ηk∣)
-
∣ ∑ k = 1 n a i k b k j ∣ ≤ ∑ k = 1 n ∣ a i k b k j ∣ = ∑ k = 1 n ∣ a i k ∣ ∣ b k j ∣ |\sum_{k=1}^na_{ik}b_{kj}|\leq\sum_{k=1}^n|a_{ik}b_{kj}|=\sum_{k=1}^n|a_{ik}||b_{kj}| ∣∑k=1naikbkj∣≤∑k=1n∣aikbkj∣=∑k=1n∣aik∣∣bkj∣
-
( a b + c d ) 2 ≤ ( a 2 + c 2 ) ( b 2 + d 2 ) (ab+cd)^2\leq(a^2+c^2)(b^2+d^2) (ab+cd)2≤(a2+c2)(b2+d2) 【柯西不等式】
在证明范数三角不等式或相容性时可以使用。
-
【三(15)】矩阵函数
1 计算 e A t e^{At} eAt
采用待定系数法求解,P71,例3.10。
Step1:求矩阵A的特征多项式 d e t ( λ I − A ) det(\lambda{I}-A) det(λI−A),并求出最小多项式 m A ( λ ) m_{A}(\lambda) mA(λ);
特征多项式计算:
( λ + a ) ( λ + b ) ( λ + c ) = λ 3 + ( a + b + c ) λ 2 + ( a b + a c + b c ) λ + a b c (\lambda+a)(\lambda+b)(\lambda+c)=\lambda^3+(a+b+c)\lambda^2+(ab+ac+bc)\lambda+abc (λ+a)(λ+b)(λ+c)=λ3+(a+b+c)λ2+(ab+ac+bc)λ+abc
最小多项式计算:试探法,P21。
Step2:根据最小多项式的最高次次数 n n n,设待定方程
f ( λ ) = e λ t = m A ( λ ) q ( λ ) + b n − 1 λ n − 1 + . . . + b 1 λ + b 0 f(\lambda)=e^{\lambda{t}}=m_{A}(\lambda)q(\lambda)+b_{n-1}\lambda^{n-1}+...+b_1\lambda+b_0 f(λ)=eλt=mA(λ)q(λ)+bn−1λn−1+...+b1λ+b0
接着,根据特征值及其重数(决定是否求导)列出方程组,求解得到系数 b 0 , b 1 , . . . , b n − 1 b_0,b_1,...,b_{n-1} b0,b1,...,bn−1;Step3:将 A A A代入上式计算 e A t e^{At} eAt。
e A t = b n − 1 A n − 1 + . . . + b 1 A + b 0 I e^{At}=b_{n-1}A^{n-1}+...+b_1A+b_0I eAt=bn−1An−1+...+b1A+b0I
本题题型只涉及3阶矩阵,且最小多项式为最高次数为2(要么两个不同特征值,要么一个2重根的特征值,后者需要求一阶导建立一个方程,注意是对 λ \lambda λ求导, t t t为参数),即待定系数只有2个。
2 求微分方程的解
Step1:计算 e − A τ b ( τ ) e^{-A\tau}b(\tau) e−Aτb(τ);
其中 e − A τ e^{-A\tau} e−Aτ便是将 e A τ e^{A\tau} eAτ矩阵中每个元素中的 τ \tau τ用 − τ -\tau −τ替换后的矩阵。
Step2:计算 ∫ 0 t e − A τ b ( τ ) d τ \int_0^t{e^{-A\tau}b(\tau)}d\tau ∫0te−Aτb(τ)dτ;
∫ a b A ( t ) d t = ( ∫ a b a i j ( t ) d t ) m × n \int_a^b{A(t)}dt=(\int_a^ba_{ij}(t)dt)_{m\times{n}} ∫abA(t)dt=(∫abaij(t)dt)m×n
Step3:计算 x ( t ) x(t) x(t)。
x ( t ) = e A t x ( 0 ) + e A t ∫ 0 t e − A τ b ( τ ) d τ x(t)=e^{At}x(0)+e^{At}\int_0^t{e^{-A\tau}b(\tau)}d\tau x(t)=eAtx(0)+eAt∫0te−Aτb(τ)dτ
历年考卷中一般本题积分项结果都很简单。
刷题总结
-
注意 λ I − A \lambda{I}-A λI−A计算时,除了将 A A A对角线元素用 λ \lambda λ减去外,其他元素应被 0 0 0减去,即取负。
-
计算 e A t e^{At} eAt时,常数项应当乘以一个单位矩阵,在计算两个矩阵元素对应相加时,注意只有对角线元素需要加上常数项值。
-
在计算过程中将 e e e的幂次项提出来,方便计算。
【四(10)】QR分解
这道题固定为使用Householder变换或Givens变换求解4阶矩阵A的QR分解。
矩阵乘法计算容易出错😢
Householder变换是镜像变换,不改变向量的长度。
【五(10)】Gerschgorim隔离特征值
本题题型:
利用盖尔圆隔离矩阵A(4阶,偶尔3阶)的特征值,并画图表示。
利用实矩阵特征值的性质改进所得结果。
实矩阵的特征值必为实数,于是A的特征值分别在区间……中。
证明矩阵可逆。
矩阵可逆的充要条件 ∣ A ∣ ≠ 0 |A|\neq0 ∣A∣=0,即证明0不是不是矩阵A的特征值(0不在盖尔圆内部)。
【六(15)】 A + A^+ A+计算及求解线性方程组
1 计算满秩分解
- Step1:初等行变换将 A A A矩阵变为Hermite标准形(行最简形) H H H;
- Step2:取 A A A的第 j 1 , j 2 , . . . j r j_1,j_2,...j_r j1,j2,...jr( j i j_i ji为第 i i i行的第一个非零元素 1 1 1位于的列)列构成矩阵 F F F(列满秩),取 H H H的前 r r r行构成矩阵 G G G(行满秩)。
2 计算 A + A^+ A+
A + = G H ( G G H ) − 1 ( F H F ) − 1 F H A^+=G^H(GG^H)^{-1}(F^HF)^{-1}F^H A+=GH(GGH)−1(FHF)−1FH
这道题也可能是已知 A + A^+ A+求 A A A,利用性质 ( A + ) + = A (A^+)^+=A (A+)+=A即可计算。
3 判断线性方程组解是否存在
线性方程组有解的充要条件为:
A
A
+
b
=
b
AA^+b=b
AA+b=b
极小范数解(有解)或极小范数最小二乘解(无解)为:
x
=
A
+
b
x=A^+b
x=A+b
先算 A + b A^+b A+b的值,在算 A A + b AA^+b AA+b。
【七(15)】线性空间与线性变换
求线性变换T在一组基下的矩阵
T ( α 1 , α 2 , . . . , α n ) ≜ ( T ( α 1 ) , T ( α 2 ) , . . . , T ( α n ) ) = ( α 1 , α 2 , . . . , α n ) A T(\alpha_1,\alpha_2,...,\alpha_n)\triangleq(T(\alpha_1),T(\alpha_2),...,T(\alpha_n))=(\alpha_1,\alpha_2,...,\alpha_n)A T(α1,α2,...,αn)≜(T(α1),T(α2),...,T(αn))=(α1,α2,...,αn)A
同一个线性变换在不同基下的矩阵关系:
T
(
α
1
,
α
2
,
.
.
.
,
α
n
)
=
(
α
1
,
α
2
,
.
.
.
,
α
n
)
A
,
T
(
β
1
,
β
2
,
.
.
.
,
β
n
)
=
(
β
1
,
β
2
,
.
.
.
,
β
n
)
B
,
(
β
1
,
β
2
,
.
.
.
,
β
n
)
=
(
α
1
,
α
2
,
.
.
.
,
α
n
)
P
,
⇒
T
(
β
1
,
β
2
,
.
.
.
,
β
n
)
=
T
[
(
α
1
,
α
2
,
.
.
.
,
α
n
)
P
]
=
T
(
α
1
,
α
2
,
.
.
.
,
α
n
)
P
=
(
α
1
,
α
2
,
.
.
.
,
α
n
)
A
P
=
(
β
1
,
β
2
,
.
.
.
,
β
n
)
P
−
1
A
P
T(\alpha_1,\alpha_2,...,\alpha_n)=(\alpha_1,\alpha_2,...,\alpha_n)A,\\ T(\beta_1,\beta_2,...,\beta_n)=(\beta_1,\beta_2,...,\beta_n)B,\\ (\beta_1,\beta_2,...,\beta_n)=(\alpha_1,\alpha_2,...,\alpha_n)P,\\ \rArr T(\beta_1,\beta_2,...,\beta_n)=T[(\alpha_1,\alpha_2,...,\alpha_n)P]=T(\alpha_1,\alpha_2,...,\alpha_n)P\\=(\alpha_1,\alpha_2,...,\alpha_n)AP=(\beta_1,\beta_2,...,\beta_n)P^{-1}AP
T(α1,α2,...,αn)=(α1,α2,...,αn)A,T(β1,β2,...,βn)=(β1,β2,...,βn)B,(β1,β2,...,βn)=(α1,α2,...,αn)P,⇒T(β1,β2,...,βn)=T[(α1,α2,...,αn)P]=T(α1,α2,...,αn)P=(α1,α2,...,αn)AP=(β1,β2,...,βn)P−1AP
生成子空间的基与维数
某组基(自然基)下坐标向量组的极大无关组与秩 ⇔ \lrArr ⇔ 元素组的极大无关组和秩 ⇔ \lrArr ⇔ 生成子空间的基与维数。
线性变换矩阵的化简(相似对角化)
证明线性变换是对称变换
即证明线性变换在线性空间的一组基下的矩阵为对称矩阵。