Mechanics 版 (精华区)

发信人: jackleo (黑马), 信区: Mechanics
标  题: [合集] 请教ODE45
发信站: 哈工大紫丁香 (Fri Jun  9 08:45:55 2006), 站内

────────────────────────────────────────
  muscleman (疏星淡月)             于  (Wed May 17 16:54:19 2006)  说道:

怎样求解如下方程组：
x1''+a*x1-b*x2=0;
x2''-c*x1+c*x2=0;
谢谢


────────────────────────────────────────
  feifeifool (学习、学习、再学习)  于  (Wed May 17 17:19:01 2006)  说道:

你为什么不到matlab版呢,不过在这个版前一阵子有过一个帖子.
【 在 muscleman (疏星淡月) 的大作中提到: 】
: 怎样求解如下方程组：
: x1''+a*x1-b*x2=0;
: x2''-c*x1+c*x2=0;
: ...................



────────────────────────────────────────
  feifeifool (学习、学习、再学习)  于  (Wed May 17 17:24:07 2006)  说道:

给个例子吧:
function Yd=DYdt2(t,Y)%其中的Y=[x1;x2;dx1;dx2]的列向量
Yd(1)=Y(3);
Yd(2)=Y(4);
Yd(3)=0.4909*Y(2)-1.5708*Y(1)+1.3695*Y(4)-72.73*Y(3);
Yd(4)=-0.1964*Y(2)+0.4909*Y(1)-64.157*Y(4)+1.369*Y(3);
Yd=Yd';
参数和命令：>> x1=0;x2=0;dx1=0;dx2=0;
>> Y0=[x1;x2;dx1;dx2];
>> [t,Y]=ode45('DYdt2',[0,10],Y0)
欢迎到matlab版来,呵呵,:)

【 在 feifeifool (学习、学习、再学习) 的大作中提到: 】
: 你为什么不到matlab版呢,不过在这个版前一阵子有过一个帖子.




────────────────────────────────────────
  zjliu (秋天的萝卜)               于  (Wed May 17 17:37:25 2006)  说道:

如下是你的方程的解:

a=1;b=2;c=3;
Df=inline('[b*x(2)-a*x(1);c*(x(1)-x(2))]',...
    't','x','flag','a','b','c');
[t,x]=ode45(Df,[0,10],[1,0],[],a,b,c);
plot(t,x(:,1),'r');hold on;
plot(t,x(:,2),'b');
legend('x1','x2');

【 在 muscleman (疏星淡月) 的大作中提到: 】
: 怎样求解如下方程组：
: x1''+a*x1-b*x2=0;
: x2''-c*x1+c*x2=0;
: ...................



────────────────────────────────────────
  jackleo (黑马)                   于  (Wed May 17 20:05:41 2006)  说道:

呵呵，两个版主，两种方法，挺全的。
【 在 zjliu (秋天的萝卜) 的大作中提到: 】
: 如下是你的方程的解:
: a=1;b=2;c=3;
: Df=inline('[b*x(2)-a*x(1);c*(x(1)-x(2))]',...
: ...................



────────────────────────────────────────
  zjliu (秋天的萝卜)               于  (Wed May 17 20:32:18 2006)  说道:

我的还不完全正确,那里是二阶导数,变为4个
一阶方程即可
【 在 jackleo (黑马) 的大作中提到: 】
: 标  题: Re: 请教ODE45
: 发信站: 哈工大紫丁香 (Wed May 17 20:05:41 2006), 转信
: 
: 呵呵，两个版主，两种方法，挺全的。
: 【 在 zjliu (秋天的萝卜) 的大作中提到: 】
: : 如下是你的方程的解:
: : a=1;b=2;c=3;
: : Df=inline('[b*x(2)-a*x(1);c*(x(1)-x(2))]',...
: : ...................
: 
: --
: 欢迎光临力学版(Mechanics)!
: 大风起兮云飞扬
: 威加海内兮归故乡
: 安得猛士兮守四方
: http://bbs.skylook.org/uploads/post-11-1111662913.jpg
: 
: 
: ※ 来源:·哈工大紫丁香 bbs.hit.edu.cn·[FROM: 202.118.247.36]




────────────────────────────────────────
  feifeifool (学习、学习、再学习)  于  (Wed May 17 21:08:44 2006)  说道:

我直接把以前人家问的问题复制粘贴过来的,省事.又能解决问题,
传说中的拿来主义,:)
【 在 zjliu (秋天的萝卜) 的大作中提到: 】
: 我的还不完全正确,那里是二阶导数,变为4个
: 一阶方程即可




────────────────────────────────────────
  zjliu (秋天的萝卜)               于  (Wed May 17 21:09:51 2006)  说道:

能解决问题就行，呵呵
【 在 feifeifool (学习、学习、再学习) 的大作中提到: 】
: 我直接把以前人家问的问题复制粘贴过来的,省事.又能解决问题,
: 传说中的拿来主义,:)