2025-09-15

Linear Algebra

Linear algebra

A fun question

In n-D space we can found at most n+1n+1 vector v1vn+1v_1\ldots v_{n+1} such that: ij,vivj<0\forall i\ne j,v_iv_j<0

An example in 3-D space(the one on the book)

The construction is Obviously(CH4CH_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 v1v_1, and it can be written as [1,0,0][1,0,\ldots 0](with some rotation)

so vi,i>1,viv1<0    vi,1<0    vi,1vj,1>0    vivjvi,1vj,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

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+2n+2 point could be divided to two convex hull with intersection.

n+2n+2 vector must be dependent: ci s.t. i=1n+2cixi=0;ixi=0\exists c_i \ s.t.\ \sum _{i=1}^{n+2} c_i x_i = 0; \sum_i x_i = 0(the second condition can be satisfied by add another all-1 dimension).

so divide the vector by sign of cic_i we got:v=iAcixi=jBcjxjv=\sum_{i\in A} c_i x_i = \sum_{j\in B} c_jx_j,so divide the eqution by iAci\sum_{i\in A} c_i, you get one point(vv) in the intersection.

We noticed 0<v2=(iAcixi)(jBcjxj)<00<v^2=(\sum_{i\in A} c_i x_i) \cdot (\sum_{j\in B} c_jx_j)<0, contradiction!

Operator's Left Inverse and Right Inverse

Existance

对算子TT来说,左逆存在等价于右逆存在.

Proof 1

注意到左逆存在等价于TT是单射,右逆存在等价于TT是满射.

又因为TT单射等价于TT满射等价于TT双射,所以左逆存在等价于右逆存在且两个逆一定相同.

Proof 2

左逆推出TT分解成初等行变换矩阵,然后一个一个取逆推右逆.

Union of finite count of proper sub-space isn't V

F is inifinite FieldU1Uk,Ui is subspace of V(F)Vi=1kUi\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}

考虑归纳,k=1k=1成立,设k1k-1个不行,反证,那么可以取vi=1k1Uiv\notin \bigcup_{i=1}^{k-1} U_i,必有vUkv\in U_k.

再取uUku\notin U_k,令 S={u+iviR}S=\{ u+iv \vert i \in R\},则对j<kj<k,每个UjU_j至多包含一个SS中的向量(否则vUjv\in U_j),而UkU_k必然没有SS中向量(否则uUku\in U_k),而SS中有无限个向量,于是SS不可能被他们包含,得证.

[think] 而这个定理甚至在有限域下有反例.

CR分解

Row Reduced Echelon Form

  • 每一行的首个非零元素是 1,这个元素称为“主元”(pivot)。
  • 每个主元所在列的其他元素都是 0,也就是说主元是该列唯一的非零元素。
  • 每个主元都在其所在行的右边位置,相对于上一行的主元。
  • 所有零行(即整行都是 0)都排在非零行的下面

AA的Row Reduced Echelon Form记为rref(A)\operatorname{rref}(A)

B=rref(A)B=\operatorname{rref}(A)的主元所在列构成集合S=ia pivot is in iS={i\vert \text{a pivot is in } i},设A=[a1an]A=[a_1\ldots a_n],则C=[aiiS]C=[a_i \vert i \in S],R=B1 rankA,1 mR=B_{1~\operatorname{rank}A,1~m}(即去除所有全00行),满足A=CRA=CR.

我们把它看成用RR组装CC的列,那么进行初等行变换不改变列之间的线性关系.

而消元之后呢,看列的话主元所在列显然是标准基,那命题是显然的了.

秩分解

An×m=Pn×nBn×mQm×mA_{n\times m}=P_{n\times n}B_{n\times m}Q_{m\times m},其中BB为只有左上角是一个 rankA×rankA\operatorname{rank} A\times \operatorname{rank} A 的单位矩阵其他位置全是00.

AA做行变换+列变换消元易得.

能不能换个视角,这个是不是在说,对 TL(V,W)T\in \mathcal L( V , W ),存在一个VV的一个基v1vnv_1\ldots v_n,WW的一个基w1wmw_1\ldots w_m使得 Tvi=[irankT]wiTv_i=[i\le \operatorname{rank} T]w_i.

那么先构造vv,我们先找一个 nullA\operatorname{null} A的基vnr+1vnv_{n-r+1}\ldots v_n,然后再任意扩充出剩下的v1vnv_1\ldots v_n.

ww,显然Tv1TvrTv_1\ldots Tv_r线性无关,再扩充wr+1wmw_{r+1}\ldots w_m得到一组基.

显然这组基满足需求.

[think] 还是对矩阵基变换理解不到位.

注意我们把上面那个再搞一搞:PP的后mrm-r列是没用的,QQ的后mrm-r行是没用的,都丢到会得到P=CRP=CR.

P=CR=[c1cr][r1TrnT]T=i=1nciriTP=CR=[c_1\ldots c_r] [r_1^T\ldots r_n^T]^T=\sum _{i = 1} ^{n} c_ir_i^T,其中每个ciriTc_ir_i^T秩为11.这就是秩分解的名字.

Ax=B有解

x,Ax=b    rankA=rank[A,b]    brangeA\begin{gathered} \exists x,Ax=b \\ \iff \operatorname{rank} A=\operatorname{rank} [A,b] \\ \iff b\in \operatorname{range} A \end{gathered}

解唯一等价于nullA=0\operatorname{null} A=0等价于n=rankAn=\operatorname{rank} A

A quiz problem

A is a real matrix,ATAu=0    Au=0\begin{gathered} A \text{ is a real matrix} , \\ A^TAu=0 \implies Au=0 \end{gathered}
A=M(T)TTu=0    v,<v,TTu>=0    v,<Tv,Tu>=0    Tu(rangeT)TurangeT<Tu,Tu>=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}

[think] 被这个题击败了,当时只想到用 nullT=(rangeT)\operatorname{null} T^*=(\operatorname{range} T)^{\perp}了,但其实是可以简单翻译过来的.

伴随和共轭转置的关系其实是显然的,内积上伴随的性质也是显然的,所以基础操作没必要用结论.做题的时候错误的感觉算子伴随和矩阵转置的距离过远(因为done right中证明是表示成规范正交基然后拆开用内积的性质,但是不看那套框架的话其实是显然的,另外对UU=VU\oplus U^\perp=V的证明掌握不好).同时左零空间.

总结就是记住了几何那边的结论但没有很好的联系到代数这边.

Several Inequations about Rank

rankA+BrankA+rankBrankABminrankA,rankBAm×nBn×s=0    rankA+rankBn\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}

Obviously

rankABrankAm×n+rankBn×sn\begin{gathered} \operatorname{rank} AB\ge \operatorname{rank} A_{m\times n}+\operatorname{rank} B_{n\times s}-n \end{gathered}

Sol 1

矩阵分解:

A=P1[Ir1,00,0]Q1B=P2[Ir2,00,0]Q2AB=P1[Ir1,00,0]Q1P2[Ir2,00,0]Q2\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}

显然P1,Q2P_1,Q_2不影响最终的秩直接扔了,而设D=Q1P2=[D1,D2D3,D4]D=Q_1P_2=\begin{bmatrix} D_1,D_2 \\ D_3,D_4 \end{bmatrix},那么你发现乘完只剩下D1D_1.

而删去矩阵一行或一列秩最多减少11,D1D_1看成DD删掉了 nrankA+nrankBn-\operatorname{rank} A + n-\operatorname{rank} B 行或列得到的.同时 rankD=n\operatorname{rank} D=n,得证.

Sol 2

考虑

C=[In,00,AB]\begin{gathered} C=\begin{bmatrix} I_n,0 \\ 0,AB \end{bmatrix} \end{gathered}

显然 rankAB+n=rankC\operatorname{rank} AB+n=\operatorname{rank} C

对它做行变换可以得到

C[In,0A,AB][In,BA,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}

而观察这个DD容易发现 rankDrankA+rankB\operatorname{rank} D\ge \operatorname{rank} A+\operatorname{rank} B,于是得证

Sol 3

考虑我们要证明 dimrangeABdimrangeBdimnullA\dim \operatorname{range} AB\ge \dim \operatorname{range} B-\dim \operatorname{null} A

考虑为什么 rangeBrangeAB\operatorname{range} B\ne \operatorname{range} AB,是因为 rangeB\operatorname{range} B中的不同元素被合成了一个,而这个合成相当于把 差是 nullA\operatorname{null} A中的元素的多个元素合成一个.所以有

dim(rangeB)/(rangeBnullA)=dimrangeA\begin{gathered} \dim (\operatorname{range} B)/(\operatorname{range} B\cap \operatorname{null} A)=\dim \operatorname{range} A \end{gathered}

显然交集小于 nullA\operatorname{null} A ,得证.

rankAC+rankCBrankC+rankACB\begin{gathered} \operatorname{rank} AC+\operatorname{rank} CB\le \operatorname{rank} C+\operatorname{rank} ACB \end{gathered}

这个结论可以直接由上一个的Sol3弄出来,考虑

[C,00,ACB]\begin{gathered} \begin{bmatrix} C,0 \\ 0,ACB \end{bmatrix} \end{gathered}

可以简单消元变成

[C,CBAC,0]\begin{gathered} \begin{bmatrix} C,CB \\ AC,0 \end{bmatrix} \end{gathered}

于是直接得证.

[think] 学会这种拼成空间再分块矩阵消元的套路.

A2=I    rank(AI)+rank(A+I)=n\begin{gathered} A^2=I \\ \implies \operatorname{rank} (A-I)+\operatorname{rank} (A+I)=n \end{gathered}
(AI)(A+I)=0    rank(AI)+rankA+In(A+I)(AI)=2I    rank(A+I)+rank(AI)nQ.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}

Eular Formula

对平面图Graph(n,m)\text{Graph}(n,m)FF个面(不含最外面),证明

首先考虑无向图的 Incidence Matrix MM,容易注意到MM中的若干行线性无关等价于这个导出子图无环.

于是看出 rankM=nc\operatorname{rank} M=n-c,cc为连通块个数.

又能看出 vnullMTv\in \operatorname{null} M^T等价于vv中的若干条边串成若干个环,会发现 dimnullMT=F\dim \operatorname{null} M^T=F

An Ex Problem

M=[A,C0,B]rankM=rankA+rankB    X,Y:AX+YB=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}

首先右推左是显然的.直接消元一下就好了.考虑左推右.

考虑

A:U1V1,B:U2V2,C:U2V1M:UV\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}

分解

rangeM=rangeBWW={[v,0][v,0]rangeM}\begin{gathered} \operatorname{range} M=\operatorname{range} B\oplus W \\ W=\{ [v,0] \vert [v,0]\in \operatorname{range} M \} \\ \end{gathered}

那么因为 [v,0]rangeM[v,0]\in \operatorname{range} M,则必然是 [v,0]=M[u1,u2]T[v,0]=M[u_1,u_2]^T,一定是Bu2=0,Au1+Cu2=vBu_2=0,Au_1+Cu_2=v.

于是 W=rangeA+C(nullB)W=\operatorname{range} A + C(\operatorname{null} B).又因为 dimW=dimrangeA\dim W=\dim \operatorname{range} A,于是有 C(nullB)rangeAC(\operatorname{null} B) \subset \operatorname{range} A.

我们再分解 U2=nullBU3U_2=\operatorname{null} B \oplus U_3,此时注意到BBU3U_3rangeB\operatorname{range} B是双射,存在Y,u,YBu=CuY',\forall u,YBu=Cu.然后通过扩充基并任意取值将YY'的定义域扩充到V2V_2得到YY.

于是(CYB)u(C-YB)u对任意uU3u\in U_300,于是 range(CYB)=(CYB)nullB=C(nullB)\operatorname{range} (C-YB)=(C-YB)\operatorname{null} B= C(\operatorname{null} B).

现在只考虑 unullBu\in \operatorname{null} B,显然v,Cu=Av\exists v,Cu=Av,那么对 nullB\operatorname{null} B的一组基u1uku_1\ldots u_k这样确定v1vkv_1\ldots v_k,就可以构造Xui=viX'u_i=v_i满足Cu=AXvCu=AX'v.再用同样的方法扩充基并任意取值将XX'的定义域扩充到U2U_2得到XX.

于是C=YB+AXC=YB+AX

[think] 感觉得到 C(nullB)=AC(\operatorname{null} B)=A这里是容易的.然后这里进行不下去,想到 nullB\operatorname{null} B去分解也是自然的. 分解后就要想办法把 nullB\operatorname{null} B之外的影响消掉,就用了CYBC-YB.而若 rangeArangeB\operatorname{range} A\subset \operatorname{range} B那么Ax=BTxAx=BTx是显然的.

投影

向一个向量投影

a,bRn<a,b><a,a>a=(aTbaaTa)=(aaTaTa)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}

向一个平面投影(平面是C(A)C(A))

p=A(ATA)1ATb\begin{gathered} p=A(A^TA)^{-1}A^Tb \end{gathered}

考虑bb的投影pC(A)p\in C(A)(bp)C(A)(b-p)\in C(A)^\perp,于是bpN(AT)b-p \in N(A^T).

于是ATb=ATpA^Tb=A^Tp,又pC(A)    x,Ax=pp\in C(A) \implies \exists x,Ax=p.

于是ATb=ATAxA^Tb=A^TAx,AA一定可以用一个满秩的,于是除过去.于是得证.

ATAx=ATb\begin{gathered} A^TAx=A^Tb \end{gathered}

一定有解

rank(ATA,ATb)=rank(AT(A,b))rankA=rankAAT\begin{gathered} \operatorname{rank} (A^TA,A^Tb) \\ =\operatorname{rank} (A^T(A,b))\le \operatorname{rank} A \\ =\operatorname{rank} AA^T \end{gathered}

所以这个证明是依赖实数的.

[think] 复数你应该把ATA^T换成AT\overline{A^T},或者说这个定理本来就应该是AT\overline{A^T}的.

P2=P,P=P    P is a projection\begin{gathered} P^2=P,P^*=P \implies P \text{ is a projection} \end{gathered}

显然那PP只能说 rangeP\operatorname{range} P 的投影.

只需证明 u,(Puu)(rangeP)\forall u,(Pu-u)\in (\operatorname{range} P)^\perp.

v,<Puu,Pv>=<PPuPu,v>=<PuPu,v>=<0,v>    Puu(rangeP)\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}

得证!

QR分解

A,Q is orthogonal matrix,R is upper triangle matrixs.t.A=QR\begin{gathered} \forall A,\exists Q \text{ is orthogonal matrix},R \text{ is upper triangle matrix} \\ s.t.\\ A=QR \end{gathered}

考虑AA可以看成把标准基变成a1ana_1\ldots a_n,那我们把a1ana_1\ldots a_n这组基用Gram-Schmidt变成b1bnb_1\ldots b_n,问题就可以变成先把标准基变成bb,再变成aa,其中第一步是等距同构,第二步中我们知道a1aia_1\ldots a_ib1bib_1\ldots b_i张成空间相同,所以第二步是上三角.

determinance

detAdetB=detAB\begin{gathered} \det A\det B=\det AB \end{gathered}

Sol1:分块矩阵.

Sol2:都拆成初等变换矩阵再乘.

Sol3:考虑定义函数 α(B)=detBAdetA\alpha(B)=\dfrac{\det BA}{\det A},容易验证它满足行列式三条公里,于是α(B)=detB\alpha(B)=\det B

Laplace Theorem

S[1,n]ZdetA=T[1,n]Z,T=SA(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}

其中A(S,T)A(S,T)表示子式,C(S,T)C(S,T)表示代数余子式

注意到你就是把SS行对应的元素钦定的时候的某个组合,我们可以先用tTti=1Ti+sSsi=1SitTt+sSs(mod2)\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次交换把这些行列顺序不变的换到前 S=T\vert S \vert=\vert T \vert 行列,则最终符号显然就是此时的符号(即两个行列式内部的符号)再乘上交换操作的符号,于是得证.

detMMi,j={a,i=jb,ij\begin{gathered} \det M \\ M_{i,j} = \begin{cases} a,i=j \\ b,i\ne j \end{cases} \end{gathered}

Sol1

考虑消元,注意到每行的和是一样的,所以我们把第一列变成所有列的和,然后用它去消后面的,模拟一下即可.

Sol2

考虑M=(ab)I+xBM=(a-b)I+xB,其中B=(b)i,jB=(b)_{i,j},应用公式

ddet(M)=trace(adjAdA)\begin{gathered} \operatorname{d} \det(M)=\operatorname{trace} (\operatorname{adj} A \operatorname{d} A) \end{gathered}

右边三项全是好算的,于是可以把导数对xx积分积回去.

Sol3

考虑M=B+(ab)IM=B+(a-b)I,若Bv=λvBv=\lambda v,则显然有(B+(ab)I)v=(λ+ab)v(B+(a-b)I)v=(\lambda+a-b)v,于是求出BB的特征值,就可以直接得到MM的特征值算行列式.

rankA=a    {T=S>a,detAS,T=0S=T=a,detAS,T0\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}

首先如果存在detAS,T=0\det A_{S,T}=\ne 0,那么SS对应的这些行必然是满秩的,所以 rankAa\operatorname{rank} A\ge a.

而如果 rankA>a\operatorname{rank} A>a,那么选出aa个线性无关行,再从这里面找到aa个线性无关列,这就是一个行列式非00的矩阵.

于是得证.

rankA=n1    rankC=1rankA<n1    rankC=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}

由上一个conclusion,第二行是显然的(全0).

对于第一行,考虑ACTAC^T=0,于是 ACT=0,rankA+rankCTnAC^T=0,\operatorname{rank} A+\operatorname{rank} C^T\le n ,又有 rankC0\operatorname{rank} C\ne 0 因为至少有一个非00.

Am×n,Bn×m,mn    detAB=S=m,S[1,n]detA[1,m]Z,SdetBS,[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}

考虑矩阵

det[InBA0]=det[InB0AB]\begin{gathered} \det \begin{bmatrix} I_n&B \\ A&0 \end{bmatrix}=\det \begin{bmatrix} I_n&B \\ 0&-AB \end{bmatrix} \end{gathered}

右边的行列式是(1)mdet(AB)(-1)^m\det (AB),考虑左边用拉普拉斯定理,则显然你的子式的列只能选AA里面的,而它的余子式就是BB左边再加上若干列,若第ii列没被子式选,则这列一定要选ii行的元素,于是BB中恰好只有SS中的行能选.

对于符号,Laplace中的是iSi+m(n+1+n+m)2\sum_{i\in S} i+\dfrac{m(n+1+n+m)}2,算代数余子式行列式的时候还有一个iSi(nm)(nm+1)2\sum_{i\notin S} i-\dfrac{(n-m)(n-m+1)}2,然后这些加起来是同余mm的.

这样符号一乘正好是00.

行列式与导数:容易注意到行列式每行不会有两个乘到一起,而最终的行列式是nn行各取一个数乘起来,于是遵循导数求导乘积的法则,或者说把每一行求导其他行不变再加起来.

Another Quiz Problem

Am×n,Bn×m s.t. ABA=A\begin{gathered} \forall A_{m\times n},\exists B_{n\times m}\ s.t.\ ABA=A \end{gathered}

Sol 1

A=P[Ik000]Q\begin{gathered} A=P\begin{bmatrix} I_k&0 \\ 0&0 \end{bmatrix}Q \end{gathered}

于是让

B=Q1[Ik000]P1\begin{gathered} B=Q^{-1}\begin{bmatrix} I_k&0 \\ 0&0 \end{bmatrix}P^{-1} \end{gathered}

Sol 2

考虑一个任意情况,对任意集合X,YX,Y和任意映射ff,你可以找到一个gg使得fgf=ff\circ g\circ f=f

显然的,因为你可以让ggff的像映射到任意一个原像.

那么这个显然的在说什么呢?它实际上在说,对任意一个ff,我们可以找到X,YX,Y各自的一个子集X,YX',Y',使得ff限制在XYX'\to Y'是双射.于是gg是这上面的逆,而这之外的随意映射就可以满足fgf=ffgf=f.

然后你再看第一个证明,那么中间那个有IkI_k的矩阵实际上就是,IkI_k对应了XYX'\to Y'的双射部分,这也是上面的做法为啥有道理.

Circulant Matrix's det

a0an1Ai,j=ai+j2modn    detA=i=1nj=1nwnijcj\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}

首先我们让PnP_n是循环移位P(x0,,xn)=(x1,,xn,x0)P(x_0,\ldots, x_n)=(x_1,\ldots, x_n,x_0),那么有

A=i=0n1aiPi\begin{gathered} A=\sum _{i = 0} ^{n-1} a_iP^i \end{gathered}

那么因为所有PP乘法的时候是保持不变子空间的,也就是说AA的特征值就是PP的特征值带入f(x)=i=0n1aixif(x)=\sum_{i=0}^{n-1}a_ix^i,所以只需要看PP的特征值.

考虑Px=λxPx=\lambda x,那么xi+1modn=λxix_{i+1\bmod n}=\lambda x_i,于是你乘一圈就有λn=1\lambda^n=1,能相应的构造出xi=wnil,λ=wnlx_i=w_n^{il},\lambda=w_n^l这样的特征向量和特征值.

于是PP的特征值就是单位根,带进去就得证.

A Determinance about Vandermonde plus 1

detA,Ai,j=xij+1\begin{gathered} \det A,A_{i,j}=x_i^j+1 \end{gathered}

容易想到用线性性去拆分,如果我们对列拆分,你发现你就是要求AA有一列换成全11之后的矩阵的行列式的和,我们设第ii列换掉之后的矩阵是BiB_i.

考虑Cramer's Rule,你注意到设Ay=bAy=b,其中bb是全11向量,则yi=detBidetAy_i=\dfrac{\det B_i}{\det A},于是我们实际上是求yiy_i,

如果我们这里用y=A1by=A^{-1}b我们可以得到一个结论:

det(A+uvT)=detA(1+vTA1u)\begin{gathered} \det(A+uv^T)=\det A(1+v^TA^{-1}u) \end{gathered}

然后对这个题,我们也不知道A1A^{-1},所以你仍然要想办法求yiy_i.

碰到范德蒙德矩阵一定要想多项式,于是你想到P(z)=i=1nyiziP(z)=\sum_{i=1}^n y_iz^i,则P(xi)=1P(x_i)=1,于是P(z)=Ai(zxi)P(z)=A\prod_i (z-x_i),于是你能通过比较常数项求出AA,然后再计算P(1)P(1)就能求出yi\sum y_i了.

[think] 碰到范德蒙德矩阵一定要想多项式

det(A+uvT)=detA(1+vTA1u)\begin{gathered} \det(A+uv^T)=\det A(1+v^TA^{-1}u) \end{gathered}

一个是上面提到的证法.

或者考虑

det(I+uvT)=(1+vTu)=(1+uTv)\begin{gathered} \det(I+uv^T)=(1+v^Tu)=(1+u^Tv) \end{gathered}

这个可以直接考虑uvTuv^T特征值是n1n-100和一个vTuv^Tu(对应特征向量uu).

然后det(A+AuvT)=detA(1+uTv)\det (A+Auv^T)=\det A(1+u^Tv)

所以让u=A1uu'=A^{-1}u代入,那么式子就成了det(A+uvT)=detA(1+vTA1u)\det (A+uv^T)=\det A(1+v^TA^{-1}u).即我们所证的.

about eigen values

相似矩阵有相同的特征多项式

detQAQ1λI=detQAQ1QλIQ1=detQ(AλI)Q1=detAλ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}

显然的吧,因为都是一个线性变换在不同的基下的.

对称矩阵的特征向量一定垂直

首先你可以谱定理启动.

否则考虑

v1Av2=v1λ2v2(v1Av2)T=λ1v1v2\begin{gathered} v_1Av_2=v_1\lambda_2 v_2 \\ (v_1Av_2)^T=\lambda_1 v_1v_2 \end{gathered}

于是就得证了.

特征多项式的nin-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}

观察AλIA-\lambda I显然.

实谱定理(矩阵版)

任意对称矩阵AA可以写成A=PDP1A=PDP^{-1},其中DD是对角矩阵,PP是正交矩阵.

首先复化,考虑AA的一组特征值特征向量Av=λvAv=\lambda v,则 vTAv=λvv^TA\overline{v}=\overline\lambda \vert\vert v \vert\vert,同时取共轭得这个还等于 λv\lambda \vert\vert v \vert\vert,于是能证明λR\lambda \in R.

然后我们说明了AA的特征值都是实的,我们取v1v_1是一个特征值,并将其扩充到一组标准正交基v1,v2vnv_1,v_2\ldots v_n,于是

A[v1vn]=[v1vn][λ1x0A1]\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}

其中xxA1A_1是任意向量/矩阵,而中间v1vnv_1\ldots v_n是一个正交矩阵设为BB,则右边那个分块矩阵是B1AB=BTABB^{-1}AB=B^TAB,于是这个也是对称的,从而x=0x=0,A1A_1是对称矩阵,从而可以归纳到n1n-1维.归纳得证.

这里x=0x=0就是TTUU^\perp不变,A1A_1对称就是TT限制仍然是自伴.

二次型

惯性定理 Inertia Theorem

任意二次型AA通过任意可逆矩阵PP合同变换到对角矩阵PDPTPDP^T,DD的对角线上00的个数,正数个数,负数个数都是确定的(与PP无关).

首先00的个数是确定的,因为这是矩阵的秩.

考虑在两组基下,假设它们对角型正数的个数分别为a,ba,b,负数是ra,rbr-a,r-b,不妨设a<ba<b.

那么我们不考虑00的那部分限制到rr维的空间上考虑,则你可以找到一个aa维的空间,上面值都是正的(要求所有负惯性指数位置都是负的),再找一个rbr-b上面都是负的.而a+rb>ra+r-b>r所以这俩空间有交,爆炸了.

[think] 注意到正惯性指数是最大的子空间满足这里面的值带进去都是正的.

所有顺序主子式大于00等价于正定

考虑对角化的过程,如果已经消除了左上角k×kk\times k的矩阵,接下来这些地方就不会动了.然后你消(k+1)×(k+1)(k+1)\times (k+1)的左上角,注意到我们前面已经进行的行列变换都是在前kk行进行的,所以前面所有的操作都不会改变当前(k+1)×(k+1)(k+1)\times (k+1)矩阵的行列式,直到你消完.于是每个顺序主子式对应的子矩阵消前和消后的行列式是一样的.

SVD

PCA(主成分分析)

你有若干个向量aia_i(若干个点)

先让aa都减去平均值得到aa'A=[a1an]A=[a_1'\ldots a_n'],令B=AATB=AA^T,你发现Bi,j=Var(xi,xj)B_{i,j}=\operatorname{Var}(x_i,x_j)(两维的协方差.)

AA进行SVD分解A=UDVA=UDV,则若s1s_1对应的方向v1v_1,直线tv1tv_1是最小化每个点到直线距离(垂线段距离)平方的直线.

设最优方向为单位向量nn.

L=(xin<xi,n>)2=ai2+n2<ai,n>22<ai,n><ai,n>=(ai2)(i<ai,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}

只要最大化第二项.考虑

i<ai,n>2=nT(aiaiT)n=nT(AAT)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}

这是二次型,把AATAA^T谱分解就看出显然最大值在它最大的特征值对应特征向量方向,也就是v1v_1方向.

你还可以看出这个值是ai2s12\sum a_i^2-s_1^2

SVD求伪逆

A=UΣVTA=U\Sigma V^T,A+=VΣ+UTA^+=V\Sigma^+ U^T,其中Σ+\Sigma^+是把所有非零对角线元素变成倒数再取转置.

用SVD逼近

A=UDVT=iuisiviT    minrankC=kAC=sk+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}

Am×n,Vn×n,Um×m,Dm×nA_{m\times n},V_{n\times n},U_{m\times m},D_{m\times n}

考虑对AA分解得A=UDVTA=UDV^T,则我们让DD只保留前kk个特征值,其余设为00得到DD',令C=UDVTC=UD'V^T,则显然ACA-C的SVD分解就是U(DD)VTU(D-D')V^T,其最大的奇异值恰为sk+1s_{k+1}(就是只用前kk个奇异值去逼近啦).

然后只需要证明这是最小值.考虑AA的前k+1k+1个右奇异向量(v1vkv_1\ldots v_k)张成的子空间和BB的零空间的交( dimnullB=nrankB=nk\dim \operatorname{null} B=n-\operatorname{rank} B=n-k,故必然有交)中的单位向量xx.

(AB)x2=Ax2=xTATAx=xTVD2VTx=(VTx)TD2(VTx)sk+12\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}

于是得证.

[think] 证明最小化的时候从分析矩阵结构改为找一个向量.因为这个向量显然一定可以在前k+1k+1维(这个每个维度都放大需求倍数的空间中任何一个向量也都被放大需求倍数).

多项式

有理域可约多项式

f(x)=i=0naixi is primitive     gcd(a0,,an)=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}

Gauss's Lemma

f,g are primitive     fg are primitive\begin{gathered} f,g \text{ are primitive } \implies fg \text{ are primitive} \end{gathered}

考虑反证,如果fgfg有公质因数pp,且p∤gp\not\vert g,证明pfp\vert f.

f,g,fgf,g,fg对应的系数列分别是an,bn,cna_n,b_n,c_n.

那么直接递推就好了,我们假设pgcd(f0,,fA,g0,,gB)p \vert \gcd(f_0,\ldots ,f_A,g_0,\ldots,g_B)(没有取1-1)

那么考虑cA+B+2c_{A+B+2},它一定有一项是fA+1gB+1f_{A+1}g_{B+1},且其他的项一定都被pp整除,于是你一定可以让AABB11.

这样走n+m+1n+m+1步一定可以把AABB弄满.

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)=f1(x)f2(x),f1,f2Q[x]f(x)=f_1(x)f_2(x),f_1,f_2\in Q[x].

那么你显然可以有fi(x)=gi(x)piqif_i(x)=g_i(x)\dfrac {p_i}{q_i}.且g1,g2 is primitiveg_1,g_2 \text{ is primitive}.

那么f(x)=p1p2q1q2g1(x)g2(x)f(x)=\dfrac{p_1p_2}{q_1q_2}g_1(x)g_2(x),其中g1g2g_1g_2还是primitive\text{primitive}的,q1q2f(x)=p1p2g1(x)g2(x)q_1q_2f(x)=p_1p_2g_1(x)g_2(x).

考虑两侧的系数的最大公因数相等,那只能是q1q2=p1p2q_1q_2=p_1p_2了.

{f(x)=xn+i=0n1aixi,aiZp s.t.0in1,paip2∤a0    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=f1f2f=f_1f_2设它们的系数分别是bn,cnb_n,c_n,根据上面高斯引理的证明归纳过程,pa0p|a_0,那么pb0c0p|b_0c_0,那么因为pp只能整除其中的一个,由上面那个由(A,B)(A,B)推到(A+1,B)(A+1,B)(A,B+1)(A,B+1)的过程,你有一边走不了,最终就直接得到pbnp\vert b_npcnp\vert c_n了.

但是显然p∤fp\not| f,它是有系数是11的!矛盾.

你看一下下面那个模pp的性质于是这个东西另一个证法是先模pp,那么f(x)xnf(x)\equiv x^n,于是它被拆了之后只能是xixnix^i x^{n-i},但这说明拆出来的两项常数项模pp都是00,和p2∤a0p^2\not|a_0矛盾.

pP,f(x)=i=0p1xi    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}
f(x)=xp1x1=y=x1(y+1)p1y=i=0p1yi(pi+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}

此时用上面定理.

若首一多项式ffZ[x]Z[x]上可约则ffFp[x]F_p[x]上可约

显然,因为你把Z[x]Z[x]上那两个因子分别模pp就好了.首一主要是防止你模pp的时候ff模成常数.

Smith Standard Form

Smith标准型是说环上元素的矩阵,相抵于一个只有主对角线上有d1dr,0,0d_1\ldots d_r,0\ldots,0的矩阵且满足di+1did_{i+1}\vert d_i.

定义DiD_i是所有ii阶子式的gcd,则di=DiDi1d_i=\dfrac{D_i}{D_{i-1}}.

如果我们消完了的dd有整除链的性质,那么最后是显然满足条件的,只需要证明DiD_i在初等行变换下不变:

  • 交换两行/列或给一行/列乘±1\pm 1(环上,其他的没逆)显然不影响
  • 把一行/列乘kk加到另一行/列,此时若一个子式同时包含或不包含涉及到的两行则它不变,如果只包含一个则用线性性相当于给他加上另一个子式,不改变gcd\gcd.

而只要证明存在一种消法消成整除链,你发现如果消完了不是整除链你总可以用辗转相减变成di=gcd(d1,,di)d_i=\gcd(d'_1,\ldots,d'_i)弄成整除链,于是得证.

这里dd叫不变因子(invariant factors),DD叫行列式因子(determinant factors).

然后如果这个环是UFD,定义初等因子elementfactorselement factors是把每个dd做素分解得到的每个picip_i^{c_i},不进行合并.

Jordan Standard Form

你可以把矩阵弄成只有分块对角矩阵,且每个块只有对角线上有相同的值,对角线上方的一条斜线都是11.

Proof1:初等变换

E(i,j,a)E(i,j,a)表示II的基础上(i,j)(i,j)位置为aa的初等变换矩阵.而容易注意到E(i,j,a)1=E(i,j,a)E(i,j,a)^{-1}=E(i,j,-a).

所以当我们做E(i,j,k)AE(i,j,k)E(i,j,k)AE(i,j,-k)这个相似变换我们是:

AA的第ii行加上kk倍第jj行,再给第jj列加上kk倍第ii列.

对一个上三角矩阵,我们发现进行这样一次操作,Ai,jA_{i,j}加上了k(Ai,iAj,j)k(A_{i,i}-A_{j,j}),且这一步对上三角矩阵只影响x<i,y>jx<i,y>j的位置(一个右上角).所以你可以以一定的顺序消掉这样的元素,得到一个分块对角矩阵,每个块对角线元素相同.

那么每个块都可以表示成λI+N\lambda I+N,NN是幂零的.而当你对这个子空间进行换基的时候λI\lambda I不变所以只要找NN的变换.

你发现对于NN,你仍然进行过上面这种操作就可以简单的用某一行的11把同行后面的所有数消掉,所以你从第一行开始从上往下消.如果遇到一行它的11不在副对角线你可以做交换把它换到前面来,这样都弄完了就得到标准的约旦标准型了.

Proof2

两矩阵A,BA,B相似等价于特征矩阵AλI,BλIA-\lambda I,B-\lambda I相抵.

首先如果A=PBP1A=PBP^{-1},那么

AλI=PBP1PλIP1=P(BλI)P1\begin{gathered} A-\lambda I=PBP^{-1}-P\lambda IP^{-1}=P(B-\lambda I)P^{-1} \end{gathered}

成立.

再看另一边:

考虑现在V(AλI)=(BλI)WV(A-\lambda I)=(B-\lambda I)W,V,WV,W是可逆的λ\lambda矩阵.

直接设 V=i=0mViλi,W=i=0mWiλiV=\sum _{i = 0} ^{m} V_i\lambda^i,W=\sum _{i = 0} ^{m} W_i\lambda^i.

那么

{Vm=WmVkAVk1=BWkWk1V0A=W0B\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}

然后你发现可以盯着一边看去做一个相消,就是给ViAV_iA的那个等式右乘AiA^i,则左边是00,右边是?

BWiAiWi1Ai=B(WiAi)(WiAi)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}

于是设P=WiAiP=\sum W_iA^i就有BP=PABP=PA.

那么你把V,WV,W取逆分别移到另一边是(AλI)W1=V1(BλI)(A-\lambda I)W^{-1}=V^{-1}(B-\lambda I).再走一遍能得到AQ=QBAQ=QB.

然后能看出AkQ=QBkA^kQ=QB^k,Q=W1(B)Q=W^{-1}(B),PAk=BkPPA^k=B^kP.

那么

PQ=W(A)Q=iWiAiQ=iWiQBi=iWi(jWj1Bj)Bi=ijWiWj1Bi+j=WW1(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}

[think] 你看那个神秘的求和是不自然的,你注意到如果我们对一开始的V(AλI)=(BλI)WV(A-\lambda I)=(B-\lambda I)W去代入λ=A\lambda =A会很妙,然后所谓的代入实际上对应了取模,所以算(BλI)Wmod(AλI)(B-\lambda I)W \bmod (A-\lambda I),算这个的时候拆开WW乘上BB算,因为没有交换律不保证余数可乘.

[think] 好吧上面那个think说的你直接代入就可以的,这种时候你要规定代入方向,然后(BλI)W=BW(λ)λW(λ)(B-\lambda I)W=BW(\lambda)-\lambda W(\lambda),此时你代入λ=A\lambda=A因为是右代入所以你要先把λW\lambda W变成WλW\lambda.

[think] 后面这个=I=I也可以直接通过代入看出来,或者再逐项系数展开.

不变因子组和初等因子组是双射

考虑Smith标准型,那么从矩阵角度,分块对角矩阵有两个块,问题是现在不满足整除链,那么只要每次,对所有素因子,取其最高次组成一个作为不变因子,就可以重排出一个大矩阵的整除链,于是得证.

V=W1W2V=W_1\oplus W_2,则TT的初等因子就是TW1T|_{W_1}的初等因子与TW2T|_{W_2}的初等因子的和(多重集的和).

考虑你有矩阵AA和矩阵BB都已经是Smith标准型,然后你要把它们拼起来.

此时考虑行列式因子DD,因为A,BA,B是对角阵拼完了还是对角阵,所以只有主子式不是00,此时kk阶行列式一定是AA里取若干个元素BB里取若干个元素的乘积,又因为整除链的存在,所以取的时候一定分别是A,BA,B里最低次的若干个.

那么因为是直接乘积再取gcd,现在只考虑所有不变因子dd某个特定素因子的次数就成了,新的大矩阵的DD对应的次数是原来两个小的min+卷积,因为整除链所以DD的次数是凸的,所以是闵和,所有作为斜率的那个dd的初等因子的次数一定都保留到了新的结果里,就完事了.

有理标准型

把矩阵的不变因子的友矩阵作为分块对角矩阵的块

由于上面说初等因子是直接拼,而初等因子确定不变因子,所以你这么拼是对的.

最大的不变因子是极小多项式

首先友矩阵的极小多项式是它的特征多项式,因为取[1,0,0,,0][1,0,0,\ldots,0],那么你发现友矩阵的前n1n-1次方作用到它身上是线性无关的,所以特征多项式的次数必须是nn次的.

然后你转乘上面的有理标准型,就显然.

于是这种证法说你考虑一个Jordan块正好对应了一个不变因子(xλi)n(x - \lambda_i)^n.

两个实矩阵在复数域上相似则也在实矩阵上相似

首先这说明复数域上,他俩不变因子相同,所以行列式因子相同.又因为行列式因子一定是实的,所以结束.

准素分解

TT的极小多项式是m(x)=ipi(x)cim(x)=\prod_i p_i(x)^{c_i},pip_i是互不相同的不可约多项式,则V=iWiV=\oplus_i W_i,且 Wi=nullpici(T)W_i=\operatorname{null} p_i^{c_i}(T).

Sol1:(from Gemini and me )

不变性是显然的,我们发现对一个WiW_i,其极小多项式是picip_i^{c_i},考虑其中一个vv,pici(T)v=0p_i^{c_i}(T)v=0,那么显然pici(T)Tv=Tpici(T)v=0p_i^{c_i}(T)Tv=Tp_i^{c_i}(T)v=0.

let fi=m(x)pici(x)if uW1W2WnuW1    p1c1(T)u=0uW2Wn    fi(T)u=0fip1c1s(x),t(x),s(x)p1c1(x)+t(x)fi(x)=1s(T)p1c1(T)u+t(T)fi(T)u=Iu=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}

这说明是直和,再说明直和是整个空间:

gcd(f1,,fn)=1    ai,iai(x)fi(x)=1    iai(T)fi(T)=Ilet Ei=ai(T)fi(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}

显然EiEj(T)=0E_iE_j(T)=0,又因为

let Ei(v)=wpici(T)w=pici(T)ai(T)fi(T)v=ai(T)m(T)=0    rangeEiWi(jEj)Ei=Ei    Ei2=EiwWi,jEj(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}

而你Ej,jiE_j,j\ne i的像空间的直和与WiW_i无交,所以一定只能是Ei(w)=wE_i(w)=w

这样实际证明了EiE_i就是像每个空间的投影,且和为II,于是得证.

Sol2: from teacher

Lemma

(f(x),g(x))=1    nullfg(T)=nullf(T)nullg(T)\begin{gathered} (f(x),g(x))=1 \implies \operatorname{null} fg(T)=\operatorname{null} f(T) \oplus \operatorname{null} g(T) \end{gathered}

首先,裴蜀定理,a(x)f(x)+b(x)y(x)=1a(x)f(x)+b(x)y(x)=1,于是如果 unullf(T)nullg(T)u\in \operatorname{null} f(T)\cap \operatorname{null} g(T),则a(T)f(T)u+b(T)y(T)u=Iu=0a(T)f(T)u+b(T)y(T)u=Iu=0,得u=0u=0.

然后考虑,左边显然包含右边,只要证明右边包含左边,注意到 f(T)unullg(T),g(T)unullf(T)f(T)u\in \operatorname{null} g(T),g(T)u\in \operatorname{null} f(T),又有a(T)f(T)u+b(T)g(T)u=ua(T)f(T)u+b(T)g(T)u=u,于是拆成u=a(T)f(T)u+b(T)g(T)uu=a(T)f(T)u+b(T)g(T)u即证.

于是直接证完了.

[think] 观察证明复杂度可知,你还是先考虑两个比较好()另外是要看出本质条件是互素.

对任意线性变换TT,v\exists v使得任意满足p(T)v=0p(T)v=0pp是极小多项式的倍式.

首先如果p=picip=p_i^{c_i},那么是显然的,否则pp就不是这个了.

否则用准素分解,则你得到若干个picip_i^{c_i},每个对应一个向量viv_i,考虑直接把它们加起来v=iviv=\sum_i v_i,那么:

0=g(T)ivi=ig(T)vi0=g(T)\sum_i v_i=\sum_i g(T)v_i,因为后面每一项分别在WiW_i中,所以一定有g(T)vi=0g(T)v_i=0,所以在每个子空间中的极小多项式picip_i^{c_i}都是g(T)g(T)的因子,于是mgm|g,又因为gg是整个的零化多项式所以gmg|m.得证.

循环分解

对线性变换TT,找到vv使得 g(T)v=0    pgg(T)v=0 \implies p|g,其中pp为极小多项式.然后找到最大的kk使得v,Tv,Tkvv,Tv,\ldots T^kv线性无关,设 W1=span(v,Tv,Tkv)W_1=\operatorname{span}( v,Tv,\ldots T^kv ),设W1W=VW_1\oplus W=V,则继续对WW重复就得到一串不变子空间,且极小多项式构成整除链.

前面我们已经证明了一定有这样的vv,同时容易看出TW1T|_{W_1}的极小多项式是pp,现在就是证明我们这么找出W1W_1后一定存在满足条件的不变子空间WW.

考虑构造线性泛函ff满足fTkv=1fT^kv=1,其他的基全设00.

W=i=0knullfTiW=\bigcap_{i=0}^k \operatorname{null} fT^i.

φ:FnFk,φ(v)=[fv,fTv,fTkv]\varphi:F^n\to F^k,\varphi(v)=[fv,fTv,\ldots fT^kv],W=nullφW=\operatorname{null} \varphi.显然是线性变换.

显然WW是不变子空间.我们要证明W1W=VW_1\oplus W=V.

考虑φW1:FkFk\varphi|_{W_1}:F^k\to F^k,且显然是双射(直接把TkvT^kv带进去就是线性无关的),于是零空间只有00,于是W1W={0}W_1\cap W=\{0\}.

接下来证W1+W=VW_1+W=V,你再考虑φ\varphi的值域既然是kk,零空间自然是nkn-k维,于是dimW+dimW1=n\dim W+\dim W_1=n,得证.

这就结束了.

[think] 显然ff就是提取那一维系数,那φ\varphi是什么玩意呢?它是一个用零空间区分W1W_1WW的映射,

复习时,以及考试后看到的题

实对称矩阵的秩等于最高阶非零主子式的阶数.

显然阶数高于秩的子式一定都是00,所以只要找到一个等于秩的非零主子式.

注意到我们把实对称矩阵分解到合同规范性后得到PDPTPDP^T,其中PP是正交阵,DD是前rr个对角线元素非00,其余元素全为00的矩阵.

注意到最终矩阵的某个子式A(S,T)A(S,T),行集合为SS,列集合为TT,实际上就是由PP只保留SS集合这些行得到XXPTP^T只保留TT集合这些列得到YY后的XDYXDY.这是矩阵乘法的规则决定的.

那么PP的前rr列中我们一定可以取出rr个线性无关行组成集合SS,这r×rr\times r的子式非00,同时在PTP^T中取对应的SS中的列,那么最后乘出来矩阵的行列式就是这三个矩阵行列式相乘非00.且恰好对应了最终矩阵中SS中行列组成的主子式.

[think] 关键是矩阵乘法中,ABAB的某个子式就是AA保留对应行BB保留对应列组成的这个性质.

D,NCn×n,D is diagonalizable, N is nilpotent, s.t. ND=DNf,gC[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}

A=D+NA=D+N.

考虑可对角化意味着可同时上三角化,此时A=D+NA=D+N,NN的对角线上全00,于是AADD特征值相同.

注意到对DD的特征空间 V1=null(DλI)V_1=\operatorname{null} (D-\lambda I)(DλI)V1(D-\lambda I)|_{V_1}是幂零的,于是V1V_1(DλI+N)v1(D-\lambda I+N)|_{v_1}也是幂零的.他是AA的广义本征空间GλG_\lambda的一部分.而AA的不同广义本征是不交的,所以这些VV是不交的.又因为这些VV的直和是全空间,所以只能是DD的特征空间就是AA的广义特征空间.这很奇妙啊!

那么每个AA的广义特征空间GλG_\lambda中有Dv=λvDv=\lambda v,就是f(A)Gλ=λf(A)|_{G_\lambda}=\lambda.

广义本征空间是 null(AλI)k\operatorname{null} (A-\lambda I)^k ,所以实际在说 fλ(mod(xλ)k)f\equiv \lambda \pmod {(x-\lambda)^k}.

而你可以拿中国剩余定理构造.

最后g=xfg=x-f是显然的.

{A,B,C,DRn×n,ABT,CDT is symmetricADTBCT=In    ATDCTB=In\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}

ABTAB^T对称这种看起来就很诡异啊,你应该想到把条件写开放在这里:

{ABTBAT=0CDTDCT=0ADTBCT=In\begin{gathered} \begin{cases} AB^T-BA^T=0 \\ CD^T-DC^T=0 \\ AD^T-BC^T=I_n \end{cases} \end{gathered}

然后集中注意力注意到

[ABCD][DTBTCTAT]=[In00In]\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}

于是由于左逆也是右逆,就有

[DTBTCTAT][ABCD]=[In00In]\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}

这样就能提取出ATDCTB=InA^TD-C^TB=I_n.

[think] 观察形式想到分块矩阵.

A2+B2=2AB    detA=detB\begin{gathered} A^2+B^2=2AB \\ \implies \det A=\det B \end{gathered}

我觉得第一步就很困难啊!注意到

A(AB)=(AB)B\begin{gathered} A(A-B)=(A-B)B \end{gathered}

X=ABX=A-B,如果XX可逆则A,BA,B相似,显然.现在考虑XX不可逆.

则存在v0,Xv=0v\ne 0,Xv=0.XBv=AXv=0XBv=AXv=0U=nullXU=\operatorname{null} XBB下不变.且这个空间中A=BA=B啊,所以XXAA下也不变.且其中detAU=detBU\det A|_U=\det B|_U.

VUV-U中,XX可逆,A,BA,B相似?这不对!这其中XX单而不一定满.正确的做法是归纳,AVU,BVUA|_{V-U},B|_{V-U}仍然满足题目中那个式子,于是就归纳下去做完了.

[think] 知乎答案说对给定某个等式的题这是套路.

Z2×2Z^{2\times 2}的可逆矩阵(行列式为11)可由下面两个生成元生成:

A=[1101]B=[0110]\begin{gathered} A=\begin{bmatrix} 1&1\\0&1 \end{bmatrix} \\ B=\begin{bmatrix} 0&-1\\1&0 \end{bmatrix} \end{gathered}

思路是对 M=[abcd]M=\begin{bmatrix} a&b\\c&d \end{bmatrix},你发现可以对a,ca,c进行辗转相减,让cc变成00.然后就a=d=1a=d=1,最后M=AbM=A^b了.