声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 896|回复: 0

[编程技巧] 为什么ode45解方程的结果不变呀?

[复制链接]
发表于 2009-5-18 15:45 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
function shiyan
promrt={'请输入起始时刻t0','请输入终止时刻t_end'};
name='请输入积分时间段.';
numlines=1;
defAns={'0','1'};
Resize='on';
R=inputdlg(promrt,name,numlines,defAns,'on');
t0=str2num(R{1})
t_end=str2num(R{2})
q10=[1;2;3;4;5;6];
[T,Y]=ode45(@shiyan01,[t0:0.2:t_end],q10)
function F=shiyan01(t,q1)
dy1=q1(1);
dy2=q1(2);
dy3=q1(3);
dy4=5;
dy5=4;
dy6=1;
aq1=[dy4;dy5;dy6];
F=[dy1;dy2;dy3;aq1];
为什么这个微风方程无论如何改变其积分区间,所得的结果都不变呀?求求那位高手帮忙解决一下呀,很着急的用。拜谢啦!!
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-14 17:46 , Processed in 0.139314 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表