计算方法复习提纲-下

第7章 解线性方程组的直接方法

Guass消去法

• 消元
  • $G^{(1)}=\left[\begin{array}{ccccc}1& 0& 0& & 0\\ -\frac{a_{21}}{a_{11}}& 1& 0& & 0\\ -\frac{a_{31}}{a_{11}}& 0& 1& & 0\\ ⋮& & & \ddots & \\ -\frac{a_{n1}}{a_{11}}& 0& 0& & 1\end{array}\right]$
  • $A^{(1)}=A$
  • $A^{(2)}=G^{(1)}{A}^{(1)}$
  • $A^{(2)}=\left[\begin{array}{cc}{a}_{1,1}^{(2)}& a_{1,2:n}^{(2)}\\ 0& A_{2:n,2:n}^{(2)}\end{array}\right]$
  • $A^{(n)}=G^{(n-1)}{A}^{(n-1)}$
• 回代
  • 进一步把$A$转化成对角矩阵
• 实现Gauss消去法时不需要交换行
  • $a_{kk}^{(k)}\ne 0⟺\det (A_{1:k,1:k})\ne 0$
    • 非奇异线性变换不改变行列式的值
    • 推论
      • $a_{11}^{(1)}=\det (A_{1:1,1:1})$
      • $a_{kk}^{(k)}=\frac{\det (A_{1:k,1:k})}{\det (A_{1:k-1,1:k-1})}$
• LU分解
  • $A^{(n)}=G^{(n)}\cdots G^{(1)}{A}^{(1)}$
  • $L^{-1}=G^{(n)}\cdots G^{(1)}$
  • $U=A^{(n)}$
  • $a_{kk}^{(k)}\ne 0$时是唯一的
    • $A$非奇异时, $L_1{U}_1=L_2{U}_2$, $I=L_1^{-1}{L}_2=U_1{U}_2^{-1}$.
• 计算量
  • 加减 $\sum _{i=1}^n{i}^2=\frac{n^3}{3}$ 次
  • 乘除 $\sum _{i=1}^n{i}^2=\frac{n^3}{3}$ 次

列主元素消去法

• 数值稳定性
  • $a_{kk}^{(k)}$ 作为分母不宜太小
• 完全主元素
  • $A_{k:n,k:n}^{(k)}$ 中绝对值最大元素交换到 $a_{kk}^{(k)}$
• 算法
  • 乘数 $\frac{a_{ki}}{a_{kk}}$ 冲掉 $a_{ki}$
  • 消元结果 $A_{k+1:n,k:n}$ 冲掉 $A$
  • 计算结果 $x$ 冲掉 $b$
  • $a_{i_{\max }{j}_{\max }}=0$ 终止
  • 行交换
    • 同时交换 $A$ 和 $b$
  • 列交换
    • 交换 $\mathrm{MapIndex}$
• 列主元素
  • 寻找完全主元素代价较大
  • 列交换比行交换麻烦
  • 可以额外计算 $\mathtt{detA}$
    • 行交换时 $\mathtt{detA}\leftarrow -\mathtt{detA}$
• PLU分解
  • $\cdots S_{x}G{S}_y\cdots =\cdots (S_{x}G{S}_x)S_x{S}_y\cdots$
  • $U=A^{(n)}={\tilde{L}}^{-1}P{A}^{(1)}$

平方根法

• 对称正定的矩阵
  • $A=LD{L}^T$
  • $D_{ii}>0$
  • $A=C{C}^T$
• 逐列计算
  • $i⩽j$
  • $a_{ij}=\sum _{k=1}^j{c}_{ik}{c}_{jk}=\sum _{k=1}^{j-1}{c}_{ik}{c}_{jk}+c_{jj}{c}_{ij}$
  • $c_{jj}$
  • $c_{ij}$ ($i>j$)
• 数值稳定性
  • $a_{jj}=\sum _{k=1}^j{l}_{jk}^2$
  • $l_{jk}^2<a_{jj}$
• 乘除 $\frac{n^3}{6}$ 次
• 避免开方运算
  • $A=(LD)L^T$
• 逐行计算
  • $i⩽j$
  • a i j = ∑ k = 1 j l i k d k l j k = ∑ k = 1 j − 1 l i