24小时热门版块排行榜    

查看: 858  |  回复: 4

huxing5066

金虫 (小有名气)

[求助] 求助1 stopt 编写代码,求出参数,模拟酶反应进程

求大神帮忙计算机模拟,实验中得出不同反应时间的各种底物和产物的浓度,列出相关系数方程,具体条件如下:


反应简化的动力学模型,酶反应中包括五步反应(1.1-1.4)和六个反应速率参数(k1、k-1、k2、k3、k4和k5),如下所示:
求助1 stopt 编写代码,求出参数,模拟酶反应进程

其中,参数k分别表示为每步反应步骤的速率常数,S表示为蔗糖的浓度,E表示酶的浓度,F表示果糖的浓度,M表示为麦芽糖的浓度,DP3表示为产物中低聚三糖的浓度,DP4表示为产物中所有低聚四糖的浓度,DP5为产物中所有低聚五糖的浓度,SE为蔗糖-酶中间体的浓度,GE为葡萄糖基-酶结合体中间体。根据质量守恒定律,参与反应的物质的速率可用一组常微分方程表示(2.1-2.9),如下所示:
求助1 stopt 编写代码,求出参数,模拟酶反应进程-1


有的公式是图片,可以看附件文档

具体数据如下:
M:S        单位        mg/ml                               
25:100        F          S                 M             DP3                DP4              DP5
0        0.00                 104.47         27.09        0.00          0.00         0.00
1        13.83         72.81         13.27        12.64         9.60         4.20
2        25.11         53.63         7.40            9.53                 18.34        9.16
3        35.06         36.00         1.90                 5.76                 17.36       16.33
4        40.19         27.06         0.00            3.53                 14.57       16.64
6        47.16         10.40         0.00            1.65            10.68       20.78
8        52.11         2.88             0.00                 2.56            8.95                22.52
10        52.36         1.34            0.00                1.90                 9.66         23.18                                        
0        0.00                 97.69         48.55       0.00                 0.00                0.00
1        16.46         69.27         28.49       22.13          8.79          2.70
2        30.61         46.43         16.58       24.59          19.12     8.42
3        40.16         27.60         11.35       22.73          27.41     11.13
4        46.47         16.94         8.00         17.60          35.03     15.25
6        50.47         3.45                 3.60                14.73          28.67     19.75
8        52.41         0.95            1.40                17.45          27.76     21.05
10        53.02         0.94                 0.00                18.60          24.75      22.72                                        
0        0.00                 93.22         70.91        0.00                  0.00               0.00
1        16.75         71.16         49.66                27.60         6.86        1.33
2        30.18         44.69         33.2           37.52        19.25       3.81
3        39.87         26.63         25.67        38.47         26.78       8.18
4        45.15         14.16         19.82        37.98         33.31       9.94
6        51.76         2.96                 17.39        33.48         38.33       10.78
8        52.76         0.66            16.47        33.41         38.33       14.02
10        53.13         0.00                 14.98        38.12         35.37       15.28                
0        0.00                 96.21         97.45         0.00            0.00                0.00
1        15.15         70.92         72.12         29.12         5.75                0.00
2        29.89         44.40         54.06         43.96         14.90        2.60
3        40.32         26.98         44.45         50.34         25.57         6.23
4        47.18         15.52         38.07         49.22         30.28         7.01
6        52.77         3.25                 32.94         49.23         38.96         8.32
8        51.66         0.56            30.00         48.99         36.02         11.80
10        48.77         0.00                 27.96         47.90         33.17         12.00
                               
0        0.00                 88.86         135.25         0.00            0.00            0.00
1        14.67         60.00         108.17         24.00         3.67           0.00
2        25.95         43.68         90.59         45.42         9.06                 0.44
3        36.55         28.01         81.00         58.00         18.47         1.20
4        43.50         17.49         74.68         64.71         24.37         3.22
6        49.00         4.30                 62.96         72.74         30.58         5.77
8        52.71         0.82                 58.82         74.87         33.33         5.12
10        52.04         0.00                 56.20         72.78         32.49         5.30                                          
0        0.00          100.00         190.00         0.00                 0.00                 0.00
1        12.69         72.71         160.02         29.91         1.99             0.00
2        22.12         46.25         132.49         46.45         5.76                 0.00
3        32.53         30.99         119.82         60.09         11.90         0.00
4        40.00         15.00         108.00         65.00         17.00         0.00
6        46.70         6.67                 94.30         71.23         21.08         0.00
8        45.39         2.82                 90.60         69.51         22.38         0.00
10        45.34         0.79            90.49         69.70         25.16         0.00
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 论文.doc
  • 2016-06-21 10:25:42, 145 K

» 猜你喜欢

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

huxing5066

金虫 (小有名气)

上面数据乱掉了。。。
time(h)        单位        mg/ml                               
0        0.00         104.47         27.09         0.00         0.00         0.00
1        13.83         72.81         13.27         12.64         9.60         4.20
2        25.11         53.63         7.40         9.53         18.34         9.16
3        35.06         36.00         1.90         5.76         17.36         16.33
4        40.19         27.06         0.00         3.53         14.57         16.64
6        47.16         10.40         0.00         1.65         10.68         20.78
8        52.11         2.88         0.00         2.56         8.95         22.52
10        52.36         1.34         0.00         1.90         9.66         23.18
0        0.00         97.69         48.55         0.00         0.00         0.00
1        16.46         69.27         28.49         22.13         8.79         2.70
2        30.61         46.43         16.58         24.59         19.12         8.42
3        40.16         27.60         11.35         22.73         27.41         11.13
4        46.47         16.94         8.00         17.60         35.03         15.25
6        50.47         3.45         3.60         14.73         28.67         19.75
8        52.41         0.95         1.40         17.45         27.76         21.05
10        53.02         0.94         0.00         18.60         24.75         22.72
0        0.00         93.22         70.91         0.00         0.00         0.00
1        16.75         71.16         49.66         27.60         6.86         1.33
2        30.18         44.69         33.25         37.52         19.25         3.81
3        39.87         26.63         25.67         38.47         26.78         8.18
4        45.15         14.16         19.82         37.98         33.31         9.94
6        51.76         2.96         17.39         33.48         38.33         10.78
8        52.76         0.66         16.47         33.41         38.33         14.02
10        53.13         0.00         14.98         38.12         35.37         15.28
0        0.00         96.21         97.45         0.00         0.00         0.00
1        15.15         70.92         72.12         29.12         5.75         0.00
2        29.89         44.40         54.06         43.96         14.90         2.60
3        40.32         26.98         44.45         50.34         25.57         6.23
4        47.18         15.52         38.07         49.22         30.28         7.01
6        52.77         3.25         32.94         49.23         38.96         8.32
8        51.66         0.56         30.00         48.99         36.02         11.80
10        48.77         0.00         27.96         47.90         33.17         12.00
0        0.00         88.86         135.25         0.00         0.00         0.00
1        14.67         60.00         108.17         24.00         3.67         0.00
2        25.95         43.68         90.59         45.42         9.06         0.44
3        36.55         28.01         81.00         58.00         18.47         1.20
4        43.50         17.49         74.68         64.71         24.37         3.22
6        49.00         4.30         62.96         72.74         30.58         5.77
8        52.71         0.82         58.82         74.87         33.33         5.12
10        52.04         0.00         56.20         72.78         32.49         5.30
0                100.00         190.00         0.00         0.00         0.00
1        12.69         72.71         160.02         29.91         1.99         0.00
2        22.12         46.25         132.49         46.45         5.76         0.00
3        32.53         30.99         119.82         60.09         11.90         0.00
4        40.00         15.00         108.00         65.00         17.00         0.00
6        46.70         6.67         94.30         71.23         21.08         0.00
8        45.39         2.82         90.60         69.51         22.38         0.00
10        45.34         0.79         90.49         69.70         25.16         0.00
2楼2016-06-21 10:49:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huxing5066

金虫 (小有名气)

我自己已解决,请版主帮我撤掉哈~~
3楼2016-06-21 17:13:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

何fox

新虫 (著名写手)

4楼2016-06-21 23:27:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

何fox

新虫 (著名写手)

5楼2016-06-22 11:58:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huxing5066 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:40 by jCd0dEvKHShX
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +14 晴天加油 2026-08-26 15/750 2026-08-29 11:06 by biily
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见