[MATLAB基础] ode45的使用 - MATLAB爱好者论坛-LabFans.com

剑八 · 2009-04-13, 17:39

请问如何使用ode45解
d2x1/dt2=x2-x1;
d2x2/dt2=30+x1-2*x2;
x1(0)=12,x2(0)=23,Dx1(0)=0,Dx2(0)=0

hrt · 2009-04-13, 20:59

首先建立函数M文件如下：
function dx=fmat(t,x) % x=[x1,x2,x3,x4],其中x3=dx1、x4=dx2；
dx=[x(3);x(4);x(2)-x(1);30+x(1)-2*x(2)];
end

再在命令行中输入：
>> t=0:100;
>> x0=[12,23,0,0];
>> [t,x]=ode45(@fmat,t,x0);
>> [t,x]

可得到结果。
~ ~

2009-04-13, 17:39	#1
剑八初级会员注册日期: 2008-11-24 住址: 福州帖子: 11 声望力: 17	ode45的使用请问如何使用ode45解 d2x1/dt2=x2-x1; d2x2/dt2=30+x1-2*x2; x1(0)=12,x2(0)=23,Dx1(0)=0,Dx2(0)=0

2009-04-13, 20:59	#2
hrt 初级会员注册日期: 2009-03-11 年龄: 38 帖子: 9 声望力: 0	回复: ode45的使用首先建立函数M文件如下： function dx=fmat(t,x) % x=[x1,x2,x3,x4],其中x3=dx1、x4=dx2； dx=[x(3);x(4);x(2)-x(1);30+x(1)-2*x(2)]; end 再在命令行中输入： >> t=0:100; >> x0=[12,23,0,0]; >> [t,x]=ode45(@fmat,t,x0); >> [t,x] 可得到结果。 ~ ~