| 查看: 981 | 回复: 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也没有显示。谢谢 |
回复此楼» 猜你喜欢
据悉今年马上要出结果了
已经有4人回复
-
2026国自然放榜时间
已经有6人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有7人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有15人回复
-
面上项目filecode邪修
已经有7人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
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
- 专业: 化学反应工程
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
20