四、网格构建与离散化：时间步长与空间步长选择、网格边界处理、均匀与非均匀网格
有限差分法的本质是将连续的偏微分方程转化为离散的代数方程组，而这一转化的核心载体就是网格。网格的构建质量直接决定了数值解的精度、稳定性和计算效率。在金融期权定价中，我们面对的是倒向随机偏微分方程（如Black-Scholes方程），网格设计不仅影响误差水平，还关系到边界条件的施加方式与计算资源的消耗。本章将系统讲解时间步长与空间步长的选取原则、边界处理策略以及均匀与非均匀网格的适用场景，帮助你在实战中构建稳健高效的定价网格。
4.1 核心原理：网格参数与误差来源
有限差分法的误差主要来源于三个方面：截断误差、边界误差和舍入误差。截断误差由差分格式替代微分算子时的高阶项忽略产生，其大小与步长直接相关。对于显式格式，时间方向采用前向差分，截断误差为O(Δt)；空间方向采用中心差分，截断误差为O(ΔS²)。因此，空间步长的缩小对精度提升更为显著，但代价是计算量呈平方级增长。
网格设计需要平衡三个核心矛盾：精度与效率（步长越小精度越高但计算越慢）、稳定性与显隐式选择（显式格式受CFL条件约束，隐式格式无条件稳定但每步计算量更大）、局部精度与全局精度（在行权价附近、障碍价格附近等关键区域需要更细的网格）。
4.2 时间步长与空间步长的选择方法
时间步长Δt和空间步长ΔS的选择必须遵循以下步骤：
- 确定空间边界Smax：通常取Smax = 3~4倍的行权价K，或Smax = K × max(1, S₀/K × e^(rT)) × 3。对于奇异期权（如障碍期权），Smax还需覆盖障碍价格B的1.5~2倍。
- 确定空间步长ΔS：先根据所需精度确定网格节点数N，令ΔS = Smax / N。一般N取200~500，美式期权或奇异期权建议N≥400。
- 确定时间步长Δt：若采用显式格式，需满足稳定性条件：
Δt ≤ ΔS² / (2σ²Smax²) （显式格式CFL条件）
实际取 Δt = 0.8 × ΔS² / (2σ²Smax²) 作为安全系数若采用隐式格式（如Crank-Nicolson），则无稳定性限制，但需根据时间精度要求确定M：Δt = T/M，通常M取100~500。Crank-Nicolson格式在时间方向具有二阶精度，是金融定价中最常用的格式。
| 格式类型 | 时间精度 | 空间精度 | 稳定性条件 | 适用场景 |
|---|---|---|---|---|
| 显式（前向欧拉） | O(Δt) | O(ΔS²) | Δt ≤ ΔS²/(2σ²S²) | 简单欧式期权、教学演示 |
| 隐式（后向欧拉） | O(Δt) | O(ΔS²) | 无条件稳定 | 美式期权、奇异期权 |
| Crank-Nicolson | O(Δt²) | O(ΔS²) | 无条件稳定（可能振荡） | 欧式/美式期权标准选择 |
4.3 网格边界处理策略
边界处理是有限差分法中最容易出错也最影响精度的环节。在空间方向，我们需要在S=0和S=Smax处施加边界条件；在时间方向，我们需要在t=T处（到期日）施加终端条件。
S=0处的边界：当标的价格为0时，期权价值满足Dirichlet边界条件。对于看涨期权，V(0,t)=0；对于看跌期权，V(0,t)=K·e^(-r(T-t))。这些条件可以直接施加在网格边界节点上。
S=Smax处的边界：当标的价格远高于行权价时，看涨期权价值近似等于S - K·e^(-r(T-t))，此时可采用线性外推或Dirichlet条件。对于看跌期权，当S=Smax时，V(Smax,t)≈0。更精细的做法是采用二阶外推：
V(Smax, t) = 2V(Smax-ΔS, t) - V(Smax-2ΔS, t) （线性外推）
或 V(Smax, t) = V(Smax-ΔS, t) + ΔS × ∂V/∂S|_{Smax-ΔS} （一阶近似）障碍期权边界：对于敲出障碍期权，当标的价格触及障碍价格B时，期权价值归零（或支付rebate）。此时需要在网格中精确设置B对应的节点位置，若B不在网格节点上，需通过插值处理。
4.4 均匀与非均匀网格的实战选择
均匀网格实现简单、代码易调试，但存在明显的效率缺陷：在远离行权价的区域，网格精度过剩；在行权价附近，精度不足。非均匀网格通过在关键区域加密、在非关键区域稀疏化，可以在相同节点数下显著提升精度。
常用的非均匀网格构造方法包括：
- 对数变换法：令x = ln(S)，在x空间使用均匀网格。由于Black-Scholes方程在对数变换后变为常系数方程，这种方法既简化了方程形式，又天然实现了S空间的不均匀分布（S小时网格密，S大时网格疏）。
- 分段均匀法：将S轴分为若干区间，每个区间内使用不同的ΔS。例如在[0, K/2]使用ΔS=2，在[K/2, 2K]使用ΔS=0.5，在[2K, Smax]使用ΔS=3。
- 正弦/余弦变换法：通过S = K + A·sinh(c·ξ)等映射函数，将均匀网格ξ映射到非均匀S网格，在K附近自然加密。
| 网格类型 | 优点 | 缺点 | 推荐场景 |
|---|---|---|---|
| 均匀网格 | 实现简单、误差分析容易、代码调试方便 | 节点利用率低、关键区域精度不足 | 教学演示、快速原型、精度要求不高的场景 |
| 对数变换网格 | 方程简化、S小区域自动加密、数值稳定性好 | 需处理变换后的边界条件、S=0处奇异 | 欧式/美式期权标准选择 |
| 分段均匀网格 | 灵活控制局部精度、实现相对简单 | 区间交界处需插值处理、可能引入额外误差 | 障碍期权、含多个奇异点的产品 |
4.5 案例拆解：障碍期权定价的网格设计实战
考虑一个向下敲出看涨期权（Down-and-Out Call）：标的价格S₀=100，行权价K=105，障碍价格B=95，无风险利率r=5%，波动率σ=30%，期限T=1年。我们需要设计网格并完成定价。
初始参数与约束：障碍价格B=95非常接近当前价格S₀=100，这意味着网格必须在B附近有足够的分辨率，否则障碍是否被触发将产生显著误差。同时，行权价K=105也需要加密。
网格设计决策：
- 设定Smax = 300（约为K的3倍），Smin = 0。
- 采用分段均匀网格：在[0, 90]区间用ΔS=5（20个节点），在[90, 110]区间用ΔS=0.5（40个节点），在[110, 300]区间用ΔS=5（38个节点），总节点数N=100。
- 时间方向采用Crank-Nicolson格式，M=200步，Δt=0.005年。
- 边界条件：S=0处V=0；S=300处采用线性外推；S=B=95处设置V=0（敲出条件）。
结果对比：使用均匀网格（ΔS=3，N=100）定价结果为4.82；使用上述非均匀网格定价结果为4.67。解析解（通过反射原理计算）为4.71。非均匀网格的误差为0.85%，而均匀网格的误差为2.3%。在相同节点数下，非均匀网格将误差降低了约60%。
风险与代价：非均匀网格在区间交界处（S=90和S=110）需要处理差分格式的修正，若处理不当会引入局部振荡。此外，当标的价格S₀恰好落在加密区间内时，定价精度最高；若S₀偏离加密区间，则优势减弱。
4.6 常见误区与总结复盘
误区一：盲目追求小步长。步长过小会导致计算量急剧增加，且当Δt和ΔS过小时，舍入误差可能超过截断误差，总误差反而增大。建议通过网格收敛性检验（halving方法）确定最优步长。
误区二：忽略时间步长与空间步长的匹配关系。即使使用隐式格式，Δt过大也会导致时间方向精度不足。经验法则：Δt应满足Δt ≤ ΔS²/(4σ²S²)以保证Crank-Nicolson格式不产生振荡。
误区三：边界条件施加不当。在Smax处使用Dirichlet条件（V=0）对于看涨期权会严重低估期权价值。必须使用线性外推或精确的渐近解作为边界条件。
误区四：非均匀网格的过度设计。对于简单欧式期权，非均匀网格带来的精度提升有限，却增加了实现复杂度。建议先用均匀网格验证模型正确性，再根据需求引入非均匀网格。