声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 999|回复: 1

[编程技巧] 【求助】微分方程与普通方程组成的方程组求解

[复制链接]
发表于 2012-3-19 18:10 | 显示全部楼层 |阅读模式

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

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

x
以下方程系数可能有些复杂,请大牛们看完后给出代码,数字可换用简单的。谢谢
dSdt=-(7.37*10^7/(1+2.48*10^-3*A)+1.46*10^7/(1+2.48*10^-3*A)/H+3570/(1+2.48*10^-3*A)/H^2)*S*O-1.09*10^15/(1+13*H)*S*B

dOdt=-2.33*10^-5*(7.37*10^7/(1+2.48*10^-3*A)+1.46*10^7/(1+2.48*10^-3*A)/H+3570/(1+2.48*10^-3*A)/H^2)*S*O

dBdt=-1.08*10^15/(1+13*H)*S*B

A=1.24*(1+0.0129/H+7.76*10^-10/(H^2))
H 满足 H+2.17*10^10/(1+62(1+1.709*10^9*H))*H=10^-14/H+1.46*10^7/(1+2.48*10^-3*A)*S/H+7140/(1+2.48*10^-3*A)*S/H^2+2*(0.2-10^6*S)/(H^2/10.2+1+H/0.0102)+(0.2-10^6*S)/(H/1000+1+0.0102/H)

t=0时, H=10^-6.17, S=2*10^-7, O=4.7*10^-11, B=7.1*10^-5
求t从0到3600, H,S,O,B的数值解
求代码以供学习
总的来说,这是三个常微分方程和一个普通方程组成的方程组,知道初值,请问如何求数值解;或者说是变系数常微分方程组,其中系数被另一个复杂的分式方程限制。

另弱问这种情况ode45好使吗?还是我太弱。。。
回复
分享到:

使用道具 举报

发表于 2012-3-22 16:36 | 显示全部楼层
搜微分代数方程组   可解
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-26 02:27 , Processed in 0.087272 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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