| 查看: 986 | 回复: 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也没有显示。谢谢 |
回复此楼» 猜你喜欢
综述论文作为代表作会不会影响评审专家的印象分？
已经有10人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有3人回复
-
国自然结果
已经有5人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有12人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有8人回复
-
据悉今年马上要出结果了
已经有8人回复
-
关于代码变化问题，想知道的进来
已经有23人回复
-
FileCode能看出啥？
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有20人回复
-
2026国自然放榜时间
已经有9人回复
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
- 专业: 化学反应工程
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
- 专业: 化学反应工程
5楼2009-10-12 19:32:54
100