| 查看: 852 | 回复: 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也没有显示。谢谢 |
回复此楼» 猜你喜欢
0856材料与化工调剂，339
已经有9人回复
-
一志愿西安交大材料学硕（英一数二）347，求调剂到高分子/材料相关专业
已经有7人回复
-
324分 085600材料与化工
已经有9人回复
-
考研调剂
已经有6人回复
-
材料专硕306英一数二
已经有3人回复
-
材料调剂
已经有5人回复
-
309分085801求调剂
已经有6人回复
-
一志愿武理材料工程302调剂环化或化工
已经有10人回复
-
311（085601）求调剂
已经有13人回复
-
335求调剂
已经有4人回复
monitor2885
主管区长
队长
- 应助: 6 (幼儿园)
- 贵宾: 0.339
- 金币: 39452
- 散金: 147
- 红花: 22
- 帖子: 4992
- 在线: 2174.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
- 专业: 化学反应工程
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
