矩阵论 - Part I

矩阵论 - Part I

文章目录

  • 矩阵论 - Part I
    • 概念索引
    • 1 线性基础
    • 2 矩阵代数
    • 3 线性方程组

概念索引

1 行列式, Laplace定理(展开), Gramer法则, 齐次线性方程组
2 同型矩阵, 矩阵转置, 矩阵共轭, 反对称矩阵, 伴随矩阵, 逆矩阵, 分块矩阵, 子式, 余子式, 秩
3 线性变换, 等价矩阵, 初等矩阵, 矩阵迭代法

1 线性基础

矩阵3大特点:

  1. 矩阵是线性的
  2. 矩阵是离散的
  3. 矩阵是代数和几何交融的

行列式

n n n行列式 D D D的值为
D = ∣ a 11 a 12 ⋯ a 1 n a 21 a 22 ⋯ a 2 n ⋮ ⋮ ⋮ a n 1 a n 2 ⋯ a n n ∣ = ∑ t = 0 n ! ( − 1 ) t a 1 p 1 a 2 p 2 ⋯ a n p n D = \left| \begin{matrix} a_{11} & a_{12} & \cdots & a_{1n} \\ a_{21} & a_{22} & \cdots & a_{2n} \\ \vdots & \vdots & &\vdots \\ a_{n1} & a_{n2} & \cdots &a_{nn} \end{matrix} \right | = \sum_{t=0}^{n!} (-1)^ta_{1p_1}a_{2p_2}\cdots a_{np_n} D=a11a21an1a12a22an2a1na2nann=t=0n!(1)ta1p1a2p2anpn
其中 ( p 1 , p 2 , ⋯   , p n ) (p_1, p_2, \cdots, p_n) (p1,p2,,pn) ( 1 , 2 , ⋯   , n ) (1, 2, \cdots, n) (1,2,,n)的一个排列, t t t为该排列的逆序数.

行列式的性质

性质1: D T = D = d e t ( a i j ) D^T = D=det(a_{ij}) DT=D=det(aij)

性质2: 更换行列式的行(或列), 行列式变号
推论: 若行列式两行(或列) 完全相同, 则行列式为 0 0 0

性质3: 行列式某行(或列)乘以一个系数 k k k, 则行列式的值变为原来的 k k k
推论: 行列式的行(或列)有相同公因子, 则可提到行列式前

性质4: 行列式的某行(或列)的元素为两数之和, 则可认为是两个行列式之和

性质5: 行列式中某两列(或行)成比例, 则行列式的值为 0 0 0

性质6: 把行列式中的一行(或一列)乘以一个系数加到另一行(或一列), 行列式的值不变

Laplace定理: 行列式的另一展开定理

Laplace定理是更加广义的行列式展开定理, 需要用到两个概念:

  1. 余子式
  2. 代数余子式

Gramer法则

对于线性方程组 D x = b Dx=b Dx=b, 若 D ≠ 0 D \neq 0 D=0, 则方程组有唯一解:
x 1 = D 1 D , x 2 = D 2 D , ⋯   , x n = D n D x_1 = \frac{D_1}{D}, x_2 = \frac{D_2}{D}, \cdots, x_n = \frac{D_n}{D} x1=DD1,x2=DD2,,xn=DDn
其中
D j = ∣ a 11 ⋯ a 1 , j − 1 b 1 a 1 , j + 1 ⋯ a 1 n a 21 ⋯ a 2 , j − 1 b 1 a 2 , j + 1 ⋯ a 2 n ⋮ ⋮ ⋮ ⋮ ⋮ a n 1 ⋯ a n , j − 1 b 1 a n , j + 1 ⋯ a n n ∣ j = 1 , 2 , ⋯   , n D_{j} = \left| \begin{matrix} a_{11} & \cdots & a_{1, j-1} &b_1& a_{1, j+1} & \cdots & a_{1n} \\ a_{21} & \cdots & a_{2, j-1} &b_1& a_{2, j+1} & \cdots & a_{2n} \\ \vdots & &\vdots & \vdots &\vdots & &\vdots \\ a_{n1} & \cdots & a_{n, j-1} &b_1& a_{n, j+1} & \cdots & a_{nn} \end{matrix} \right | \qquad j = 1, 2, \cdots, n Dj=a11a21an1a1,j1a2,j1an,j1b1b1b1a1,j+1a2,j+1an,j+1a1na2nannj=1,2,,n
如果线性方程组无解或有两个不同的解, 则起系数行列式 D = 0 D=0 D=0.

齐次线性方程组的解: 非零解的存在性

  • 如果线性方程组 D x = b Dx=b Dx=b中, b = 0 b=0 b=0, 则称该方程组为齐次线性方程组, 此时 x = 0 \bm{x}=0 x=0为该方程组的一个解, 称零解; 若齐次线性方程组存在非零解, 则必有 D = 0 D=0 D=0. 即 D = 0 D=0 D=0是齐次线性方程组有非零解的充要条件.
  • 此外, 齐次线性方程组的非零解必定是无穷多解

2 矩阵代数

基本概念:

  • 矩阵
  • n阶方阵
  • 列矩阵(1维向量)
  • 1行n列的矩阵称为行阶方阵
  • 零矩阵
  • 同型矩阵, 行列数相等的矩阵
  • 相等矩阵, 对应元素相等的同型矩阵

广义上, 矩阵是一个线性变换

  • 线性变换(矩阵是一个变换, 且是线性变换)
  • 恒等变换: 单位矩阵
  • 对角变换: 对角矩阵

矩阵运算

  • 矩阵加法 (必须为同型矩阵)
  • 数乘运算 (和行列式区别开, 矩阵乘以一个数, 等于矩阵中的每个元素乘以该数)
  • 矩阵乘法: 矩阵 A A A乘以矩阵 B B B, 必须满足 A A A列数= B B B行数; 一般矩阵乘法不满足交换律; 矩阵乘法满足结合律
  • 幂运算
  • 矩阵转置
    - 性质1: 主对角元素不变
    - 性质2: ( A T ) T = A (A^T)^T=A (AT)T=A
    - 性质3: ( A + B ) T = A T + B T (A+B)^T = A^T + B^T (A+B)T=AT+BT
    - 性质4: : ( k A ) T = k A T (kA)^T = kA^T (kA)T=kAT
    - 性质5: ( A B ) T = B T A T (AB)^T = B^TA^T (AB)T=BTAT
  • 共轭矩阵: 若矩阵 A = ( a i j ) A=(a_{ij}) A=(aij)元素为复数, A ‾ = ( a i j ‾ ) \overline{A} = (\overline{a_{ij}}) A=(aij)称为矩阵 A A A的共轭矩阵
    - 性质1: A ‾ ‾ = A \overline{\overline{A} } = A A=A
    - 性质2: A + B ‾ = A ‾ + B ‾ \overline{A+B} = \overline{A} + \overline{B} A+B=A+B
    - 性质3: k A ‾ = k ‾ A ‾ \overline{kA} = \overline{k}\overline{A} kA=kA
    - 性质4: A B ‾ = A ‾ B ‾ \overline{AB} = \overline{A}\overline{B} AB=AB

对称矩阵和反对称矩阵

A T = A A^T = A AT=A, 则称 A A A对称矩阵, 若 A T = − A A^T = - A AT=A, 则称 A A A反对称矩阵
推论: 任意方阵 A A A可以写成对称矩阵和反对称矩阵之和, 具体有 A = 1 2 ( A + A T ) + 1 2 ( A − A T ) A = \frac{1}{2}(A + A^T) + \frac{1}{2}(A-A^T) A=21(A+AT)+21(AAT)
显然, 右边的两项分别是对称矩阵反对称矩阵.

方阵行列式

方阵元素构成的行列式, 记为 ∣ A ∣ |A| A d e t A detA detA
性质:

  • d e t A T = d e t A det A^T = det A detAT=detA
  • d e t ( k A ) = k n d e t A det(kA) = k^n det A det(kA)=kndetA
  • d e t ( A B ) = d e t A d e t B det(AB) = detA detB det(AB)=detAdetB

伴随矩阵

由行列式 d e t A detA detA各元素组成的代数余子式 A i j A_{ij} Aij构成的矩阵的转置 A ∗ = ( A i j ) T A^* = (A_{ij})^T A=(Aij)T称为 A A A伴随矩阵, A ∗ = [ A 11 A 21 ⋯ A n 1 A 12 A 22 ⋯ A n 2 ⋮ ⋮ ⋮ A 1 n A n 2 ⋯ A n n ] = ( A i j ) T A^* = \left[ \begin{matrix} A_{11} & A_{21} & \cdots & A_{n1} \\ A_{12} & A_{22} & \cdots & A_{n2} \\ \vdots & \vdots & &\vdots \\ A_{1n} & A_{n2} & \cdots &A_{nn} \end{matrix} \right] = (A_{ij})^T A=A11A12A1nA21A22An2An1An2Ann=(Aij)T
- 性质1: A ∗ A = A ∗ A = ∣ A ∣ E A^*A = A^*A=|A|E AA=AA=AE
- 性质2: A A A可逆, 当且仅当 A ∗ A^* A可逆
- 性质3: 若 A A A可逆, 则 ( A − 1 ) ∗ = ( A ∗ ) − 1 (A^{-1})^* = (A^*)^{-1} (A1)=(A)1

逆矩阵: 对于 n n n阶方阵

逆矩阵存在则唯一

矩阵 A A A可逆, 当且仅当 ∣ A ∣ ≠ 0 |A|\neq 0 A=0, 且 A − 1 = A ∗ ∣ A ∣ A^{-1} = \frac{A^*}{|A|} A1=AA 式中 A ∗ A^* A为伴随矩阵

∣ A ∣ = 0 |A|=0 A=0时, A A A奇异矩阵, ∣ A ∣ ≠ 0 |A|\neq 0 A=0时, 称非奇异矩阵, 非奇异矩阵存在逆矩阵.

逆矩阵的性质

  • 可逆矩阵的逆矩阵可逆, 且为原矩阵
  • A A A可逆且 k ≠ 0 k\neq 0 k=0, 则 ( k A ) − 1 = 1 k A − 1 (kA)^{-1} = \frac{1}{k}A^{-1} (kA)1=k1A1
  • A A A可逆则 A T A^T AT可逆, 且有 ( A T ) − 1 = ( A − 1 ) T (A^T)^{-1} = (A^{-1})^T (AT)1=(A1)T
  • A , B A, B A,B为同阶矩阵且可逆, 则 A B AB AB可逆, 且有 ( A B ) − 1 = B − 1 A − 1 (AB)^{-1} = B ^{-1}A^{-1} (AB)1=B1A1
  • 矩阵和逆矩阵的行列式互为倒数
  • 可逆对称矩阵的逆矩阵是对称矩阵, 可逆反对称矩阵的逆矩阵也是反对称矩阵
  • ∣ A ∣ ≠ 0 |A|\neq 0 A=0, 定义 { A 0 = E , A − k = ( A − 1 ) k \left \{ \begin{aligned} &A^0 = E, \\ & A^{-k}=(A^{-1})^k\end{aligned}\right. {A0=E,Ak=(A1)k

分块矩阵

若分块矩阵呈对角化, 则可以通过求对角线上的分块矩阵的逆, 来求矩阵的逆.

k k k阶子式

k k k阶子式: 任取 m × n m\times n m×n阶矩阵中的 k k k行, k k k ( k ≤ m , k ≤ n ) (k\leq m, k \leq n) (km,kn)位于行列式交叉处的 k 2 k^2 k2个元素, 不改变相对位置而得到的 k k k阶行列式, 称为 A A A k k k阶子式, 共有 C m k C ˙ n k C_m^k\dot C_n^k CmkC˙nk k k k阶子式

矩阵的秩

设矩阵 A A A中有一个不等于0的 r r r阶子式 D D D, 且所有 r + 1 r+1 r+1阶子式全为0, 则称 r r r矩阵 A A A的秩, 记为 R ( A ) R(A) R(A)
矩阵在初等变换中, 秩保持不变, 即如果 A ∼ B A \sim B AB, 则 R ( A ) = R ( B ) R(A) = R(B) R(A)=R(B)

求矩阵的秩:

对矩阵进行初等行变换为行阶梯矩阵, 则行阶梯矩阵中非0元素的行数即为次矩阵的秩

矩阵方程: 讨论 A x = b Ax=b Ax=b的系数矩阵和增广矩阵和方程解的关系
齐次矩阵方程

n元齐次矩阵方程 A x = 0 Ax = 0 Ax=0有非零解的充要条件是矩阵 A A A的秩 R ( A ) < n R(A) < n R(A)<n,其中 A A A m × n m \times n m×n矩阵. 这是对Gramer法则的推广, 因为Gramer法则仅适用于 m = n m=n m=n的情况.

非齐次矩阵方程

n元非齐次矩阵方程 A x = b Ax = b Ax=b 有非零解的充要条件是矩阵 A A A增广矩阵 [ A    b ] [A~~b] [A  b]的秩相等, 其中 A A A m × n m \times n m×n矩阵.

n n n维向量

线性组合, 基, 线性相关, 线性无关,

3 线性方程组

  可以通过Gramer法则和求矩阵的逆两种方式求解线性方程组, 但矩阵求逆计算量过大, Gramer法则无法求解 D = 0 D=0 D=0时的方程, 因此需要讨论其他方法.

等价矩阵

两个 m × n m \times n m×n阶矩阵 A A A B B B, 如果存在可逆矩阵 P , Q P, Q P,Q, 满足 B = Q A P B= QAP B=QAP, 那么称矩阵 A , B A, B A,B等价, 记为 A ∼ B A\sim B AB, 即 A A A经过有限次初等变换得到 B B B.

初等矩阵

由单位矩阵 E E E经过一次初等变换得到的矩阵称为初等矩阵, 三种基本初等变换对应三种不同的初等矩阵:
(1) 两行对调, 或两列对调 – 换法矩阵
(2) 用非零数c乘以矩阵的某一行 – 倍法矩阵, 特点是 ( i , i ) (i, i) (i,i)元为 c c c
(3) 将矩阵的第 j j j行乘非零数 k k k, 加到第 i i i行 – 消法矩阵, 特点是 ( i , j ) (i, j) (i,j)元为 k k k

初等行变换和初等列变换

A A A m × n m \times n m×n阶矩阵, 对 A A A初等行变换等于左乘 m m m阶初等矩阵, 对 A A A初等列变换相当于右乘 n n n阶初等矩阵.

一切线性变换都可以用矩阵表示

A A A为可逆矩阵, 则存在有限个初等矩阵 P 1 , P 2 , ⋯   , P l P_1, P_2, \cdots, P_l P1,P2,,Pl使得 A = P 1 P 2 ⋯ P l A=P_1P_2\cdots P_l A=P1P2Pl

m × n m \times n m×n 阶矩阵 A ∼ B A \sim B AB的充要条件是 P A Q = B PAQ=B PAQ=B式中, P P P m m m阶矩阵, Q Q Q n n n阶矩阵.

利用初等变换求逆

将可逆矩阵 n n n A A A n n n阶单位矩阵 E E E组合为 n × 2 n n\times 2n n×2n阶矩阵 [ A    E ] [A~~E] [A  E], 对该矩阵进行初等行变换, 将 A A A变为 E E E, 则原来 E E E的位置就变为了 A − 1 A^{-1} A1.

上述方法可以推广到求 A − 1 B A^{-1}B A1B, 即只要对 [ A    B ] [A ~~ B] [A  B]进行初等行变换, 当 A A A变为 E E E时, B B B的位置即变为 A − 1 B A^{-1}B A1B

齐次线性方程组的基础解系

解向量, 满足 A x = 0 Ax=0 Ax=0的向量 ξ \xi ξ

解向量构成的集合, 对于加法和数乘运算封闭, 构成解空间

解空间的维度: n − R ( A ) n - R(A) nR(A)

解空间的一组, 称为齐次方程的基础解系, 基础解系的线性组合形式称为方程的通解

( A T A ) x = 0 (A^TA)x=0 (ATA)x=0 A x = 0 Ax=0 Ax=0方程同解 → \rightarrow R ( A T A ) = R ( A ) R(A^TA)=R(A) R(ATA)=R(A)

非齐次线性方程组的解

η 1 \eta_1 η1 η 2 \eta_2 η2 A x = b Ax=b Ax=b的解, 则 η 1 − η 2 \eta_1 - \eta_2 η1η2 A x = 0 Ax=0 Ax=0的解

η \eta η是非齐次线性方程组的解, ξ \xi ξ是齐次线性方程组的解, 则 η + ξ \eta+\xi η+ξ是非齐次线性方程组的解

矩阵迭代法

考虑矩阵方程 A x = b Ax=b Ax=b A = B − C A=B-C A=BC 其中 B B B为非奇异矩阵, 则原矩阵方程可以写成 B x = b + C x Bx=b+Cx Bx=b+Cx 迭代形式可以写成 B x k + 1 = C x k + b Bx^{k+1} = Cx^k + b Bxk+1=Cxk+b 更进一步 x k + 1 = L x k + b x^{k+1} = Lx^k + b xk+1=Lxk+b 其中 L = B − 1 C L=B^{-1}C L=B1C, 上述迭代法收敛的条件是 lim ⁡ k → ∞ x k = x \lim_{k\rightarrow \infty} x^k = x klimxk=x 根据迭代公式有 x ( k ) − x = L k ( x ( 0 ) − x ) → 0 x^{(k)} - x = L^k (x^{(0)} - x) \rightarrow 0 x(k)x=Lk(x(0)x)0 因为 x ( 0 ) x^{(0)} x(0)可以使任意初值, 因此, 必有 lim ⁡ k → ∞ L k = 0 \lim_{k\rightarrow \infty} L^k = 0 klimLk=0

Gauss-Seidel迭代

B = D − E , C = F B = D-E, \quad C=F B=DE,C=F, 其中 D D D为对角矩阵 E E E为下三角矩阵, F F F为上三角矩阵

Jacobi迭代

B = D , C = E + F B = D, \quad C=E+F B=D,C=E+F, 其中 D D D为对角矩阵 E E E为下三角矩阵, F F F为上三角矩阵

Newton迭代

一种迭代求逆矩阵的迭代算法, 迭代公式为 x k + 1 = x k ( 2 E − A x k ) x_{k+1} = x_k(2E - Ax_k) xk+1=xk(2EAxk)

你可能感兴趣的:(理学)