Math ClubMath Club
v1 · padrão canônico

第35课 — 通过矩阵求解方程组

Cramer, escalonamento de Gauss, matriz inversa. Quando cada método é melhor.

Used in: 1.º ano EM (15 anos) · Equiv. Math II japonês · Equiv. Klasse 11 alemã

Ax=bx=A1bA\mathbf{x} = \mathbf{b} \quad \Rightarrow \quad \mathbf{x} = A^{-1}\mathbf{b}
Choose your door

Rigorous notation, full derivation, hypotheses

求解方法

矩阵形式

{a1x+b1y+c1z=d1a2x+b2y+c2z=d2a3x+b3y+c3z=d3\begin{cases} a_1 x + b_1 y + c_1 z = d_1 \\ a_2 x + b_2 y + c_2 z = d_2 \\ a_3 x + b_3 y + c_3 z = d_3 \end{cases}Ax=bA\mathbf{x} = \mathbf{b}A3×3A_{3 \times 3}x,bR3\mathbf{x}, \mathbf{b} \in \mathbb{R}^3

方法1——高斯消元

初等变换(不改变解):

  1. 交换两行。
  2. 一行乘以非零标量。
  3. 一行倍数加到另一行。

目标:将 [Ab][A | \mathbf{b}] 三角化为阶梯形。然后回代。

方法2——克拉默

Ax=bA\mathbf{x} = \mathbf{b}detA0\det A \neq 0xi=detAidetAx_i = \frac{\det A_i}{\det A}

其中 AiA_iAA 的第 ii 列被 b\mathbf{b} 替换。

方法3——逆

x=A1b\mathbf{x} = A^{-1} \mathbf{b}。可通过 [AI][IA1][A | I] \to [I | A^{-1}] 行变换计算 A1A^{-1}

何时使用

方法代价何时使用
克拉默O(n4)O(n^4)n3n \leq 3,教学
高斯O(n3)O(n^3)实际标准
O(n3)O(n^3) + 每方程组 O(n2)O(n^2)多个 b\mathbf bAA
LUO(n3)O(n^3) + 每方程组 O(n2)O(n^2)多方程组,比逆好

方程组分类

  • 确定的(SPD):唯一解(detA0\det A \neq 0,或 rank(A)=rank([Ab])=n\text{rank}(A) = \text{rank}([A|\mathbf b]) = n)。
  • 不定的(SPI):无穷多解(rank(A)=rank([Ab])<n\text{rank}(A) = \text{rank}([A|\mathbf b]) < n)。
  • 不可解(SI):无解(rank(A)<rank([Ab])\text{rank}(A) < \text{rank}([A|\mathbf b]))。

鲁谢-卡佩利定理

Ax=bA\mathbf x = \mathbf b 有解 ⟺ rank(A)=rank([Ab])\text{rank}(A) = \text{rank}([A|\mathbf b])。两者 = 未知数个数时解唯一。

Exercise list

46 exercises · 11 with worked solution (25%)

Application 32Understanding 2Modeling 10Challenge 1Proof 1
  1. Ex. 35.1Application
    用克拉默解:{x+y=52xy=1\begin{cases} x + y = 5 \\ 2x - y = 1 \end{cases}。(答:x=2,y=3x = 2, y = 3。)
  2. Ex. 35.2ApplicationAnswer key
    用消元法解:{3x+2y=11xy=2\begin{cases} 3x + 2y = 11 \\ x - y = 2 \end{cases}。(答:x=3,y=1x = 3, y = 1。)
  3. Ex. 35.3Application
    用消元法解 {x+2yz=42x+y+z=6xy+2z=3\begin{cases} x + 2y - z = 4 \\ 2x + y + z = 6 \\ x - y + 2z = 3 \end{cases}
  4. Ex. 35.4ApplicationAnswer key
    齐次方程组 Ax=0A\mathbf{x} = \mathbf{0}detA=50\det A = 5 \neq 0。解?(答:平凡 x=0\mathbf x = \mathbf 0。)
  5. Ex. 35.5Application
    kk 取何值时方程组 {x+2y=3kx+4y=6\begin{cases} x + 2y = 3 \\ kx + 4y = 6 \end{cases} 有无穷多解?(答:k=2k = 2。)
  6. Ex. 35.6Application
    上述方程组在 kk 取何值时无解?
  7. Ex. 35.7Application
    {2x+y=5x3y=1\begin{cases} 2x + y = 5 \\ x - 3y = 1 \end{cases} 的矩阵形式。计算 A1bA^{-1}\mathbf{b}
  8. Ex. 35.8Application
    用克拉默解 {xy+z=12x+yz=4x+2y+z=0\begin{cases} x - y + z = 1 \\ 2x + y - z = 4 \\ -x + 2y + z = 0 \end{cases}
  9. Ex. 35.9Application
    用消元法解 {2xy=0x+3y=7\begin{cases} 2x - y = 0 \\ x + 3y = 7 \end{cases}。(答:x=1,y=2x = 1, y = 2。)
  10. Ex. 35.10ApplicationAnswer key
    用消元法验证 {x+y+z=32x+2y+2z=63x+3y+3z=9\begin{cases} x + y + z = 3 \\ 2x + 2y + 2z = 6 \\ 3x + 3y + 3z = 9 \end{cases} 有无穷多解。
  11. Ex. 35.11Application
    用逆解:{4x+3y=112x+y=5\begin{cases} 4x + 3y = 11 \\ 2x + y = 5 \end{cases}。(答:x=2,y=1x = 2, y = 1。)
  12. Ex. 35.12Application
    通过 [AI][A|I] 消元法求 (121011001)\begin{pmatrix} 1 & 2 & 1 \\ 0 & 1 & 1 \\ 0 & 0 & 1 \end{pmatrix}A1A^{-1}
  13. Ex. 35.13Application
    方程组 {x+y+z=1x+y+z=2\begin{cases} x + y + z = 1 \\ x + y + z = 2 \end{cases}——解?(答:无解——不相容。)
  14. Ex. 35.14Application
    4 个方程 2 个未知数——通常超定,无精确解。实际中用最小二乘。
  15. Ex. 35.15Application
    2 个方程 4 个未知数——欠定,无穷多解。
  16. Ex. 35.16ApplicationAnswer key
    {0,1x+0,2y=0,30,4x0,5y=0,1\begin{cases} 0{,}1 x + 0{,}2 y = 0{,}3 \\ 0{,}4 x - 0{,}5 y = 0{,}1 \end{cases}——先乘以10。
  17. Ex. 35.17Application
    {x+yz=02xy+z=0\begin{cases} x + y - z = 0 \\ 2x - y + z = 0 \end{cases}(2x3 齐次)的通解。
  18. Ex. 35.18Application
    证明齐次解 + 非齐次特解 = 通解。
  19. Ex. 35.19Application
    检查相容性:{x+y=32x+2y=7\begin{cases} x + y = 3 \\ 2x + 2y = 7 \end{cases}。(答:不相容。)
  20. Ex. 35.20Application
    克拉默给 x=Dx/Dx = D_x/DDD 取何值时方法失败?(答:D=0D = 0。)
  21. Ex. 35.21Application
    用消元法解 {x+2y+3z=64x+5y+6z=157x+8y+10z=25\begin{cases} x + 2y + 3z = 6 \\ 4x + 5y + 6z = 15 \\ 7x + 8y + 10z = 25 \end{cases}
  22. Ex. 35.22Application
    用克拉默解 {2xy+3z=9x+yz=03x+2y+z=5\begin{cases} 2x - y + 3z = 9 \\ x + y - z = 0 \\ 3x + 2y + z = 5 \end{cases}
  23. Ex. 35.23ApplicationAnswer key
    kk 取何值时方程组 {x+y+z=1x+2y+kz=2x+4y+k2z=4\begin{cases} x + y + z = 1 \\ x + 2y + kz = 2 \\ x + 4y + k^2 z = 4 \end{cases} 有唯一解?
  24. Ex. 35.24Application
    求所有满足 (123456)x=(00)\begin{pmatrix} 1 & 2 & 3 \\ 4 & 5 & 6 \end{pmatrix} \mathbf x = \begin{pmatrix} 0 \\ 0 \end{pmatrix}x\mathbf x
  25. Ex. 35.25Application
    用回代解阶梯形 {x+2yz=33y+2z=15z=10\begin{cases} x + 2y - z = 3 \\ 3y + 2z = 1 \\ 5z = 10 \end{cases}
  26. Ex. 35.26Application
    对方程组 {x+ky=1kx+y=1\begin{cases} x + ky = 1 \\ kx + y = 1 \end{cases},按 kk 分类。
  27. Ex. 35.27ApplicationAnswer key
    {x+2yz=12x+4y2z=2x+y+z=3\begin{cases} x + 2y - z = 1 \\ 2x + 4y - 2z = 2 \\ x + y + z = 3 \end{cases}——注意冗余。
  28. Ex. 35.28ApplicationAnswer key
    通过高斯-若尔当法求 A=(123014001)A = \begin{pmatrix} 1 & 2 & 3 \\ 0 & 1 & 4 \\ 0 & 0 & 1 \end{pmatrix}A1A^{-1}
  29. Ex. 35.29ApplicationAnswer key
    用 35.28 得到的 A1A^{-1}Ax=(6,5,1)TA\mathbf x = (6, 5, 1)^T
  30. Ex. 35.30Application
    4×44 \times 4{x1+x2=1x2+x3=2x3+x4=3x1+x4=4\begin{cases} x_1 + x_2 = 1 \\ x_2 + x_3 = 2 \\ x_3 + x_4 = 3 \\ x_1 + x_4 = 4 \end{cases}
  31. Ex. 35.31Application
    识别矩阵 (123246111)\begin{pmatrix} 1 & 2 & 3 \\ 2 & 4 & 6 \\ 1 & 1 & 1 \end{pmatrix} 的秩。
  32. Ex. 35.32Application
    将鲁谢-卡佩利应用于 {x+y=12x+2y=3\begin{cases} x + y = 1 \\ 2x + 2y = 3 \end{cases}——分类。
  33. Ex. 35.33Modeling
    3 网孔电路中,基尔霍夫定律给出 3x3 方程组。建模:3 个网孔电流,电阻,电源。
  34. Ex. 35.34ModelingAnswer key
    在经济学中,IS-LM 模型生成 2x2 方程组:产出 YY 和利率 rr 同时。
  35. Ex. 35.35Modeling
    混合:100 ml 含 3 种物质的溶液。已知浓度——3×33 \times 3 方程组。
  36. Ex. 35.36ModelingAnswer key
    4 节点 3 个未知力的桁架——消元法。
  37. Ex. 35.37Modeling
    在统计中,最小二乘 XTXβ=XTyX^TX\beta = X^Ty 是线性方程组。XMn×pX \in M_{n\times p}npn \gg p,方程组维度?
  38. Ex. 35.38ModelingAnswer key
    在控制中,求 u\mathbf u 使 y=C(sIA)1Bu\mathbf y = C(sI - A)^{-1}B \mathbf u 接近参考。
  39. Ex. 35.39Modeling
    在 CG 中,光的辐射度求解导致方程组 (IF)B=E(I - F)\mathbf B = \mathbf E——稀疏。
  40. Ex. 35.40Modeling
    在二次优化 min12xTHx+cTx\min \frac{1}{2}\mathbf x^T H \mathbf x + \mathbf c^T \mathbf x 中,临界点:Hx=cH\mathbf x = -\mathbf c
  41. Ex. 35.41Modeling
    在 ML 中,岭回归:(XTX+λI)β=XTy(X^TX + \lambda I)\beta = X^T\mathbf y。为何正则化?
  42. Ex. 35.42Modeling
    在 PDE 中,1D 热方程离散化导致三对角方程组——Thomas 算法 O(n)O(n)
  43. Ex. 35.43Understanding
    证明方程组 Ax=0A\mathbf{x} = \mathbf{0} 始终有 x=0\mathbf{x} = \mathbf{0}
  44. Ex. 35.44Understanding
    证明若 AA 可逆,Ax=0A\mathbf{x} = \mathbf{0} 只有 x=0\mathbf{x} = \mathbf{0}
  45. Ex. 35.45Challenge
    用克拉默和高斯解同一 3x3 方程组——比较运算次数。
  46. Ex. 35.46Proof
    证明消元法保持解集——逐次操作。

本课参考来源

Updated on 2026-04-30 · Author(s): Clube da Matemática

Found an error? Open an issue on GitHub or submit a PR — open source forever.