zjadms.m

来自「基于matlab开发软件 的选择输入类型的数值计算仿真程序 迭代方法包括 」· M 代码 · 共 13 行

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;