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

来源:学生作业帮助网 编辑:作业帮 时间:2024/05/01 00:05:37
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

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
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);
报错.请问是怎么回事?

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

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 求大神支招,我想用ode45来求解一个微分方程的数值解,但是MATLAB老是说有错,怎么办?程序是clear all;options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]);[T,Y]=ode45(@guanxi,[0 12],[0 1 1],options);对应 guanxi.m 是funct 用MATLAB求微分方程分别用Euler法和ode45解下列常微分方程并与解析解比较 matlab解微分方程组求微分方程组的数值解并作图.要求是用ode45的方法来解 matlab ode45 可以解二阶微分方程吗 Matlab 已知y''-3y'+2y=0 y(0)=1 y'(0)=0 ,求y(0.5) 会的帮忙写下程序 用ode45实现最好 我在编matlab程序的过程中,需要用ode45解一个一阶微分方程组,我在编matlab程序的过程中,需要用ode45解一个一阶微分方程组,关键问题是这个方程组的系数需要前面的程序传进去,所以我不会用, 大神求解动力学微分方程,用matlab中ode45编程实现!其中的参数可以随便设置, 用MATLAB如何用ODE45求解如下微分方程已知J,Tb,Fx和R 用matlab求解常微分这个程序怎么编的啊?我想编一个matlab程序,使用ode45求解:dy/dx=(y^2+5)*(x+6).x范围是0到5,y0=0.5 matlab解微分方程组ode45用ode45 解微分方程组 以前基本都没碰过matlab,自己试着编了一点儿,搞不定啊, 只剩这么多分了,好人一生平安,雅安加油!@liuliangsxd 您好,下面的是我根据您给的程序,完 matlab ode45 求方程(1+x^2)y''=2xy',y(0)=1,y'(0)=3 的数值解,我的运行不行,我编的函数文件:function dy=sdre(x,y) dy=[y(2);2*x/(1+x^2)*y(2)]; 程序 [x,y]=ode45('sdre',[0 10],[1 3]);plot(x,y) matlab用ode45 解微分方程如何用ODE45解微分方程 其中初值给出方程如下,请帮忙写出方程在程序中如何表达,我用function dy=func(x,y,r)dy = zeros(size(y));dy(1) = y(1)dy(2) = sqrt((1+r^2)/r^2*log((1+r^2)/(1+y(1)^2))+y(1 matlab中ode45,4和5分别代表什么? 急求MATLAB解常微分方程组数值解的解决方案我用MATLAB7.4的ODE45解四元常微分方程组,程序如下:function xdot = eqx(t,x);m = 0.5;alf = 0.15;K=[0.0265,0.1,0.15];Q=20;a=0.8;xdot = zeros(4,1);xdot(1) = Q - alf * x(1) * x(2) + matlab ode45 求解二阶常微分方程新手想用ode45求解上述二阶三元常微分方程组,几次尝试均未成功,请大神帮忙编写matlab代码,最好有图, 用matlab中ode45解方程组 x1'=-2*x1+x2 x2'=-3*x1+6*x2 初始值自定 用matlab绘制等量异种电荷的的三维电场线,初始条件怎么确定使用ode45