| 查看: 1179 | 回复: 10 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
monitor2885至尊木虫 (职业作家)
队长
|
[交流]
【求助】简单的ode45程序修改
|
||
|
function abc clear all;clc global y u u=[1;2;3;4;5]; t=0:4; [t,y]=ode45(@fun,t,[0 2]); plot(y) function dydt=fun(t,y) global u v=u*2; dydt=zeros(2,1); dydt(1)=y(2); dydt(2)=v*y(1)-2*v*y(2); 哪里错了呢？怎么改？谢谢 |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有27人回复
-
所感
已经有3人回复
-
要不要辞职读博？
已经有7人回复
-
不自信的我
已经有11人回复
-
北核录用
已经有3人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
7楼2010-03-18 16:14:52
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
★
monitor2885(金币+5):帮忙看下面一个 2010-03-17 21:33
kuhailangyu(金币+1):欢迎积极参与 2010-03-18 16:24
monitor2885(金币+5):帮忙看下面一个 2010-03-17 21:33
kuhailangyu(金币+1):欢迎积极参与 2010-03-18 16:24
|
function abc clear all clc u=[1;2;3;4;5]; tspan=0:1:4; tt=[];yy=[]; y0=[0 2]; for i=1:length(tspan)-1 [t,y]=ode45(@fun,[tspan(i),tspan(i+1)],y0,[],u(i)); tt=[tt;t];yy=[yy;y]; y0=y(end,: ); end figure(1) plot(tt,yy) function dydt=fun(t,y,u) v=u*2; dydt=[y(2); v*y(1)-2*v*y(2)]; |
2楼2010-03-17 19:45:45
monitor2885
至尊木虫 (职业作家)
队长
- 应助: 6 (幼儿园)
- 贵宾: 0.339
- 金币: 38693.5
- 散金: 147
- 红花: 22
- 帖子: 4942
- 在线: 2166小时
- 虫号: 701577
- 注册: 2009-02-15
- 性别: GG
- 专业: 生态系统生态学
|
比如，现在2个参数，即u和v，下面哪里不对呢？ clear all;clc global yy u=[1;5;15;20;25]; v=[6;12;18;24;30]; tspan=0:1:4; tt=[];yy=[]; y0=[0 2]; for i=1:length(tspan)-1 [t,y]=ode45(@fun,[tspan(i),tspan(i+1)],y0,[],(u(i),v(i))); tt=[tt;t]; yy=[yy;y]; y0=y(end,  ;end figure;plot(tt,yy) function dydt=fun(t,y,u,v) dydt=zeros(2,1); dydt(1)=u+y(2); dydt(2)=v*y(1)-2*v*y(2); [ Last edited by monitor2885 on 2010-3-17 at 21:34 ] |
3楼2010-03-17 21:33:30
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
4楼2010-03-17 23:04:14