【复杂海洋环境下的抛物方程高级求解器：原理与算法体系】第2章 算法内核与数值实现：从理论公式到工程代码

目录2.1 宽角近似的数学基础与Padé系数体系2.2 能量守恒修正的严格推导2.3 自启动器的算子构造2.4 弹性海底的抛物方程矢量形式2.5 海底衰减的复波数建模2.6 三对角系统的快速求解与变深度处理2.7 标准验证案例与精度控制2.8 旋转算子Padé与数值稳定性控制2.9 求解器代码级数据流与模块架构2.1 宽角近似的数学基础与Padé系数体系2.1.1 平方根算子的函数论性质抛物方程的核心困难在于处理伪微分算子 $Q=\sqrt{1+X}$ 。在形式解 $u(r+\Delta r)=\exp[ik_0\Delta r(Q-1)]u(r)$ 中，$X$ 为深度方向的微分-乘法复合算子，其谱分布于复平面的负实轴附近。直接对 $X$ 进行Taylor展开仅在 $|X| \ll 1$ 时收敛，对应小角度传播。为扩展有效角域，需采用有理函数逼近。设 $R_{N,N}(X)$ 为 $N/N$ 型Padé有理逼近，满足：$$\sqrt{1+X} = R_{N,N}(X) + \mathcal{O}(X^{2N+1}) \quad (2.1)$$其中：$$R_{N,N}(X) = 1 + \sum_{j=1}^N \frac{a_j X}{1 + b_