| 查看: 2906 | 回复: 21 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题 已有7人参与
|
|||
|
向各位大侠求助matlab求解微分方程组遇到的一个问题： 假设有以下的微分方程组 u'(t)=-a*f1(t)*y(t); v'(t)=-b*v(t)+a*f1(t)*y(t); w'(t)=b*v(t); x'(t)=-a*f2(t)*v(t)+a*f1(t)*y(t)+a*f3(t)*y(t)+c*y(t); y'(t)= a*f2(t)*v(t) - a*f1(t)*y(t) - a*f3(t)*y(t) - c*y(t); 其中 f1(t)=0 if u(t)=0 =1 if u(t)>0; f2(t)=0 if x(t)=0 =1 if x(t)>0; f3(t)=0 if v(t)=0 =1 if v(t)>0; a=1;b=0.1;c=0.5; 且u(0)=1000;v(0)=1;w(0)=0;x(0)=200;y(0)=0.现在要得到各函数随t的变化曲线，如果微分方程组中没有f1(t)，f2(t)，f2(t)这三个方程，那么在matlab中可以直接用ode来得到，但是现在有这三个方程，我就没办法了，试了好几种方法都没用,想请教各位高手，在matlab下怎么才能求解这个方程？ [ Last edited by ybnav on 2011-3-11 at 22:27 ] |
回复此楼» 猜你喜欢
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
要不要辞职读博？
已经有3人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
自荐读博
已经有3人回复
-
不自信的我
已经有5人回复
|
我在matlab里是这样做的： 如果不考虑f1(t)、f2(t)、f3(t)，令它们的值都等于1 函数 function dydt=xyz(t,y) dydt=[ -y(5) -0.1*y(2)+y(5) 0.1*y(2) -y(2)+y(5)+y(5)+0.5*y(5) y(2)-y(5)-y(5)-0.5*y(5) ]; 脚本 clear all close all clc tspan = [0 60]; y0 = [1000;1;0;200;0]; [t,Y]=ode45(@xyz,tspan,y0); plot(t,Y(:,1),'r-','LineWidth',1.5); hold on; plot(t,Y(:,2),'g-','LineWidth',1.5); hold on; plot(t,Y(:,3),'b.','LineWidth',1.5); hold on; plot(t,Y(:,4),'g.','LineWidth',1.5); hold on; plot(t,Y(:,5),'r-.','LineWidth',1.5); axis([0 60 0 1000]) legend('Y(1)','Y(2)','Y(3)','Y(4)','Y(5)') grid on 这是没有加上f1(t)、f2(t)、f3(t)这三个的情况，如果加上这三个函数，应该怎么改呢？ |
13楼2011-03-11 22:37:59
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小木虫(金币+0.5):给个红包，谢谢回帖交流
|
把你的函数贴出来，我稍改改就可以求解了。 或者参考：http://muchong.com/bbs/viewthread.php?tid=2764636### 把微分方程写成分段函数的形式即可。 |
2楼2011-03-11 15:42:54
|
就是这个微分方程组 u'(t)=-*f1(t)*y(t); v'(t)=-0.1*v(t)+*f1(t)*y(t); w'(t)=0.1*v(t); x'(t)=-f2(t)*v(t)+f1(t)*y(t)+f3(t)*y(t)+0.5*y(t); y'(t)= f2(t)*v(t) -f1(t)*y(t) -f3(t)*y(t) -0.5*y(t); 其中 f1(t)=0 if u(t)=0 =1 if u(t)>0; f2(t)=0 if x(t)=0 =1 if x(t)>0; f3(t)=0 if v(t)=0 =1 if v(t)>0; 且u(0)=1000;v(0)=1;w(0)=0;x(0)=200;y(0)=0. 只要求在matlab中画出u(t)、v(t)、w(t)、x(t)、y(t)在t的某一段区间内，比如[0,60]的曲线图，还是没想出来改怎么办，同学说可能需要迭代，我也不知道怎么办，万分感谢版主同志！ [ Last edited by ybnav on 2011-3-11 at 22:28 ] |
3楼2011-03-11 17:31:40
|
刚方程写错了，应该是 就是这个微分方程组 u'(t)=-f1(t)*y(t); v'(t)=-0.1*v(t)+f1(t)*y(t); w'(t)=0.1*v(t); x'(t)=-f2(t)*v(t)+f1(t)*y(t)+f3(t)*y(t)+0.5*y(t); y'(t)= f2(t)*v(t) -f1(t)*y(t) -f3(t)*y(t) -0.5*y(t); 其中 f1(t)=0 if u(t)=0 =1 if u(t)>0; f2(t)=0 if x(t)=0 =1 if x(t)>0; f3(t)=0 if v(t)=0 =1 if v(t)>0; 且u(0)=1000;v(0)=1;w(0)=0;x(0)=200;y(0)=0. 只要求在matlab中画出u(t)、v(t)、w(t)、x(t)、y(t)在t的某一段区间内，比如[0,60]的曲线图，还是没想出来改怎么办，同学说可能需要迭代，我也不知道怎么办，万分感谢版主同志！ [ Last edited by ybnav on 2011-3-11 at 22:31 ] |
4楼2011-03-11 17:33:08