55、带状结构线性系统的直接求解方法

带状结构线性系统的直接求解方法

1. 引言

在许多实际问题中，我们常常会遇到线性方程组的求解问题。某些算法产生的线性方程组具有特定的结构，例如由泊松方程离散化得到的线性方程组具有带状结构。这种结构的线性方程组可以利用其稀疏性来设计更高效的求解方法。下面我们将介绍几种针对带状结构线性方程组的求解方法，特别是三对角系统的求解。

2. 带状矩阵与三对角系统

• 带状矩阵定义 ：对于矩阵 (A = (a_{ij}) {i,j = 1,\cdots,n} \in R^{n\times n})，如果存在 (r \in N)，(r \leq n)，使得当 (|i - j| > r) 时，(a {ij} = 0)，则称矩阵 (A) 为带状矩阵，(r) 称为矩阵 (A) 的半带宽。
• 三对角矩阵 ：当 (r = 1) 时，带状矩阵称为三对角矩阵。三对角系统是指系数矩阵为三对角矩阵的线性方程组，例如一维泊松方程离散化后就会得到三对角系统。

3. 三对角系统的高斯消元法

对于三对角矩阵 (A) 的线性方程组 (Ax = y)，可以使用高斯消元法求解。
- 步骤 ：
1. 正向消元（无主元选择）的第 (k) 步（(k < n)）：
- 计算 (l_{ik} := \frac{a_{ik}^{(k)}}{a_{kk}^{(k)}})，其中 (i = k + 1,\cdots,n)。
- 从第 (i) 行（(i = k + 1,\