33、数值方法求解边界值与初边值问题

数值方法求解边界值与初边值问题

1. 有限差分逼近

在区间 $[0, 1]$ 上，我们引入网格点 $x_j = jh$，其中 $n \geq 2$ 为整数，$h = 1/n$ 是网格间距。通过用二阶中心有限差分替换 $u’‘(x_j)$，我们可以将问题近似为一个有限序列 ${u_j} {j = 0}^{n}$，满足：
[
-\frac{u {j + 1} - 2u_j + u_{j - 1}}{h^2} = f(x_j), \quad j = 1, \cdots, n - 1
]
且 $u_0 = u_n = 0$。这可以写成矩阵形式 $A_{fd}u = f$，其中 $A_{fd}$ 是对称的 $(n - 1) \times (n - 1)$ 有限差分矩阵，定义为 $A_{fd} = h^{-2}\text{tridiag}_{n - 1}(-1, 2, -1)$。该矩阵是行对角占优的，并且是正定的，这意味着上述线性系统有唯一解。

为了分析稳定性，我们定义离散内积：
[
(w_h, v_h) h = h\sum {k = 0}^{n} c_kw_kv_k
]
其中 $c_0 = c_n = 1/2$，$c_k = 1$ 对于 $k = 1, \cdots, n - 1$。基于此，我们可以证明算子 $L_h$ 是对称且正定的。

对于收敛性分析，我们引入局部截断误差 $\tau_h(x_j) = (L_hu)(x_j) - f(x_j)$。通过泰勒级数展开，我们可以得到 $\tau_h(x_j)$ 的表达式，并证明当 $h \to 0$ 时，$|\tau_h|