示例:杜芬方程
無阻尼、非強迫運動的杜芬方程為
-
其中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.