24小时热门版块排行榜    

查看: 1226  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +10 Tide man 2026-08-11 10/500 2026-08-12 20:59 by 陈晨晨陈啊
[基金申请] 不应该看fileCode +5 且听虎啸 2026-08-12 6/300 2026-08-12 16:26 by Tide man
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] filecode +13 documentary 2026-08-10 15/750 2026-08-12 15:56 by 云上清扬
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
请填处理意见