24小时热门版块排行榜    

查看: 3360  |  回复: 20

kingfisher8288

铜虫 (小有名气)

[交流] 【求助】matlab 求解微分方程中的未知参数已有4人参与

求助,本人刚学matlab,不太熟悉。现在碰到一个问题。
一个反应体系会同时进行n个化学反应,彼此反应物和产物互相关联,我列了一系列的反应速率方程,当所有的反应速率常数已知时,用ode15s解微分方程可以得到任何反应物浓度随时间的变化。
现在其中有几个速率常数是未知的,但是我知道一种反应物的浓度随时间变化的曲线(反应时间及其对应的反应物浓度),用什么方法可以求出这几个未知的速率常数?
ps:已知的反应曲线数量是大于未知数的,我觉得理论上应该可以求出那几个速率常数。
各位大侠帮帮忙,感激不尽!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

forcal

禁虫 (初入文坛)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
xiegangmai(金币+1): 谢谢应助!可否请虫友介绍一下forcal啊? 2011-04-11 21:05:29
本帖内容被屏蔽

3楼2011-04-11 16:47:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forcal

禁虫 (初入文坛)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
xiegangmai(金币+2): 谢谢应助! 2011-04-15 21:57:50
本帖内容被屏蔽

9楼2011-04-15 08:24:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forcal

禁虫 (初入文坛)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzuwangshilei(金币+1): 寻优很不方便 2011-04-15 09:49:53
本帖内容被屏蔽

10楼2011-04-15 08:37:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by kingfisher8288 at 2011-04-10 15:14:13:
求助,本人刚学matlab,不太熟悉。现在碰到一个问题。
一个反应体系会同时进行n个化学反应,彼此反应物和产物互相关联,我列了一系列的反应速率方程,当所有的反应速率常数已知时,用ode15s解微分方程可以得到任 ...

感觉这个问题有些难度。

微分方程组中,已知一种反应物的浓度随时间变化的曲线,而有几个未知系数,需把系数拟合出来,好像还没看到过类似的做法。

鼓励虫友们讨论。
明德厚学、求是创新
2楼2011-04-10 16:21:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kingfisher8288

铜虫 (小有名气)

谢谢

在上面的系统中,还有个问题,就是如果未知的反应系数k=a*k1+b*k2+c*k3,k1,k2,k3是三个常数,a,b,c是三个物性参数,当给定一组a,b,c的值时,会有对应的一组一种反应物随浓度的变化曲线。
也就是比如有n组a,b,c的值,及其对应的n组一种反应物随浓度的变化曲线,如何确定k1,k2,k3的值呢?
4楼2011-04-11 22:31:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forcal

禁虫 (初入文坛)


小木虫(金币+0.5):给个红包,谢谢回帖交流
本帖内容被屏蔽

5楼2011-04-14 03:11:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kingfisher8288

铜虫 (小有名气)


xiegangmai(金币+1): 鼓励自己动手,积极讨论交流! 2011-04-14 23:24:32
引用回帖:
Originally posted by forcal at 2011-04-14 03:11:13:
属于优化拟合问题,请给出实际例子(公式模型及数据)以方便讨论。

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forcal

禁虫 (初入文坛)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
xiegangmai(金币+1): 谢谢建议!对照着数学模型,方便理解。 2011-04-14 23:25:24
本帖内容被屏蔽

7楼2011-04-14 19:46:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kingfisher8288

铜虫 (小有名气)

引用回帖:
Originally posted by forcal at 2011-04-14 19:46:12:
matlab代码看不太懂,能否用普通的数学语言描述一下?

另外,优化拟合问题,可以求得全局最优解,但不一定满足用户要求,例如公式模型不合适,或者解不符合物理意义等。故还是完整描述问题比较好。

谢谢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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kingfisher8288 的主题更新
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 职称论文 +8 xiaoqi_66 2024-11-09 9/450 2024-11-14 14:10 by cstiger555
[找工作] 咨询一下江西的高校待遇,人文氛围怎么样? +5 akslis2024 2024-11-09 5/250 2024-11-14 13:53 by 啄木鸟、
[硕博家园] 研究生的生活该是什么样 +3 lqy0719 2024-11-14 3/150 2024-11-14 13:47 by 破戒刀A
[教师之家] 大学老师 +7 考研一路顺风 2024-11-13 9/450 2024-11-14 13:31 by 考研一路顺风
[基金申请] 广东省自然科学基金-面上项目消息有了么? +3 pachang 2024-11-14 3/150 2024-11-14 12:29 by zb789
[基金申请] 76批博后基金 +3 feiyi3986 2024-11-14 3/150 2024-11-14 11:50 by puly
[有机交流] 硫谜用双氧水氧化出现的困惑(老板说一定要用双氧水氧化做好)求指教 10+4 luowenwu 2024-11-10 15/750 2024-11-14 10:52 by luolinfeng
[论文投稿] 要不要撤稿另投 10+4 wangzhesd 2024-11-09 8/400 2024-11-14 09:09 by LIU_V
[硕博家园] 实验进度 +8 麻辣烫。 2024-11-11 13/650 2024-11-13 23:25 by 麻辣烫。
[论文投稿] 找不到审稿人 +4 吉祥如意141913 2024-11-13 4/200 2024-11-13 16:33 by 13333105938
[论文投稿] 核心期刊那个阶段最难 +3 公共茅斯 2024-11-11 3/150 2024-11-13 15:27 by 13333105938
[论文投稿] 投稿系统中的通讯作者和文章中的通讯作者不一样,文章目前被录用了? +3 babybabygo 2024-11-12 4/200 2024-11-13 12:40 by babybabygo
[论文投稿] 请问有EI期刊生物工程类的比较靠谱的代理推荐吗 20+3 Chen2000c 2024-11-07 4/200 2024-11-13 11:01 by 国际期刊小编
[教师之家] 北大教授何怀宏曾如此描述他的同行 +9 zju2000 2024-11-09 9/450 2024-11-13 10:11 by tfang
[论文投稿] 求JACS和ACIE的最新文章模板 10+4 Jane060390 2024-11-10 4/200 2024-11-13 09:13 by 北京莱茵润色
[基金申请] 变性人申请基金咋算性别? +19 leicahfs 2024-11-11 19/950 2024-11-13 00:09 by cw277
[论文投稿] 二审审一年的佛系编辑 +10 thefoxrain 2024-11-09 15/750 2024-11-12 19:27 by lide966
[论文投稿] 投稿求助 40+3 liukun11 2024-11-12 4/200 2024-11-12 19:13 by liukun11
[有机交流] 求助NMR +5 苯巴比妥! 2024-11-09 5/250 2024-11-12 11:20 by 88817753
[论文投稿] renewable energy 添加作者 15+3 Hebauwww12 2024-11-09 4/200 2024-11-11 08:42 by 北京莱茵润色
信息提示
请填处理意见