示例:杜芬方程
无阻尼、非强迫运动的杜芬方程为
-
其中t > 0,0 < ε ≪ 1。[4]
假设初值为
-
使用摄动法,假设级数解为x(t) = x0(t) + ε x1(t) + … 。可以得到,级数的前两项为
-
此近似解会随着时间无限地增大,与该方程所描述的物理系统不符。而导致这一原因的是其中的长期项t sin t 随时间而不断增大。为使近似解随时间变化仍然有效,可以采用如下的庞加莱-林德斯泰特方法。
此方法中,不仅近似解本身表示为渐近展开,时间t也表示为级数形式
- 其中
由于解的角频率的领头项为1,取ω0 = 1。于是,原方程变为
-
初值则不变。假设解的形式为 x(τ) = x0(τ) + ε x1(τ) + … ,通过ε的零阶与一阶项可以得到
-
取ω1 = 3⁄8便可消除长期项。按此继续进行分析,便可得到更高阶的精度。以下为精确到ε一阶精度的近似解为
-
参考文献
- ^ Drazin, P.G., Nonlinear systems, Cambridge University Press, 1992, ISBN 0-521-40668-4 , pp. 181–186.
- ^ Poincaré, H., Les Méthodes Nouvelles de la Mécanique Célèste II, New York: Dover Publ., 1957 [1893] , §123–§128.
- ^ A. Lindstedt, Abh. K. Akad. Wiss. St. Petersburg 31, No. 4 (1882)
- ^ J. David Logan. Applied Mathematics, Second Edition, John Wiley & Sons, 1997. ISBN 0-471-16513-1.