解常微分方程的多步法(有Adams外插法,Adams内插法,一般多步法(Milne法))是数值求解常微分方程初值问题的重要方法。
这类方法通过利用前几步的计算结果来构造高阶差商,从而获得更高的精度和效率。
Adams外插法(如Adams-Bashforth方法)是一类显式多步法,其公式基于函数在前几个节点处的导数近似,计算简单且易于实现,但稳定性相对较弱。Adams内插法(如Adams-Moulton方法)则是隐式多步法,通常作为Adams外插法的校正步骤,通过迭代求解来提高精度和稳定性。
一般多步法中的Milne法是一种典型的显式四阶方法,它通过对积分区间进行分段线性或更高次插值,将微分方程转化为代数方程组求解。
这些方法在实际应用中,常结合预测-校正策略使用,即先用外插法(如Milne法)进行预测,再用内插法(如Adams-Moulton法)进行校正,以平衡计算效率与数值稳定性。
选择合适的多步法需根据具体问题的特性(如刚性、精度要求等)以及计算资源进行综合考量。

admas外插法
%admas外插法,具有k-1阶的精度,下面编写k等于3的外插公式
clear
clc
syms x y
f=3*x^2;
a=0;
b=10;
y0=0;
y1=0.001;
y2=0.008;
y3=0.027;
h=0.1;
n=(b-a)/h+1;
xx=zeros(n,1);
yy=zeros(n,1);
yy=zeros(n,1);
yy(1)=y0;
yy(2)=y1;
yy(3)=y2;
yy(4)=y3;
for i=1:n
xx(i)=a+h*(i-1);
end
for i=5:n
k1=55/24*h*(subs(f,{
x,y},{
xx(i-1),yy(i-1)}));
k2=-
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删