解方程

线性方程组:m∗nm*n 的系数矩阵 AA ,长度为 nn 的解向量 xx ,和长度为 mm 的偏置为 bb 。需要解 Ax=bAx=b 。

矩阵初等变换:

  1. 交换两行

  2. 某行倍加到另一行上

  3. 某行乘非 00 倍数

增广矩阵:把 AA 和 bb 拼到一起

阶梯型矩阵:从上到下,00 个数逐行严格增加,直到下方的全 00 行

最简阶梯型:在阶梯型的基础上,有两个条件:

  1. 每行第一个非 00 元是 11

  2. 每行第一个非 00 元所在的列,其他行全 00

(让我想到高斯-约旦消元法)

在最简阶梯型的基础上,可以调整列的位置关系。使得非全 00 行的个数为 rr ,且对于 i≤ri\le r ,第 ii 行的 00 只存在于长度为 i−1i-1 的前缀。

然后讨论解的状态,先看最后的全 00 行,如果存在一行的偏置非 00 ,则无解。

否则一定有解。设 rr 是非全 00 行的行数,如果 r=nr=n ,通过回代可以确定唯一解。如果 r<nr<n ,则解无限。

具体的,可以发现确定 xr+1,xr+2,…,xnx_{r+1},x_{r+2},\dots,x_n 之后,x1∼rx_{1\sim r} 都可以回代得到。我们称之为,有 n−rn-r 个自由变量。

齐次线性方程组: 不存在无解情况,其余讨论基本相同。

我们的讨论都是在数域 FF 上进行的,即可以进行加减乘除运算的代数结构。

向量空间

向量空间:包括向量的集合 AA 和域 BB ,它需要满足:

  1. AA 上的加法构成了阿贝尔群。

  2. AA 和 BB 之间需要定义数乘运算 ⋅:B×A→A\cdot:B×A\rightarrow A ,需要满足以下性质:

封闭性:∀b∈B,a∈A,ba∈A\forall b\in B,a\in A,ba\in A

结合律:∀b,b′∈B,a∈A\forall b,b'\in B,a\in A ,有 bb′a=b(b′a)bb'a=b(b'a) 。

分配律:∀b∈B,a,a′∈A\forall b\in B,a,a'\in A ,有 b(a+a′)=ba+ba′b(a+a')=ba+ba' 。∀b,b′∈B,a∈A\forall b,b'\in B,a\in A ,有 (b+b′)a=ba+b′a(b+b')a=ba+b'a 。

单位元作用:∀a∈A,1B⋅a=a\forall a\in A,1_B\cdot a=a 。

一些记号:设系数矩阵为 A∈Fm∗nA\in F^{m*n} ,解向量为 x∈Fnx\in F^{n} ,我们要解 Ax=βAx=\beta 。

设 AA 的每一列是 α1,α2,…,αn∈Fm\alpha_1,\alpha_2,\dots,\alpha_n\in F^m ,

可以把 AA 记作 (α1,α2,…,αn)(\alpha_1,\alpha_2,\dots,\alpha_n) 。

方程组可以看做 x1α1+x2α2+⋯+xnαn=βx_1\alpha_1+x_2\alpha_2+\dots+x_n\alpha_n=\beta 。

若 x1=k1,x2=k2,…,xn=knx_1=k_1,x_2=k_2,\dots,x_n=k_n 为解,则 β=k1α1+k2α2+⋯+knαn\beta=k_1\alpha_1+k_2\alpha_2+\dots+k_n\alpha_n

线性表述

β∈Fm\beta\in F^m 可由 α1,α2,…αs∈Fm\alpha_1,\alpha_2,\dots\alpha_s\in F^m 线性表述 ↔\leftrightarrow β=x1α1+x2α2+⋯+xsαs\beta=x_1\alpha_1+x_2\alpha_2+\dots+x_s\alpha_s 有解

线性子空间

对于线性空间 (A,B,+,⋅)(A,B,+,\cdot) ,若 A′⊆AA'\subseteq A 且满足以下条件:

∀α,β∈A′,αβ∈A′\forall \alpha,\beta\in A',\alpha\beta\in A' 。

∀α∈A′,k∈B,kα∈A′\forall \alpha \in A',k\in B,k\alpha\in A' 。

则 A′A' 是 AA 的子空间。

例子:齐次方程的解空间就是一个子空间。我们称一般方程的解空间为仿射空间。

线性无关

一组向量 α1,α2,…,αn\alpha_1,\alpha_2,\dots,\alpha_n 线性相关,当且仅当 α1x1+⋯+αnxn=0\alpha_1x_1+\dots +\alpha_nx_n=0 存在非 00 解。

不是线性相关,那就是线性无关。

注:α1,α2,…,αn\alpha_1,\alpha_2,\dots,\alpha_n 能表出 β\beta 只是 α1,α2,…,αn,β\alpha_1,\alpha_2,\dots,\alpha_n,\beta 线性相关的充分条件。

例子:α1=(1,0,0),α2=(0,0,0),β=(0,1,0)\alpha_1=(1,0,0),\alpha_2=(0,0,0),\beta =(0,1,0)

性质:

  1. 如果这些向量的一个子集线性相关,则所有向量都线性相关

  2. 若 s>ns>n ,则 α1,α2,…,αs∈Fn\alpha_1,\alpha_2,\dots,\alpha_s\in F^n 必线性相关

  3. 向量组相关,则其缩短组相关,向量组无关,则其延伸组无关

极大无关组

对于向量组 {α1,α2,…,αm}\{\alpha_1,\alpha_2,\dots,\alpha_m\} ,它的一个极大无关组 αi1,αi2,…,αir\alpha_{i_1},\alpha_{i_2},\dots,\alpha_{i_r} 满足:

  1. αi1,αi2,…,αir\alpha_{i_1},\alpha_{i_2},\dots,\alpha_{i_r} 无关

  2. ∀1≤i≤m\forall 1\le i\le m ,αi1,αi2,…,αir,αi\alpha_{i_1},\alpha_{i_2},\dots,\alpha_{i_r},\alpha_i 线性相关

条件 22 和 α1,α2,…,αm\alpha_1,\alpha_2,\dots,\alpha_{m} 均可被 αi1,αi2,…,αir\alpha_{i_1},\alpha_{i_2},\dots,\alpha_{i_r} 线性表出等价。

性质:含有非 00 的向量组,极大无关组必存在

向量组的等价关系

对于向量组 I1,I2I_1,I_2 ,称 I1I_1 能被 I2I_2 线性表出,当且仅当 ∀α∈I1\forall \alpha\in I_1 ,α\alpha 能被 I2I_2 线性表出。

I1,I2I_1,I_2 能互相表出,则称 I1,I2I_1,I_2 等价。显然这个关系满足自反性,传递性,对称性。

性质:一个向量组和它的任意极大无关组等价

性质:若 I1I_1 能表出 I2I_2 且 ∣I1∣<∣I2∣|I_1|<|I_2| ,则 I2I_2 线性相关(逆否:若 I1I_1 能表出 I2I_2 且 I2I_2 线性无关,则 ∣I1∣≥∣I2∣|I_1|\geq |I_2| )

性质:一个向量组的所有极大无关组包含的向量数目都是相等的。

向量组的秩

一个向量组的极大无关组包含的向量数目称为该向量组的秩。特例: r(0,0,…,0)=0r(0,0,\dots,0)=0 。

若两个向量组等价,则它们的秩相等。

性质:若 I1I_1 能表出 I2I_2 ,则 rank(I1)≤rank(I2)rank(I_1)\le rank(I_2) 。

矩阵的秩

设 A=(α1,…,αm)A=(\alpha_1,\dots,\alpha_m) ,则 α1,…,αm\alpha_1,\dots,\alpha_m 的秩即为 AA 的列秩,其中 A∈Fn×mA\in F^{n×m} ,同理可以定义行秩。

性质:行秩等于列秩。因为初等变换不会改变行秩列秩,消成最简阶梯型后又易得行秩=列秩。

所以我们将行秩列秩统称为矩阵的秩。

满秩:r(Am,n)=min⁡(m,n)r(A_{m,n})=\min(m,n) 。若 m=nm=n ,则称为满秩方阵。

我们可以通过化成最简阶梯型找到极大无关组。

向量组生成的子空间

对向量组 α1,…,αs∈Fn\alpha_1,\dots,\alpha_s\in F^n , {k1α1+⋯+ksαs∣k1,k2,…,ks∈F}⊆Fn\{k_1\alpha_1+\dots+k_s\alpha_s|k_1,k_2,\dots,k_s\in F\}\subseteq F^n 满足数乘,加法封闭,它是 FnF^n 的子空间。

我们将其称为该向量组生成的子空间,记作 L(α1,…,αs)L(\alpha_1,\dots,\alpha_s) 性质:它的极大无关组生成的子空间是相同的。

子空间的基

设 UU 是 KnK^n 的一个子空间,若向量组 I=(α1,α2,…,αr)I=(\alpha_1,\alpha_2,\dots,\alpha_r) 满足以下条件:

  1. II 线性无关

  2. UU 中每一个向量都能由 II 表出

则 II 是 UU 的一个基。

性质:任意两个基都是等价的,我们把一个向量空间的基的大小记为 dim⁡FU\dim_F U .

∀β∈U\forall \beta\in U ,可以把 β\beta 表示成 k1α1+⋯+krαrk_1\alpha_1+\dots+k_r\alpha_r ,且表示方式唯一。

称 (k1,…,kr)(k_1,\dots,k_r) 为 β\beta 在基 α1,…,αr\alpha_1,\dots,\alpha_r 下的坐标。

性质:U=L(α1,…,αr)U=L(\alpha_1,\dots,\alpha_r) 。显然 dim⁡L(α1,…,αr)=rank(α1,…,αr)\dim L(\alpha_1,\dots,\alpha_r)=rank(\alpha_1,\dots,\alpha_r)

性质:dd 维空间中大小 >d>d 的向量组必相关

性质:dd 维空间中任意 dd 个无关的向量都构成一组基

性质:若 U⊆VU\subseteq V 且 dim⁡U=dim⁡V\dim U=\dim V ,则 U=VU=V

方程解空间

现在我们来分析一下齐次方程 Ax=0Ax=0 的解空间 UU,令 A∈Fm×n,x∈FnA\in F^{m×n},x\in F^n 。定理: dimU=n−r(A)dim_U=n-r(A)

原因是你发现把 AA 消元之后自由变量的个数就是 n−r(A)n-r(A) ,它们的任意一个取法就对应了 UU 中的一个元素。

写出齐次方程的解空间:求出通解,找到 n−rn-r 个基(取一个自由变量为 11 ,其他自由变量为 00),从而确定基

非齐次方程怎么分析呢?

对于非齐次方程 α1x1+α2x2+⋯+αnxn=β\alpha_1x_1+\alpha_2x_2+\dots+\alpha_nx_n=\beta ,考虑 α\alpha 构成的矩阵 AA 和增广矩阵 A′A' 。

性质:有解 ↔\leftrightarrow r(A)=r(A′)r(A)=r(A') 。

定理:该方程的解集是 r0+Ur_0+U ,其中 r0r_0 为任一特解,UU 是 β=0\beta=0 时的解空间。

行列式

首先我们阐释一下映射中的像和原像:对于 f:A→Bf:A\to B ,S⊆AS\subseteq A ,称 SS 的像是 f(S)f(S) ,即 {f(x)∣x∈S}\{f(x)|x\in S\}

对于 T⊆BT\subseteq B ,称 TT 的原像是 f−1(T)f^{-1}(T) ,即 {x∣f(x)∈T}\{x|f(x)\in T\}

定义线性函数 f:Fn→Ff:F^n\to F 满足 f(kα+lβ)=kf(α)+lf(β)f(k\alpha+l\beta)=kf(\alpha)+lf(\beta)

然后我们称 f:Fn×Fn⋯×Fn→Ff:F^n×F^n\dots×F^n\to F (kk 个 FnF^n)为 FnF^n 上的 kk 重线性函数,若固定其他 k−1k-1 个向量,剩下的 Fn→FF^n\to F 是线性的。

k=1k=1 时即线性函数,k=2k=2 时称为双线性函数。

例:设 fi,j(α,β)=αiβjf_{i,j}(\alpha,\beta)=\alpha_i\beta_j ,那么 fi,jf_{i,j} 是双线性函数。

我们记 Lk(Fn)L^k(F^n) 为 FnF^n 上的 kk 重线性函数集,这个集合关于函数 +,⋅+,\cdot 构成了向量空间,且 fi1,i2,…,ikf_{i_1,i_2,\dots,i_k} 这 nkn^k 个函数是它的一组基,

因为能发现 f(α1,…,αk)=∑i1,i2,…,ik≤nf(ϵi1,ϵi2,…,ϵik)∏j=1kαj,ijf(\alpha_1,\dots,\alpha_k)=\sum\limits_{i_1,i_2,\dots,i_k\le n}f(\epsilon_{i_1},\epsilon_{i_2},\dots,\epsilon_{i_k})\prod\limits_{j=1}^k\alpha_{j,i_j} ,这个性质后面有用。

反对称: f∈Lk(Fn)f\in L^k(F^n) 满足若 αi=αj(i≠j)\alpha_i=\alpha_j(i\neq j) ,则 f(α1…αi…αj…αk)=0f(\alpha_1\dots\alpha_i\dots\alpha_j\dots\alpha_k)=0

其实这等价于 f(α1…αi…αj…αk)+f(α1…αj…αi…αk)=0f(\alpha_1\dots\alpha_i\dots\alpha_j\dots\alpha_k)+f(\alpha_1\dots\alpha_j\dots\alpha_i\dots\alpha_k)=0

怎么证明呢?先证条件 11 →\to 条件 22。考察 f(…αi+αj…αi+αj… )=f(…αj…αi… )+f(…αi…αj… )+f(…αi…αi… )+f(…αj…αj… )f(\dots \alpha_i+\alpha_j\dots\alpha_i+\alpha_j\dots)=f(\dots\alpha_j\dots\alpha_i\dots)+f(\dots\alpha_i\dots\alpha_j\dots)+f(\dots\alpha_i\dots\alpha_i\dots)+f(\dots\alpha_j\dots\alpha_j\dots)

把等于 00 的部分代入即可。条件 22 →\to 条件 11 是容易的。

定义 D:Mn(F)→FD:M_n(F)\to F 为 FnF^n 上的 nn 阶行列式函数,当且仅当: 它是 n-重线性函数,它反对称,且 D(En)=1D(E_n)=1

显然 DD 是唯一确定的,但我们尝试严谨的说明这个事情。

设 MijM_{ij} 是去掉 ai,ja_{i,j} 所在的行和列之后构成的 n−1n-1 阶方阵,称 Aij=(−1)i+jDn−1(Mij)A_{ij}=(-1)^{i+j}D_{n-1}(M_{ij}) 为 aija_{ij} 的代数余子式。

接下来归纳构造行列式函数,任取 1≤j≤n1\le j\le n ,令 Dn(A)=∑i=1nai,jAi,jD_n(A)=\sum\limits_{i=1}^n a_{i,j}A_{i,j} 可以证明它确实满足行列式函数的三个条件。

我们换一个方式写出来这个函数。考虑 f(α1,…,αn)=∑i1,i2,…,in≤nf(ϵi1,ϵi2,…,ϵin)∏j=1nαj,ijf(\alpha_1,\dots,\alpha_n)=\sum\limits_{i_1,i_2,\dots,i_n\le n}f(\epsilon_{i_1},\epsilon_{i_2},\dots,\epsilon_{i_n})\prod\limits_{j=1}^n\alpha_{j,i_j} ,而如果 ij=iki_j=i_k 则根据反对称 f(ϵi1,ϵi2,…,ϵin)=0f(\epsilon_{i_1},\epsilon_{i_2},\dots,\epsilon_{i_n})=0 。看来 {in}\{i_n\} 一定是一个排列。而且 f(ϵi1,ϵi2,…,ϵin)f(\epsilon_{i_1},\epsilon_{i_2},\dots,\epsilon_{i_n}) 也容易计算出来,因为它可以从 {1,2,…,n}\{1,2,\dots,n\} 开始不断交换相邻项得到,根据反对称性,我们直接看它的逆序数的奇偶性即可。

所以 f(α1,…,αn)=∑p∈Symn(−1)τ(p)∏i=1nαi,pif(\alpha_1,\dots,\alpha_n)=\sum\limits_{p\in Sym_n}(-1)^{\tau(p)}\prod\limits_{i=1}^n \alpha_{i,p_i} ,其中 τ(p)\tau(p) 表示逆序对数。

求解行列式

最朴素的方法是用高斯消元消成上三角矩阵。

然后对于一些特殊的行列式,我们可以进行一些手玩以进行消元。

另一个方法是用代数余子式。我们称 ∣A∣|A| 按第 ii 行展开为 ∣A∣=∑j=1nai,jAi,j|A|=\sum\limits_{j=1}^n a_{i,j}A_{i,j} (对于一些特殊的稀疏矩阵很有效)

另一个 trick 是把原矩阵在最后添加一行一列,使得 ∀i≤n,ai,n+1=0\forall i\le n,a_{i,n+1}=0 ,而 an+1,n+1=1a_{n+1,n+1}=1 。

还有一个事情:det⁡(α1,…,αk+β,…,αn)=det⁡(α1,…,αk,…,αn)+det⁡(α1,…,β,…,αn)\det(\alpha_1,\dots,\alpha_k+\beta,\dots,\alpha_n)=\det(\alpha_1,\dots,\alpha_k,\dots,\alpha_n)+\det(\alpha_1,\dots,\beta,\dots,\alpha_n)

Laplace 展开

我们先定义 kk 阶子式的代数余子式:

设这个子式下标是 i1,i2,…,ik,j1,j2,…,jki_1,i_2,\dots,i_k,j_1,j_2,\dots,j_k ,那么它的代数余子式就是 (−1)i1+i2+⋯+ik+j1+j2+⋯+jk(-1)^{i_1+i_2+\dots+i_k+j_1+j_2+\dots+j_k} 乘以余下的矩阵的行列式。

Laplace 展开(多行/多列展开):

考虑 A=(ai,j)nA=(a_{i,j})_n 中的 kk 个行 i1<i2<⋯<iki_1<i_2<\dots<i_k ,这 kk 个行对应了 t=(nk)t=\tbinom{n}{k} 个 kk 阶子式,我们记为 D1,D2,…,DtD_1,D_2,\dots,D_t ,对应的代数余子式为 A1,A2,…,AtA_1,A_2,\dots,A_t ,则 ∣A∣=∑i=1tDiAi|A|=\sum\limits_{i=1}^t D_iA_i

例:我们可以证明 ∣A0BC∣=∣A∣⋅∣C∣\left|\begin{array}{}A&0\\B&C\end{array}\right|=|A|\cdot |C| ,其中 A,CA,C 都是方阵。理由就是设 AA 的行数为 nn ,对前 nn 行展开即可。

范德蒙德行列式

设 D(a1,a2,…,an)D(a_1,a_2,\dots,a_n) 是满足 Ai,j=aij−1A_{i,j}=a_i^{j-1} 的矩阵的行列式。

考虑以下消元:从第 n−1n-1 列到第 11 列 ,依次把第 ii 行乘上 −a1-a_1 加到第 i+1i+1 行上。那除了 a1,1a_{1,1} 第 11 列变成全 00 ,而 ∀i≥2,j≥2\forall i\geq 2,j\geq 2 ,ai,j′=aij−a1aij−1=ai,hj−1(ai−a1)a'_{i,j}=a_i^{j}-a_1a_i^{j-1}=a_i^{ ,hj-1}(a_i-a_1) 。

于是按第 11 列展开,并对每一行提取公因子 ai−a1a_i-a_1 ,可得 D(a1,a2,…,an)∏i=2n(ai−a1)D(a2,…,an)D(a_1,a_2,\dots,a_n)\prod\limits_{i=2}^n (a_i-a_1)D(a_2,\dots,a_n)

于是有 D(a1,a2,…,an)=∏1≤i<j≤n(aj−ai)D(a_1,a_2,\dots,a_n)=\prod\limits_{1\le i<j\le n}(a_j-a_i) 。

克莱姆法则

考虑方程 x1α1+x2α2+⋯+xnαn=βx_1\alpha_1+x_2\alpha_2+\dots+x_n\alpha_n=\beta ,矩阵 A=(α1,α2,…,αn)∈Fn×nA=(\alpha_1,\alpha_2,\dots,\alpha_n)\in F^{n×n} ,则方程组有唯一解 ↔det⁡(A)≠0↔rank(A)=n\leftrightarrow \det(A)\neq 0\leftrightarrow rank(A)=n ,且这个解为 xi=∣Ai∣∣A∣x_i=\frac{|A_i|}{|A|} ,其中 Ai=(α1,…,αi−1,β,αi+1,…,αn)A_i=(\alpha_1,\dots,\alpha_{i-1},\beta,\alpha_{i+1},\dots,\alpha_n) 。

这个东西是能直接解出来的,考虑 det⁡(Ai)=det⁡(α1,…,αi−1,∑xjαj,αi+1,…,αn)=∑jxjdet⁡(α1,…,αi−1,αj,αi+1,…,αn)=xidet⁡(A)\det(A_i)=\det(\alpha_1,\dots,\alpha_{i-1},\sum x_j\alpha_j,\alpha_{i+1},\dots,\alpha_n)=\sum\limits_{j}x_j\det(\alpha_1,\dots,\alpha_{i-1},\alpha_{j},\alpha_{i+1},\dots,\alpha_n)=x_i\det(A) 。于是 xi=∣Ai∣∣A∣x_i=\frac{|A_i|}{|A|} 。

这里可以再讨论一下行列式和秩的联系。

引理:若 AA 有一个 kk 阶子式非 00 ,则 rank(A)≥krank(A)\geq k 。

定理:r(A)=k↔r(A)=k\leftrightarrow 存在一个非 00 的 kk 阶子式且所有 k+1k+1 阶子式皆为 00 。

一个记号:A(i1,…,ikj1,…,jk)A\begin{pmatrix}i_1,\dots,i_k\\j_1,\dots,j_k\end{pmatrix} 表示 AA 的一个 kk 阶子式

矩阵运算

加法,数乘,零元,负元,乘法皆显然

还可以对矩阵做分块,比如写成行/列向量的拼接,这样往往方便我们分析乘法的性质。比如说有 (α1⋮αn)B=(α1B⋮αnB)\begin{pmatrix}\alpha_1\\\vdots \\\alpha_n\end{pmatrix}B=\begin{pmatrix}\alpha_1B\\\vdots \\\alpha_nB\end{pmatrix} 。

性质:r(A+B)≤r(A)+r(B),r(AB)≤min⁡(r(A),r(B))r(A+B)\le r(A)+r(B),r(AB)\le \min(r(A),r(B))。若 An∗mBm∗p=0A_{n*m}B_{m*p}=0 ,则 r(A)+r(B)≤mr(A)+r(B)\le m。r(ATA)=r(AAT)=r(A)r(A^TA)=r(AA^T)=r(A) 。

我们来证为何 r(AAT)=r(A)r(AA^T)=r(A) ,只需要证 AATX=0AA^TX=0 和 AX=0AX=0 同解。

一个方向显然。另一个方向:AATX=0→XTAATX=0→∣∣ATX∣∣2=0→ATX=0AA^TX=0\to X^TAA^TX=0\to ||A^TX||^2=0\to A^TX=0

同理还有 r(AATA)=r(A)r(AA^TA)=r(A)

向量 α1,α2,…,αn\alpha_1,\alpha_2,\dots,\alpha_n 能表出 β1,β2,…,βm\beta_1,\beta_2,\dots,\beta_m 当且仅当存在矩阵 Xn∗mX_{n*m} 使得 (α1,α2,…,αn)X=(β1,β2,…,βm)(\alpha_1,\alpha_2,\dots,\alpha_n)X=(\beta_1,\beta_2,\dots,\beta_m) 。因为 βj=∑iαiXi,j\beta_j=\sum\limits_{i} \alpha_i X_{i,j} 。

观察:对矩阵的行初等变换全部可以看做左乘了一个方阵,同理列变换是右乘一个方阵。

第一类变换:交换两行 第二类变换:一行乘 kk 倍 第三类变换:一行乘 kk 倍加到另一行

只用第三类变换就能把矩阵消成:一个子矩阵是 r∗rr*r 的对角矩阵,其余位置全 00 。而且第三类变换不改变矩阵行列式。

这样就能证明 ∣A∣∣B∣=∣AB∣|A||B|=|AB| :

设 A=x1x2…xpDy1y2…yqA=x_1x_2\dots x_pDy_1y_2\dots y_q 。则 ∣AB∣=∣x1x2…xpDy1y2…yqB∣=∣Dy1y2…yqB∣|AB|=|x_1x_2\dots x_pDy_1y_2\dots y_qB|=|Dy_1y_2\dots y_qB| 。注意到对角矩阵乘上 y1y2…yqBy_1y_2\dots y_qB 其实是在做第二类变换,于是它就等于 ∣D∣⋅∣y1y2…yqB∣=∣D∣⋅∣B∣=∣A∣∣B∣|D|\cdot |y_1y_2\dots y_qB|=|D|\cdot |B|=|A||B| 。

现在考虑矩阵的逆 A−1A^{-1} 。我们有 AA 可逆当且仅当 ∣A∣≠0|A|\neq 0 。且逆元唯一:如果 B,CB,C 都是 AA 的逆,那么 B=BE=BAC=(BA)C=CB=BE=BAC=(BA)C=C 。

A−1A^{-1} 怎么计算:考虑伴随矩阵 A∗A^{*} 满足 Ai,j∗=Aj,iA^{*}_{i,j}=A_{j,i} ,其中 Ai,jA_{i,j} 表示代数余子式。则根据行列式的展开/异行展开,可得 A∗A=AA∗=∣A∣IA^{*}A=AA^{*}=|A|I ,于是 A−1=A∗∣A∣A^{-1}=\frac{A^{*}}{|A|} 。

更高效的方法:对 (AI)\begin{pmatrix}A&I\end{pmatrix} 做行初等变换得到 (IB)\begin{pmatrix}I&B\end{pmatrix} ,有 B=A−1B=A^{-1} ,因为初等变换可以看做乘初等矩阵。

性质:(AB)−1=B−1A−1(AB)^{-1}=B^{-1}A^{-1}

降阶定理:∣A+BC∣=∣A∣⋅∣I+CA−1B∣|A+BC|=|A|\cdot |I+CA^{-1}B| 这可以通过分块矩阵推导得到

有些地方也叫做 矩阵行列式引理 ,根据调查,在 OI 中竟然是有用的

正交

对于某线性子空间 V⊆RnV\subseteq R^n ,我们如此定义正交基:

V=span(β1,…,βs)V=span(\beta_1,\dots,\beta_s) ,且 ∀i≠j\forall i\neq j ,βi⋅βj=0\beta_i\cdot \beta_j=0 。

进一步的,如果还满足 ∀i,βi⋅βi=1\forall i,\beta_i\cdot \beta_i=1 ,我们称为标准正交基。

怎么求出一组正交基呢?设 V=span(α1,…,αs)V=span(\alpha_1,\dots,\alpha_s) ,我们可以如此构造:

β1=α1,β2=α2−β1⋅α2β1⋅β1β1\beta_1=\alpha_1,\beta_2=\alpha_2-\frac{\beta_1\cdot \alpha_2}{\beta_1\cdot \beta_1}\beta_1 ,以此类推,βs=αs−∑i=1s−1αs⋅βiβi⋅βiβi\beta_s=\alpha_s-\sum\limits_{i=1}^{s-1}\frac{\alpha_s\cdot \beta_i}{\beta_i\cdot \beta_i}\beta_i 。

由正交基求标准正交基是 trivial 的。

现在我们定义 正交阵 :

Q∈Rn×nQ\in R^{n×n} 是正交阵,当且仅当 QTQ=IQ^TQ=I 。

你发现如果设 Q=(η1,…,ηn)Q=(\eta_1,\dots,\eta_n) ,那 η1,…,ηn\eta_1,\dots,\eta_n 其实是 RnR^n 的一组标准正交基。

线性映射

接下来,对于 V1,V2V_1,V_2 是 RnR^n 的线性子空间,我们定义线性映射 f:V1→V2f:V_1\to V_2 满足:f(α1+α2)=f(α1)+f(α2)f(\alpha_1+\alpha_2)=f(\alpha_1)+f(\alpha_2) 以及 f(kα)=kf(α)f(k\alpha)=kf(\alpha) ,也就是说它保加法,保数乘。

线性映射有一些良好性质。若 ff 是单射,则 V1V_1 的无关组映射后还是无关组。若 ff 是同构,则 ff 将 V1V_1 的基映射到 V2V_2 的基。

若存在 V1V_1 到 V2V_2 的同构映射,则称 V1,V2V_1,V_2 同构,记作 V1≅V2V_1\cong V_2 。

可以证明 V1≅V2⇔dim⁡V1=dim⁡V2V_1\cong V_2\Leftrightarrow \dim V_1=\dim V_2

现在我们定义线性映射的核和像。定义 ker⁡f={α∈V1∣f(α)=0}⊆V1\ker f=\{\alpha\in V_1|f(\alpha)=0\}\subseteq V_1 定义 Imf={f(α)∣α∈V1}⊆V2\text{Im} f=\{f(\alpha)|\alpha\in V_1\}\subseteq V_2

ker 和 Im 有良好的性质。考虑 m∗nm*n 的矩阵以及 f:Fn→Fmf:F^n\to F^m 满足 f(X)=AXf(X)=AX ,我们有 dim⁡ker⁡f=n−r(A),dim⁡Imf=r(A)\dim \ker f=n-r(A),\dim \text{Im} f=r(A) 。

这让我想起期中考试的最后一题:证明dim⁡{BX∣ABX=0}=r(B)−r(AB)\dim \{BX|ABX=0\}=r(B)-r(AB) 。其实很简单,设 V1={BX}V_1=\{BX\} ,设 f:V1→Fnf:V_1\to F^n 满足 f(P)=APf(P)=AP ,那么 dim⁡{BX∣ABX=0}=ker⁡f=dim⁡V1−Imf=r(B)−r(AB)\dim \{BX|ABX=0\}=\ker f=\dim V_1-\text{Im} f=r(B)-r(AB) 。

我们可以观察线性映射的整体结构。对于线性空间 V1,V2V_1,V_2 ,考察 V1,V2V_1,V_2 间的所有线性映射,即 HomF(V1,V2)\text{Hom}_F(V_1,V_2) 。发现每个线性映射都可以用一个矩阵表示,因为我们能取出 V1V_1 的基,确定 f(α1),…,f(αn)f(\alpha_1),\dots,f(\alpha_n) 即可。所以我们称其矩阵表示为 (f(α1),…,f(αn))=(β1,…,βm)Am×n(f(\alpha_1),\dots,f(\alpha_n))=(\beta_1,\dots,\beta_m)A_{m×n} 。这让我们感受到研究线性映射其实就是研究矩阵。

我们可以在 HomF(V,V)Hom_F(V,V) 定义加法和乘法,和矩阵的加法乘法类似。有 (f+g)(x)=f(x)+g(x),(fg)(x)=f(g(x))(f+g)(x)=f(x)+g(x),(fg)(x)=f(g(x)) 。

接下来我们可以观察同一线性映射的不同表示间的关系:

考虑 HomF(V1,V2)Hom_F(V_1,V_2) 的某映射 ff ,设 V1V_1 的两组基是 (α1,…,αm),(α1′,…,αm′)(\alpha_1,\dots,\alpha_m),(\alpha'_1,\dots,\alpha'_m) ,V2V_2 的两组基是 (β1,…,βn),(β1′,…,βn′)(\beta_1,\dots,\beta_n),(\beta'_1,\dots,\beta'_n) 。

我们设 (f(α1),…,f(αm))=(β1,…,βn)A,(f(α1′),…,f(αm′))=(β1′,…,βn′)B,(α1′,…,αm′)=(α1,…,αm)S,(β1′,…,βn′)=(β1,…,βn)T(f(\alpha_1),\dots,f(\alpha_m))=(\beta_1,\dots,\beta_n)A,(f(\alpha'_1),\dots,f(\alpha'_m))=(\beta'_1,\dots,\beta'_n)B,(\alpha'_1,\dots,\alpha'_m)=(\alpha_1,\dots,\alpha_m)S,(\beta'_1,\dots,\beta'_n)=(\beta_1,\dots,\beta_n)T 。我们有 TB=ASTB=AS ,于是 T−1AS=BT^{-1}AS=B 。

反之,对任意可逆矩阵 T,ST,S ,若 AA 是 ff 的某个矩阵表示,那 B=TASB=TAS 也是 ff 的某个矩阵表示。

同理,考虑 HomF(V,V)Hom_F(V,V) 中的某映射 ff ,设 VV 的两组基是 (α1,…,αm)(\alpha_1,\dots,\alpha_m) 和 (β1,…,βm)(\beta_1,\dots,\beta_m)。令 (f(α1),…,f(αm))=(α1,…,αm)A,(f(β1),…,f(βm))=(β1,…,βm)B,(β1,…,βm)=(α1,…,αm)U(f(\alpha_1),\dots,f(\alpha_m))=(\alpha_1,\dots,\alpha_m)A,(f(\beta_1),\dots,f(\beta_m))=(\beta_1,\dots,\beta_m)B,(\beta_1,\dots,\beta_m)=(\alpha_1,\dots,\alpha_m)U ,我们有 U−1AU=BU^{-1}AU=B 。

反之,对线性映射 ff 和某矩阵表示 AA 和可逆矩阵 TT ,都有 B=T−1ATB=T^{-1}AT 亦为 ff 在某矩阵表示。

上面的事情启发了我们什么呢?

相抵和相似

相抵

在 Mm×n(F)M_{m×n}(F) 上,我们定义等价关系 ∼\sim 满足 $A\sim B\Leftrightarrow $ AA 能通过初等变换得到 BB ,也就是存在可逆矩阵 S,TS,T 使得 SAT=BSAT=B 。你发现这相当于 A∼B⇔r(A)=r(B)A\sim B\Leftrightarrow r(A)=r(B) 。这被称为 A,BA,B 相抵 。

通过之前的分析我们知道:同一线性映射的不同表示是某个等价类。所以对任一线性映射 f:V1→V2f:V_1\to V_2 ,我们都能找到 V1V_1 的一组基 (α1,…,αm)(\alpha_1,\dots,\alpha_m) ,V2V_2 的一组基 β1,…,βn\beta_1,\dots,\beta_n ,使得 (f(α1),…,f(αm))=(β1,…,βn)(Er000)(f(\alpha_1),\dots,f(\alpha_m))=(\beta_1,\dots,\beta_n)\begin{pmatrix}E_r&0\\0&0\end{pmatrix} 。

趣题:证明任一矩阵 Am×nA_{m×n} 都能拆成一个列满秩矩阵和行满秩矩阵的乘积。

设其秩为 rr ,有 A=S(Er000)T=S(Er0)(Er0)TA=S\begin{pmatrix}E_r&0\\0&0\end{pmatrix}T=S\begin{pmatrix}E_r\\0\end{pmatrix}\begin{pmatrix}E_r&0\end{pmatrix}T ,其中 S,TS,T 是可逆矩阵。那么 S(Er0)S\begin{pmatrix}E_r\\0\end{pmatrix} 列满秩,(Er0)T\begin{pmatrix}E_r&0\end{pmatrix}T 行满秩。

相似

我们再在 Mn∗n(F)M_{n*n}(F) 上定义等价关系:A∼B⇔∃A\sim B\Leftrightarrow \exists 可逆阵 UU 使得 U−1AU=BU^{-1}AU=B 。我们称这样的 A,BA,B 是相似的。

相似有更多良好的性质。

首先它也是等价关系。其次,如果 A∼BA\sim B ,那么 det(A)=det(B)det(A)=det(B) ,tr(A)=tr(B)tr(A)=tr(B) 。注:对于任意矩阵 A,BA,B ,都有 tr(AB)=tr(BA)tr(AB)=tr(BA) 。

接下来,类似于相抵关系下的标准形,我们希望找到相似关系的标准形,那其实就是希望找到 UU 使得 D=U−1AUD=U^{-1}AU 是对角阵。如果能找到这样的 UU ,就称 AA 可对角化,DD 是 AA 的相似标准形。

经过一些计算,发现其实就是要找到线性无关的 nn 个向量 α1,…,αn\alpha_1,\dots,\alpha_n ,使得 Aαi=λiαiA\alpha_i=\lambda_i\alpha_i 。

这启发我们研究 Aα=λ0αA\alpha=\lambda_0\alpha 的性质。我们称 λ0\lambda_0 是 AA 的特征值,α\alpha 是 AA 的属于特征值 λ0\lambda_0 的一个特征向量。

不难发现 λ0\lambda_0 是特征值 ⇔(A−λ0I)α=0\Leftrightarrow (A-\lambda_0I)\alpha=0 有非零解 ⇔\Leftrightarrow det⁡(A−λ0I)=0\det(A-\lambda_0I)=0 。

这引出了 特征多项式 :设 f(λ)=∣λI−A∣f(\lambda)=|\lambda I-A| ,那么特征值就是这个函数的所有零点。我们把 ker⁡(λiI−A)\ker (\lambda_i I-A) 称为 AA 的属于 λi\lambda_i 的特征子空间 ,记为 ViV_i 。

然后可以证明:若 λi≠λj\lambda_i\neq \lambda_j ,则 Vi∩Vj={0}V_i\cap V_j=\{0\} 。

于是 AnA_n 可对角化 ⇔\Leftrightarrow ∑i=1rdim⁡Vi=n\sum\limits_{i=1}^r \dim V_i=n

我们设 f(λ)=∣λE−A∣=∏i=1r(λ−λi)αif(\lambda)=|\lambda E-A|=\prod\limits_{i=1}^r (\lambda-\lambda_i)^{\alpha_i} 。

称 αi\alpha_i 是 λi\lambda_i 的代数重数,dim⁡Vi\dim V_i 是几何重数。可以证明几何重数 ≤\le 代数重数。

特征值的函数保持性 设 λ0\lambda_0 为 AA 的特征值,α\alpha 是对应特征向量,有 g(λ0)g(\lambda_0) 是 g(A)g(A) 的特征向量,且 α\alpha 是 g(A)g(A) 的属于 g(λ0)g(\lambda_0) 的特征向量。证明很简单:可以归纳证明 Akα=λ0kαA^k\alpha=\lambda_0^k\alpha 。

可以导出:如果非零矩阵 AA 满足 Ak=0A^k=0 ,那么 AA 不可对角化。因为 Ak=0→λk=0→λ=0A^k=0\to \lambda^k=0\to \lambda=0 ,特征值只有 00 且可对角化,就只能导出 00 了。

例题:若 A=A2A=A^2 ,则 AA 必能对角化。因为 AA 特征值满足 λ=λ2\lambda=\lambda^2 ,于是特征值只有 0,10,1 。往证:dim⁡ker⁡A+dim⁡ker⁡(A−I)=n\dim \ker A+\dim \ker (A-I)=n 。这是经典题

实对称矩阵对角化

命题 1:实对称矩阵的特征多项式在复数域下每一个根都是实数。

命题 2:实对称矩阵属于不同特征值的特征向量是正交的。

命题 3:实对称矩阵一定正交相似于对角矩阵。

二次型和合同

我们称 nn 元多项式 f(x1,…,xn)=∑i,jai,jxixjf(x_1,\dots,x_n)=\sum\limits_{i,j}a_{i,j}x_ix_j 为二次型 ,称 A=(ai,j)n×nA=(a_{i,j})_{n×n} 为该二次型的矩阵,其中 ai,j=aj,ia_{i,j}=a_{j,i} 。设 X=(x1,…,xn)TX=(x_1,\dots,x_n)^T ,其实就有 f(x1,…,xn)=XTAXf(x_1,\dots,x_n)=X^TAX 。

Extra:抽象版本的二次型

VV 上的二次型指一个函数 g:V→Fg:V\to F 满足

  1. g(kα)=k2g(α)∀k∈F,α∈Vg(k\alpha)=k^2g(\alpha)\forall k\in F,\alpha\in V
  2. φ(α,β)=12(g(α+β)−g(α)−g(β))\varphi(\alpha,\beta)=\frac{1}{2}(g(\alpha+\beta)-g(\alpha)-g(\beta)) 为 VV 上的双线性函数。

可以证明在这个定义下, g(∑i=1nxiαi)=∑i=1ng(αi)xi2+2∑i<jxixjφ(αi,αj)g(\sum\limits_{i=1}^n x_i\alpha_i)=\sum\limits_{i=1}^n g(\alpha_i)x_i^2+2\sum\limits_{i<j}x_ix_j\varphi(\alpha_i,\alpha_j)。

于是对称双线性函数 φ\varphi 和二次型相互决定,我们有 φ(α,β)=12(g(α+β)−g(α)−g(β))\varphi(\alpha,\beta)=\frac{1}{2}(g(\alpha+\beta)-g(\alpha)-g(\beta)) 和 g(α)=φ(α,α)g(\alpha)=\varphi(\alpha,\alpha) 。

然后我们可以研究二次型之间的等价关系,对于两个二次型 XTAXX^TAX 和 YTBYY^TBY ,如果存在可逆矩阵 CC 使得:若 X=CYX=CY (我们称这个过程是 非退化线性替换)则 XTAX=YTBYX^TAX=Y^TBY ,则称这两个二次型等价。你发现这个条件其实就是说存在可逆矩阵 CC 使得 CTAC=BC^TAC=B 。我们称这样的矩阵 AA 与 BB 合同 ,记作 A≃BA\simeq B 。

容易证明合同也是等价关系,我们记合同关系下的某个等价类为 合同类

现在我们想研究某个二次型 XTAXX^TAX 是否能等价于只含平方项的二次型,我们称这个是 XTAXX^TAX 的 标准形。

如果对称矩阵 AA 合同于一个对角矩阵,就称其为 AA 的 合同标准形

惯性定理 一个实对称矩阵一定合同于唯一的规范形。规范形如:即对角矩阵且只有 1,−1,01,-1,0 ,且 11 都写在前面,其次是 −1-1 。我们称这里的 1 的个数为 正惯性指数。

两个对称矩阵合同 ⇔\Leftrightarrow 它们秩相等,正惯性指数相等

正定

接下来考虑研究 正定矩阵 :对任意非零向量 α\alpha ,都有 αTAα>0\alpha^TA\alpha>0 。

同理可以定义半正定,即任意向量都满足 αTAα≥0\alpha^TA\alpha\geq 0

定理:以下四个条件等价:

\textcircled{1} AA 正定 \textcircled{2} A≃EA\simeq E \textcircled{3} ∀X∈RnXTAX≥0\forall X\in R^n\quad X^TAX\geq 0 且取等当且仅当 X=0X=0 \textcircled{4} AA 的特征值均为正实数

同时半正定等价于 A≃(Ir000)A\simeq \begin{pmatrix}I_r&0\\0&0\end{pmatrix}

正定阵的性质:

AA 正定则存在可逆 CC 使得 A=CTCA=C^TC

AA 正定则 A−1,ATA^{-1},A^T 均正定

AA 正定则 AkA^{k} 正定

AA 正定 BB 半正定则 A+BA+B 正定

AnA_n 正定 且 Pn×mP_{n×m} 列满秩则 PTAPP^TAP 正定

一个有趣的例题:考虑 A=12(E+αTα)A=\frac{1}{2}(E+\alpha^T\alpha) 是否正定,其中 α=(1,1,…,1)T\alpha=(1,1,\dots,1)^T 。

发现设 B=αTαB=\alpha^T\alpha ,有 B2=nBB^2=nB ,所以 BB 的特征值只会有 00 或 nn ,那么 AA 的特征值有 n+12,12\frac{n+1}{2},\frac{1}{2} ,自然是正定的

正定之主子式判别:

AA 正定 ⇔\Leftrightarrow AA 的各级主子式均正实数 ⇔\Leftrightarrow AA 的各级顺序主子式均正实数

证明过程蕴含了 AA 正定 ⇒\Rightarrow ∣A∣≤∏i=1nai,i|A|\le \prod\limits_{i=1}^n a_{i,i} 这档事

定理: AA 正定/半正定 ⇒\Rightarrow 存在唯一正定/半正定 SS 使得 A=S2A=S^2

酉空间

考虑共轭双线性函数 φ:V×V→C\varphi: V×V \to C 满足 φ(∑xiαi,∑yjαj)=∑i∑jxi‾yjφ(αi,αj)\varphi(\sum\limits x_i\alpha_i,\sum\limits y_j\alpha_j)=\sum\limits_{i}\sum\limits_{j}\overline{x_i}y_j\varphi(\alpha_i,\alpha_j) ,可以记作 φ(α,β)=X‾TAY\varphi(\alpha,\beta)=\overline{X}^TAY 。

复合同

$X\simeq Y\Leftrightarrow $ 存在可逆矩阵 CC 使得 C‾TAC=B\overline{C}^TAC=B(类比合同)

Hermite

定义 Hermite 双线性函数:是共轭双线性函数且 φ(α,β)=φ(β‾,α)\varphi(\alpha,\beta)=\varphi(\overline{\beta},\alpha) 。其矩阵表示满足 A=A‾TA=\overline{A}^T 。反 Hermite 满足 A=−A‾TA=-\overline{A}^T 。(类比实对称矩阵)

正定 酉空间 内积

称 VV 上的 hermite 双线性函数 hh 是正定的当且仅当 ∀\forall 非零向量 α\alpha 都有 φ(α,α)>0\varphi(\alpha,\alpha)> 0 。称 (V,h)(V,h) 为酉空间。hh 称为酉空间 VV 的内积。

标注内积:(X,Y)=X‾TY(X,Y)=\overline{X}^TY

正交化

记忆:βn=αn−∑i=1n−1(βi,αn)(βi,βi)βi\beta_n=\alpha_n-\sum\limits_{i=1}^{n-1}\frac{(\beta_i,\alpha_n)}{(\beta_i,\beta_i)}\beta_i

酉矩阵

U‾TU=E\overline{U}^TU=E (类比正交矩阵)

定理:nn 阶 Hermite 矩阵的 nn 个特征值都是实数

定理:Hermite 矩阵酉相似于对角形,也即存在酉矩阵 UU 使得 U−1AUU^{-1}AU 是对角矩阵。(类比:实对称矩阵正交相似于对角矩阵)

定理:Hermite 矩阵复合同于规范矩阵 (Ep−Eq0)\begin{pmatrix}E_p&&\\&-E_q&\\&&0\end{pmatrix} (类比:规范形,惯性原理)

定理:设 HH 是 nn 阶 Hermite 矩阵,则

HH 正定 ⇔\Leftrightarrow HH 复合同于 EnE_n ⇔\Leftrightarrow ∀X∈Cn,X‾THX≥0\forall X\in \mathbb{C}^n,\overline{X}^THX\geq 0 ,等号成立当且仅当 X=0X=0 ⇔\Leftrightarrow HH 的特征值皆为正实数

定理:设 HH 是 nn 阶 Hermite 矩阵,则

HH 正定 ⇔\Leftrightarrow HH 各阶主子式皆为正实数 ⇔\Leftrightarrow HH 各阶顺序主子式皆为实数