应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】简单的ode45程序修改
北京石油化工学院2026年研究生招生接收调剂公告
111/2返回列表12下一页
查看: 1228  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

monitor2885

至尊木虫 (职业作家)

队长

[交流] 【求助】简单的ode45程序修改

function abc
clear all;clc
global y u
u=[1;2;3;4;5];
t=0:4;
[t,y]=ode45(@fun,t,[0 2]);
plot(y)

function dydt=fun(t,y)
global u
v=u*2;
dydt=zeros(2,1);
dydt(1)=y(2);
dydt(2)=v*y(1)-2*v*y(2);

哪里错了呢？怎么改？谢谢
回复此楼

» 猜你喜欢
Retirement
1楼 2010-03-17 19:23:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师
★
monitor2885(金币+5):帮忙看下面一个 2010-03-17 21:33
kuhailangyu(金币+1):欢迎积极参与 2010-03-18 16:24
function abc
clear all
clc
u=[1;2;3;4;5];
tspan=0:1:4;
tt=[];yy=[];
y0=[0 2];
for i=1:length(tspan)-1
    [t,y]=ode45(@fun,[tspan(i),tspan(i+1)],y0,[],u(i));
    tt=[tt;t];yy=[yy;y];
    y0=y(end,: );
end
figure(1)
plot(tt,yy)

function dydt=fun(t,y,u)
v=u*2;
dydt=[y(2);
      v*y(1)-2*v*y(2)];
赞一下(1人) 回复此楼
2楼2010-03-17 19:45:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

monitor2885

至尊木虫 (职业作家)

队长
引用回帖:
Originally posted by change0618 at 2010-03-17 19:45:45:
function abc
clear all
clc
u=[1;2;3;4;5];
tspan=0:1:4;
tt=[];yy=[];
y0=[0 2];
for i=1:length(tspan)-1
    [t,y]=ode45(@fun,[tspan(i),tspan(i+1)],y0,[],u(i));
    tt=[tt;t];yy=[yy;y];
     ...

比如，现在2个参数，即u和v，下面哪里不对呢？
clear all;clc
global yy
u=[1;5;15;20;25];
v=[6;12;18;24;30];
tspan=0:1:4;
tt=[];yy=[];
y0=[0 2];
for i=1:length(tspan)-1
    [t,y]=ode45(@fun,[tspan(i),tspan(i+1)],y0,[],(u(i),v(i)));
    tt=[tt;t];
    yy=[yy;y];
    y0=y(end,;
end
figure;plot(tt,yy)

function dydt=fun(t,y,u,v)
dydt=zeros(2,1);
dydt(1)=u+y(2);
dydt(2)=v*y(1)-2*v*y(2);

[ Last edited by monitor2885 on 2010-3-17 at 21:34 ]
赞一下 回复此楼
Retirement
3楼2010-03-17 21:33:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师
★ ★
monitor2885(金币+5):继续发金币，谢谢哥们帮忙解决问题 2010-03-18 14:29
kuhailangyu(金币+2):欢迎积极参与 2010-03-18 16:24
[t,y]=ode45(@fun,[tspan(i),tspan(i+1)],y0,[],      (u(i),v(i))      );
改成
[t,y]=ode45(@fun,[tspan(i),tspan(i+1)],y0,[],      u(i),v(i)        );


global yy  可以去掉
赞一下(1人) 回复此楼
4楼2010-03-17 23:04:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

monitor2885

至尊木虫 (职业作家)

队长
引用回帖:
Originally posted by change0618 at 2010-03-17 23:04:14:
[t,y]=ode45(@fun,[tspan(i),tspan(i+1)],y0,[],      (u(i),v(i))      );
改成
[t,y]=ode45(@fun,[tspan(i),tspan(i+1)],y0,[],      u(i),v(i)        );


global yy  可以去掉

global yy u v
u=[1;5;15;20;25];
v=[6;12;18;24;30];
tspan=0:1:4;
tt=[];yy=[];
y0=[0 2];
for i=1:length(tspan)-1
    [t,y]=ode45(@fun,[tspan(i),tspan(i+1)],y0,[],u(i),v(i));
    tt=[tt;t];
    yy=[yy;y];
    y0=y(end,;
end
figure;plot(tt,yy)

function dydt=fun(t,y)
global u v
r=u+2;
s=v-2;
dydt=zeros(2,1);
dydt(1)=r+y(2);
dydt(2)=s*y(1)-2*s*y(2);

还是不对，再改改谢谢
赞一下 回复此楼
Retirement
5楼2010-03-18 14:28:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师
★
kuhailangyu(金币+1):欢迎积极参与 2010-03-18 16:24
monitor2885(金币+5):global yy之后，我发现yy有208个，为何不是5个呢？ 2010-03-18 19:43
function abc
clear all
clc
u=[1;5;15;20;25];
v=[6;12;18;24;30];

tspan=0:1:4;
tt=[];yy=[];
y0=[0 2];
for i=1:length(tspan)-1
    [t,y]=ode45(@fun,[tspan(i),tspan(i+1)],y0,[],u(i),v(i));
    tt=[tt;t];
    yy=[yy;y];
    y0=y(end,: );
end

plot(tt,yy)


function dydt=fun(t,y,u,v)
r=u+2;
s=v-2;
dydt=[r+y(2);
      s*y(1)-2*s*y(2)];
赞一下(1人) 回复此楼
6楼2010-03-18 15:38:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

张轩中

铜虫 (小有名气)
??? function abc
    |
Error: Function definitions are not permitted at the prompt or in scripts.
赞一下 回复此楼
7楼2010-03-18 16:14:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

monitor2885

至尊木虫 (职业作家)

队长
引用回帖:
Originally posted by change0618 at 2010-03-18 15:38:18:
function abc
clear all
clc
u=[1;5;15;20;25];
v=[6;12;18;24;30];

tspan=0:1:4;
tt=[];yy=[];
y0=[0 2];
for i=1:length(tspan)-1
    [t,y]=ode45(@fun,[tspan(i),tspan(i+1)],y0,[],u(i),v(i));
...

global yy之后，我发现yy有208个，为何不是5个呢？我知道ode45是变步长，但是如果写成ode45(@fun,tspan,y0),就会有5个yy，而不是208个。为什么
比如说，下面的简单程序：y就是5个
global y
tspan=0:1:4;
y0=[0 2];
[t,y]=ode45(@fun,tspan,y0);

function dydt=fun(t,y)
dydt=zeros(2,1);
dydt(1)=y(2);
dydt(2)=y(1)-2*y(2);

[ Last edited by monitor2885 on 2010-3-18 at 19:48 ]
赞一下 回复此楼
Retirement
8楼2010-03-18 19:45:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师
★
adu886886(金币+1):谢谢提供意见 2010-03-19 08:29
monitor2885(金币+5):[t,y]=ode45(@fun,[tspan(i):1:tspan(i+1)],y0,[],u(i),v(i)); 还是208个点呀 2010-03-19 13:09
tspan=0:1:4;
[t,y]=ode45(@fun,tspan,y0);
这两句会令t=[0 1 2 3 4],并将其对应的值赋给y，所以y只有5个；

而[t,y]=ode45(@fun,[tspan(i),tspan(i+1)],y0,[],u(i),v(i));
此句中的[tspan(i),tspan(i+1)],只给定了一个区间，
t的元素个数一般取其计算过程中的个数
如果你只想取五个值，那就用tspan(i):tspan(i+1),
你画图的话会看到五个点，但实际上是8个点，中间三个有重复
如果想彻底只有5个，你可以自己抠掉。
当然还有其他方法抽取那五个点的y值。
程序肯定是计算了很多很多个点，但是具体你要取几个，那在于你的意图了。
赞一下(1人) 回复此楼
9楼2010-03-18 21:46:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

monitor2885

至尊木虫 (职业作家)

队长
引用回帖:
Originally posted by change0618 at 2010-03-18 21:46:37:
tspan=0:1:4;
[t,y]=ode45(@fun,tspan,y0);
这两句会令t=[0 1 2 3 4],并将其对应的值赋给y，所以y只有5个；

而[t,y]=ode45(@fun,[tspan(i),tspan(i+1)],y0,[],u(i),v(i));
此句中的[tspan(i),tspan(i+1)],只 ...

[t,y]=ode45(@fun,[tspan(i):1:tspan(i+1)],y0,[],u(i),v(i));
还是 208 个点呀
赞一下 回复此楼
Retirement
10楼2010-03-19 13:09:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 monitor2885 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求材料调剂，一志愿郑州大学289分 +13 硕星赴 2026-04-03 13/650 2026-04-03 11:41 by 纯净水2025
[考研] 环境285分，过六级，求调剂 +9 xhr12 2026-04-02 9/450 2026-04-03 10:23 by 化学化工硕士招�
[考研] 325分化学调剂 +5 15771691647 2026-04-02 5/250 2026-04-03 09:58 by ChemPharm
[考研] 理学07化学 303求调剂 +17 睿08 2026-03-27 18/900 2026-04-03 09:30 by zhouxiaoyu
[考研] 282求调剂 +3 aaa车辆 2026-04-02 3/150 2026-04-02 21:55 by zllcz
[考研] 343求调剂085601 +5 要努力学习x 2026-03-29 6/300 2026-04-02 21:18 by 1104338198
[考研] 070305高分子化学与物理 304分求调剂 +16 c297914 2026-03-28 16/800 2026-04-02 20:50 by dongzh2009
[考研] 293求调剂 +4 珂珂乐 2026-04-02 4/200 2026-04-02 20:10 by 6781022
[考研] 270求调剂 +8 小杰pp 2026-03-31 10/500 2026-04-02 12:57 by yulian1987
[考研] 085900土木水利336分求调剂 +4 Zhangjiangj 2026-03-31 6/300 2026-04-02 11:40 by 1753564080
[考研] 0805求调剂 +8 是水分 2026-03-31 8/400 2026-04-02 10:46 by guanxin1001
[考研] 一志愿346上海大学生物学 +3 上海大学346调剂 2026-04-01 3/150 2026-04-02 08:36 by w虫虫123
[考研] 求调剂，一志愿北林食品与营养095500，301分，已过六级，有科研经历 +4 快乐储蓄罐 2026-03-31 4/200 2026-04-01 09:26 by JourneyLucky
[考研] 070300化学专业279调剂 +10 哈哈哈^_^ 2026-03-31 10/500 2026-03-31 23:13 by liu823948201
[考研] 求调剂 生物学 377分 +6 zzll03 2026-03-31 6/300 2026-03-31 17:33 by 唐沐儿
[考研] 生物考研337分求调剂 +4 cgxin 2026-03-30 6/300 2026-03-31 14:18 by 记事本2026
[考研] 0703化学 +20 妮妮ninicgb 2026-03-27 20/1000 2026-03-31 13:33 by 无际的草原
[考研] 313求调剂 +6 卖个关子吧 2026-03-31 6/300 2026-03-31 10:58 by Jaylen.
[考研] 293分求调剂，外语为俄语 +5 加一一九 2026-03-31 5/250 2026-03-31 09:39 by zhshch
[考研] 356求调剂 +4 gysy?s?a 2026-03-28 4/200 2026-03-29 10:32 by 唐沐儿
信息提示
关闭
请填处理意见
关闭