matlab程序ode45

来源：学生作业帮编辑：作业帮分类：综合作业时间：2024/05/26 13:26:06

matlab程序ode45
用ode45求解微分方程组dx=20/7*x-yz,dy=-10y+xz,dz=-4z+xy,如何编写,还要画图.初值(3,-4,2).对初值很敏感,t=0.001:0.001:20.
function dx=myfun(t,x)
dx(1)=20/7*x(1)-x(2)*x(3);
dx(2)=-10*x(2)+x(1)*x(3);
dx(3)=-4*x(3)+x(1)*x(2);
dx=dx(:);
x0=[3,-4,2];
t0=0.001:0.001:20;
[t,x]=ode45('myfun',[0.001,20],x0); %ode45会自动调整步长
plot(t,x)
legend('x','y','z')
程序中运行出来
Input argument "x" is undefined.
Error in ==> myfun at 2
dx(1)=20/7*x(1)-x(2)*x(3);
报错.请问是怎么回事?

function dx=myfun(t,x)
dx(1)=20/7*x(1)-x(2)*x(3);
dx(2)=-10*x(2)+x(1)*x(3);
dx(3)=-4*x(3)+x(1)*x(2);
dx=dx(:);这部分保存为m函数文件
命令行运行>> x0=[3,-4,2];t0=0.001:0.001:20;[t,x]=ode45('myfun',[0.001,20],x0); %ode45会自动调整步长plot(t,x) legend('x','y','z')>> 结果

matlab ode45 用matlab的ode45求解如下微分方程的程序 matlab ode45 解方程 MATLAB 中的ODE45 matlab函数ode45出错,如果能调好整个完整程序,愿意加100分我在编matlab程序的过程中,需要用ode45解一个一阶微分方程组, matlab解微分方程组ode45 matlab ode45 求解二阶常微分方程 matlab用ode45 解微分方程 matlab ode45解微分方程组 matlab ode45求解微分方程组 matlab ode45 求解常微分方程组