应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab拟合酶促反应机制
81/1返回列表
查看: 2172  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

晴朗1989

金虫 (小有名气)

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

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
走遍天下
1楼 2013-08-26 23:15:27
已阅   回复此楼   关注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 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭