| 查看: 4584 | 回复: 20 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
kingfisher8288铜虫 (小有名气)
|
[交流]
【求助】matlab 求解微分方程中的未知参数 已有4人参与
|
||||
|
求助，本人刚学matlab，不太熟悉。现在碰到一个问题。 一个反应体系会同时进行n个化学反应，彼此反应物和产物互相关联，我列了一系列的反应速率方程，当所有的反应速率常数已知时，用ode15s解微分方程可以得到任何反应物浓度随时间的变化。 现在其中有几个速率常数是未知的，但是我知道一种反应物的浓度随时间变化的曲线（反应时间及其对应的反应物浓度），用什么方法可以求出这几个未知的速率常数？ ps：已知的反应曲线数量是大于未知数的，我觉得理论上应该可以求出那几个速率常数。 各位大侠帮帮忙，感激不尽！ |
回复此楼» 收录本帖的淘帖专辑推荐
 matlab | DFT模拟 |
» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有10人回复
-
2026国自然函评费到账
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
范进中举一文的中心思想
已经有4人回复
-
跳槽后在研项目怎么办？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab解微分方程
已经有10人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 帮帮忙看看我的难题
已经有15人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
kingfisher8288
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 197.2
- 散金: 47
- 帖子: 56
- 在线: 39.5小时
- 虫号: 708266
- 注册: 2009-02-25
- 性别: MM
- 专业: 污染控制化学
|
谢谢forcal关注，我不知道什么是数学语言，我就大概解释一下吧，k1，k2，k3等是反应速率常数；r1，r2，r3 等是基元反应速率；y（1），y（2），y（3）等代表反应物的浓度；基元反应速率等于反应速率常数乘以相应的反应物浓度；z1，z2，z3 分别代表三种反应物的反应速率，它是由各个基元反应加和而成的，因为对于每一种物质来说有的基元反应是消耗该物质，有的是生成该物质。对于【y（1），y（2），y（3），y(4)】 有一组初始值 [30.0*10^-6;10.2*10^-6;4.1*10^-6;6.7]。k3是未知的，但是它和a，b，c 三个物性常数有关，k3 = k=a*ka+b*kb+c*kc; 当给定一组物性常数【a，b，c 】=[100， 5， 20]时，获得一组y（1）（一种反应物的浓度）对时间的数据t=[0, 10, 20, 30, 40] y(1)=[30.0*10^-6, 17.5*10^-6, 15*10^-6, 14*10^-6, 13.5*10^-6 ]； 当改变物性常数的值为[a, b, c]=[120, 15, 25]时， t=[0 10 20 30 40] y(1)=[30.0*10^-6, 12.5*10^-6, 10*10^-6, 8*10^-6, 7*10^-6 ]；当然还可以通过改变物性常数获得n组数据。我需要确定的是反应速率常数对于物性常数的相关系数 ka，kb，kc。ka，kb，kc 有可能为正值，也可能为负值，也就是说a，b，c可能正作用于k3，也可能负作用于k3. 不知道我的描述是不是能有助于你的理解，谢谢！ |
8楼2011-04-14 23:26:48
2楼2011-04-10 16:21:49
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
xiegangmai(金币+1): 谢谢应助！可否请虫友介绍一下forcal啊？ 2011-04-11 21:05:29
小木虫(金币+0.5):给个红包，谢谢回帖交流
xiegangmai(金币+1): 谢谢应助！可否请虫友介绍一下forcal啊？ 2011-04-11 21:05:29
|
本帖内容被屏蔽 |
3楼2011-04-11 16:47:03
kingfisher8288
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 197.2
- 散金: 47
- 帖子: 56
- 在线: 39.5小时
- 虫号: 708266
- 注册: 2009-02-25
- 性别: MM
- 专业: 污染控制化学
4楼2011-04-11 22:31:34
40