24小时热门版块排行榜    

查看: 3541  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 帮忙看看fileCode +4 wwncly 2026-08-10 8/400 2026-08-10 22:07 by demonstreets
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +3 Tide man 2026-08-10 3/150 2026-08-10 21:48 by GZUZM
[基金申请] filecode +8 documentary 2026-08-10 8/400 2026-08-10 19:20 by loufangrui
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 11/550 2026-08-10 17:09 by 六两废铜
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +14 Lanmanbaby 2026-08-09 22/1100 2026-08-10 14:57 by wwncly
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +4 布布和一二 2026-08-10 4/200 2026-08-10 14:50 by 医学老男孩
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见