zjadms.m

基于matlab开发软件 的选择输入类型的数值计算仿真程序 迭代方法包括 1可自启动的4阶龙格库塔法 2需要使用其他方法启动的阿达姆斯法

function [x2,y2]=zjadms(a,b,c,d,u,h,n,x0)
x=x0;
y=c*x;
for i=1:n
    if (i==1)
        x(:,i+1)=x(:,i)+h*(a*x(:,i)+b*u);
    else
        x(:,i+1)=x(:,i)+h/2*(3*(a*x(:,i)+b*u)+(a*x(:,i-1)+b*u));
    end
    y(i+1)=c*x(:,i+1);
end
x2=x;
y2=y;