| 查看: 3419 | 回复: 20 | |||||
kingfisher8288铜虫 (小有名气)
|
[交流]
【求助】matlab 求解微分方程中的未知参数已有4人参与
|
||||
|
求助，本人刚学matlab，不太熟悉。现在碰到一个问题。 一个反应体系会同时进行n个化学反应，彼此反应物和产物互相关联，我列了一系列的反应速率方程，当所有的反应速率常数已知时，用ode15s解微分方程可以得到任何反应物浓度随时间的变化。 现在其中有几个速率常数是未知的，但是我知道一种反应物的浓度随时间变化的曲线（反应时间及其对应的反应物浓度），用什么方法可以求出这几个未知的速率常数？ ps：已知的反应曲线数量是大于未知数的，我觉得理论上应该可以求出那几个速率常数。 各位大侠帮帮忙，感激不尽！ |
回复此楼» 收录本帖的淘帖专辑推荐
 matlab | DFT模拟 |
» 猜你喜欢
研一新生（无科研经历）求助
已经有13人回复
-
[吃瓜]网传江苏某高校人事处处长被打
已经有12人回复
-
多所高校已采取教师退出机制或教师末位淘汰机制。大家怎么看？
已经有18人回复
-
是博后，有国自然青年基金，省青年基金，文章5篇左右，能找到个什么样的岗位
已经有17人回复
-
请问有在北京市属高校工作的老师吗？待遇一般怎么样的
已经有11人回复
-
美国的华人副教授可以用中文交流吗
已经有3人回复
-
年末总结—基金挂蛋
已经有16人回复
-
浙大女生被MIT开除，女生称自己在浙江大学本科阶段就一直用这种方式收集数据发论文
已经有11人回复
-
nature comm, 浙大今年293篇，北大，288篇
已经有13人回复
-
自己找博士可真难啊
已经有22人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab解微分方程
已经有10人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 帮帮忙看看我的难题
已经有15人回复
- matlab解既有微分方程又有普通方程的方程组
已经有14人回复
- 【求助】用matlab计算一个常微分方程的解
已经有11人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】Matlab解分数阶微分方程问题
已经有18人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
★ ★
小木虫(金币+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
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
xiegangmai(金币+2): 谢谢应助！ 2011-04-15 21:57:50
小木虫(金币+0.5):给个红包，谢谢回帖交流
xiegangmai(金币+2): 谢谢应助！ 2011-04-15 21:57:50
|
本帖内容被屏蔽 |
9楼2011-04-15 08:24:52
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzuwangshilei(金币+1): 寻优很不方便 2011-04-15 09:49:53
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzuwangshilei(金币+1): 寻优很不方便 2011-04-15 09:49:53
|
本帖内容被屏蔽 |
10楼2011-04-15 08:37:47
xiegangmai
版主 (职业作家)
我没头衔
- 仿真EPI: 3
- 应助: 157 (高中生)
- 贵宾: 5.735
- 金币: 41393.5
- 散金: 2394
- 红花: 63
- 沙发: 4
- 帖子: 4044
- 在线: 1329.9小时
- 虫号: 653607
- 注册: 2008-11-13
- 性别: GG
- 专业: 机械动力学
- 管辖: 仿真模拟
2楼2011-04-10 16:21:49
kingfisher8288
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 197.2
- 散金: 47
- 帖子: 56
- 在线: 39.5小时
- 虫号: 708266
- 注册: 2009-02-25
- 性别: MM
- 专业: 污染控制化学
4楼2011-04-11 22:31:34
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小木虫(金币+0.5):给个红包，谢谢回帖交流
|
本帖内容被屏蔽 |
5楼2011-04-14 03:11:13
kingfisher8288
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 197.2
- 散金: 47
- 帖子: 56
- 在线: 39.5小时
- 虫号: 708266
- 注册: 2009-02-25
- 性别: MM
- 专业: 污染控制化学
★
xiegangmai(金币+1): 鼓励自己动手，积极讨论交流！ 2011-04-14 23:24:32
xiegangmai(金币+1): 鼓励自己动手，积极讨论交流！ 2011-04-14 23:24:32
|
function ydot = f(t, y) k1 = 3.5*10^4; k2 = 1.0*10^3; k3 = k=a*ka+b*kb+c*kc; r1 = k1*y(1)*y(2); r2 = k2*y(1)*y(3); r3 = k3*y(1)*y(4); z1 = -r1-r2-r3; z2 = -r1; z3 = -r2+r1; z4 = -r3+r2; ydot = [z1; z2; z3; z4;]; y0 = [30.0*10^-6;10.2*10^-6;4.1*10^-6;6.7] 当[a b c]=[100 5 20]时 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 ] 求ka，kb，kc 我自己的模型东西太多了，看着费劲。上面的例子和数据是我自己编的，但基本能描述我所要的东西，请帮忙解答，谢谢。 [ Last edited by kingfisher8288 on 2011-4-14 at 04:26 ] |
6楼2011-04-14 04:22:10
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
xiegangmai(金币+1): 谢谢建议！对照着数学模型，方便理解。 2011-04-14 23:25:24
小木虫(金币+0.5):给个红包，谢谢回帖交流
xiegangmai(金币+1): 谢谢建议！对照着数学模型，方便理解。 2011-04-14 23:25:24
|
本帖内容被屏蔽 |
7楼2011-04-14 19:46:12
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
5