19、边界值问题的数值解法

边界值问题的数值解法

1. 边界值问题概述

在许多实际应用中，除了初始值问题（n 阶微分方程伴随着在自变量同一值处指定的 n 个初始条件），还会遇到边界值问题（BVP）。在边界值问题里，辅助条件通常在自变量的不同值处指定，一般是在系统的端点处。

边界值问题可以通过打靶法或有限差分法进行数值求解：
- 打靶法 ：先通过猜测缺失的初始条件，将边界值问题转化为初始值问题，然后用数值方法（如四阶龙格 - 库塔法，RK4）求解该初始值问题，并检验解是否满足原边界条件。
- 有限差分法 ：把系统区间划分为若干子区间，用有限差分近似代替微分方程中的导数，从而将微分方程转化为代数方程组，再用相关方法求解。

2. 不同类型的边界值问题

2.1 二阶边界值问题

二阶微分方程的一般形式为 $y’’ = f(x, y, y’)$，其中 $a \leq x \leq b$，并受两个边界条件约束，边界条件通常在端点 $a$ 和 $b$ 处指定。由于边界条件可以是 $y$ 的给定值或 $y’$ 的给定值，因此会遇到不同形式的边界条件。

2.2 常见边界条件

• 狄利克雷边界条件 ：给定端点处 $y$ 的值，即 $y(a) = y_a$，$y(b) = y_b$。
• 诺伊曼边界条件 ：给定端点处 $y’$ 的值，即 $y’(a) = y’_a$，$y’(b) = y’_b$。
• 混