| 查看: 1252 | 回复: 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); 哪里错了呢？怎么改？谢谢 |
回复此楼» 猜你喜欢
一个化合物的合成路线：CAS:367929-02-0 名称：8β-乙烯基雌二醇
已经有3人回复
-
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业
已经有3人回复
-
河北省自然基金
已经有8人回复
-
收到国自然专家邀请后几年才会有本子送过来评
已经有3人回复
-
考博
已经有3人回复
-
有没有快的中文核心比较快录用的，纳米材料光催化
已经有4人回复
-
本人42，博士刚毕业，现在找不到工作，怎么办？:(
已经有21人回复
-
有人投过CCC中国控制会议吗？
已经有3人回复
-
3,4-二羟基苯乙酮如何纯化？
已经有5人回复
-
国基评审
已经有10人回复
monitor2885
至尊木虫 (知名作家)
队长
- 应助: 6 (幼儿园)
- 贵宾: 0.339
- 金币: 39623
- 散金: 147
- 红花: 22
- 帖子: 5014
- 在线: 2175.5小时
- 虫号: 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
- 专业: 化学反应工程
★
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
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
4楼2010-03-17 23:04:14
monitor2885
至尊木虫 (知名作家)
队长
- 应助: 6 (幼儿园)
- 贵宾: 0.339
- 金币: 39623
- 散金: 147
- 红花: 22
- 帖子: 5014
- 在线: 2175.5小时
- 虫号: 701577
- 注册: 2009-02-15
- 性别: GG
- 专业: 生态系统生态学
|
global yy u v 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) global u v r=u+2; s=v-2; dydt=zeros(2,1); dydt(1)=r+y(2); dydt(2)=s*y(1)-2*s*y(2); 还是不对，再改改谢谢 |
5楼2010-03-18 14:28:35