24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1506  |  回复: 8

小虫6686

铜虫 (初入文坛)


[交流] 【求助】急求各位大侠,关于用ode45解微分方程

本人想解一组耦合的微分方程,大概形式为
dy/dT=a*y*x;
dx/dT=b*y/x;
应该可以用ode45来解,但是遇到了个问题,请教各位大侠。
其中的a 和b 不是常数,而是由x 和y 组成的,而且有循环行为,大概形式为
a=1-exp(x*a);
求各位帮忙解答,这种情况应该如何求解


我尝试把a写成变量,在写m函数是把a 作为一个方程,但不知道ode45是否可以包含非微分方程的情况,望各位牛人给指条明路。

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小虫6686

铜虫 (初入文坛)


追加条件

也给个b的具体形式

b=1-exp(y*b^2)

多谢各位了
2楼2011-01-20 13:52:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
小虫6686(金币+3): 希望给出进一步解答,多谢 2011-01-22 13:08:48
引用回帖:
Originally posted by 小虫6686 at 2011-01-20 13:36:11:
本人想解一组耦合的微分方程,大概形式为
dy/dT=a*y*x;
dx/dT=b*y/x;
应该可以用ode45来解,但是遇到了个问题,请教各位大侠。
其中的a 和b 不是常数,而是由x 和y 组成的,而且有循环行为,大概形式为
a=1- ...

初值条件呢?
a = 0
b = 0
x = 1
y = 1
算不算一组解呢?

[ Last edited by xiegangmai on 2011-1-20 at 18:48 ]
3楼2011-01-20 18:46:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)



xiegangmai(金币+1): 谢谢应助,Simulink应该能解这类问题。 2011-01-20 21:04:36
小虫6686(金币+1): 感谢 2011-01-21 12:35:25
干脆simulink就是了
4楼2011-01-20 19:17:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)


好像不太行,貌似出现数据环
5楼2011-01-21 00:47:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小虫6686

铜虫 (初入文坛)


引用回帖:
Originally posted by xiegangmai at 2011-01-20 18:46:01:


初值条件呢?
a = 0
b = 0
x = 1
y = 1
算不算一组解呢?

[ Last edited by xiegangmai on 2011-1-20 at 18:48 ]

初始条件类似,
在T=0时,满足a = 1
b = 0
x = 1
y = 0
这样的话能解?
6楼2011-01-21 12:33:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)


引用回帖:
Originally posted by 小虫6686 at 2011-01-21 12:33:05:



初始条件类似,
在T=0时,满足a = 1
b = 0
x = 1
y = 0
这样的话能解?

这样是不是就有解析解了
7楼2011-01-22 14:13:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WF1987

铜虫 (小有名气)


我试试看
8楼2011-12-23 17:22:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WF1987

铜虫 (小有名气)



小木虫(金币+0.5):给个红包,谢谢回帖
没有成功,盼高手出现。
[t y]=ode45('weifenzu',[0 20],[1 0]);
plot(t,y(:,1),'-')
%………………………………………………
function dydx=weifenzu(t,y)
    a=1-exp(y(1)*a);%如果能把a、b都化到左边就好办了
   b=1-exp(y(2)*b^2);
    dy(1)=b*y(2)/y(1);
    dy(2)=a*y(2)*y(1);
9楼2011-12-23 20:44:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小虫6686 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 找调剂 +10 楚乔乔 2026-04-01 10/500 2026-04-05 22:19 by syh9288
[考研] 复试调剂 +13 呼呼?~+123456 2026-04-05 13/650 2026-04-05 22:07 by 醉翁wl
[考研] 0703化学321分求调剂 +17 三dd. 2026-03-30 18/900 2026-04-05 18:07 by 蓝云思雨
[考研] 复试调剂 +3 asdasdassda 2026-04-05 3/150 2026-04-05 17:26 by zhousanduo
[考研] 328分调剂 +6 门men 2026-04-04 6/300 2026-04-05 13:40 by imissbao
[考研] 一志愿江南大学085501机械工程专硕326分,本科佳木斯大学 +5 顾若浮生 2026-04-03 9/450 2026-04-05 09:57 by 1753564080
[考研] 283求调剂 +10 A child 2026-04-04 10/500 2026-04-05 08:22 by qlm5820
[考研] 能动调剂326专硕 +4 wan112233 2026-04-04 4/200 2026-04-04 22:47 by yu221
[考研] 278求调剂 +3 依旧! 2026-04-02 4/200 2026-04-04 20:27 by 蓝云思雨
[考研] 一志愿上海大学生物学346 +3 上海大学346调剂 2026-04-03 3/150 2026-04-04 20:20 by dongzh2009
[考研] 一志愿华北电力大学(北京),材料科学与工程学硕265,求调剂 +11 yelck 2026-04-03 12/600 2026-04-04 19:52 by dongzh2009
[考研] 11408 一志愿西电,277分求调剂 +4 zhouzhen654 2026-04-03 4/200 2026-04-04 18:10 by 猪会飞
[考研] 一志愿南昌大学324求调剂 +9 hanamiko 2026-03-30 9/450 2026-04-04 11:04 by 猪会飞
[考研] 350一志愿北京航空航天大学08500材料科学与工程求调剂 +5 kjnasfss 2026-04-03 5/250 2026-04-03 22:29 by 无际的草原
[考研] 295求调剂 +3 尚偌呀 2026-04-03 4/200 2026-04-03 21:23 by zhq0425
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-03 16:46 by wxiongid
[考研] 085600 295分求调剂 +19 W55j 2026-03-30 23/1150 2026-04-03 09:53 by 千千运气
[考研] 材料340分调剂 +7 夏夜晚风_long 2026-04-02 9/450 2026-04-02 21:20 by dongzh2009
[考研] 085410 一志愿211 22408分数359求调剂 +3 123456789qw 2026-03-31 4/200 2026-04-02 00:06 by 义文wang
[考研] 生物与医药考研调剂 +5 铁憨憨123425 2026-03-31 5/250 2026-04-01 18:01 by syh9288
信息提示
请填处理意见