Linear algebra
A fun question
In n-D space we can found at most n + 1 n+1 n + 1 vector v 1 … v n + 1 v_1\ldots v_{n+1} v 1 … v n + 1 such that: ∀ i ≠ j , v i v j < 0 \forall i\ne j,v_iv_j<0 ∀ i = j , v i v j < 0
An example in 3-D space(the one on the book)
The construction is Obviously(C H 4 CH_4 C H 4 )
Chosen a vector, the others must be in a semisphere. We say two semisphere is seperated by plane A.
And we notice that: if two vector's shadow on plane A construct a acute angle, their dot product must be positive, transformed it into a 2-D problems which is easy to solve.
Generalize to n-D
We choose a vector v 1 v_1 v 1 , and it can be written as [ 1 , 0 , … 0 ] [1,0,\ldots 0] [ 1 , 0 , … 0 ] (with some rotation)
so ∀ v i , i > 1 , v i v 1 < 0 ⟹ v i , 1 < 0 ⟹ v i , 1 v j , 1 > 0 ⟹ v i v j − v i , 1 v j , 1 < 0 \forall v_i,i>1, v_iv_1<0 \implies v_{i,1} < 0 \implies v_{i,1}v_{j,1}>0 \implies v_iv_j-v_{i,1}v_{j,1} < 0 ∀ v i , i > 1 , v i v 1 < 0 ⟹ v i , 1 < 0 ⟹ v i , 1 v j , 1 > 0 ⟹ v i v j − v i , 1 v j , 1 < 0
then transformed it into the (n-1)-D situation.
The construction can be easily give during the induction
Another Proof
given by Bing!
Lemma: Radon Partition
In n-D space, n + 2 n+2 n + 2 point could be divided to two convex hull with intersection.
n + 2 n+2 n + 2 vector must be dependent: ∃ c i s . t . ∑ i = 1 n + 2 c i x i = 0 ; ∑ i x i = 0 \exists c_i \ s.t.\ \sum _{i=1}^{n+2} c_i x_i = 0; \sum_i x_i = 0 ∃ c i s . t . ∑ i = 1 n + 2 c i x i = 0 ; ∑ i x i = 0 (the second condition can be satisfied by add another all-1 dimension).
so divide the vector by sign of c i c_i c i we got:v = ∑ i ∈ A c i x i = ∑ j ∈ B c j x j v=\sum_{i\in A} c_i x_i = \sum_{j\in B} c_jx_j v = ∑ i ∈ A c i x i = ∑ j ∈ B c j x j ,so divide the eqution by ∑ i ∈ A c i \sum_{i\in A} c_i ∑ i ∈ A c i , you get one point(v v v ) in the intersection.
We noticed 0 < v 2 = ( ∑ i ∈ A c i x i ) ⋅ ( ∑ j ∈ B c j x j ) < 0 0<v^2=(\sum_{i\in A} c_i x_i) \cdot (\sum_{j\in B} c_jx_j)<0 0 < v 2 = ( ∑ i ∈ A c i x i ) ⋅ ( ∑ j ∈ B c j x j ) < 0 , contradiction!
Operator's Left Inverse and Right Inverse
Existance
Proof 1
注意到左逆存在等价于T T T 是单射,右逆存在等价于T T T 是满射.
又因为T T T 单射等价于T T T 满射等价于T T T 双射,所以左逆存在等价于右逆存在且两个逆一定相同.
Proof 2
左逆推出T T T 分解成初等行变换矩阵,然后一个一个取逆推右逆.
Union of finite count of proper sub-space isn't V
F is inifinite Field ∀ U 1 … U k , U i is subspace of V ( F ) V ≠ ⋃ i = 1 k U i \begin{gathered}
F \text{ is inifinite Field} \\
\forall U_1\ldots U_k,U_i \text{ is subspace of } V(F) \\
V\ne \bigcup_{i=1}^k U_i
\end{gathered} F is inifinite Field ∀ U 1 … U k , U i is subspace of V ( F ) V = i = 1 ⋃ k U i
考虑归纳,k = 1 k=1 k = 1 成立,设k − 1 k-1 k − 1 个不行,反证,那么可以取v ∉ ⋃ i = 1 k − 1 U i v\notin \bigcup_{i=1}^{k-1} U_i v ∈ / ⋃ i = 1 k − 1 U i ,必有v ∈ U k v\in U_k v ∈ U k .
再取u ∉ U k u\notin U_k u ∈ / U k ,令 S = { u + i v ∣ i ∈ R } S=\{ u+iv \vert i \in R\} S = { u + i v ∣ i ∈ R } ,则对j < k j<k j < k ,每个U j U_j U j 至多包含一个S S S 中的向量(否则v ∈ U j v\in U_j v ∈ U j ),而U k U_k U k 必然没有S S S 中向量(否则u ∈ U k u\in U_k u ∈ U k ),而S S S 中有无限个向量,于是S S S 不可能被他们包含,得证.
[think] 而这个定理甚至在有限域下有反例.
CR分解
Row Reduced Echelon Form
每一行的首个非零元素是 1,这个元素称为“主元”(pivot)。
每个主元所在列的其他元素都是 0,也就是说主元是该列唯一的非零元素。
每个主元都在其所在行的右边位置,相对于上一行的主元。
所有零行(即整行都是 0)都排在非零行的下面
A A A 的Row Reduced Echelon Form记为rref ( A ) \operatorname{rref}(A) rref ( A )
设B = rref ( A ) B=\operatorname{rref}(A) B = rref ( A ) 的主元所在列构成集合S = i ∣ a pivot is in i S={i\vert \text{a pivot is in } i} S = i ∣ a pivot is in i ,设A = [ a 1 … a n ] A=[a_1\ldots a_n] A = [ a 1 … a n ] ,则C = [ a i ∣ i ∈ S ] C=[a_i \vert i \in S] C = [ a i ∣ i ∈ S ] ,R = B 1 rank A , 1 m R=B_{1~\operatorname{rank}A,1~m} R = B 1 rank A , 1 m (即去除所有全0 0 0 行),满足A = C R A=CR A = C R .
我们把它看成用R R R 组装C C C 的列,那么进行初等行变换不改变列之间的线性关系.
而消元之后呢,看列的话主元所在列显然是标准基,那命题是显然的了.
秩分解
A n × m = P n × n B n × m Q m × m A_{n\times m}=P_{n\times n}B_{n\times m}Q_{m\times m} A n × m = P n × n B n × m Q m × m ,其中B B B 为只有左上角是一个 rank A × rank A \operatorname{rank} A\times \operatorname{rank} A rank A × rank A 的单位矩阵其他位置全是0 0 0 .
A A A 做行变换+列变换消元易得.
能不能换个视角,这个是不是在说,对 T ∈ L ( V , W ) T\in \mathcal L( V , W ) T ∈ L ( V , W ) ,存在一个V V V 的一个基v 1 … v n v_1\ldots v_n v 1 … v n ,W W W 的一个基w 1 … w m w_1\ldots w_m w 1 … w m 使得 T v i = [ i ≤ rank T ] w i Tv_i=[i\le \operatorname{rank} T]w_i T v i = [ i ≤ rank T ] w i .
那么先构造v v v ,我们先找一个 null A \operatorname{null} A null A 的基v n − r + 1 … v n v_{n-r+1}\ldots v_n v n − r + 1 … v n ,然后再任意扩充出剩下的v 1 … v n v_1\ldots v_n v 1 … v n .
对w w w ,显然T v 1 … T v r Tv_1\ldots Tv_r T v 1 … T v r 线性无关,再扩充w r + 1 … w m w_{r+1}\ldots w_m w r + 1 … w m 得到一组基.
显然这组基满足需求.
[think] 还是对矩阵基变换理解不到位.
注意我们把上面那个再搞一搞:P P P 的后m − r m-r m − r 列是没用的,Q Q Q 的后m − r m-r m − r 行是没用的,都丢到会得到P = C R P=CR P = C R .
而 P = C R = [ c 1 … c r ] [ r 1 T … r n T ] T = ∑ i = 1 n c i r i T P=CR=[c_1\ldots c_r] [r_1^T\ldots r_n^T]^T=\sum _{i = 1} ^{n} c_ir_i^T P = C R = [ c 1 … c r ] [ r 1 T … r n T ] T = ∑ i = 1 n c i r i T ,其中每个c i r i T c_ir_i^T c i r i T 秩为1 1 1 .这就是秩分解的名字.
Ax=B有解
∃ x , A x = b ⟺ rank A = rank [ A , b ] ⟺ b ∈ range A \begin{gathered}
\exists x,Ax=b \\
\iff \operatorname{rank} A=\operatorname{rank} [A,b] \\
\iff b\in \operatorname{range} A
\end{gathered} ∃ x , A x = b ⟺ rank A = rank [ A , b ] ⟺ b ∈ range A
解唯一等价于null A = 0 \operatorname{null} A=0 null A = 0 等价于n = rank A n=\operatorname{rank} A n = rank A
A quiz problem
A is a real matrix , A T A u = 0 ⟹ A u = 0 \begin{gathered}
A \text{ is a real matrix} , \\
A^TAu=0 \implies Au=0
\end{gathered} A is a real matrix , A T A u = 0 ⟹ A u = 0
A = M ( T ) T ∗ T u = 0 ⟺ ∀ v , < v , T ∗ T u > = 0 ⟺ ∀ v , < T v , T u > = 0 ⟺ T u ∈ ( range T ) ⊥ ∵ T u ∈ range T ∴ < T u , T u > = 0 , u = 0 \begin{gathered}
A=\mathcal M( T ) \\
T^*Tu=0 \\
\iff \forall v,<v,T^*Tu>=0 \\
\iff \forall v,<Tv,Tu>=0 \\
\iff Tu\in (\operatorname{range} T)^\perp \\
\because Tu\in \operatorname{range} T \\
\therefore <Tu,Tu>=0,u=0
\end{gathered} A = M ( T ) T ∗ T u = 0 ⟺ ∀ v , < v , T ∗ T u >= 0 ⟺ ∀ v , < T v , T u >= 0 ⟺ T u ∈ ( range T ) ⊥ ∵ T u ∈ range T ∴ < T u , T u >= 0 , u = 0
[think] 被这个题击败了,当时只想到用 null T ∗ = ( range T ) ⊥ \operatorname{null} T^*=(\operatorname{range} T)^{\perp} null T ∗ = ( range T ) ⊥ 了,但其实是可以简单翻译过来的.
伴随和共轭转置的关系其实是显然的,内积上伴随的性质也是显然的,所以基础操作没必要用结论.做题的时候错误的感觉算子伴随和矩阵转置的距离过远(因为done right中证明是表示成规范正交基然后拆开用内积的性质,但是不看那套框架的话其实是显然的,另外对U ⊕ U ⊥ = V U\oplus U^\perp=V U ⊕ U ⊥ = V 的证明掌握不好).同时左零空间.
总结就是记住了几何那边的结论但没有很好的联系到代数这边.
Several Inequations about Rank
rank A + B ≤ rank A + rank B rank A B ≤ min rank A , rank B A m × n B n × s = 0 ⟹ rank A + rank B ≤ n \begin{gathered}
\operatorname{rank} A+B \le \operatorname{rank} A+\operatorname{rank} B \\
\operatorname{rank} AB \le \min \operatorname{rank} A,\operatorname{rank} B \\
A_{m\times n}B_{n\times s}=0 \implies \operatorname{rank} A+\operatorname{rank} B\le n
\end{gathered} rank A + B ≤ rank A + rank B rank A B ≤ min rank A , rank B A m × n B n × s = 0 ⟹ rank A + rank B ≤ n
rank A B ≥ rank A m × n + rank B n × s − n \begin{gathered}
\operatorname{rank} AB\ge \operatorname{rank} A_{m\times n}+\operatorname{rank} B_{n\times s}-n
\end{gathered} rank A B ≥ rank A m × n + rank B n × s − n
Sol 1
矩阵分解:
A = P 1 [ I r 1 , 0 0 , 0 ] Q 1 B = P 2 [ I r 2 , 0 0 , 0 ] Q 2 A B = P 1 [ I r 1 , 0 0 , 0 ] Q 1 P 2 [ I r 2 , 0 0 , 0 ] Q 2 \begin{gathered}
A=P_1 \begin{bmatrix}
I_{r_1},0 \\
0,0
\end{bmatrix}Q_1 \\
B=P_2 \begin{bmatrix}
I_{r_2},0 \\
0,0
\end{bmatrix}Q_2 \\
AB=P_1\begin{bmatrix}
I_{r_1},0 \\
0,0
\end{bmatrix}Q_1P_2\begin{bmatrix}
I_{r_2},0 \\
0,0
\end{bmatrix}Q_2
\end{gathered} A = P 1 [ I r 1 , 0 0 , 0 ] Q 1 B = P 2 [ I r 2 , 0 0 , 0 ] Q 2 A B = P 1 [ I r 1 , 0 0 , 0 ] Q 1 P 2 [ I r 2 , 0 0 , 0 ] Q 2
显然P 1 , Q 2 P_1,Q_2 P 1 , Q 2 不影响最终的秩直接扔了,而设D = Q 1 P 2 = [ D 1 , D 2 D 3 , D 4 ] D=Q_1P_2=\begin{bmatrix}
D_1,D_2 \\
D_3,D_4
\end{bmatrix} D = Q 1 P 2 = [ D 1 , D 2 D 3 , D 4 ] ,那么你发现乘完只剩下D 1 D_1 D 1 .
而删去矩阵一行或一列秩最多减少1 1 1 ,D 1 D_1 D 1 看成D D D 删掉了 n − rank A + n − rank B n-\operatorname{rank} A + n-\operatorname{rank} B n − rank A + n − rank B 行或列得到的.同时 rank D = n \operatorname{rank} D=n rank D = n ,得证.
Sol 2
考虑
C = [ I n , 0 0 , A B ] \begin{gathered}
C=\begin{bmatrix}
I_n,0 \\
0,AB
\end{bmatrix}
\end{gathered} C = [ I n , 0 0 , A B ]
显然 rank A B + n = rank C \operatorname{rank} AB+n=\operatorname{rank} C rank A B + n = rank C
对它做行变换可以得到
C → [ I n , 0 A , A B ] → [ I n , − B A , 0 ] = D \begin{gathered}
C \to \begin{bmatrix}
I_n,0 \\
A,AB
\end{bmatrix} \\
\to \begin{bmatrix}
I_n,-B \\
A,0
\end{bmatrix}=D \\
\end{gathered} C → [ I n , 0 A , A B ] → [ I n , − B A , 0 ] = D
而观察这个D D D 容易发现 rank D ≥ rank A + rank B \operatorname{rank} D\ge \operatorname{rank} A+\operatorname{rank} B rank D ≥ rank A + rank B ,于是得证
Sol 3
考虑我们要证明 dim range A B ≥ dim range B − dim null A \dim \operatorname{range} AB\ge \dim \operatorname{range} B-\dim \operatorname{null} A dim range A B ≥ dim range B − dim null A
考虑为什么 range B ≠ range A B \operatorname{range} B\ne \operatorname{range} AB range B = range A B ,是因为 range B \operatorname{range} B range B 中的不同元素被合成了一个,而这个合成相当于把 差是 null A \operatorname{null} A null A 中的元素的多个元素合成一个.所以有
dim ( range B ) / ( range B ∩ null A ) = dim range A \begin{gathered}
\dim (\operatorname{range} B)/(\operatorname{range} B\cap \operatorname{null} A)=\dim \operatorname{range} A
\end{gathered} dim ( range B ) / ( range B ∩ null A ) = dim range A
显然交集小于 null A \operatorname{null} A null A ,得证.
rank A C + rank C B ≤ rank C + rank A C B \begin{gathered}
\operatorname{rank} AC+\operatorname{rank} CB\le \operatorname{rank} C+\operatorname{rank} ACB
\end{gathered} rank A C + rank C B ≤ rank C + rank A C B
这个结论可以直接由上一个的Sol3弄出来,考虑
[ C , 0 0 , A C B ] \begin{gathered}
\begin{bmatrix} C,0 \\
0,ACB \end{bmatrix}
\end{gathered} [ C , 0 0 , A C B ]
可以简单消元变成
[ C , C B A C , 0 ] \begin{gathered}
\begin{bmatrix} C,CB \\
AC,0 \end{bmatrix}
\end{gathered} [ C , C B A C , 0 ]
于是直接得证.
[think] 学会这种拼成空间再分块矩阵消元的套路.
A 2 = I ⟹ rank ( A − I ) + rank ( A + I ) = n \begin{gathered}
A^2=I \\
\implies \operatorname{rank} (A-I)+\operatorname{rank} (A+I)=n
\end{gathered} A 2 = I ⟹ rank ( A − I ) + rank ( A + I ) = n
( A − I ) ( A + I ) = 0 ⟹ rank ( A − I ) + rank A + I ≤ n ( A + I ) − ( A − I ) = 2 I ⟹ rank ( A + I ) + rank ( A − I ) ≥ n Q.E.D \begin{gathered}
(A-I)(A+I)=0 \\
\implies \operatorname{rank} (A-I)+\operatorname{rank} A+I \le n \\
(A+I) - (A-I)=2I \\
\implies \operatorname{rank} (A+I)+\operatorname{rank} (A-I)\ge n
\\
\text{Q.E.D}
\end{gathered} ( A − I ) ( A + I ) = 0 ⟹ rank ( A − I ) + rank A + I ≤ n ( A + I ) − ( A − I ) = 2 I ⟹ rank ( A + I ) + rank ( A − I ) ≥ n Q.E.D
对平面图Graph ( n , m ) \text{Graph}(n,m) Graph ( n , m ) 有F F F 个面(不含最外面),证明
首先考虑无向图的 Incidence Matrix M M M ,容易注意到M M M 中的若干行线性无关等价于这个导出子图无环.
于是看出 rank M = n − c \operatorname{rank} M=n-c rank M = n − c ,c c c 为连通块个数.
又能看出 v ∈ null M T v\in \operatorname{null} M^T v ∈ null M T 等价于v v v 中的若干条边串成若干个环,会发现 dim null M T = F \dim \operatorname{null} M^T=F dim null M T = F
An Ex Problem
M = [ A , C 0 , B ] rank M = rank A + rank B ⟺ ∃ X , Y : A X + Y B = C \begin{gathered}
M=\begin{bmatrix} A,C \\0,B \end{bmatrix} \\
\operatorname{rank} M=\operatorname{rank} A+\operatorname{rank} B \iff \exists X,Y:AX+YB=C
\end{gathered} M = [ A , C 0 , B ] rank M = rank A + rank B ⟺ ∃ X , Y : A X + Y B = C
首先右推左是显然的.直接消元一下就好了.考虑左推右.
考虑
A : U 1 → V 1 , B : U 2 → V 2 , C : U 2 → V 1 M : U → V \begin{gathered}
A:U_1\to V_1,B:U_2\to V_2,C:U_2\to V_1 \\
M:U\to V
\end{gathered} A : U 1 → V 1 , B : U 2 → V 2 , C : U 2 → V 1 M : U → V
分解
range M = range B ⊕ W W = { [ v , 0 ] ∣ [ v , 0 ] ∈ range M } \begin{gathered}
\operatorname{range} M=\operatorname{range} B\oplus W \\
W=\{ [v,0] \vert [v,0]\in \operatorname{range} M \} \\
\end{gathered} range M = range B ⊕ W W = {[ v , 0 ] ∣ [ v , 0 ] ∈ range M }
那么因为 [ v , 0 ] ∈ range M [v,0]\in \operatorname{range} M [ v , 0 ] ∈ range M ,则必然是 [ v , 0 ] = M [ u 1 , u 2 ] T [v,0]=M[u_1,u_2]^T [ v , 0 ] = M [ u 1 , u 2 ] T ,一定是B u 2 = 0 , A u 1 + C u 2 = v Bu_2=0,Au_1+Cu_2=v B u 2 = 0 , A u 1 + C u 2 = v .
于是 W = range A + C ( null B ) W=\operatorname{range} A + C(\operatorname{null} B) W = range A + C ( null B ) .又因为 dim W = dim range A \dim W=\dim \operatorname{range} A dim W = dim range A ,于是有 C ( null B ) ⊂ range A C(\operatorname{null} B) \subset \operatorname{range} A C ( null B ) ⊂ range A .
我们再分解 U 2 = null B ⊕ U 3 U_2=\operatorname{null} B \oplus U_3 U 2 = null B ⊕ U 3 ,此时注意到B B B 在U 3 U_3 U 3 到 range B \operatorname{range} B range B 是双射,存在Y ′ , ∀ u , Y B u = C u Y',\forall u,YBu=Cu Y ′ , ∀ u , Y B u = C u .然后通过扩充基并任意取值将Y ′ Y' Y ′ 的定义域扩充到V 2 V_2 V 2 得到Y Y Y .
于是( C − Y B ) u (C-YB)u ( C − Y B ) u 对任意u ∈ U 3 u\in U_3 u ∈ U 3 为0 0 0 ,于是 range ( C − Y B ) = ( C − Y B ) null B = C ( null B ) \operatorname{range} (C-YB)=(C-YB)\operatorname{null} B= C(\operatorname{null} B) range ( C − Y B ) = ( C − Y B ) null B = C ( null B ) .
现在只考虑 u ∈ null B u\in \operatorname{null} B u ∈ null B ,显然∃ v , C u = A v \exists v,Cu=Av ∃ v , C u = A v ,那么对 null B \operatorname{null} B null B 的一组基u 1 … u k u_1\ldots u_k u 1 … u k 这样确定v 1 … v k v_1\ldots v_k v 1 … v k ,就可以构造X ′ u i = v i X'u_i=v_i X ′ u i = v i 满足C u = A X ′ v Cu=AX'v C u = A X ′ v .再用同样的方法扩充基并任意取值将X ′ X' X ′ 的定义域扩充到U 2 U_2 U 2 得到X X X .
于是C = Y B + A X C=YB+AX C = Y B + A X
[think] 感觉得到 C ( null B ) = A C(\operatorname{null} B)=A C ( null B ) = A 这里是容易的.然后这里进行不下去,想到 null B \operatorname{null} B null B 去分解也是自然的. 分解后就要想办法把 null B \operatorname{null} B null B 之外的影响消掉,就用了C − Y B C-YB C − Y B .而若 range A ⊂ range B \operatorname{range} A\subset \operatorname{range} B range A ⊂ range B 那么A x = B T x Ax=BTx A x = B T x 是显然的.
投影
向一个向量投影
a , b ∈ R n < a , b > < a , a > a = ( a T b a a T a ) = ( a ⋅ a T a T a ) b \begin{gathered}
a,b\in R^n \\
\dfrac{<a,b>}{<a,a>} a = (\dfrac{a^Tba}{a^T a})= (\dfrac{a\cdot a^T}{a^T a}) b
\end{gathered} a , b ∈ R n < a , a > < a , b > a = ( a T a a T ba ) = ( a T a a ⋅ a T ) b
向一个平面投影(平面是C ( A ) C(A) C ( A ) )
p = A ( A T A ) − 1 A T b \begin{gathered}
p=A(A^TA)^{-1}A^Tb
\end{gathered} p = A ( A T A ) − 1 A T b
考虑b b b 的投影p ∈ C ( A ) p\in C(A) p ∈ C ( A ) 有( b − p ) ∈ C ( A ) ⊥ (b-p)\in C(A)^\perp ( b − p ) ∈ C ( A ) ⊥ ,于是b − p ∈ N ( A T ) b-p \in N(A^T) b − p ∈ N ( A T ) .
于是A T b = A T p A^Tb=A^Tp A T b = A T p ,又p ∈ C ( A ) ⟹ ∃ x , A x = p p\in C(A) \implies \exists x,Ax=p p ∈ C ( A ) ⟹ ∃ x , A x = p .
于是A T b = A T A x A^Tb=A^TAx A T b = A T A x ,A A A 一定可以用一个满秩的,于是除过去.于是得证.
A T A x = A T b \begin{gathered}
A^TAx=A^Tb
\end{gathered} A T A x = A T b
一定有解
rank ( A T A , A T b ) = rank ( A T ( A , b ) ) ≤ rank A = rank A A T \begin{gathered}
\operatorname{rank} (A^TA,A^Tb) \\
=\operatorname{rank} (A^T(A,b))\le \operatorname{rank} A \\
=\operatorname{rank} AA^T
\end{gathered} rank ( A T A , A T b ) = rank ( A T ( A , b )) ≤ rank A = rank A A T
所以这个证明是依赖实数的.
[think] 复数你应该把A T A^T A T 换成A T ‾ \overline{A^T} A T ,或者说这个定理本来就应该是A T ‾ \overline{A^T} A T 的.
P 2 = P , P ∗ = P ⟹ P is a projection \begin{gathered}
P^2=P,P^*=P \implies P \text{ is a projection}
\end{gathered} P 2 = P , P ∗ = P ⟹ P is a projection
显然那P P P 只能说 range P \operatorname{range} P range P 的投影.
只需证明 ∀ u , ( P u − u ) ∈ ( range P ) ⊥ \forall u,(Pu-u)\in (\operatorname{range} P)^\perp ∀ u , ( P u − u ) ∈ ( range P ) ⊥ .
∀ v , < P u − u , P v > = < P P u − P u , v > = < P u − P u , v > = < 0 , v > ⟹ P u − u ∈ ( range P ) ⊥ \begin{gathered}
\forall v,<Pu-u,Pv>=<PPu-Pu,v>=<Pu-Pu,v>=<0,v> \\
\implies Pu-u\in (\operatorname{range} P)^\perp
\end{gathered} ∀ v , < P u − u , P v >=< P P u − P u , v >=< P u − P u , v >=< 0 , v > ⟹ P u − u ∈ ( range P ) ⊥
得证!
QR分解
∀ A , ∃ Q is orthogonal matrix , R is upper triangle matrix s . t . A = Q R \begin{gathered}
\forall A,\exists Q \text{ is orthogonal matrix},R \text{ is upper triangle matrix} \\ s.t.\\
A=QR
\end{gathered} ∀ A , ∃ Q is orthogonal matrix , R is upper triangle matrix s . t . A = QR
考虑A A A 可以看成把标准基变成a 1 … a n a_1\ldots a_n a 1 … a n ,那我们把a 1 … a n a_1\ldots a_n a 1 … a n 这组基用Gram-Schmidt变成b 1 … b n b_1\ldots b_n b 1 … b n ,问题就可以变成先把标准基变成b b b ,再变成a a a ,其中第一步是等距同构,第二步中我们知道a 1 … a i a_1\ldots a_i a 1 … a i 和b 1 … b i b_1\ldots b_i b 1 … b i 张成空间相同,所以第二步是上三角.
determinance
det A det B = det A B \begin{gathered}
\det A\det B=\det AB
\end{gathered} det A det B = det A B
Sol1:分块矩阵.
Sol2:都拆成初等变换矩阵再乘.
Sol3:考虑定义函数 α ( B ) = det B A det A \alpha(B)=\dfrac{\det BA}{\det A} α ( B ) = det A det B A ,容易验证它满足行列式三条公里,于是α ( B ) = det B \alpha(B)=\det B α ( B ) = det B
Laplace Theorem
∀ S ⊂ [ 1 , n ] ∩ Z det A = ∑ T ⊂ [ 1 , n ] ∩ Z , ∣ T ∣ = ∣ S ∣ A ( S , T ) C ( S , T ) \begin{gathered}
\forall S\subset [1,n]\cap Z \\
\det A=\sum _{T\subset [1,n]\cap Z,\vert T \vert =\vert S \vert } A(S,T)C(S,T)
\end{gathered} ∀ S ⊂ [ 1 , n ] ∩ Z det A = T ⊂ [ 1 , n ] ∩ Z , ∣ T ∣ = ∣ S ∣ ∑ A ( S , T ) C ( S , T )
其中A ( S , T ) A(S,T) A ( S , T ) 表示子式,C ( S , T ) C(S,T) C ( S , T ) 表示代数余子式
注意到你就是把S S S 行对应的元素钦定的时候的某个组合,我们可以先用∑ t ∈ T t − ∑ i = 1 ∣ T ∣ i + ∑ s ∈ S s − ∑ i = 1 ∣ S ∣ i ≡ ∑ t ∈ T t + ∑ s ∈ S s ( m o d 2 ) \sum_{t\in T} t-\sum_{i=1}^{\vert T\vert}i+\sum_{s\in S} s-\sum_{i=1}^{\vert S\vert}i\equiv \sum_{t\in T}t+\sum_{s\in S}s \pmod 2 ∑ t ∈ T t − ∑ i = 1 ∣ T ∣ i + ∑ s ∈ S s − ∑ i = 1 ∣ S ∣ i ≡ ∑ t ∈ T t + ∑ s ∈ S s ( mod 2 ) 次交换把这些行列顺序不变的换到前 ∣ S ∣ = ∣ T ∣ \vert S \vert=\vert T \vert ∣ S ∣ = ∣ T ∣ 行列,则最终符号显然就是此时的符号(即两个行列式内部的符号)再乘上交换操作的符号,于是得证.
求
det M M i , j = { a , i = j b , i ≠ j \begin{gathered}
\det M \\
M_{i,j} = \begin{cases}
a,i=j \\
b,i\ne j
\end{cases}
\end{gathered} det M M i , j = { a , i = j b , i = j
Sol1
考虑消元,注意到每行的和是一样的,所以我们把第一列变成所有列的和,然后用它去消后面的,模拟一下即可.
Sol2
考虑M = ( a − b ) I + x B M=(a-b)I+xB M = ( a − b ) I + x B ,其中B = ( b ) i , j B=(b)_{i,j} B = ( b ) i , j ,应用公式
d det ( M ) = trace ( adj A d A ) \begin{gathered}
\operatorname{d} \det(M)=\operatorname{trace} (\operatorname{adj} A \operatorname{d} A)
\end{gathered} d det ( M ) = trace ( adj A d A )
右边三项全是好算的,于是可以把导数对x x x 积分积回去.
Sol3
考虑M = B + ( a − b ) I M=B+(a-b)I M = B + ( a − b ) I ,若B v = λ v Bv=\lambda v B v = λ v ,则显然有( B + ( a − b ) I ) v = ( λ + a − b ) v (B+(a-b)I)v=(\lambda+a-b)v ( B + ( a − b ) I ) v = ( λ + a − b ) v ,于是求出B B B 的特征值,就可以直接得到M M M 的特征值算行列式.
rank A = a ⟺ { ∣ T ∣ = ∀ ∣ S ∣ > a , det A S , T = 0 ∃ ∣ S ∣ = ∣ T ∣ = a , det A S , T ≠ 0 \begin{gathered}
\operatorname{rank} A=a \\
\iff \begin{cases}
\vert T \vert =\forall \vert S \vert > a,\det A_{S,T}=0 \\
\exists \vert S \vert =\vert T \vert =a,\det A_{S,T}\ne 0
\end{cases}
\end{gathered} rank A = a ⟺ { ∣ T ∣ = ∀∣ S ∣ > a , det A S , T = 0 ∃∣ S ∣ = ∣ T ∣ = a , det A S , T = 0
首先如果存在det A S , T = ≠ 0 \det A_{S,T}=\ne 0 det A S , T = = 0 ,那么S S S 对应的这些行必然是满秩的,所以 rank A ≥ a \operatorname{rank} A\ge a rank A ≥ a .
而如果 rank A > a \operatorname{rank} A>a rank A > a ,那么选出a a a 个线性无关行,再从这里面找到a a a 个线性无关列,这就是一个行列式非0 0 0 的矩阵.
于是得证.
rank A = n − 1 ⟺ rank C = 1 rank A < n − 1 ⟺ rank C = 0 \begin{gathered}
\operatorname{rank} A=n-1 \iff \operatorname{rank} C=1 \\
\operatorname{rank} A<n-1 \iff \operatorname{rank} C=0
\end{gathered} rank A = n − 1 ⟺ rank C = 1 rank A < n − 1 ⟺ rank C = 0
由上一个conclusion,第二行是显然的(全0).
对于第一行,考虑A C T AC^T A C T =0,于是 A C T = 0 , rank A + rank C T ≤ n AC^T=0,\operatorname{rank} A+\operatorname{rank} C^T\le n A C T = 0 , rank A + rank C T ≤ n ,又有 rank C ≠ 0 \operatorname{rank} C\ne 0 rank C = 0 因为至少有一个非0 0 0 .
A m × n , B n × m , m ≤ n ⟹ det A B = ∑ ∣ S ∣ = m , S ⊂ [ 1 , n ] det A [ 1 , m ] ∩ Z , S det B S , [ 1 , m ] ∩ Z \begin{gathered}
A_{m\times n},B_{n\times m},m\le n \\
\implies \det AB=\sum _{\vert S \vert =m,S\subset [1,n]} \det A_{[1,m]\cap Z,S}\det B_{S,[1,m]\cap Z}
\end{gathered} A m × n , B n × m , m ≤ n ⟹ det A B = ∣ S ∣ = m , S ⊂ [ 1 , n ] ∑ det A [ 1 , m ] ∩ Z , S det B S , [ 1 , m ] ∩ Z
考虑矩阵
det [ I n B A 0 ] = det [ I n B 0 − A B ] \begin{gathered}
\det \begin{bmatrix}
I_n&B \\
A&0
\end{bmatrix}=\det
\begin{bmatrix}
I_n&B \\
0&-AB
\end{bmatrix}
\end{gathered} det [ I n A B 0 ] = det [ I n 0 B − A B ]
右边的行列式是( − 1 ) m det ( A B ) (-1)^m\det (AB) ( − 1 ) m det ( A B ) ,考虑左边用拉普拉斯定理,则显然你的子式的列只能选A A A 里面的,而它的余子式就是B B B 左边再加上若干列,若第i i i 列没被子式选,则这列一定要选i i i 行的元素,于是B B B 中恰好只有S S S 中的行能选.
对于符号,Laplace中的是∑ i ∈ S i + m ( n + 1 + n + m ) 2 \sum_{i\in S} i+\dfrac{m(n+1+n+m)}2 ∑ i ∈ S i + 2 m ( n + 1 + n + m ) ,算代数余子式行列式的时候还有一个∑ i ∉ S i − ( n − m ) ( n − m + 1 ) 2 \sum_{i\notin S} i-\dfrac{(n-m)(n-m+1)}2 ∑ i ∈ / S i − 2 ( n − m ) ( n − m + 1 ) ,然后这些加起来是同余m m m 的.
这样符号一乘正好是0 0 0 .
行列式与导数:容易注意到行列式每行不会有两个乘到一起,而最终的行列式是n n n 行各取一个数乘起来,于是遵循导数求导乘积的法则,或者说把每一行求导其他行不变再加起来.
Another Quiz Problem
∀ A m × n , ∃ B n × m s . t . A B A = A \begin{gathered}
\forall A_{m\times n},\exists B_{n\times m}\ s.t.\
ABA=A
\end{gathered} ∀ A m × n , ∃ B n × m s . t . A B A = A
Sol 1
A = P [ I k 0 0 0 ] Q \begin{gathered}
A=P\begin{bmatrix}
I_k&0 \\
0&0
\end{bmatrix}Q
\end{gathered} A = P [ I k 0 0 0 ] Q
于是让
B = Q − 1 [ I k 0 0 0 ] P − 1 \begin{gathered}
B=Q^{-1}\begin{bmatrix}
I_k&0 \\
0&0
\end{bmatrix}P^{-1}
\end{gathered} B = Q − 1 [ I k 0 0 0 ] P − 1
Sol 2
考虑一个任意情况,对任意集合X , Y X,Y X , Y 和任意映射f f f ,你可以找到一个g g g 使得f ∘ g ∘ f = f f\circ g\circ f=f f ∘ g ∘ f = f
显然的,因为你可以让g g g 把f f f 的像映射到任意一个原像.
那么这个显然的在说什么呢?它实际上在说,对任意一个f f f ,我们可以找到X , Y X,Y X , Y 各自的一个子集X ′ , Y ′ X',Y' X ′ , Y ′ ,使得f f f 限制在X ′ → Y ′ X'\to Y' X ′ → Y ′ 是双射.于是g g g 是这上面的逆,而这之外的随意映射就可以满足f g f = f fgf=f f g f = f .
然后你再看第一个证明,那么中间那个有I k I_k I k 的矩阵实际上就是,I k I_k I k 对应了X ′ → Y ′ X'\to Y' X ′ → Y ′ 的双射部分,这也是上面的做法为啥有道理.
Circulant Matrix's det
∀ a 0 … a n − 1 A i , j = a i + j − 2 m o d n ⟹ det A = ∏ i = 1 n ∑ j = 1 n w n i j c j \begin{gathered}
\forall a_0\ldots a_{n-1} \\
A_{i,j}=a_{i+j-2\bmod n} \\
\implies \det A=\prod _{i = 1} ^{n} \sum _{j = 1} ^{n} w_n^{ij}c_j
\end{gathered} ∀ a 0 … a n − 1 A i , j = a i + j − 2 mod n ⟹ det A = i = 1 ∏ n j = 1 ∑ n w n ij c j
首先我们让P n P_n P n 是循环移位P ( x 0 , … , x n ) = ( x 1 , … , x n , x 0 ) P(x_0,\ldots, x_n)=(x_1,\ldots, x_n,x_0) P ( x 0 , … , x n ) = ( x 1 , … , x n , x 0 ) ,那么有
A = ∑ i = 0 n − 1 a i P i \begin{gathered}
A=\sum _{i = 0} ^{n-1} a_iP^i
\end{gathered} A = i = 0 ∑ n − 1 a i P i
那么因为所有P P P 乘法的时候是保持不变子空间的,也就是说A A A 的特征值就是P P P 的特征值带入f ( x ) = ∑ i = 0 n − 1 a i x i f(x)=\sum_{i=0}^{n-1}a_ix^i f ( x ) = ∑ i = 0 n − 1 a i x i ,所以只需要看P P P 的特征值.
考虑P x = λ x Px=\lambda x P x = λ x ,那么x i + 1 m o d n = λ x i x_{i+1\bmod n}=\lambda x_i x i + 1 mod n = λ x i ,于是你乘一圈就有λ n = 1 \lambda^n=1 λ n = 1 ,能相应的构造出x i = w n i l , λ = w n l x_i=w_n^{il},\lambda=w_n^l x i = w n i l , λ = w n l 这样的特征向量和特征值.
于是P P P 的特征值就是单位根,带进去就得证.
A Determinance about Vandermonde plus 1
求
det A , A i , j = x i j + 1 \begin{gathered}
\det A,A_{i,j}=x_i^j+1
\end{gathered} det A , A i , j = x i j + 1
容易想到用线性性去拆分,如果我们对列拆分,你发现你就是要求A A A 有一列换成全1 1 1 之后的矩阵的行列式的和,我们设第i i i 列换掉之后的矩阵是B i B_i B i .
考虑Cramer's Rule,你注意到设A y = b Ay=b A y = b ,其中b b b 是全1 1 1 向量,则y i = det B i det A y_i=\dfrac{\det B_i}{\det A} y i = det A det B i ,于是我们实际上是求y i y_i y i ,
如果我们这里用y = A − 1 b y=A^{-1}b y = A − 1 b 我们可以得到一个结论:
det ( A + u v T ) = det A ( 1 + v T A − 1 u ) \begin{gathered}
\det(A+uv^T)=\det A(1+v^TA^{-1}u)
\end{gathered} det ( A + u v T ) = det A ( 1 + v T A − 1 u )
然后对这个题,我们也不知道A − 1 A^{-1} A − 1 ,所以你仍然要想办法求y i y_i y i .
碰到范德蒙德矩阵一定要想多项式,于是你想到P ( z ) = ∑ i = 1 n y i z i P(z)=\sum_{i=1}^n y_iz^i P ( z ) = ∑ i = 1 n y i z i ,则P ( x i ) = 1 P(x_i)=1 P ( x i ) = 1 ,于是P ( z ) = A ∏ i ( z − x i ) P(z)=A\prod_i (z-x_i) P ( z ) = A ∏ i ( z − x i ) ,于是你能通过比较常数项求出A A A ,然后再计算P ( 1 ) P(1) P ( 1 ) 就能求出∑ y i \sum y_i ∑ y i 了.
[think] 碰到范德蒙德矩阵一定要想多项式
det ( A + u v T ) = det A ( 1 + v T A − 1 u ) \begin{gathered}
\det(A+uv^T)=\det A(1+v^TA^{-1}u)
\end{gathered} det ( A + u v T ) = det A ( 1 + v T A − 1 u )
一个是上面提到的证法.
或者考虑
det ( I + u v T ) = ( 1 + v T u ) = ( 1 + u T v ) \begin{gathered}
\det(I+uv^T)=(1+v^Tu)=(1+u^Tv)
\end{gathered} det ( I + u v T ) = ( 1 + v T u ) = ( 1 + u T v )
这个可以直接考虑u v T uv^T u v T 特征值是n − 1 n-1 n − 1 个0 0 0 和一个v T u v^Tu v T u (对应特征向量u u u ).
然后det ( A + A u v T ) = det A ( 1 + u T v ) \det (A+Auv^T)=\det A(1+u^Tv) det ( A + A u v T ) = det A ( 1 + u T v )
所以让u ′ = A − 1 u u'=A^{-1}u u ′ = A − 1 u 代入,那么式子就成了det ( A + u v T ) = det A ( 1 + v T A − 1 u ) \det (A+uv^T)=\det A(1+v^TA^{-1}u) det ( A + u v T ) = det A ( 1 + v T A − 1 u ) .即我们所证的.
about eigen values
det Q A Q − 1 − λ I = det Q A Q − 1 − Q λ I Q − 1 = det Q ( A − λ I ) Q − 1 = det A − λ I \begin{gathered}
\det QAQ^{-1}-\lambda I \\
=\det QAQ^{-1}-Q\lambda I Q^{-1} \\
= \det Q(A-\lambda I)Q^{-1} \\
=\det A-\lambda I
\end{gathered} det Q A Q − 1 − λ I = det Q A Q − 1 − Q λ I Q − 1 = det Q ( A − λ I ) Q − 1 = det A − λ I
显然的吧,因为都是一个线性变换在不同的基下的.
首先你可以谱定理启动.
否则考虑
v 1 A v 2 = v 1 λ 2 v 2 ( v 1 A v 2 ) T = λ 1 v 1 v 2 \begin{gathered}
v_1Av_2=v_1\lambda_2 v_2 \\
(v_1Av_2)^T=\lambda_1 v_1v_2
\end{gathered} v 1 A v 2 = v 1 λ 2 v 2 ( v 1 A v 2 ) T = λ 1 v 1 v 2
于是就得证了.
特征多项式的n − i n-i n − i 项系数等于:( − 1 ) i × sum of all A’s principle minor of order i (-1)^i\times \text{sum of all A's principle minor of order i} ( − 1 ) i × sum of all A’s principle minor of order i
实谱定理(矩阵版)
任意对称矩阵A A A 可以写成A = P D P − 1 A=PDP^{-1} A = P D P − 1 ,其中D D D 是对角矩阵,P P P 是正交矩阵.
首先复化,考虑A A A 的一组特征值特征向量A v = λ v Av=\lambda v A v = λ v ,则 v T A v ‾ = λ ‾ ∣ ∣ v ∣ ∣ v^TA\overline{v}=\overline\lambda \vert\vert v \vert\vert v T A v = λ ∣∣ v ∣∣ ,同时取共轭得这个还等于 λ ∣ ∣ v ∣ ∣ \lambda \vert\vert v \vert\vert λ ∣∣ v ∣∣ ,于是能证明λ ∈ R \lambda \in R λ ∈ R .
然后我们说明了A A A 的特征值都是实的,我们取v 1 v_1 v 1 是一个特征值,并将其扩充到一组标准正交基v 1 , v 2 … v n v_1,v_2\ldots v_n v 1 , v 2 … v n ,于是
A [ v 1 … v n ] = [ v 1 … v n ] [ λ 1 x 0 A 1 ] \begin{gathered}
A[v_1\ldots v_n]=[v_1\ldots v_n]
\begin{bmatrix}
\lambda_1&x \\
0&A_1
\end{bmatrix}
\end{gathered} A [ v 1 … v n ] = [ v 1 … v n ] [ λ 1 0 x A 1 ]
其中x x x 和A 1 A_1 A 1 是任意向量/矩阵,而中间v 1 … v n v_1\ldots v_n v 1 … v n 是一个正交矩阵设为B B B ,则右边那个分块矩阵是B − 1 A B = B T A B B^{-1}AB=B^TAB B − 1 A B = B T A B ,于是这个也是对称的,从而x = 0 x=0 x = 0 ,A 1 A_1 A 1 是对称矩阵,从而可以归纳到n − 1 n-1 n − 1 维.归纳得证.
这里x = 0 x=0 x = 0 就是T T T 在U ⊥ U^\perp U ⊥ 不变,A 1 A_1 A 1 对称就是T T T 限制仍然是自伴.
二次型
惯性定理 Inertia Theorem
任意二次型A A A 通过任意可逆矩阵P P P 合同变换到对角矩阵P D P T PDP^T P D P T ,D D D 的对角线上0 0 0 的个数,正数个数,负数个数都是确定的(与P P P 无关).
首先0 0 0 的个数是确定的,因为这是矩阵的秩.
考虑在两组基下,假设它们对角型正数的个数分别为a , b a,b a , b ,负数是r − a , r − b r-a,r-b r − a , r − b ,不妨设a < b a<b a < b .
那么我们不考虑0 0 0 的那部分限制到r r r 维的空间上考虑,则你可以找到一个a a a 维的空间,上面值都是正的(要求所有负惯性指数位置都是负的),再找一个r − b r-b r − b 上面都是负的.而a + r − b > r a+r-b>r a + r − b > r 所以这俩空间有交,爆炸了.
[think] 注意到正惯性指数是最大的子空间满足这里面的值带进去都是正的.
考虑对角化的过程,如果已经消除了左上角k × k k\times k k × k 的矩阵,接下来这些地方就不会动了.然后你消( k + 1 ) × ( k + 1 ) (k+1)\times (k+1) ( k + 1 ) × ( k + 1 ) 的左上角,注意到我们前面已经进行的行列变换都是在前k k k 行进行的,所以前面所有的操作都不会改变当前( k + 1 ) × ( k + 1 ) (k+1)\times (k+1) ( k + 1 ) × ( k + 1 ) 矩阵的行列式,直到你消完.于是每个顺序主子式对应的子矩阵消前和消后的行列式是一样的.
SVD
PCA(主成分分析)
你有若干个向量a i a_i a i (若干个点)
先让a a a 都减去平均值得到a ′ a' a ′ 令A = [ a 1 ′ … a n ′ ] A=[a_1'\ldots a_n'] A = [ a 1 ′ … a n ′ ] ,令B = A A T B=AA^T B = A A T ,你发现B i , j = Var ( x i , x j ) B_{i,j}=\operatorname{Var}(x_i,x_j) B i , j = Var ( x i , x j ) (两维的协方差.)
对A A A 进行SVD分解A = U D V A=UDV A = U D V ,则若s 1 s_1 s 1 对应的方向v 1 v_1 v 1 ,直线t v 1 tv_1 t v 1 是最小化每个点到直线距离(垂线段距离)平方的直线.
设最优方向为单位向量n n n .
L = ∑ ( x i − n < x i , n > ) 2 = ∑ a i 2 + n 2 < a i , n > 2 − 2 < a i , n > < a i , n > = ( ∑ a i 2 ) − ( ∑ i < a i , n > 2 ) \begin{gathered}
L=\sum (x_i-n<x_i,n>)^2 \\
=\sum a_i^2+n^2<a_i,n>^2-2<a_i,n><a_i,n> \\
=(\sum a_i^2)-(\sum_i<a_i,n>^2)
\end{gathered} L = ∑ ( x i − n < x i , n > ) 2 = ∑ a i 2 + n 2 < a i , n > 2 − 2 < a i , n >< a i , n > = ( ∑ a i 2 ) − ( i ∑ < a i , n > 2 )
只要最大化第二项.考虑
∑ i < a i , n > 2 = n T ( ∑ a i a i T ) n = n T ( A A T ) n \begin{gathered}
\sum_i <a_i,n>^2 \\
=n^T(\sum a_i a_i^T)n \\
=n^T(AA^T)n
\end{gathered} i ∑ < a i , n > 2 = n T ( ∑ a i a i T ) n = n T ( A A T ) n
这是二次型,把A A T AA^T A A T 谱分解就看出显然最大值在它最大的特征值对应特征向量方向,也就是v 1 v_1 v 1 方向.
你还可以看出这个值是∑ a i 2 − s 1 2 \sum a_i^2-s_1^2 ∑ a i 2 − s 1 2
SVD求伪逆
A = U Σ V T A=U\Sigma V^T A = U Σ V T ,A + = V Σ + U T A^+=V\Sigma^+ U^T A + = V Σ + U T ,其中Σ + \Sigma^+ Σ + 是把所有非零对角线元素变成倒数再取转置.
用SVD逼近
A = U D V T = ∑ i u i s i v i T ⟹ min rank C = k ∥ A − C ∥ = s k + 1 \begin{gathered}
A=UDV^T=\sum_i u_is_iv_i^T \\
\implies \min_{\operatorname{rank} C=k} \Vert A-C \Vert = s_{k+1}
\end{gathered} A = U D V T = i ∑ u i s i v i T ⟹ rank C = k min ∥ A − C ∥ = s k + 1
设A m × n , V n × n , U m × m , D m × n A_{m\times n},V_{n\times n},U_{m\times m},D_{m\times n} A m × n , V n × n , U m × m , D m × n
考虑对A A A 分解得A = U D V T A=UDV^T A = U D V T ,则我们让D D D 只保留前k k k 个特征值,其余设为0 0 0 得到D ′ D' D ′ ,令C = U D ′ V T C=UD'V^T C = U D ′ V T ,则显然A − C A-C A − C 的SVD分解就是U ( D − D ′ ) V T U(D-D')V^T U ( D − D ′ ) V T ,其最大的奇异值恰为s k + 1 s_{k+1} s k + 1 (就是只用前k k k 个奇异值去逼近啦).
然后只需要证明这是最小值.考虑A A A 的前k + 1 k+1 k + 1 个右奇异向量(v 1 … v k v_1\ldots v_k v 1 … v k )张成的子空间和B B B 的零空间的交( dim null B = n − rank B = n − k \dim \operatorname{null} B=n-\operatorname{rank} B=n-k dim null B = n − rank B = n − k ,故必然有交)中的单位向量x x x .
而
∥ ( A − B ) x ∥ 2 = ∥ A x ∥ 2 = x T A T A x = x T V D 2 V T x = ( V T x ) T D 2 ( V T x ) ≥ s k + 1 2 \begin{gathered}
\Vert (A-B)x \Vert^2 \\
=\Vert Ax \Vert^2 \\
=x^TA^TAx \\
=x^TVD^2V^Tx \\
=(V^Tx)^TD^2(V^Tx) \\
\ge s_{k+1}^2
\end{gathered} ∥ ( A − B ) x ∥ 2 = ∥ A x ∥ 2 = x T A T A x = x T V D 2 V T x = ( V T x ) T D 2 ( V T x ) ≥ s k + 1 2
于是得证.
[think] 证明最小化的时候从分析矩阵结构改为找一个向量.因为这个向量显然一定可以在前k + 1 k+1 k + 1 维(这个每个维度都放大需求倍数的空间中任何一个向量也都被放大需求倍数).
多项式
有理域可约多项式
f ( x ) = ∑ i = 0 n a i x i is primitive ⟺ gcd ( a 0 , … , a n ) = 1 \begin{gathered}
f(x)= \sum _{i = 0} ^{n} a_ix^i \text{ is primitive } \\
\iff \gcd(a_0,\ldots ,a_n)=1
\end{gathered} f ( x ) = i = 0 ∑ n a i x i is primitive ⟺ g cd( a 0 , … , a n ) = 1
Gauss's Lemma
f , g are primitive ⟹ f g are primitive \begin{gathered}
f,g \text{ are primitive } \implies fg \text{ are primitive}
\end{gathered} f , g are primitive ⟹ f g are primitive
考虑反证,如果f g fg f g 有公质因数p p p ,且p ∤ g p\not\vert g p ∣ g ,证明p ∣ f p\vert f p ∣ f .
令f , g , f g f,g,fg f , g , f g 对应的系数列分别是a n , b n , c n a_n,b_n,c_n a n , b n , c n .
那么直接递推就好了,我们假设p ∣ gcd ( f 0 , … , f A , g 0 , … , g B ) p \vert \gcd(f_0,\ldots ,f_A,g_0,\ldots,g_B) p ∣ g cd( f 0 , … , f A , g 0 , … , g B ) (没有取− 1 -1 − 1 )
那么考虑c A + B + 2 c_{A+B+2} c A + B + 2 ,它一定有一项是f A + 1 g B + 1 f_{A+1}g_{B+1} f A + 1 g B + 1 ,且其他的项一定都被p p p 整除,于是你一定可以让A A A 或B B B 加1 1 1 .
这样走n + m + 1 n+m+1 n + m + 1 步一定可以把A A A 或B B B 弄满.
f ( x ) ∈ Z [ x ] , f is primitive ⟹ ( f ( x ) is inreducible in Z[x] ⟺ f ( x ) is inreducible in Q[x] ) \begin{gathered}
f(x)\in Z[x],f \text{ is primitive} \\
\implies
(f(x) \text{ is inreducible in Z[x]} \iff f(x) \text{ is inreducible in Q[x]} )
\end{gathered} f ( x ) ∈ Z [ x ] , f is primitive ⟹ ( f ( x ) is inreducible in Z[x] ⟺ f ( x ) is inreducible in Q[x] )
考虑如果f ( x ) = f 1 ( x ) f 2 ( x ) , f 1 , f 2 ∈ Q [ x ] f(x)=f_1(x)f_2(x),f_1,f_2\in Q[x] f ( x ) = f 1 ( x ) f 2 ( x ) , f 1 , f 2 ∈ Q [ x ] .
那么你显然可以有f i ( x ) = g i ( x ) p i q i f_i(x)=g_i(x)\dfrac {p_i}{q_i} f i ( x ) = g i ( x ) q i p i .且g 1 , g 2 is primitive g_1,g_2 \text{ is primitive} g 1 , g 2 is primitive .
那么f ( x ) = p 1 p 2 q 1 q 2 g 1 ( x ) g 2 ( x ) f(x)=\dfrac{p_1p_2}{q_1q_2}g_1(x)g_2(x) f ( x ) = q 1 q 2 p 1 p 2 g 1 ( x ) g 2 ( x ) ,其中g 1 g 2 g_1g_2 g 1 g 2 还是primitive \text{primitive} primitive 的,q 1 q 2 f ( x ) = p 1 p 2 g 1 ( x ) g 2 ( x ) q_1q_2f(x)=p_1p_2g_1(x)g_2(x) q 1 q 2 f ( x ) = p 1 p 2 g 1 ( x ) g 2 ( x ) .
考虑两侧的系数的最大公因数相等,那只能是q 1 q 2 = p 1 p 2 q_1q_2=p_1p_2 q 1 q 2 = p 1 p 2 了.
{ f ( x ) = x n + ∑ i = 0 n − 1 a i x i , a i ∈ Z ∃ p s . t . ∀ 0 ≤ i ≤ n − 1 , p ∣ a i p 2 ∤ a 0 ⟹ f ( x ) is irreducible in Q[x] \begin{gathered}
\begin{cases}
f(x)=x^n+\sum _{i = 0} ^{n-1} a_ix^i,a_i\in Z \\
\exists p \ s.t.\\
\forall 0\le i\le n-1 ,p | a_i \\
p^2 \not| a_0
\end{cases} \\
\implies f(x) \text{ is irreducible in Q[x]}
\end{gathered} ⎩ ⎨ ⎧ f ( x ) = x n + ∑ i = 0 n − 1 a i x i , a i ∈ Z ∃ p s . t . ∀0 ≤ i ≤ n − 1 , p ∣ a i p 2 ∣ a 0 ⟹ f ( x ) is irreducible in Q[x]
反证,考虑f = f 1 f 2 f=f_1f_2 f = f 1 f 2 设它们的系数分别是b n , c n b_n,c_n b n , c n ,根据上面高斯引理的证明归纳过程,p ∣ a 0 p|a_0 p ∣ a 0 ,那么p ∣ b 0 c 0 p|b_0c_0 p ∣ b 0 c 0 ,那么因为p p p 只能整除其中的一个,由上面那个由( A , B ) (A,B) ( A , B ) 推到( A + 1 , B ) (A+1,B) ( A + 1 , B ) 或( A , B + 1 ) (A,B+1) ( A , B + 1 ) 的过程,你有一边走不了,最终就直接得到p ∣ b n p\vert b_n p ∣ b n 或p ∣ c n p\vert c_n p ∣ c n 了.
但是显然p ∤ f p\not| f p ∣ f ,它是有系数是1 1 1 的!矛盾.
你看一下下面那个模p p p 的性质于是这个东西另一个证法是先模p p p ,那么f ( x ) ≡ x n f(x)\equiv x^n f ( x ) ≡ x n ,于是它被拆了之后只能是x i x n − i x^i x^{n-i} x i x n − i ,但这说明拆出来的两项常数项模p p p 都是0 0 0 ,和p 2 ∤ a 0 p^2\not|a_0 p 2 ∣ a 0 矛盾.
p ∈ P , f ( x ) = ∑ i = 0 p − 1 x i ⟹ f ( x ) is irreducible \begin{gathered}
p\in P,f(x)=\sum _{i = 0} ^{p-1} x^i \\
\implies f(x) \text{ is irreducible}
\end{gathered} p ∈ P , f ( x ) = i = 0 ∑ p − 1 x i ⟹ f ( x ) is irreducible
f ( x ) = x p − 1 x − 1 = y = x − 1 ( y + 1 ) p − 1 y = ∑ i = 0 p − 1 y i ( p i + 1 ) \begin{gathered}
f(x)=\dfrac{x^p-1}{x-1} \\
\xlongequal{ y=x-1 } \dfrac{(y+1)^p-1}{y} \\
=\sum _{i = 0} ^{p-1} y^i \binom{p}{i+1}
\end{gathered} f ( x ) = x − 1 x p − 1 y = x − 1 y ( y + 1 ) p − 1 = i = 0 ∑ p − 1 y i ( i + 1 p )
此时用上面定理.
若首一多项式f f f 在Z [ x ] Z[x] Z [ x ] 上可约则f f f 在F p [ x ] F_p[x] F p [ x ] 上可约
显然,因为你把Z [ x ] Z[x] Z [ x ] 上那两个因子分别模p p p 就好了.首一主要是防止你模p p p 的时候f f f 模成常数.
Smith标准型是说环上元素的矩阵,相抵于一个只有主对角线上有d 1 … d r , 0 … , 0 d_1\ldots d_r,0\ldots,0 d 1 … d r , 0 … , 0 的矩阵且满足d i + 1 ∣ d i d_{i+1}\vert d_i d i + 1 ∣ d i .
定义D i D_i D i 是所有i i i 阶子式的gcd,则d i = D i D i − 1 d_i=\dfrac{D_i}{D_{i-1}} d i = D i − 1 D i .
如果我们消完了的d d d 有整除链的性质,那么最后是显然满足条件的,只需要证明D i D_i D i 在初等行变换下不变:
交换两行/列或给一行/列乘± 1 \pm 1 ± 1 (环上,其他的没逆)显然不影响
把一行/列乘k k k 加到另一行/列,此时若一个子式同时包含或不包含涉及到的两行则它不变,如果只包含一个则用线性性相当于给他加上另一个子式,不改变gcd \gcd g cd .
而只要证明存在一种消法消成整除链,你发现如果消完了不是整除链你总可以用辗转相减变成d i = gcd ( d 1 ′ , … , d i ′ ) d_i=\gcd(d'_1,\ldots,d'_i) d i = g cd( d 1 ′ , … , d i ′ ) 弄成整除链,于是得证.
这里d d d 叫不变因子(invariant factors),D D D 叫行列式因子(determinant factors).
然后如果这个环是UFD,定义初等因子e l e m e n t f a c t o r s element factors e l e m e n t f a c t or s 是把每个d d d 做素分解得到的每个p i c i p_i^{c_i} p i c i ,不进行合并.
你可以把矩阵弄成只有分块对角矩阵,且每个块只有对角线上有相同的值,对角线上方的一条斜线都是1 1 1 .
Proof1:初等变换
记E ( i , j , a ) E(i,j,a) E ( i , j , a ) 表示I I I 的基础上( i , j ) (i,j) ( i , j ) 位置为a a a 的初等变换矩阵.而容易注意到E ( i , j , a ) − 1 = E ( i , j , − a ) E(i,j,a)^{-1}=E(i,j,-a) E ( i , j , a ) − 1 = E ( i , j , − a ) .
所以当我们做E ( i , j , k ) A E ( i , j , − k ) E(i,j,k)AE(i,j,-k) E ( i , j , k ) A E ( i , j , − k ) 这个相似变换我们是:
给A A A 的第i i i 行加上k k k 倍第j j j 行,再给第j j j 列加上k k k 倍第i i i 列.
对一个上三角矩阵,我们发现进行这样一次操作,A i , j A_{i,j} A i , j 加上了k ( A i , i − A j , j ) k(A_{i,i}-A_{j,j}) k ( A i , i − A j , j ) ,且这一步对上三角矩阵只影响x < i , y > j x<i,y>j x < i , y > j 的位置(一个右上角).所以你可以以一定的顺序消掉这样的元素,得到一个分块对角矩阵,每个块对角线元素相同.
那么每个块都可以表示成λ I + N \lambda I+N λ I + N ,N N N 是幂零的.而当你对这个子空间进行换基的时候λ I \lambda I λ I 不变所以只要找N N N 的变换.
你发现对于N N N ,你仍然进行过上面这种操作就可以简单的用某一行的1 1 1 把同行后面的所有数消掉,所以你从第一行开始从上往下消.如果遇到一行它的1 1 1 不在副对角线你可以做交换把它换到前面来,这样都弄完了就得到标准的约旦标准型了.
Proof2
两矩阵A , B A,B A , B 相似等价于特征矩阵A − λ I , B − λ I A-\lambda I,B-\lambda I A − λ I , B − λ I 相抵.
首先如果A = P B P − 1 A=PBP^{-1} A = P B P − 1 ,那么
A − λ I = P B P − 1 − P λ I P − 1 = P ( B − λ I ) P − 1 \begin{gathered}
A-\lambda I=PBP^{-1}-P\lambda IP^{-1}=P(B-\lambda I)P^{-1}
\end{gathered} A − λ I = P B P − 1 − P λ I P − 1 = P ( B − λ I ) P − 1
成立.
再看另一边:
考虑现在V ( A − λ I ) = ( B − λ I ) W V(A-\lambda I)=(B-\lambda I)W V ( A − λ I ) = ( B − λ I ) W ,V , W V,W V , W 是可逆的λ \lambda λ 矩阵.
直接设 V = ∑ i = 0 m V i λ i , W = ∑ i = 0 m W i λ i V=\sum _{i = 0} ^{m} V_i\lambda^i,W=\sum _{i = 0} ^{m} W_i\lambda^i V = ∑ i = 0 m V i λ i , W = ∑ i = 0 m W i λ i .
那么
{ − V m = − W m V k A − V k − 1 = B W k − W k − 1 V 0 A = W 0 B \begin{gathered}
\begin{cases}
-V_m=-W_m \\
V_kA-V_{k-1}=BW_k-W_{k-1} \\
V_0A=W_0B
\end{cases}
\end{gathered} ⎩ ⎨ ⎧ − V m = − W m V k A − V k − 1 = B W k − W k − 1 V 0 A = W 0 B
然后你发现可以盯着一边看去做一个相消,就是给V i A V_iA V i A 的那个等式右乘A i A^i A i ,则左边是0 0 0 ,右边是?
∑ B W i A i − W i − 1 A i = B ( ∑ W i A i ) − ∑ ( W i A i ) A = 0 \begin{gathered}
\sum BW_iA^i-W_{i-1}A^i \\
=B(\sum W_iA^i)-\sum (W_iA^i)A \\
=0
\end{gathered} ∑ B W i A i − W i − 1 A i = B ( ∑ W i A i ) − ∑ ( W i A i ) A = 0
于是设P = ∑ W i A i P=\sum W_iA^i P = ∑ W i A i 就有B P = P A BP=PA B P = P A .
那么你把V , W V,W V , W 取逆分别移到另一边是( A − λ I ) W − 1 = V − 1 ( B − λ I ) (A-\lambda I)W^{-1}=V^{-1}(B-\lambda I) ( A − λ I ) W − 1 = V − 1 ( B − λ I ) .再走一遍能得到A Q = Q B AQ=QB A Q = QB .
然后能看出A k Q = Q B k A^kQ=QB^k A k Q = Q B k ,Q = W − 1 ( B ) Q=W^{-1}(B) Q = W − 1 ( B ) ,P A k = B k P PA^k=B^kP P A k = B k P .
那么
P Q = W ( A ) Q = ∑ i W i A i Q = ∑ i W i Q B i = ∑ i W i ( ∑ j W j − 1 B j ) B i = ∑ i ∑ j W i W j − 1 B i + j = W W − 1 ( B ) = I \begin{gathered}
PQ \\
=W(A)Q \\
=\sum_i W_iA^iQ \\
=\sum_i W_iQB^i \\
=\sum_i W_i(\sum_j W^{-1}_jB^j)B^i \\
=\sum_i \sum_j W_iW^{-1}_jB^{i+j} \\
=WW^{-1}(B) \\
=I
\end{gathered} P Q = W ( A ) Q = i ∑ W i A i Q = i ∑ W i Q B i = i ∑ W i ( j ∑ W j − 1 B j ) B i = i ∑ j ∑ W i W j − 1 B i + j = W W − 1 ( B ) = I
[think] 你看那个神秘的求和是不自然的,你注意到如果我们对一开始的V ( A − λ I ) = ( B − λ I ) W V(A-\lambda I)=(B-\lambda I)W V ( A − λ I ) = ( B − λ I ) W 去代入λ = A \lambda =A λ = A 会很妙,然后所谓的代入实际上对应了取模,所以算( B − λ I ) W m o d ( A − λ I ) (B-\lambda I)W \bmod (A-\lambda I) ( B − λ I ) W mod ( A − λ I ) ,算这个的时候拆开W W W 乘上B B B 算,因为没有交换律不保证余数可乘.
[think] 好吧上面那个think说的你直接代入就可以的,这种时候你要规定代入方向,然后( B − λ I ) W = B W ( λ ) − λ W ( λ ) (B-\lambda I)W=BW(\lambda)-\lambda W(\lambda) ( B − λ I ) W = B W ( λ ) − λW ( λ ) ,此时你代入λ = A \lambda=A λ = A 因为是右代入所以你要先把λ W \lambda W λW 变成W λ W\lambda W λ .
[think] 后面这个= I =I = I 也可以直接通过代入看出来,或者再逐项系数展开.
考虑Smith标准型,那么从矩阵角度,分块对角矩阵有两个块,问题是现在不满足整除链,那么只要每次,对所有素因子,取其最高次组成一个作为不变因子,就可以重排出一个大矩阵的整除链,于是得证.
若V = W 1 ⊕ W 2 V=W_1\oplus W_2 V = W 1 ⊕ W 2 ,则T T T 的初等因子就是T ∣ W 1 T|_{W_1} T ∣ W 1 的初等因子与T ∣ W 2 T|_{W_2} T ∣ W 2 的初等因子的和(多重集的和).
考虑你有矩阵A A A 和矩阵B B B 都已经是Smith标准型,然后你要把它们拼起来.
此时考虑行列式因子D D D ,因为A , B A,B A , B 是对角阵拼完了还是对角阵,所以只有主子式不是0 0 0 ,此时k k k 阶行列式一定是A A A 里取若干个元素B B B 里取若干个元素的乘积,又因为整除链的存在,所以取的时候一定分别是A , B A,B A , B 里最低次的若干个.
那么因为是直接乘积再取gcd,现在只考虑所有不变因子d d d 某个特定素因子的次数就成了,新的大矩阵的D D D 对应的次数是原来两个小的min+卷积,因为整除链所以D D D 的次数是凸的,所以是闵和,所有作为斜率的那个d d d 的初等因子的次数一定都保留到了新的结果里,就完事了.
有理标准型
把矩阵的不变因子的友矩阵作为分块对角矩阵的块
由于上面说初等因子是直接拼,而初等因子确定不变因子,所以你这么拼是对的.
首先友矩阵的极小多项式是它的特征多项式,因为取[ 1 , 0 , 0 , … , 0 ] [1,0,0,\ldots,0] [ 1 , 0 , 0 , … , 0 ] ,那么你发现友矩阵的前n − 1 n-1 n − 1 次方作用到它身上是线性无关的,所以特征多项式的次数必须是n n n 次的.
然后你转乘上面的有理标准型,就显然.
于是这种证法说你考虑一个Jordan块正好对应了一个不变因子( x − λ i ) n (x - \lambda_i)^n ( x − λ i ) n .
首先这说明复数域上,他俩不变因子相同,所以行列式因子相同.又因为行列式因子一定是实的,所以结束.
准素分解
T T T 的极小多项式是m ( x ) = ∏ i p i ( x ) c i m(x)=\prod_i p_i(x)^{c_i} m ( x ) = ∏ i p i ( x ) c i ,p i p_i p i 是互不相同的不可约多项式,则V = ⊕ i W i V=\oplus_i W_i V = ⊕ i W i ,且 W i = null p i c i ( T ) W_i=\operatorname{null} p_i^{c_i}(T) W i = null p i c i ( T ) .
Sol1:(from Gemini and me )
不变性是显然的,我们发现对一个W i W_i W i ,其极小多项式是p i c i p_i^{c_i} p i c i ,考虑其中一个v v v ,p i c i ( T ) v = 0 p_i^{c_i}(T)v=0 p i c i ( T ) v = 0 ,那么显然p i c i ( T ) T v = T p i c i ( T ) v = 0 p_i^{c_i}(T)Tv=Tp_i^{c_i}(T)v=0 p i c i ( T ) T v = T p i c i ( T ) v = 0 .
let f i = m ( x ) p i c i ( x ) if u ∈ W 1 ∩ W 2 ⊕ … ⊕ W n ≠ ∅ u ∈ W 1 ⟹ p 1 c 1 ( T ) u = 0 u ∈ W 2 ⊕ … ⊕ W n ⟹ f i ( T ) u = 0 ∵ f i ⊥ p 1 c 1 ∴ ∃ s ( x ) , t ( x ) , s ( x ) p 1 c 1 ( x ) + t ( x ) f i ( x ) = 1 s ( T ) p 1 c 1 ( T ) u + t ( T ) f i ( T ) u = I u = u ⟹ u = 0 \begin{gathered}
\text{let } f_i=\dfrac{m(x)}{p_i^{c_i}(x)} \\
\text{if } u\in W_1\cap W_2\oplus \ldots\oplus W_n\ne \emptyset \\
u\in W_1 \implies p_1^{c_1}(T)u=0 \\
u\in W_2\oplus \ldots\oplus W_n \implies f_i(T)u=0 \\
\because f_i \perp p_1^{c_1} \\
\therefore \exists s(x),t(x),s(x)p_1^{c_1}(x)+t(x)f_i(x)=1 \\
s(T)p_1^{c_1}(T)u+t(T)f_i(T)u=Iu=u \\
\implies u=0
\end{gathered} let f i = p i c i ( x ) m ( x ) if u ∈ W 1 ∩ W 2 ⊕ … ⊕ W n = ∅ u ∈ W 1 ⟹ p 1 c 1 ( T ) u = 0 u ∈ W 2 ⊕ … ⊕ W n ⟹ f i ( T ) u = 0 ∵ f i ⊥ p 1 c 1 ∴ ∃ s ( x ) , t ( x ) , s ( x ) p 1 c 1 ( x ) + t ( x ) f i ( x ) = 1 s ( T ) p 1 c 1 ( T ) u + t ( T ) f i ( T ) u = I u = u ⟹ u = 0
这说明是直和,再说明直和是整个空间:
gcd ( f 1 , … , f n ) = 1 ⟹ ∃ a i , ∑ i a i ( x ) f i ( x ) = 1 ⟹ ∑ i a i ( T ) f i ( T ) = I let E i = a i ( T ) f i ( T ) \begin{gathered}
\gcd(f_1,\ldots,f_n)=1 \\
\implies \exists a_i,\sum_i a_i(x)f_i(x)=1 \\
\implies \sum_i a_i(T)f_i(T)=I \\
\text{let } E_i=a_i(T)f_i(T) \\
\end{gathered} g cd( f 1 , … , f n ) = 1 ⟹ ∃ a i , i ∑ a i ( x ) f i ( x ) = 1 ⟹ i ∑ a i ( T ) f i ( T ) = I let E i = a i ( T ) f i ( T )
显然E i E j ( T ) = 0 E_iE_j(T)=0 E i E j ( T ) = 0 ,又因为
let E i ( v ) = w p i c i ( T ) w = p i c i ( T ) a i ( T ) f i ( T ) v = a i ( T ) m ( T ) = 0 ⟹ range E i ⊂ W i ( ∑ j E j ) E i = E i ⟹ E i 2 = E i ∀ w ∈ W i , ∑ j E j ( T ) w = w \begin{gathered}
\text{let } E_i(v)=w \\
p_i^{c_i}(T)w=p_i^{c_i}(T)a_i(T)f_i(T)v \\
=a_i(T)m(T)=0 \\
\implies \operatorname{range} E_i \subset W_i \\
(\sum_j E_j)E_i=E_i \implies E_i^2=E_i \\
\forall w\in W_i,
\sum_j E_j(T)w=w \\
\end{gathered} let E i ( v ) = w p i c i ( T ) w = p i c i ( T ) a i ( T ) f i ( T ) v = a i ( T ) m ( T ) = 0 ⟹ range E i ⊂ W i ( j ∑ E j ) E i = E i ⟹ E i 2 = E i ∀ w ∈ W i , j ∑ E j ( T ) w = w
而你E j , j ≠ i E_j,j\ne i E j , j = i 的像空间的直和与W i W_i W i 无交,所以一定只能是E i ( w ) = w E_i(w)=w E i ( w ) = w
这样实际证明了E i E_i E i 就是像每个空间的投影,且和为I I I ,于是得证.
Sol2: from teacher
Lemma
( f ( x ) , g ( x ) ) = 1 ⟹ null f g ( T ) = null f ( T ) ⊕ null g ( T ) \begin{gathered}
(f(x),g(x))=1 \implies \operatorname{null} fg(T)=\operatorname{null} f(T) \oplus \operatorname{null} g(T)
\end{gathered} ( f ( x ) , g ( x )) = 1 ⟹ null f g ( T ) = null f ( T ) ⊕ null g ( T )
首先,裴蜀定理,a ( x ) f ( x ) + b ( x ) y ( x ) = 1 a(x)f(x)+b(x)y(x)=1 a ( x ) f ( x ) + b ( x ) y ( x ) = 1 ,于是如果 u ∈ null f ( T ) ∩ null g ( T ) u\in \operatorname{null} f(T)\cap \operatorname{null} g(T) u ∈ null f ( T ) ∩ null g ( T ) ,则a ( T ) f ( T ) u + b ( T ) y ( T ) u = I u = 0 a(T)f(T)u+b(T)y(T)u=Iu=0 a ( T ) f ( T ) u + b ( T ) y ( T ) u = I u = 0 ,得u = 0 u=0 u = 0 .
然后考虑,左边显然包含右边,只要证明右边包含左边,注意到 f ( T ) u ∈ null g ( T ) , g ( T ) u ∈ null f ( T ) f(T)u\in \operatorname{null} g(T),g(T)u\in \operatorname{null} f(T) f ( T ) u ∈ null g ( T ) , g ( T ) u ∈ null f ( T ) ,又有a ( T ) f ( T ) u + b ( T ) g ( T ) u = u a(T)f(T)u+b(T)g(T)u=u a ( T ) f ( T ) u + b ( T ) g ( T ) u = u ,于是拆成u = a ( T ) f ( T ) u + b ( T ) g ( T ) u u=a(T)f(T)u+b(T)g(T)u u = a ( T ) f ( T ) u + b ( T ) g ( T ) u 即证.
于是直接证完了.
[think] 观察证明复杂度可知,你还是先考虑两个比较好()另外是要看出本质条件是互素.
对任意线性变换T T T ,∃ v \exists v ∃ v 使得任意满足p ( T ) v = 0 p(T)v=0 p ( T ) v = 0 的p p p 是极小多项式的倍式.
首先如果p = p i c i p=p_i^{c_i} p = p i c i ,那么是显然的,否则p p p 就不是这个了.
否则用准素分解,则你得到若干个p i c i p_i^{c_i} p i c i ,每个对应一个向量v i v_i v i ,考虑直接把它们加起来v = ∑ i v i v=\sum_i v_i v = ∑ i v i ,那么:
0 = g ( T ) ∑ i v i = ∑ i g ( T ) v i 0=g(T)\sum_i v_i=\sum_i g(T)v_i 0 = g ( T ) ∑ i v i = ∑ i g ( T ) v i ,因为后面每一项分别在W i W_i W i 中,所以一定有g ( T ) v i = 0 g(T)v_i=0 g ( T ) v i = 0 ,所以在每个子空间中的极小多项式p i c i p_i^{c_i} p i c i 都是g ( T ) g(T) g ( T ) 的因子,于是m ∣ g m|g m ∣ g ,又因为g g g 是整个的零化多项式所以g ∣ m g|m g ∣ m .得证.
循环分解
对线性变换T T T ,找到v v v 使得 g ( T ) v = 0 ⟹ p ∣ g g(T)v=0 \implies p|g g ( T ) v = 0 ⟹ p ∣ g ,其中p p p 为极小多项式.然后找到最大的k k k 使得v , T v , … T k v v,Tv,\ldots T^kv v , T v , … T k v 线性无关,设 W 1 = span ( v , T v , … T k v ) W_1=\operatorname{span}( v,Tv,\ldots T^kv ) W 1 = span ( v , T v , … T k v ) ,设W 1 ⊕ W = V W_1\oplus W=V W 1 ⊕ W = V ,则继续对W W W 重复就得到一串不变子空间,且极小多项式构成整除链.
前面我们已经证明了一定有这样的v v v ,同时容易看出T ∣ W 1 T|_{W_1} T ∣ W 1 的极小多项式是p p p ,现在就是证明我们这么找出W 1 W_1 W 1 后一定存在满足条件的不变子空间W W W .
考虑构造线性泛函f f f 满足f T k v = 1 fT^kv=1 f T k v = 1 ,其他的基全设0 0 0 .
则 W = ⋂ i = 0 k null f T i W=\bigcap_{i=0}^k \operatorname{null} fT^i W = ⋂ i = 0 k null f T i .
设φ : F n → F k , φ ( v ) = [ f v , f T v , … f T k v ] \varphi:F^n\to F^k,\varphi(v)=[fv,fTv,\ldots fT^kv] φ : F n → F k , φ ( v ) = [ f v , f T v , … f T k v ] ,W = null φ W=\operatorname{null} \varphi W = null φ .显然是线性变换.
显然W W W 是不变子空间.我们要证明W 1 ⊕ W = V W_1\oplus W=V W 1 ⊕ W = V .
考虑φ ∣ W 1 : F k → F k \varphi|_{W_1}:F^k\to F^k φ ∣ W 1 : F k → F k ,且显然是双射(直接把T k v T^kv T k v 带进去就是线性无关的),于是零空间只有0 0 0 ,于是W 1 ∩ W = { 0 } W_1\cap W=\{0\} W 1 ∩ W = { 0 } .
接下来证W 1 + W = V W_1+W=V W 1 + W = V ,你再考虑φ \varphi φ 的值域既然是k k k ,零空间自然是n − k n-k n − k 维,于是dim W + dim W 1 = n \dim W+\dim W_1=n dim W + dim W 1 = n ,得证.
这就结束了.
[think] 显然f f f 就是提取那一维系数,那φ \varphi φ 是什么玩意呢?它是一个用零空间区分W 1 W_1 W 1 和W W W 的映射,
复习时,以及考试后看到的题
显然阶数高于秩的子式一定都是0 0 0 ,所以只要找到一个等于秩的非零主子式.
注意到我们把实对称矩阵分解到合同规范性后得到P D P T PDP^T P D P T ,其中P P P 是正交阵,D D D 是前r r r 个对角线元素非0 0 0 ,其余元素全为0 0 0 的矩阵.
注意到最终矩阵的某个子式A ( S , T ) A(S,T) A ( S , T ) ,行集合为S S S ,列集合为T T T ,实际上就是由P P P 只保留S S S 集合这些行得到X X X 和P T P^T P T 只保留T T T 集合这些列得到Y Y Y 后的X D Y XDY X D Y .这是矩阵乘法的规则决定的.
那么P P P 的前r r r 列中我们一定可以取出r r r 个线性无关行组成集合S S S ,这r × r r\times r r × r 的子式非0 0 0 ,同时在P T P^T P T 中取对应的S S S 中的列,那么最后乘出来矩阵的行列式就是这三个矩阵行列式相乘非0 0 0 .且恰好对应了最终矩阵中S S S 中行列组成的主子式.
[think] 关键是矩阵乘法中,A B AB A B 的某个子式就是A A A 保留对应行B B B 保留对应列组成的这个性质.
∀ D , N ∈ C n × n , D is diagonalizable, N is nilpotent , s . t . N D = D N ∃ f , g ∈ C [ x ] , f ( D + N ) = D , g ( D + N ) = N \begin{gathered}
\forall D,N\in C^{n\times n},D \text{ is diagonalizable, }N \text{ is nilpotent}, \\
\ s.t.\ ND=DN \\
\exists f,g\in C[x],f(D+N)=D,g(D+N)=N
\end{gathered} ∀ D , N ∈ C n × n , D is diagonalizable, N is nilpotent , s . t . N D = D N ∃ f , g ∈ C [ x ] , f ( D + N ) = D , g ( D + N ) = N
设A = D + N A=D+N A = D + N .
考虑可对角化意味着可同时上三角化,此时A = D + N A=D+N A = D + N ,N N N 的对角线上全0 0 0 ,于是A A A 和D D D 特征值相同.
注意到对D D D 的特征空间 V 1 = null ( D − λ I ) V_1=\operatorname{null} (D-\lambda I) V 1 = null ( D − λ I ) 上( D − λ I ) ∣ V 1 (D-\lambda I)|_{V_1} ( D − λ I ) ∣ V 1 是幂零的,于是V 1 V_1 V 1 中( D − λ I + N ) ∣ v 1 (D-\lambda I+N)|_{v_1} ( D − λ I + N ) ∣ v 1 也是幂零的.他是A A A 的广义本征空间G λ G_\lambda G λ 的一部分.而A A A 的不同广义本征是不交的,所以这些V V V 是不交的.又因为这些V V V 的直和是全空间,所以只能是D D D 的特征空间就是A A A 的广义特征空间.这很奇妙啊!
那么每个A A A 的广义特征空间G λ G_\lambda G λ 中有D v = λ v Dv=\lambda v D v = λ v ,就是f ( A ) ∣ G λ = λ f(A)|_{G_\lambda}=\lambda f ( A ) ∣ G λ = λ .
广义本征空间是 null ( A − λ I ) k \operatorname{null} (A-\lambda I)^k null ( A − λ I ) k ,所以实际在说 f ≡ λ ( m o d ( x − λ ) k ) f\equiv \lambda \pmod {(x-\lambda)^k} f ≡ λ ( mod ( x − λ ) k ) .
而你可以拿中国剩余定理构造.
最后g = x − f g=x-f g = x − f 是显然的.
{ A , B , C , D ∈ R n × n , A B T , C D T is symmetric A D T − B C T = I n ⟹ A T D − C T B = I n \begin{gathered}
\begin{cases}
A,B,C,D\in R^{n \times n}, \\
AB^T,CD^T \text{ is symmetric} \\
AD^T-BC^T=I_n \\
\end{cases} \\
\implies A^TD-C^TB=I_n
\end{gathered} ⎩ ⎨ ⎧ A , B , C , D ∈ R n × n , A B T , C D T is symmetric A D T − B C T = I n ⟹ A T D − C T B = I n
A B T AB^T A B T 对称这种看起来就很诡异啊,你应该想到把条件写开放在这里:
{ A B T − B A T = 0 C D T − D C T = 0 A D T − B C T = I n \begin{gathered}
\begin{cases}
AB^T-BA^T=0 \\
CD^T-DC^T=0 \\
AD^T-BC^T=I_n
\end{cases}
\end{gathered} ⎩ ⎨ ⎧ A B T − B A T = 0 C D T − D C T = 0 A D T − B C T = I n
然后集中注意力注意到
[ A B C D ] [ D T − B T − C T A T ] = [ I n 0 0 I n ] \begin{gathered}
\begin{bmatrix} A&B\\C&D \end{bmatrix}
\begin{bmatrix} D^T&-B^T\\-C^T&A^T \end{bmatrix}
=\begin{bmatrix} I_n&0\\0&I_n \end{bmatrix}
\end{gathered} [ A C B D ] [ D T − C T − B T A T ] = [ I n 0 0 I n ]
于是由于左逆也是右逆,就有
[ D T − B T − C T A T ] [ A B C D ] = [ I n 0 0 I n ] \begin{gathered}
\begin{bmatrix} D^T&-B^T\\-C^T&A^T \end{bmatrix}
\begin{bmatrix} A&B\\C&D \end{bmatrix}
=\begin{bmatrix} I_n&0\\0&I_n \end{bmatrix}
\end{gathered} [ D T − C T − B T A T ] [ A C B D ] = [ I n 0 0 I n ]
这样就能提取出A T D − C T B = I n A^TD-C^TB=I_n A T D − C T B = I n .
[think] 观察形式想到分块矩阵.
A 2 + B 2 = 2 A B ⟹ det A = det B \begin{gathered}
A^2+B^2=2AB \\
\implies \det A=\det B
\end{gathered} A 2 + B 2 = 2 A B ⟹ det A = det B
我觉得第一步就很困难啊!注意到
A ( A − B ) = ( A − B ) B \begin{gathered}
A(A-B)=(A-B)B
\end{gathered} A ( A − B ) = ( A − B ) B
设X = A − B X=A-B X = A − B ,如果X X X 可逆则A , B A,B A , B 相似,显然.现在考虑X X X 不可逆.
则存在v ≠ 0 , X v = 0 v\ne 0,Xv=0 v = 0 , X v = 0 .X B v = A X v = 0 XBv=AXv=0 X B v = A X v = 0 则 U = null X U=\operatorname{null} X U = null X 在B B B 下不变.且这个空间中A = B A=B A = B 啊,所以X X X 在A A A 下也不变.且其中det A ∣ U = det B ∣ U \det A|_U=\det B|_U det A ∣ U = det B ∣ U .
而 V − U V-U V − U 中,X X X 可逆,A , B A,B A , B 相似?这不对!这其中X X X 单而不一定满.正确的做法是归纳,A ∣ V − U , B ∣ V − U A|_{V-U},B|_{V-U} A ∣ V − U , B ∣ V − U 仍然满足题目中那个式子,于是就归纳下去做完了.
[think] 知乎答案说对给定某个等式的题这是套路.
Z 2 × 2 Z^{2\times 2} Z 2 × 2 的可逆矩阵(行列式为1 1 1 )可由下面两个生成元生成:
A = [ 1 1 0 1 ] B = [ 0 − 1 1 0 ] \begin{gathered}
A=\begin{bmatrix} 1&1\\0&1 \end{bmatrix} \\
B=\begin{bmatrix} 0&-1\\1&0 \end{bmatrix}
\end{gathered} A = [ 1 0 1 1 ] B = [ 0 1 − 1 0 ]
思路是对 M = [ a b c d ] M=\begin{bmatrix} a&b\\c&d \end{bmatrix} M = [ a c b d ] ,你发现可以对a , c a,c a , c 进行辗转相减,让c c c 变成0 0 0 .然后就a = d = 1 a=d=1 a = d = 1 ,最后M = A b M=A^b M = A b 了.