24小时热门版块排行榜    

查看: 2400  |  回复: 7

晴朗1989

金虫 (小有名气)

[求助] matlab拟合酶促反应机制

酶促酯化反应动力学符合乒乓模型,请问怎样用matlab来拟合实验数值,求得动力学参数啊
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与 2013-08-27 17:10:57
y有具体模型和数据吗?详细点一起放上来。
2楼2013-08-27 09:20:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

晴朗1989

金虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2013-08-27 09:20:14
y有具体模型和数据吗?详细点一起放上来。

具体的数值都写在word里面了 帮忙求解下 谢谢
走遍天下
3楼2013-08-28 22:45:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

你的word文档在哪?
4楼2013-08-29 09:04:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

晴朗1989

金虫 (小有名气)

引用回帖:
4楼: Originally posted by dingd at 2013-08-29 09:04:12
你的word文档在哪?

不好意思 word没传上去

公式



数值

A/B   80                160           240          320         480           
5   0.469263 0.360785 0.2500821 0.242265 0.121735233
15  0.621502        0.528029        0.411572        0.363633        0.312755
25  0.706035        0.788566        1.222889        1.07098        0.674924
30  0.674924        1.027269        1.531503        0.989152        0.85732
40  1.265887        1.620876        1.640166   1.453717        1.352839

第一列为A 第一行为B 其余的为V 要求出Kma kmb kib  Vmax
走遍天下
5楼2013-08-29 10:18:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

晴朗1989

金虫 (小有名气)

引用回帖:
4楼: Originally posted by dingd at 2013-08-29 09:04:12
你的word文档在哪?

公式是这个 上面的有点问题
matlab拟合酶促反应机制
公式.jpg

走遍天下
6楼2013-08-29 10:21:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

1stOpt求解代码(请查下公式是否正确):
CODE:
Variable A,B,v;
Function v=Vmax*A*B/(Kmb*A+Kma*B*(1+B/Kib)+A*B);
Data;
//A, B, V
5        80        0.469263
15        80        0.621502
25        80        0.706035
30        80        0.674924
40        80        1.265887
5        160        0.360785
15        160        0.528029
25        160        0.788566
30        160        1.027269
40        160        1.620876
5        240        0.2500821
15        240        0.411572
25        240        1.222889
30        240        1.531503
40        240        1.640166
5        320        0.242265
15        320        0.363633
25        320        1.07098
30        320        0.989152
40        320        1.453717
5        480        0.121735233
15        480        0.312755
25        480        0.674924
30        480        0.85732
40        480        1.352839

结果:
均方差(RMSE):0.151984472391896
残差平方和(SSE):0.577481996206077
相关系数(R): 0.944715339166438
相关系数之平方(R^2): 0.892487072056358
确定系数(DC): 0.890957450482387
卡方系数(Chi-Square): 0.418122466961142
F统计(F-Statistic): 64.1951241140902

参数        最佳估算
----------        -------------
vmax        -2.54549931799406
kmb        -174.912897327456
kma        -39.6638546619627
kib        294.760121529861
matlab拟合酶促反应机制-1
c2.jpg


matlab拟合酶促反应机制-2
c1.jpg

7楼2013-08-29 11:25:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

晴朗1989

金虫 (小有名气)

引用回帖:
7楼: Originally posted by dingd at 2013-08-29 11:25:00
1stOpt求解代码(请查下公式是否正确):


Variable A,B,v;
Function v=Vmax*A*B/(Kmb*A+Kma*B*(1+B/Kib)+A*B);
Data;
//A, B, V
5        80        0.469263
15        80        0.621502
25        80        0.706035
30        80        0.674924
40        80        1 ...

非常感谢
公式是对的 但是值怎么求出来是负的呢 还有我看人家做的图都是那样的 能帮忙解释下做出来的图不
matlab拟合酶促反应机制-3
1.jpg

走遍天下
8楼2013-08-29 12:19:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 晴朗1989 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳又变了 +7 wuchongjun 2026-08-20 8/400 2026-08-20 15:49 by 开心的小狮子
[基金申请] 言之凿凿,8月21号(明天)放榜 +14 医学老男孩 2026-08-20 16/800 2026-08-20 15:47 by LNP@mRNA
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 21/1050 2026-08-20 10:37 by Vivilian
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见