A study of approximate analytical solutions of a kind of typical second-order nonlinear different equation
0 引言
在非惯性转动参照系中研究力学体系的运动, 常常会遇到一类分子分母均含非线性项的二阶非线性微分方程[1], 此类方程可以通过Mathematica软件直接得到它的数值解, 但无法得到其近似解析解.文献[2]用同伦渐近法得到了三阶近似周期解, 但求解过程过于繁复, 求解思路比较抽象. Adomian分解法是求解非线性微分方程的一有效方法, 已广泛应用于求非线性常微分方程、非线性偏微分方程的近似解析解中, 并取得了一系列的成果[3-17].此方法的主要思路是先将非线性微分方程分解成线性最高阶常微分部分、一般的线性部分、非线性部分和非齐次部分, 把非线性微分方程的解也分解成无穷个解分量, 然后利用Adomian多项式和逆算符的性质, 由低阶解分量递推求得高阶解分量, 最后将全部解分量相加得到非线性微分方程的近似解析解.用Adomian分解法求非线性微分方程近似解的关键和难点在于求得与非线性项相对应的Adomian多项式, 对于分母不含非线性项的情况比较容易求得其Adomian多项式, 但对于分子分母中都含有非线性项的情况, 其相应的Adomian多项式较难求得.本文首先给出用Adomian分解法求二阶非线性微分方程近似解的基本步骤.其次, 通过研究在非惯性转动参照系中力学体系的运动, 建立一类分子分母均含非线性项的二阶非线性微分方程.第三, 计算与二阶非线性微分方程的非线性项相应的前5项Adomian多项式.第四, 得到在给定初始条件和参数下的近似解析解, 并与直接用Mathematica软件作出的解曲线以及与同伦渐近法[2]所得解曲线进行比较.结果表明, 在第一个1/4周期时间内, 近似解的解曲线与直接用Mathematica软件作出的数值解解曲线十分吻合, 并且其误差比由文献[2]所得的解曲线更小; 但当时间较长时, 由Adomian分解法得到的解曲线与其他方法得到的解曲线出现了较明显的偏差, 主要原因是由于Adomian分解法要在$n\to \infty
$时其解才会接近数值解, 而本文只取了前面5项.另外, 本文对取3项、4项、5项截断的解曲线进行了分析, 发现5项截断的解曲线比3项、4项截断的解曲线更接近数值解解曲线.
1 Adomian分解法的基本理论
Adomian分解法是求解非线性微分方程的一有效方法[3], 其基本的思路是, 设非线性微分方程可表示成
$
\begin{equation}
Lx+Nx+Rx=g(t),
\end{equation}
$
|
(1) |
且有初始条件
$
\begin{equation}
x_0 =\alpha , \quad \dot {x}_0 =\beta ,
\end{equation}
$
|
(2) |
这里$\alpha , \beta $是实数,
$L$代表方程中最高阶求导的线性可逆算符, $R$是其余线性算符,
$N$是作用在独立变量$x$及其高阶导数上的非线性算符,
$g(t)$是方程的非齐次部分.对于二阶非线性微分方程, $L=\frac{\text
d^2}{\text dt^2}$.
方程(1)可以写成
$
\begin{equation}
\label{eq3} Lx=-Nx-Rx+g(t).
\end{equation}
$
|
(3) |
因$L$是可逆的, 将$L^-$作用于式(3), 得
$
\begin{equation}
\label{eq4} x=\alpha +\beta t-L^-Nx-L^-Rx+L^-g(t).
\end{equation}
$
|
(4) |
令
$
\begin{equation}
\label{eq5} Nx=\sum\limits_{n=0}^\infty {\lambda ^n} a_n ,
\end{equation}
$
|
(5) |
其中$\lambda $是一参数, $a_n $为Adomian多项式, 其计算公式为[4]
$
\begin{equation}
\label{eq6} a_n =\frac{1}{n!}\frac{\partial ^n}{\partial \lambda
^n}\Big[N\Big(\sum\limits_{k=0}^\infty {\lambda ^kx_k \Big)}
\Big]_{\lambda =0} , \quad n=0, 1, 2, 3, \cdots ,
\end{equation}
$
|
(6) |
其中$x_k $是解分量, 且
$
\begin{equation}
\label{eq7} x=\sum\limits_{k=0}^\infty {x_k } .
\end{equation}
$
|
(7) |
将式(5)和式(7)代入式(4), 得各解分量的递推公式
$
\begin{equation}
\label{eq8} \left( {{\begin{array}{*{20}c}
{x_0 } \hfill \\
{x_1 } \hfill \\
{x_2 } \hfill \\
\;\vdots \hfill \\
{x_n } \hfill \\
\end{array} }} \right)=\left( {{\begin{array}{*{20}c}
\quad\quad{\alpha +\beta t+L^-g(t)} \hfill \\
\quad {-L^-Na_0 -L^-Rx_0 } \hfill \\
\quad {-L^-Na_1 -L^-Rx_1 } \hfill \\
\quad\qquad\qquad\vdots \hfill \\
{-L^-Na_{n-1} -L^-Rx_{n-1} } \hfill \\
\end{array} }} \right).
\end{equation}
$
|
(8) |
则非线性微分方程的$n$项近似解为
$
\begin{equation}
\varphi =\sum\limits_{i=0}^{n-1} {x_i } .
\end{equation}
$
|
(9) |
2 二阶非线性微分方程的构建及其近似解析解的求解方法
如图 1所示, 一质量为$m$的光滑小环, 套在一光滑的抛物线形金属丝上, 并可沿着金属丝滑动, 抛物线形金属丝以角速度$\omega $绕轴匀速转动.设抛物线的方程为$x^2=4py$, 则小环在$x$方向的运动微分方程为[1]
$
\begin{equation}
\label{eq10} \ddot {x}+\frac{x\dot {x}^2}{x^2+4p^2}-\frac{4p^2\omega
^2x}{x^2+4p^2}+\frac{2pgx}{x^2+4p^2}=0.
\end{equation}
$
|
(10) |
方程(10)可以写成
$
\begin{equation}
\label{eq11} \ddot {x}=-\frac{x\dot
{x}^2}{x^2+4p^2}-\frac{bx}{x^2+4p^2},
\end{equation}
$
|
(11) |
其中,
$
\begin{equation}
\label{eq12} b=4p^2(\frac{g}{2p}-\omega ^2)=4p^2(\omega _e^2 -\omega
^2),
\end{equation}
$
|
(12) |
$\omega _e =\sqrt {\frac{g}{2p}}
$为小环保持相对平衡时金属丝的转动角速度, 当$\omega <\omega _e $时,
$b>0$; 当$\omega > \omega _e$时, $b<0$; 当$\omega =\omega _e
$时, $b=0$.
方程(11)可以写成
$
\begin{equation}
\label{eq13} Lx=-N_1 x-N_2 x,
\end{equation}
$
|
(13) |
其中
$
N_1 x=\frac{x\dot
{x}^2}{x^2+4p^2}=\sum\limits_{n=0}^\infty {\lambda ^n} a_{1, n} ,
$
|
(14a) |
$
N_2 x=\frac{bx}{x^2+4p^2}=\sum\limits_{n=0}^\infty {\lambda ^n}
a_{2, n} .
$
|
(14b) |
与式(13)相对应的解分量的递推公式为
$
\begin{equation}
\label{eq15} \left( {{\begin{array}{*{20}c}
{x_0 } \hfill \\
{x_1 } \hfill \\
{x_2 } \hfill \\
\;\vdots \hfill \\
{x_n } \hfill \\
\end{array} }} \right)=\left( {{\begin{array}{*{20}c}
\quad\quad\quad\quad\quad{\alpha +\beta t} \hfill \\
\quad\; {-L^-N_1 a_{1, 0} -L^-N_2 a_{2, 0} } \hfill \\
\quad\; {-L^-N_1 a_{1, 1} -L^-N_2 a_{2, 1} } \hfill \\
\qquad\qquad\qquad\vdots \hfill \\
{-L^-N_1 a_{1, n-1} -L^-N_2 a_{2, n-1} } \hfill \\
\end{array} }} \right)\quad ,
\end{equation}
$
|
(15) |
只要算出式(14a)、式(14b)中的Adomian多项式$a_{1, n} , a_{2, n} $, 代入式(15)就可得到非线性微分方程(11)的各解分量, 各解分量求和便得近似解析解.
根据式(6)可计算与$N_1 x=\frac{x\dot {x}^2}{x^2+4p^2}$对应的Adomian多项式$a_{1, n} $.其前5项Adomian多项式为
$
a_{1, 0} =\, \frac{x_0 \dot {x}_0^2 }{x_0^2 +4p^2},
$
|
(16a) |
$
a_{1, 1} =\, \frac{x_1 \dot {x}_0^2 }{x_0^2 +4p^2}+\frac{2x_0 \dot
{x}_0 \dot {x}_1 }{x_0^2 +4p^2}-\frac{2x_0^2 x_1 \dot {x}_0^2
}{(x_0^2 +4p^2)^2},
$
|
(16b) |
$
a_{1, 2} =\, \frac{x_2 \dot {x}_0^2 }{x_0^2 +4p^2}+\frac{2x_1 \dot
{x}_0 \dot {x}_1 }{x_0^2 +4p^2}+\frac{x_0 \dot {x}_1^2 }{x_0^2
+4p^2}+\frac{2x_0 \dot
{x}_0 \dot {x}_2 }{x_0^2 +4p^2} \notag\\
\ \ \ \ \ \ \ \ \ \ -\frac{3x_0 x_1^2 \dot {x}_0^2 }{(x_0^2 +4p^2)^2}-\frac{4x_0^2 x_1
\dot {x}_0 \dot {x}_1 }{(x_0^2 +4p^2)^2}-\frac{2x_0^2 x_2 \dot
{x}_0^2 }{(x_0^2 +4p^2)^2}+\frac{4x_0^3 x_1^2 \dot {x}_0^2 }{(x_0^2
+4p^2)^3} ,
$
|
(16c) |
$
a_{1, 3} =\, \frac{x_3 \dot {x}_0^2 }{x_0^2 +4p^2}+\frac{2x_2 \dot
{x}_0 \dot {x}_1 }{x_0^2 +4p^2}+\frac{x_1 \dot {x}_1^2 }{x_0^2
+4p^2}+\frac{2x_1 \dot {x}_0 \dot {x}_2 }{x_0^2 +4p^2}+\frac{2x_0
\dot {x}_1 \dot {x}_2 }{x_0^2
+4p^2}+\frac{2x_0 \dot {x}_0 \dot {x}_3 }{x_0^2 +4p^2} \notag\\
\ \ \ \ \ \ \ \ \ \ -\frac{x_1^3 \dot {x}_0^2 }{(x_0^2 +4p^2)^2}-\frac{2x_0^2 x_1 \dot {x}_1^2
}{(x_0^2 +4p^2)^2}-\frac{2x_0^2 x_2 \dot {x}_0^2 }{(x_0^2
+4p^2)^2}+\frac{8x_1^3 x_1^2 \dot {x}_0^2 }{(x_0^2
+4p^2)^3}-\frac{8x_0^4 x_1^3 \dot {x}_0^2 }{(x_0^2 +4p^2)^4} ,
$
|
(16d) |
$
a_{1, 4} =\, \frac{x_0 \dot {x}_2^2 }{x_0^2 +4p^2}+\frac{x_2 \dot
{x}_1^2 }{x_0^2 +4p^2}+\frac{2x_3 \dot {x}_0 \dot {x}_1 }{x_0^2
+4p^2}+\frac{2x_1 \dot {x}_1 \dot {x}_2 }{x_0^2 +4p^2}+\frac{2x_0
\dot {x}_1 \dot {x}_3
}{x_0^2 +4p^2}\notag \\
\ \ \ \ \ \ \ \ \ \ +\frac{2x_1 \dot {x}_0 \dot {x}_3 }{x_0^2
+4p^2}+\frac{2x_0 \dot {x}_0 \dot {x}_4 }{x_0^2 +4p^2}+\frac{2x_2
\dot {x}_0 \dot {x}_2 }{x_0^2 +4p^2}+\frac{x_4 \dot {x}_0^2 }{x_0^2
+4p^2}{\kern 1pt}{\kern 1pt}{\kern 1pt}-\frac{2x_1^3 \dot {x}_0 \dot
{x}_1 }{(x_0^2
+4p^2)^2} \notag \\
\ \ \ \ \ \ \ \ \ \ -\frac{3x_0 x_1^2 \dot {x}_1^2 }{(x_0^2 +4p^2)^2}-\frac{2x_0^2 x_2 \dot
{x}_1^2 }{(x_0^2 +4p^2)^2}-\frac{2x_0^2 x_4 \dot {x}_0^2 }{(x_0^2
+4p^2)^2}-\frac{3x_1^2 x_2 \dot {x}_0^2 }{(x_0^2 +4p^2)^2} \notag \\
\ \ \ \ \ \ \ \ \ \ -\frac{3x_0 x_2^2 \dot {x}_0^2 }{(x_0^2 +4p^2)^2}-\frac{12x_0 x_1 x_2 \dot
{x}_0 \dot {x}_1 }{(x_0^2 +4p^2)^2}-\frac{6x_0 x_1 x_3 \dot {x}_0^2
}{(x_0^2
+4p^2)^2}-\frac{4x_0^2 x_1 \dot {x}_1 \dot {x}_2 }{(x_0^2 +4p^2)^2} \notag \\
\ \ \ \ \ \ \ \ \ \ -\frac{4x_0^2 x_1 \dot {x}_0 \dot {x}_3 }{(x_0^2 +4p^2)^2}-\frac{4x_0^2 x_3
\dot {x}_0 \dot {x}_1 }{(x_0^2 +4p^2)^2}-\frac{6x_0 x_1^2 \dot {x}_0
\dot {x}_2 }{(x_0^2 +4p^2)^2}-\frac{4x_0^2 x_2 \dot {x}_0 \dot {x}_2
}{(x_0^2
+4p^2)^2}\notag \\
\ \ \ \ \ \ \ \ \ \ +\frac{4x_0^3 x_1^2 \dot {x}_1^2 }{(x_0^2 +4p^2)^3}+\frac{5x_0 x_1^4
\dot {x}_0^2 }{(x_0^2 +4p^2)^3}+\frac{4x_0^3 x_2^2 \dot {x}_0^2
}{(x_0^2
+4p^2)^3}+\frac{16x_0^3 x_1 x_2 \dot {x}_0 \dot {x}_1 }{(x_0^2 +4p^2)^3} \notag \\
\ \ \ \ \ \ \ \ \ \ +\frac{16x_0^2 x_1^3 \dot {x}_0 \dot {x}_1 }{(x_0^2 +4p^2)^3}+\frac{8x_0^3 x_1^2
\dot {x}_0 \dot {x}_2 }{(x_0^2 +4p^2)^3}+\frac{24x_0^2 x_1^2 x_2
\dot {x}_0^2
}{(x_0^2 +4p^2)^3}+\frac{8x_0^3 x_1 x_3 \dot {x}_0^2 }{(x_0^2 +4p^2)^3}\notag \\
\ \ \ \ \ \ \ \ \ \ -\frac{20x_0^3 x_1^4 \dot {x}_0^2 }{(x_0^2
+4p^2)^4}-\frac{16x_0^4 x_1^3 \dot {x}_0 \dot {x}_1 }{(x_0^2
+4p^2)^4}-\frac{24x_0^4 x_1^2 x_2 \dot {x}_0^2 }{(x_0^2
+4p^2)^4}+\frac{16x_0^5 x_1^4 \dot {x}_0^2 }{(x_0^2 +4p^2)^5}.
$
|
(16e) |
与$N_2 x=\frac{bx}{x^2+4p^2}$对应的前5项Adomian多项式为
$
a_{2, 0} =\frac{bx_0 }{x_0^2 +4p^2},
$
|
(17a) |
$
a_{2, 1} =\frac{bx_1 }{x_0^2 +4p^2}-\frac{2bx_0^2 x_1 }{(x_0^2
+4p^2)^2},
$
|
(17b) |
$
a_{2, 2} =\frac{bx_2 }{x_0^2 +4p^2}-\frac{3bx_0 x_1^2 }{(x_0^2
+4p^2)^2}-\frac{2bx_0^2 x_2 }{(x_0^2 +4p^2)^2}+\frac{4bx_0^3 x_1^2
}{(x_0^2 +4p^2)^3},
$
|
(17c) |
$a_{2, 3} =\frac{bx_3 }{x_0^2 +4p^2}-\frac{6bx_0 x_1 x_2 }{(x_0^2
+4p^2)^2}-\frac{bx_1^3 }{(x_0^2 +4p^2)^2}-\frac{2bx_0^2 x_3 }{(x_0^2
+4p^2)^2} \notag
\\
\qquad\;\;\;+\frac{8bx_0^2 x_1^3 }{(x_0^2 +4p^2)^3}+\frac{8bx_0^3
x_1 x_2 }{(x_0^2 +4p^2)^3}-\frac{8bx_0^4 x_1^3
}{(x_0^2 +4p^2)^4}
$
|
(17d) |
$
a_{2, 4} =\frac{bx_4 }{x_0^2 +4p^2}-\frac{6bx_0 x_1 x_3 }{(x_0^2
+4p^2)^2}-\frac{3bx_1^2 x_2 }{(x_0^2 +4p^2)^2}-\frac{3bx_0 x_2^2
}{(x_0^2
+4p^2)^2} \notag\\
\qquad\;\;\;-\frac{2bx_0^2 x_4 }{(x_0^2 +4p^2)^2}+\frac{24bx_0^2
x_1^2 x_2 }{(x_0^2 +4p^2)^3}+\frac{5bx_0 x_1^4
}{(x_0^2 +4p^2)^3}+\frac{4bx_0^3 x_2^2 }{(x_0^2 +4p^2)^3}\notag \\
\qquad\;\;\;+\frac{8bx_0^3 x_1 x_3 }{(x_0^2
+4p^2)^3}-\frac{20bx_0^3 x_1^4 }{(x_0^2 +4p^2)^4}-\frac{24bx_0^4
x_1^2 x_2 }{(x_0^2 +4p^2)^4}+\frac{16bx_0^5 x_1^4 }{(x_0^2 +4p^2)^5}
$
|
(17e) |
3 近似解析解与结果比较
将式(16), 式(17)代入式(15)能得到各解分量, 求和后可得非线性微分方程近似解析解的一般表达式.但由于式(16), 式(17)过于复杂且与初始条件$x_0 , \dot {x}_0 $和参数$p, b$有关, 难以从所得解析解的一般表达式中讨论其准确性.因此, 我们给定一组初始条件$x_0 , \dot {x}_0 $和参数$p, b$, 得到解析解的表达式并作出其解曲线, 并将解曲线与直接用Mathematica软件作出的解曲线以及与同伦渐近法[2]所得解曲线进行比较, 来说明结果的准确性.
设初始条件和参数为$x_0 =1$、$\dot {x}_0 =0$、$p=b=0. 25$, 将它们代入式(16)和式(17)可得$a_{1, n} , a_{2, n} (n=0, 1, 2, 3, 4)$, 具体为
$
a_{1, 2} =\frac{4\dot {x}_1^2 }{5},
$
|
(18c) |
$
a_{1, 3} =-\frac{12x_1 \dot {x}_1^2 }{25}+\frac{8}{5}\dot {x}_1
\dot {x}_2 \quad ,
$
|
(18d) |
$
a_{1, 4} =-\frac{12}{25}x_2 \dot {x}_1^2 -\frac{24}{25}x_1 \dot
{x}_1 \dot {x}_2 +\frac{4}{5}\dot {x}_2^2 +\frac{8}{5}\dot {x}_1
\dot {x}_3 +\frac{16}{125}x_1^2 \dot {x}_1^2 ,
$
|
(18e) |
$
a_{2, 0} =\frac{1}{5},
$
|
(18f) |
$
a_{2, 1} =-\frac{3x_1 }{25},
$
|
(18g) |
$
a_{2, 2} =\frac{4}{125}x_1^2 -\frac{3}{25}x_2 ,
$
|
(18h) |
$
a_{2, 3} =\frac{28}{625}x_1^3 +\frac{8}{125}x_1 x_2
-\frac{3}{25}x_3 ,
$
|
(18i) |
$
a_{2, 4} ={\kern 1pt}{\kern 1pt}-\frac{3x_4 }{25}+\frac{8x_1 x_3
}{125}+\frac{4x_2^2 }{125}+\frac{84x_1^2 x_2 }{625}-\frac{304x_1^4
}{3125}.
$
|
(18j) |
将初始条件和式(18)依次代入式(15), 可得如下解分量, 为
$
x_1 =-\frac{1}{10}t^2,
$
|
(19b) |
$
x_2 =-\frac{1}{1 000}t^4,
$
|
(19c) |
$
x_3 =-\frac{1}{375}t^4-\frac{11}{750 000}t^6,
$
|
(19d) |
$
x_4 =-\frac{11}{93 750}t^6+\frac{229}{350}\times 10^{-6}t^8,
$
|
(19e) |
$
x_5 =-\frac{16}{140 625}t^6-\frac{257}{656 25 000}t^8+\frac{289
397}{236 250}\times 10^{-7}t^{10}.
$
|
(19f) |
则在3项、4项、5项截断的近似解析解分别为
$
\varphi ^3=\sum\limits_{i=0}^3 {x_i }
=1-\frac{1}{10}t^2-\frac{11}{3 000}t^4-\frac{11}{750 000}t^6,
$
|
(20) |
$
\varphi ^4=\sum\limits_{i=0}^4 {x_i }
=1-\frac{1}{10}t^2-\frac{11}{3 000}t^4-\frac{33}{250
000}t^6+\frac{229}{350}\times 10^{-6}t^8,
$
|
(21) |
$
\begin{align}
&\varphi ^5=\sum\limits_{i=0}^5 {x_i }
=1-\frac{1}{10}t^2-\frac{11}{3 000}t^4-\frac{553}{225}\times 10^{-4}t^6 \notag\\
&\qquad\qquad\qquad-\frac{137}{42}\times 10^{-6}t^8+\frac{289
397}{236 250}\times 10^{-7}t^{10}.
\end{align}
$
|
(22) |
式(22)中第二个等号后的第一项"1"表示初位置, 第二项"$-\frac{1}{10}t^2$"反映初加速度的大小, 且初加速度大小由Adomian多项式$a_{1, 0}, a_{2, 0} $之和决定, 即由初始条件$x_0, \dot {x}_0 $和微分方程的参数$p, b$确定.
在相同的初始条件和参数下, 文献[2]中式(58)表示的近似解为
$
\begin{equation}
\varphi =1.092~937~297\cos (\Omega t)-0.123~160~203\cos (3\Omega
t)+0.030~222~906\cos (5\Omega t),
\end{equation}
$
|
(23) |
其中$\Omega =0.596~211~722$.
如图 2所示, 红色曲线是根据式(22)作出的解曲线, 蓝色曲线是根据式(23)作出的解曲线, 黑色曲线是用Mathematica软件作出的解曲线.从图 2中可以看出, 在第一个1/4周期内, 由近似解式(22)作出的解曲线与直接用Mathematica软件作出的解曲线十分吻合, 且其误差比用同伦渐近法得到的解曲线还小.但是在时间较大区域, 近似解式(22)作出的解曲线会出现较大的误差.产生误差的主要原因是, 用Adomian分解法解非线性微分方程, 当$n\to \infty
$时才达到较高的近似度.因文中方程过于复杂, 本文只计算到前5项, 从而在时间较长时出现较大误差.
如图 3所示, 蓝色曲线、绿色曲线和红色曲线分别表示在3项、4项、5项截断时的近似解析解曲线, 黑色曲线是用Mathematica软件作出的解曲线.从图 3可见, 在第一个1/4周期时间内, 5项截断比3项、4项截断更接近数值解.
另外, Adomian多项式(16)、(17)是由初始条件$x_0 , \dot {x}_0
$、非线性微分方程的参数$p, b$及解分量$x_i^n , \dot {x}_i^m (i, n,
m=1, 2, \cdots, 4)$表示的多项式, 其中的初始条件$x_0 , \dot {x}_0
$和参数$p, b$可以任意选取.由递推公式(15)得到的解分量具有普遍性, 任取一组初始条件$x_0 , \dot {x}_0 $和参数$p, b$, 就可得到一组相应的解析解表达式及其解曲线, 并与直接用Mathematica软件作出的解曲线进行比较, 仍可得出在第一个1/4周期时间内两者十分吻合的结论.
设初始条件和参数变为$x_0 =1, \dot {x}_0 =0, p=0. 25, b=1$, 按上述相同方法可得其在5项截断的近似解析解
$
\begin{align}
\varphi ^5=\sum\limits_{i=0}^5 {x_i }
=\, &1-\frac{2}{5}t^2-\frac{22}{375}t^4-\frac{6~558}{140~625}t^6 \notag\\
&-\frac{7~541}{8~203~125}t^8+\frac{1~145~066}{9~228~515~625}t^{10}.
\end{align}
$
|
(24) |
如图 4所示.红色曲线是根据式(24)作出的解曲线, 黑色曲线是用Mathematica软件作出的解曲线.从图 4中可以看出, 在第一个1/4周期内, 由近似解式(24)作出的解曲线与用Mathematica软件作出的解曲线比较吻合, 说明在不同的参数下, 用Adomian分解法得到的近似解的准确性具有普遍性, 与参数的选取无关.
4 结论
研究力学系统在非惯性转动参照系中运动时常常会遇到分子分母中均含非线性项的二阶微分方程, 因方程复杂而无法直接求得其解析解, 从而无法了解力学体系的运动规律.本文用Adomian分解法得到了在给定初始条件和参数下的3项、4项和5项截断的近似解析解, 并通过作图与直接用Mathematica软件作出的解曲线以及与同伦渐近法所得解曲线进行了比较.结果表明, 在第一个1/4周期时间内(如图 2所示), 5项截断的近似解的解曲线与直接用Mathematica软件作出的解曲线十分吻合, 并且其误差比文献[2]用同伦渐近法得到的解曲线还小.由于Adomian分解法要在$n\to \infty $时其解才会接近数值解, 因此在只计算有限项近似解, 且在时间取值较大时会存在较大的误差, 这也是Adomian分解法的不足之处.另外, 对于数值解具有周期性的运动, 用多项式表示近似解存在收敛性的问题, 这一问题我们将另文研究.