24小时热门版块排行榜    

查看: 1030  |  回复: 5

priest

木虫 (正式写手)

[求助] 求助:Matlab求解微分代数方程 已有1人参与

大家好,我现在想重复一个已经发表的模型的计算结果,模型包括16个方程,5个微分方程,11个代数方程,变量用荧光笔标记出来了。想请问大家,怎么用Matlab ode15s解呢,或者提供一些思路? 具体的方程列在上传的图片。多谢大家!可以追加金币!

求助:Matlab求解微分代数方程
IMG_1695.jpg
回复此楼

» 猜你喜欢

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

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

木头608

金虫 (著名写手)

帮你up一下
论成败,人生豪迈,只不过需要等待!
2楼2014-10-05 23:23:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

应该很简单,你把没有荧光的物理量的值给定,并给出初始条件,就可以解出来。
3楼2014-10-06 21:47:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyang6

金虫 (小有名气)

你的问题解决了吗?我也遇到这个问题
4楼2017-12-13 10:49:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ecnufsy

金虫 (正式写手)

【答案】应助回帖

微分代数方程 DAE求解和 常微分方程的隐式解法一样,代数方程和差分后的微分方程相同对待即可。
看图中的代数方程式显式表达式,貌似可以带入微分方程里,这样将代数方程消除。然后求5个解纯微分方程。
CODE:
matlab

的话写一个
CODE:
@dfun(t, y)

函数,其中
CODE:
y(1) = [K]i

CODE:
y(2) = [Ca]i

CODE:
y(3) = [Na]i

CODE:
y(4) = x

CODE:
y(5) = h


然后
CODE:
dy

的表达式是可以用11个显式的代数方程代进去的,仔细写
CODE:
dy

就行了。
5楼2017-12-13 13:44:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ecnufsy

金虫 (正式写手)

引用回帖:
5楼: Originally posted by ecnufsy at 2017-12-13 13:44:47
微分代数方程 DAE求解和 常微分方程的隐式解法一样,代数方程和差分后的微分方程相同对待即可。
看图中的代数方程式显式表达式,貌似可以带入微分方程里,这样将代数方程消除。然后求5个解纯微分方程。
用matlab的 ...

没有preview功能,也不能编辑,排版好丑
6楼2017-12-13 13:46:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 priest 的主题更新
信息提示
请填处理意见