| 查看: 1036 | 回复: 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也没有显示。谢谢 |
回复此楼» 猜你喜欢
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
梦想
已经有4人回复
-
中青基了要发朋友圈吗？
已经有6人回复
-
基金不中，共勉
已经有12人回复
-
能否申诉？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上意见出来了
已经有19人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
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
- 金币: 40132
- 散金: 147
- 红花: 22
- 帖子: 5076
- 在线: 2190.1小时
- 虫号: 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
- 金币: 40132
- 散金: 147
- 红花: 22
- 帖子: 5076
- 在线: 2190.1小时
- 虫号: 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
投票:
30