| 查看: 833 | 回复: 14 | |||
| 当前主题已经存档。 | |||
monitor2885至尊木虫 (职业作家)
队长
|
[交流]
【求助】简单的ode45函数参数
|
||
|
function dydt=ivpodefun(t,y,u) a=[1;2;3;4;5];b=[6;7;8;9;10]; for i=1:1:5 u1((i-1)*2<=t & t<=2*i)=a(i); u2((i-1)*2<=t & t<=2*i)=b(i); end plot(t,u1);hold on plot(t,u2);hold off dydt=zeros(2,1); dydt(1)=u2*y(2); dydt(2)=u1*(1-y(1)^2)*y(2)-y(1); 命令窗口输入 t=0:10; [t,y]=ode45(@ivpodefun,t,[0 2]); plot(t,y) 似乎结果没错，u1和u2在计算过程中的确取到了不同值，但是为什么它俩打印不出来呢？只有一个点。我需要这两个参数显示出来，在workspace也没有显示。谢谢 |
回复此楼» 猜你喜欢
能源材料化学课题组招收硕士研究生8-10名
已经有11人回复
-
求材料调剂
已经有9人回复
-
材料专硕英一数二306
已经有5人回复
-
085600材料与化工调剂 324分
已经有9人回复
-
0703化学调剂
已经有12人回复
-
0703化学 305求调剂
已经有4人回复
-
0703化学调剂，求各位老师收留
已经有10人回复
-
271材料工程求调剂
已经有5人回复
-
281求调剂（0805）
已经有16人回复
-
304求调剂
已经有6人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
2楼2009-10-12 18:40:05
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
刚弄明白你啥意思3楼2009-10-12 18:52:20
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
★ ★ ★ ★ ★ ★
monitor2885(金币+6,VIP+0):分成2个图就死机了，怎么回事 10-12 20:16
monitor2885(金币+6,VIP+0):分成2个图就死机了，怎么回事 10-12 20:16
|
function test t=0:10; [t,y]=ode45(@ivpodefun,t,[0 2]); figure(2) plot(t,y,'bo-') function dydt=ivpodefun(t,y,u) a=[1;2;3;4;5];b=[6;7;8;9;10]; for i=1:1:5 u1((i-1)*2<=t & t<=2*i)=a(i); u2((i-1)*2<=t & t<=2*i)=b(i); end figure(1) plot(t,u1,'ro');hold on plot(t,u2,'b^'); xlim([0 10]) ylim([0 10]) dydt=zeros(2,1); dydt(1)=u2*y(2); dydt(2)=u1*(1-y(1)^2)*y(2)-y(1); [ Last edited by change0618 on 2009-10-12 at 19:18 ] |
4楼2009-10-12 19:15:10
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
5楼2009-10-12 19:32:54
monitor2885
至尊木虫 (职业作家)
队长
- 应助: 6 (幼儿园)
- 贵宾: 0.339
- 金币: 39404.5
- 散金: 147
- 红花: 22
- 帖子: 4990
- 在线: 2171.8小时
- 虫号: 701577
- 注册: 2009-02-15
- 性别: GG
- 专业: 生态系统生态学
6楼2009-10-12 20:16:52
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
7楼2009-10-12 20:29:19
monitor2885
至尊木虫 (职业作家)
队长
- 应助: 6 (幼儿园)
- 贵宾: 0.339
- 金币: 39404.5
- 散金: 147
- 红花: 22
- 帖子: 4990
- 在线: 2171.8小时
- 虫号: 701577
- 注册: 2009-02-15
- 性别: GG
- 专业: 生态系统生态学
8楼2009-10-12 21:20:44
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
9楼2009-10-12 21:45:15
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
10楼2009-10-12 22:30:36
