用小波方法求解Ermakov-Pinney方程的数值解

《Applications in Engineering Science》：A wavelet approach to the solution of Ermakov-Pinney equation

　　

摘要：Ermakov–Pinney（EP）方程因其精确可解性及与时变谐振子(time-dependent harmonic oscillator)的联系而广为人知，但对于含复杂或时变参数之情形，数值方法是必需的。本研究采用Haar小波方法(Haar Wavelet Method, WHM)求解多种形式的EP方程，包括初值问题(initial value problem)、边值问题(boundary value problem)及含时变频率(time-dependent ω(t))与耦合边界条件(coupled boundary conditions)的情形。研究人员通过求解一系列初边值问题验证该方法之有效性，并将数值结果与精确解及四阶Runge–Kutta(Runge–Kutta, RK)方法所得结果进行比较。据研究人员所知，这是首次将Haar小波方法系统性地应用于含时变频率及耦合边界条件的EP方程。研究结果表明，Haar小波方法不仅易于实现，且在精度上优于Runge–Kutta方法。

论文解读：用Haar小波方法求解Ermakov–Pinney（EP）方程的数值解

一、研究背景与意义

Ermakov–Pinney（EP）方程 ?″+ω2(t)y = k/y3（其中y(t)为未知函数，ω(t)为时变频率，k为常数）源于1880年Ermakov的工作并经Pinney(1950)推广，其与线性时变谐振子之解存在显式联系——若x?(t)、x?(t)为对应线性齐次方程?+ω2(t)x=0之两个线性无关解且W为其Wronskian，则EP方程之通解为 y(t)=√[Ax?2(t)+Bx?2(t)+2Cx?(t)x?(t)]，常数满足 C2?AB=k/W2。该方程在量子场论、宇宙学Friedmann–Robertson–Walker模型中之尺度因子演化、非线性光学及弹性力学中均有重要应用。

尽管特殊情形下EP方程可精确求解，但当ω(t)为任意时变函数或涉及耦合边界条件时，解析解通常不存在，须借助数值方法。传统数值方法如有限差分、Runge–Kutta（RK）、打靶法(shooting method)在处理长时间积分或快速振荡解时易积累误差。小波方法尤其Haar小波因具紧支撑(compact support)、正交性及稀疏表示能力，在ODE/PDE数值求解中展现优势，但尚未被系统应用于EP方程。本文首次将Haar小波配置法(Haar Wavelet Collocation Method, HWCM)统一框架应用于EP方程之初值问题(IVP)、边值问题(BVP)及含时变频率和耦合边界条件之情形，并与RK方法对比验证精度。

二、主要关键技术方法

研究人员建立Haar小波方法数值框架：定义区间[0,1]上Haar尺度函数h?(x)=φ(x)（单位方波）及小波函数h?(x)(i=2,…,2^J+1)，其中J为分辨率水平(resolution level)，M=2^J为最大小波指标相关参数；构造Haar函数之n重积分p_i,n(x)=∫?? p_i,n?1(t)dt（p_i,0=h?(x)）。对于二阶ODE y″=f(x,y,y′)，设y″(x)=∑_i=1^2Ma?h?(x)，逐次积分代入初值条件得 y′(x)=η?+∑a?p_i,1(x) 及 y(x)=η?+η?x+∑a?p_i,2(x)；在配置点(collocation points) x_r=(r?0.5)/(2M) (r=1,…,2M)处离散原方程得到关于Haar系数a?之代数方程组（非线性情形用迭代求解）；边值问题通过引入未知初值y′(0)并利用边界条件确定之。误差评估采用L_∞范数（最大绝对误差）及最大相对误差(Max. Rel. Err.)。所有计算使用MATLAB 2019a实现。

三、研究结果

2. The Haar wavelet method（Haar小波方法理论框架）

研究人员详细阐述了L2(R)空间上多分辨分析(multiresolution analysis, MRA)、Haar尺度函数φ(x)与小波母函数ψ(x)之定义，给出区间[0,1]上归一化Haar函数族{h?(x)}之分段表达式及索引关系 i=2^j+k+1（j为分辨率级，k为平移参数，j=0,1,…,J；k=0,…,2^j?1），推导了一重与二重积分p_i,1(x)、p_i,2(x)之闭式分段形式，并说明如何将二阶ODE之二阶导数展开为Haar级数后通过双重积分结合初/边值条件转化为原函数近似表达式，再经配置点离散获得待求Haar系数之代数方程组。以非线性初值问题 y″?(y′)2+y′=0, y(0)=?ln2, y′(0)=1/2 验证方法收敛性，结果显示随J增大L_∞与最大相对误差迅速下降（J=8时L_∞≈7.15×10??），证实Haar小波方法对非线性ODE有效且较RK更精确。

3. Solution of the EP equation（EP方程之求解与算例）

(A) Initial value problems for the EP equation（EP方程初值问题）

• •
  Example 3.1：常频ω(t)=1，c=0.6，初值y(0)=1, y′(0)=0，精确解 y2=cos2t+0.6sin2t。WHM与精确解吻合良好，绝对误差小于RK方法。
• •
  Example 3.2：时变频率 ω(t)=1?3/(4t2)（实际计算域避开奇点），c=0.6，相同初值。此为含时变参数之非平凡情形，WHM与RK均给出合理数值解，WHM误差更小。
• •
  Example 3.3（阻尼Pinney方程 Damped Pinney Equation）：?″+2ε?′+ω2(t)y=c/y3，ε=0.1, ω=1, c=0.6，初值同前。WHM成功处理一阶导项，数值解光滑，绝对误差低于RK。

(B) Boundary value problems for the EP equation（EP方程边值问题）

• •
  Example 3.4：ω2(t)=t2, c=0.6，边值y(0)=1, y(1)=1.5，无精确解。WHM给出稳定数值解。
• •
  Example 3.5：ω2(t)=e^t, c=0.2，边值y(0)=1, y(1)=1.1，无精确解。WHM给出数值解。
• •
  Example 3.6：常频ω2=4, c=4，边值y(0)=1, y(π/4)=√2，精确解 y(t)=√[1+sin2(2t)+sin(4t)/2]^1/2（原文表述）。WHM数值解与精确解符合良好，验证BVP框架下方法准确性。
• •
  Example 3.7（Coupled Pinney Equations 耦合Pinney方程组）：
  y?″+y?=1/y?3+y?, y?(0)=1, y?(π/4)=1；
  y?″+y?=1/y?3+y?, y?(0)=1, y?(π/4)=1。
  无精确解，WHM获稳定耦合数值解。

四、讨论与结论

研究人员指出，Haar小波配置法可统一处理EP方程之初值问题与边值问题（含时变频率及耦合边界条件），无需线性化或变换，保留了原方程物理结构。数值对比表明WHM在实现简便性与计算精度上均优于经典Runge–Kutta方法，尤其在BVP及无精确解情形下仍给出可靠数值解；误差分析显示提升分辨率J可使L_∞与最大相对误差按数量级下降，证实方法收敛。据研究人员所知，此为首次将Haar小波方法系统性应用于含时变频率及耦合边界条件之EP方程。

结论(Conclusions)： 研究人员将Haar小波方法应用于各类Ermakov–Pinney方程（含时变ω(t)、阻尼项及耦合Pinney方程组），通过初值、边值算例验证了该方法之简便性与高精度，数值结果优于Runge–Kutta方法。此工作为EP方程及广义非线性ODE之小波数值求解提供了新途径，未来可扩展至分数阶EP方程(fractional-order EP equation)或采用高阶小波及其它小波基（如ultraspherical、Bernoulli、biorthogonal spline wavelets）进行对比研究。