Matlab 版 (精华区)

发信人: hahn (有奇~~Not Over But another Start...), 信区: Matlab
标  题: [合集] 紧急求助！！
发信站: 哈工大紫丁香 (Fri Jan 19 12:38:58 2007), 站内

────────────────────────────────────────
  wofei (ZZ)                       于  (Thu Jun 22 22:32:42 2006)  说道:

请高手帮助解一下方程,据说mapple和matlab 可以解
Dx1=x2
Dx2=(1.5*x1^2)/x4
Dx3=-4*x1-0.5*x2
Dx4=-0.5*x1-x2-x3
求出近似解即可，论文急需！！！！先谢谢好心的兄弟姐妹啦
为了论文我已经两天没有好好睡觉了



────────────────────────────────────────
  zjliu (秋天的萝卜)               于  (Thu Jun 22 22:45:11 2006)  说道:

ode45 in matlab
【 在 wofei (ZZ) 的大作中提到: 】
: 请高手帮助解一下方程,据说mapple和matlab 可以解
: Dx1=x2
: Dx2=(1.5*x1^2)/x4
: ...................



────────────────────────────────────────
  feifeifool (爱芳)                于  (Fri Jun 23 07:32:47 2006)  说道:

function xd=qiujie(t,x)
xd=[x(2);1.5*x(1)^2/x(4);-4*x(1)-0.5*x(2);-0.5*x(1)-x(2)-x(3)];
%%%%%%%%%%%%%%%%%%%
clc;clear;
tspan=[1,10];
x0=[1;1;1;10];
[t,x]=ode45('qiujie',tspan,x0);
【 在 zjliu (秋天的萝卜) 的大作中提到: 】
: ode45 in matlab




────────────────────────────────────────
  wofei (ZZ)                       于  (Fri Jun 23 09:59:46 2006)  说道:

还是不会啊，有没有好心人啊，快急死了，详细一点的

【 在 zjliu (秋天的萝卜) 的大作中提到: 】
: ode45 in matlab





────────────────────────────────────────
  feifeifool (爱芳)                于  (Fri Jun 23 10:22:01 2006)  说道:

不是给你写程序了么？自己看。
【 在 wofei (ZZ) 的大作中提到: 】
: 还是不会啊，有没有好心人啊，快急死了，详细一点的




────────────────────────────────────────
  wofei (ZZ)                       于  (Fri Jun 23 12:28:33 2006)  说道:

    谢谢老大，帮了我大忙啦！！！感激!!!!
【 在 feifeifool (爱芳) 的大作中提到: 】
: function xd=qiujie(t,x)
: xd=[x(2);1.5*x(1)^2/x(4);-4*x(1)-0.5*x(2);-0.5*x(1)-x(2)-x(3)];
: %%%%%%%%%%%%%%%%%%%
: clc;clear;
: tspan=[1,10];
: x0=[1;1;1;10];
: [t,x]=ode45('qiujie',tspan,x0);





────────────────────────────────────────
  wofei (ZZ)                       于  (Fri Jun 23 13:18:38 2006)  说道:

还想问一下，我的是末端约束，T=无穷大，x3（T）=0，x4（T）=0
怎么在程序中体现出来。
再次感谢！！！！
【 在 wofei (ZZ) 的大作中提到: 】
:     谢谢老大，帮了我大忙啦！！！感激!!!!





────────────────────────────────────────
  feifeifool (爱芳)                于  (Fri Jun 23 13:24:29 2006)  说道:

这个是边值问题,用ode45时解决不了的,无穷大用1/eps可以近似的代替吧.
可以考虑打靶法解决.你还可以看看bvp4c这个函数试一试能不能直接用，没用
过这个函数
【 在 wofei (ZZ) 的大作中提到: 】
: 还想问一下，我的是末端约束，T=无穷大，x3（T）=0，x4（T）=0
: 怎么在程序中体现出来。
: 再次感谢！！！！
: ...................



────────────────────────────────────────
  wofei (ZZ)                       于  (Fri Jun 23 15:25:14 2006)  说道:

    好，我试试　，你觉得maple能不能解这个非线性微分方程组？
【 在 feifeifool (爱芳) 的大作中提到: 】
: 这个是边值问题,用ode45时解决不了的,无穷大用1/eps可以近似的代替吧.
: 可以考虑打靶法解决.你还可以看看bvp4c这个函数试一试能不能直接用，没用
: 过这个函数





────────────────────────────────────────
  feifeifool (爱芳)                于  (Fri Jun 23 21:00:55 2006)  说道:

应该可以，我没用过。
【 在 wofei (ZZ) 的大作中提到: 】
:     好，我试试　，你觉得maple能不能解这个非线性微分方程组？