24小时热门版块排行榜    

查看: 4310  |  回复: 4
【悬赏金币】回答本帖问题,作者zyugiec将赠送您 40 个金币

zyugiec

金虫 (正式写手)

[求助] 求人帮用软件进行数据拟合 已有1人参与

有谁能帮我分别运行下面2段代码吗?用1stOPt软件撰写的,可能我的电脑配置太差,一运行软件就自动关闭
,非常感谢
Title "FF Xylulose kinetic-50";
Parameter A1=[0,], A3=[0,], E1=[0,], k2=[0,], E3=[0,],k4=[0,];
ConstStr k1=A1*exp(-E1/(8.314*(t2+273.15))), k3=A3*exp(-E3/(8.314*(t2+273.15)));
//InitialODEValue t1=0,x=20,M=0,F=0;
Variable t1,t2,X,M,F;
ODEFunction X'=-(k1+k3)*X;
            M'=k1*X-k4*M;
            F'=k1*X+k4*M-k2*F;

Data;
// 0    120    20    0    0
10    120    18.7515    NAN    0.087
20    120    17.8565    NAN    0.251
40    120    16.518    NAN    0.454
60    120    15.426    NAN    0.7585
90    120    15.356    NAN    0.965
120    120    14.6095    NAN    1.398
Data;
0    130    20    0    0
10    130    16.526    NAN    0.53241
20    130    15.3045    NAN    0.83307
40    130    13.951    NAN    1.31745
60    130    11.94    NAN    1.96916
90    130    10.396    NAN    2.97429
120    130    9.4005    NAN    3.46468
Data;
0    140    20    0    0
10    140    12.8645    NAN    0.793
20    140    11.1065    NAN    1.7735
40    140    9.863    NAN    3.1285
60    140    6.0725    NAN    3.7515
90    140    2.413    NAN    3.783
120    140    2.3695    NAN    4.847
Data;
0    150    20    0    0
10    150    8.6215    NAN    2.86759
20    150    6.054    NAN    4.01073
40    150    4.5935    NAN    4.45106
60    150    3.5335    NAN    5.14721
90    150    2.7405    NAN    5.84066
120    150    1.686    NAN    5.89082
Data;
0    160    20    0    0
10    160    6.825    NAN    4.3936
20    160    3.408    NAN    6.1605
40    160    1.5125    NAN    7.33934
60    160    0.565    NAN    7.56042
90    160    0.196    NAN    7.36719
120    160    0.1695    NAN    7.45462






Title "FF Xylulose kinetic-50";
Parameter A1=[0,], A2=[0,], A3=[0,], A4=[0,], E1=[0,], E2=[0,], E3=[0,],E4=[0,];
ConstStr k1=A1*exp(-E1/(8.314*(t2+273.15))), k2=A2*exp(-E2/(8.314*(t2+273.15))), k3=A3*exp(-E3/(8.314*(t2+273.15))), k4=A4*exp(-E4/(8.314*(t2+273.15)));
InitialODEValue t1=0,x=20,M=0,F=0;
Variable t1,t2,X,M,F;
ODEFunction X'=-(k1+k3)*X;
            M'=k1*X-k4*M;
            F'=k1*X+k4*M-k2*F;

Data;
// 0        120        20        0        0
10        120        18.7515        NAN        0.087
20        120        17.8565        NAN        0.251
40        120        16.518        NAN        0.454
60        120        15.426        NAN        0.7585
90        120        15.356        NAN        0.965
120        120        14.6095        NAN        1.398
0        130        20        0        0
10        130        16.526        NAN        0.53241
20        130        15.3045        NAN        0.83307
40        130        13.951        NAN        1.31745
60        130        11.94        NAN        1.96916
90        130        10.396        NAN        2.97429
120        130        9.4005        NAN        3.46468
0        140        20        0        0
10        140        12.8645        NAN        0.793
20        140        11.1065        NAN        1.7735
40        140        9.863        NAN        3.1285
60        140        6.0725        NAN        3.7515
90        140        2.413        NAN        3.783
120        140        2.3695        NAN        4.847
0        150        20        0        0
10        150        8.6215        NAN        2.86759
20        150        6.054        NAN        4.01073
40        150        4.5935        NAN        4.45106
60        150        3.5335        NAN        5.14721
90        150        2.7405        NAN        5.84066
120        150        1.686        NAN        5.89082
0        160        20        0        0
10        160        6.825        NAN        4.3936
20        160        3.408        NAN        6.1605
40        160        1.5125        NAN        7.33934
60        160        0.565        NAN        7.56042
90        160        0.196        NAN        7.36719
120        160        0.1695        NAN        7.45462
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
可以Matlab实现吗?
QQ:2120156492
2楼2025-07-03 14:37:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyugiec

金虫 (正式写手)

引用回帖:
2楼: Originally posted by hzlhm at 2025-07-03 14:37:06
可以Matlab实现吗?

我还没试过MATLAB
3楼2025-07-03 14:52:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

6参数 需要正版才能运行,试用版最多只能拟合4参数
发自小木虫手机客户端
数值计算
4楼2025-07-06 10:25:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jungfraujoch

木虫 (正式写手)

直接使用MATLAB吧,而且有带界面的工具箱可以手动调整,只需要确保输入正确公式即可。Curve Fitting ToolBox,后续也可以把运行过程导出代码。
Labview、Matlab、C#、Python等,QQ:一00九3八067三,注明:小木虫。
5楼2025-07-15 11:55:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zyugiec 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[基金申请] 申请科研项目纯靠实力能成功吗? 200+5 David1982 2026-10-13 8/400 2026-10-14 21:38 by lhw_GPA4.0
[硕博家园] 中科院海洋所招收海洋材料相关领域博士研究生及研究助理 +3 sujiner 2026-10-08 4/200 2026-10-14 20:40 by qicaiyun
[职场人生] 40+中年人失业后,还能找到好工作吗? +3 xqy2009 2026-10-13 3/150 2026-10-14 20:17 by 43358052
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
信息提示
请填处理意见