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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 理论与计算方法 » 用1stopt计算非线性方程组的未知参数
北京石油化工学院2026年研究生招生接收调剂公告
151/2返回列表12下一页
查看: 2246  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

l1003785517

新虫 (小有名气)

[求助] 用1stopt计算非线性方程组的未知参数

各位虫友能否帮忙看看
        代码如下，其中a1和a2都可以写成k=k0.exp^(-Ea/RT)，即反应速率常数，T是变量。我把参数合并是因为会出现过拟合现象，但即使是这样，我算出来还是会出现参数结果不一样，而且我是希望能够求出k0和Ea来的，不知道有没有什么解决的办法，下面的代码里面只放了一个温度下的实验数据。我把另外几个温度的数据也放到代码的后面了，各位虫友看能帮我解决一下这个难题吗，若有可能帮忙分别计算一下a1和a2中的k0和Ea？
        因为后面还要对不同的动力学方程进行计算，所以想问一下这个是不是得用高阶付费的1stopt软件才能做呢？我用的网上注册版1.5的无法得到正确的结果。谢谢！如果不需要的话还望虫友在解答的过程中能发一下原代码，万分感谢，这里我把剩余的金币都奉上！还望相助，谢谢！
//用1stopt计算非线性方程组的未知参数
Title "kinetic values";
Variables x(1:3),y(1:2);
Parameters a(1:5);
ShareModel;
Function y1=a1*x1/(1+a3*x1+a4*x2+a5*x3);
          y2=(a2*x2-a1*x1)/(1+a3*x1+a4*x2+a5*x2);
Data;
//T=275度 x1,x2,x3;y1,y2;
6.31317        0.029302 0.007644 0.0574 0.0345
6.29412        0.044590 0.011466 0.0407 0.0298
6.27507        0.059878 0.015288 0.0241 0.0251
6.26237        0.069433 0.018473 0.0278 0.0203
6.24713        0.080899 0.022295 0.0315 0.0208
6.23316        0.091091 0.026117 0.0268 0.0194
6.22109 0.099372 0.029939 0.0220 0.0160
6.21093        0.107016 0.032487 0.0172 0.0127
//另外几个实验数据
//T=282度
6.14807        0.133133        0.069433        0.2894        0.1089
6.062345        0.186004        0.102557        0.1519        0.0714
6.021705        0.202566        0.126763        0.0144        0.0338
5.991225        0.212758        0.147147        0.0425        -0.0038
5.968365        0.220402        0.162435        0.0705        0.0217
5.937885        0.234416        0.178997        0.0575        0.0297
5.91566        0.244608        0.1911        0.0445        0.0201
5.892165        0.254163        0.205114        0.0316        0.0105
//T=287度
6.0452        0.187278        0.118482        0.4615        0.1918
5.89026        0.277095        0.184093        0.273        0.1397
5.80136        0.326781        0.223587        0.0845        0.0876
5.74675        0.356083        0.249067        0.0829        0.0356
5.69976        0.377741        0.274547        0.0813        0.0243
5.67182        0.386022        0.294294        0.0677        0.0176
5.648325        0.390481        0.313404        0.0541        0.0152
5.614035        0.400673        0.33761        0.0404        0.0129
//T=292度
5.90677        0.220402        0.224224        0.6236        0.1761
5.741035        0.303212        0.307671        0.4216        0.1584
5.559425        0.383474        0.409591        0.2197        0.1408
5.434965        0.437619        0.480298        0.2059        0.1231
5.311775        0.497497        0.543998        0.1921        0.0677
5.26669        0.514696        0.572026        0.1219        0.0332
5.229225        0.527436        0.596869        0.0817        0.0195
5.201285        0.53508        0.617253        0.0585        0.0058@月只蓝@beefly
回复此楼

» 猜你喜欢
1楼 2018-07-03 16:50:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

l1003785517

新虫 (小有名气)
引用回帖:
3楼: Originally posted by 月只蓝 at 2018-07-04 07:36:35
如2楼所说，公式k=k0.exp^(-Ea/RT)与a的关系说明不清楚，而且此公式的写法本身就是错误的，正确的是“k=k0*exp(-Ea/R/T)，作为求助者，像已知的常数R=8.314应该给出。

...

是的，作为一个求助者真的不应该出现这种错误，很抱歉。
我在代码中给出的公式，Function y1=a1*x1/(1+a3*x1+a4*x2+a5*x3);
                                                    y2=(a2*x2-a1*x1)/(1+a3*x1+a4*x2+a5*x2);
     a1=k1*exp^(-Ea1/R/T);
     a2=k2*exp^(-Ea2/R/T);
     这是一个阿伦尼乌斯公式，k1,k2都是指前因子，Ea1，Ea2是活化能，这几个参数都是我希望求出来的，但是若果放到公式进去会出现过拟合，所以我才把它们合并了放在代码里。R是常数8.314，T是反应温度，代码中我给出了4个温度下的实验数据。
    希望能够得到您的帮助？万分感谢！
赞一下(1人) 回复此楼
6楼2018-07-04 12:55:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
独孤神宇: 金币+1 2018-07-05 14:30:10
l1003785517: 金币+20, ★★★很有帮助 2018-07-05 15:28:26
CODE:
Title "kinetic values";
Constant R=8.314;
VarConstant T=[275,282,287,292];
ConstStr a1=k1*exp(-Ea1/R/T),a2=k2*exp(-Ea2/R/T);
Variables x(1:3),y(1:2);
Parameters a(3:5);
ParameterDomain = [0,];
ShareModel;
Function y1=a1*x1/(1+a3*x1+a4*x2+a5*x3);
         y2=(a2*x2-a1*x1)/(1+a3*x1+a4*x2+a5*x2);
Data; //T=275
6.31317        0.029302 0.007644 0.0574 0.0345
6.29412        0.044590 0.011466 0.0407 0.0298
6.27507        0.059878 0.015288 0.0241 0.0251
6.26237        0.069433 0.018473 0.0278 0.0203
6.24713        0.080899 0.022295 0.0315 0.0208
6.23316        0.091091 0.026117 0.0268 0.0194
6.22109 0.099372 0.029939 0.0220 0.0160
6.21093        0.107016 0.032487 0.0172 0.0127
Data; //T=282
6.14807        0.133133        0.069433        0.2894        0.1089
6.062345        0.186004        0.102557        0.1519        0.0714
6.021705        0.202566        0.126763        0.0144        0.0338
5.991225        0.212758        0.147147        0.0425        -0.0038
5.968365        0.220402        0.162435        0.0705        0.0217
5.937885        0.234416        0.178997        0.0575        0.0297
5.91566        0.244608        0.1911        0.0445        0.0201
5.892165        0.254163        0.205114        0.0316        0.0105
Data; //T=287
6.0452        0.187278        0.118482        0.4615        0.1918
5.89026        0.277095        0.184093        0.273        0.1397
5.80136        0.326781        0.223587        0.0845        0.0876
5.74675        0.356083        0.249067        0.0829        0.0356
5.69976        0.377741        0.274547        0.0813        0.0243
5.67182        0.386022        0.294294        0.0677        0.0176
5.648325        0.390481        0.313404        0.0541        0.0152
5.614035        0.400673        0.33761        0.0404        0.0129
Data; //T=292
5.90677        0.220402        0.224224        0.6236        0.1761
5.741035        0.303212        0.307671        0.4216        0.1584
5.559425        0.383474        0.409591        0.2197        0.1408
5.434965        0.437619        0.480298        0.2059        0.1231
5.311775        0.497497        0.543998        0.1921        0.0677
5.26669        0.514696        0.572026        0.1219        0.0332
5.229225        0.527436        0.596869        0.0817        0.0195
5.201285        0.53508        0.617253        0.0585        0.0058

结果不好，自己好好检查下公式和数据：

Root of Mean Square Error (RMSE): 0.0829320200375007
Sum of Squared Residual: 0.440174076640027
Correlation Coef. (R): 0.0189691245411682
R-Square: 0.000359827685858348
Adjusted R-Square: 0.51886698386634
Determination Coef. (DC): -2.88389955816035
F-Statistic: 0.226197729277559

Parameter                  Best Estimate
--------------------        -------------
k1        4.60347494177944E35
ea1        166179.343552407
k2        9.73280643547253E23
ea2        93748.2902053327
a3        65737.9703181681
a4        4.29079677286443E-8
a5        49386056.7270752
赞一下(1人) 回复此楼
11楼2018-07-04 23:05:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
月只蓝: 金币+5, 感谢指导！ 2018-07-04 07:39:00
月只蓝: 金币+5, 感谢指导！ 2018-07-04 07:39:33
最好把公式写清楚：
Function y1=a1*x1/(1+a3*x1+a4*x2+a5*x3);
          y2=(a2*x2-a1*x1)/(1+a3*x1+a4*x2+a5*x2);

与：
k=k0.exp^(-Ea/RT)
都是什么关系？公式“k=k0.exp^(-Ea/RT)”应该是“k=k0.exp^(-Ea/(R*T))”吧？R又是多少？
赞一下 回复此楼
2楼2018-07-03 22:40:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
如2楼所说，公式k=k0.exp^(-Ea/RT)与a的关系说明不清楚，而且此公式的写法本身就是错误的，正确的是“k=k0*exp(-Ea/R/T)，作为求助者，像已知的常数R=8.314应该给出。

发自小木虫Android客户端
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
3楼2018-07-04 07:36:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

l1003785517

新虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2018-07-03 22:40:05
最好把公式写清楚：
Function y1=a1*x1/(1+a3*x1+a4*x2+a5*x3);
          y2=(a2*x2-a1*x1)/(1+a3*x1+a4*x2+a5*x2);

与：
k=k0.exp^(-Ea/RT)
都是什么关系？公式“k=k0.exp^(-Ea/RT)”应该是“k=k0.exp^(-E ...

首先感谢您的回复，对我的描述不清楚表示抱歉。
      代码中的公式，Function y1=a1*x1/(1+a3*x1+a4*x2+a5*x3);
                                             y2=(a2*x2-a1*x1)/(1+a3*x1+a4*x2+a5*x2);
     a1=k1.exp^(-Ea1/RT);
     a2=k2.exp^(-Ea2/RT);
     这是一个阿伦尼乌斯公式，k1,k2都是指前因子，Ea1，Ea2是活化能，这几个参数都是我希望求出来的，但是若果放到公式进去会出现过拟合，所以我才把它们合并了放在代码里。R是常数8.314，T是反应温度，代码中我给出了4个温度下的实验数据。
     不知道您是否能帮我解决这个问题呢？万分感谢！
赞一下 回复此楼
4楼2018-07-04 12:48:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

l1003785517

新虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2018-07-03 22:40:05
最好把公式写清楚：
Function y1=a1*x1/(1+a3*x1+a4*x2+a5*x3);
          y2=(a2*x2-a1*x1)/(1+a3*x1+a4*x2+a5*x2);

与：
k=k0.exp^(-Ea/RT)
都是什么关系？公式“k=k0.exp^(-Ea/RT)”应该是“k=k0.exp^(-E ...

不好意思这两个公式应该是这样，我写错了。
a1=k1*exp^(-Ea1/R/T)；
a2=k2*exp^(-Ea2/R/T);
赞一下 回复此楼
5楼2018-07-04 12:52:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
l1003785517: 金币+20, ★有帮助 2018-07-04 20:21:40
还是说的不明不白：
1：a3,a4和a5与a1，a2采用相同的处理方式？还是视为3个独立的待求参数？
2：对应4个T有四组数据，是想分别求出这四组数据各自对应的一组参数还是这四组数据拥有一组共同的参数？
3：待求参数有范围限制吗？比如k1，k2必须在[0,1]之间。
赞一下 回复此楼
7楼2018-07-04 17:52:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
6楼: Originally posted by l1003785517 at 2018-07-04 12:55:03
是的，作为一个求助者真的不应该出现这种错误，很抱歉。
我在代码中给出的公式，Function y1=a1*x1/(1+a3*x1+a4*x2+a5*x3);
                                                    y2=(a2*x2-a1*x1)/(1+a3*x1+a4 ...

建议你求助于dingd，以我所知，他是做拟合的大神。

发自小木虫Android客户端
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
8楼2018-07-04 18:31:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

l1003785517

新虫 (小有名气)
引用回帖:
7楼: Originally posted by dingd at 2018-07-04 17:52:32
还是说的不明不白：
1：a3,a4和a5与a1，a2采用相同的处理方式？还是视为3个独立的待求参数？
2：对应4个T有四组数据，是想分别求出这四组数据各自对应的一组参数还是这四组数据拥有一组共同的参数？
3：待求参数 ...

额，只有a1和a2才是这样，a3,a4,a5是三个独立的待求参数，我是希望用这四组数据求出同一组参数来，但是我不知道用1stopt该怎么写，范围限制的话就是所有参数都要是非负，麻烦您了！
赞一下 回复此楼
9楼2018-07-04 20:14:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

l1003785517

新虫 (小有名气)
引用回帖:
8楼: Originally posted by 月只蓝 at 2018-07-04 18:31:35
建议你求助于dingd，以我所知，他是做拟合的大神。
...

好的，非常感谢您的帮助！谢谢
赞一下 回复此楼
10楼2018-07-04 20:15:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 l1003785517 的主题更新
151/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂求助（生物与医药） +3 @6952 2026-04-06 3/150 2026-04-06 21:48 by Ecowxq666！
[考研] 材料与化工363求推荐 +11 zh096 2026-04-04 11/550 2026-04-06 19:14 by guanxin1001
[考研] 材料调剂 +12 一样YWY 2026-04-05 13/650 2026-04-06 15:38 by lin-da
[考研] 324求调剂 +3 k可乐 2026-04-05 4/200 2026-04-06 09:54 by 蓝云思雨
[考研] 377求调剂 +6 by.ovo 2026-04-05 6/300 2026-04-05 22:18 by dongzh2009
[考研] 085600，320分求调剂 +7 大馋小子 2026-04-01 8/400 2026-04-05 21:19 by 学员8dgXkO
[考研] 296求调剂 +3 汪！？！ 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 284求调剂 +7 徐同学_001 2026-04-04 13/650 2026-04-05 17:19 by yulian1987
[考研] 306分材料与化工求调剂 +7 黎吧啦啦你很有� 2026-04-03 7/350 2026-04-05 17:18 by Hdyxbekcb
[考研] 考研调剂 +3 mcbbc 2026-04-04 3/150 2026-04-05 10:03 by barlinike
[考研] 一志愿北京化工大学，初试成绩350求调剂 +9 沿岸?贝壳 2026-04-04 14/700 2026-04-05 01:09 by 沿岸?贝壳
[考研] 材料工程085601数二英一335求调剂 +6 双马尾痞老板2 2026-03-31 6/300 2026-04-04 22:29 by hemengdong
[考研] 323求调剂 +8 李佳乐1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
[考研] 359求调剂 +7 hhhhaaaa$ 2026-04-04 7/350 2026-04-04 18:49 by imissbao
[考研] 求材料调剂，一志愿郑州大学289分 +15 硕星赴 2026-04-03 15/750 2026-04-04 01:01 by userper
[考研] 326求调剂 +3 9ahye 2026-04-02 4/200 2026-04-03 08:43 by Jaylen.
[考研] 302求调剂 +9 zyx上岸！ 2026-04-02 9/450 2026-04-02 23:07 by 马儿快快地跑
[考研] 314求调剂 +11 1xiaojun23 2026-03-31 12/600 2026-04-02 12:31 by 1xiaojun23
[考研] 一志愿北交材料工程总分358 +5 cs0106 2026-04-01 7/350 2026-04-01 11:45 by wangjy2002
[考研] 一志愿 南京航空航天大学 ，080500材料科学与工程学硕 +10 @taotao 2026-03-31 11/550 2026-04-01 09:43 by xiayizhi
信息提示
关闭
请填处理意见
关闭