【有獎交流】積極回復本帖子，參與交流，就有機會分得作者 liufeisheng 的 18 個金幣，回帖就立即獲得 1 個金幣，每人有 1 次機會

liufeisheng

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 134.3
帖子: 47
在線: 10.3小時
蟲號: 1379323

[交流] 【求助】參數(shù)擬合；用lsqnonlin擬合氣固催化動力學模型如何得到最優(yōu)初值、結(jié)果？

【求助】參數(shù)擬合；用lsqnonlin擬合氣固催化動力學模型時不同參數(shù)初值得到不同的結(jié)果，請問如何得到最優(yōu)初值→最優(yōu)結(jié)果？
若用其他軟件，有無類似教程（案例）推薦？希望案例詳細直接些，小白不勝感激

學校無1stOpt\Lingo, 有Matlab\Origin Pro\Aspen\Maple

我的Matlab程序：
clc, clear;
a=textread('CH4.txt');    %讀取txt的數(shù)據(jù)
y0=a(:,1);                      %提出因變量反應(yīng)速率r的數(shù)據(jù)
x0=a(:,[2:5]);                   %自變量數(shù)據(jù)矩陣
k0=[5 7.8 4.6 69.4];       %待擬合k(1)-k(4)參數(shù)初值
lb=[-inf -inf -inf -inf];    %k(1)-k(4)下邊界
ub=[+inf +inf +inf +inf];%k(1)-k(4)上邊界
fun=@(k)(k(1).*(a(:,2)).^0.5.*(a(:,3)).^0.5).*(1+k(2).*(a(:,2)).^0.5.*(a(:,3)).^0.5+k(3).*(a(:,2)).^0.5.*(a(:,3)).^(-0.5)+k(4).*a(:,5)).^(-2)-y0;
options = optimoptions(@lsqnonlin,'Algorithm','trust-region-reflective');
[k,resnorm,residual,exitflag,output,lambda,jacobian]=lsqnonlin(fun,k0,lb,ub,options);
ci=nlparci(k,residual,'jacobian',jacobian);

輸出結(jié)果：
k =
0.006998532596973 8.341744654161287 5.000558918360854  78.381545090556088
ci=
1.0e+03 *
  -0.000233641132182 0.000247638197376
  -0.135382312161568 0.152065801469891
  -0.115657322370384 0.125658440207105
  -1.746709223913331 1.903472314094444

CH4.txt 數(shù)據(jù)：
3.37298E-05 0.045483955 0.192267913 0.003064331 0.000940799
3.21878E-05 0.046131365 0.192383275 0.002923178 0.001038222
3.1274E-05 0.046183108 0.192564555 0.002841214 0.001038421
3.14155E-05 0.046821577 0.192744839 0.002852714 0.000778737
3.84628E-05 0.047128941 0.194625771 0.054185236 0.000912819
3.52457E-05 0.046876849 0.194669018 0.053872792 0.000782523
3.51931E-05 0.046017958 0.194717184 0.053899781 0.00101082
3.44983E-05 0.046459181 0.194844229 0.053840691 0.001011045
3.88481E-05 0.047733607 0.195747671 0.105901047 0.001046088
4.38977E-05 0.046653568 0.196054833 0.106435268 0.000915325
4.09879E-05 0.046313497 0.195876834 0.106032215 0.000980822
3.9251E-05 0.047132934 0.19592966 0.105862742 0.001078872
2.8978E-05 0.037387191 0.155197021 0.002659132 0.00105032
2.95116E-05 0.037516347 0.155353051 0.00270676 0.000951828
2.79537E-05 0.037481776 0.155383917 0.002564274 0.000952052
2.86171E-05 0.038352003 0.155459011 0.002624085 0.000886453
2.77704E-05 0.036271974 0.155810267 0.043244686 0.00105194
2.81628E-05 0.037183107 0.155909698 0.043267296 0.000723187
2.61699E-05 0.036920894 0.155947996 0.043074371 0.000920381
2.8119E-05 0.036925417 0.156023334 0.043234911 0.001117742
3.58298E-05 0.037319526 0.156426627 0.085190713 0.001283718
3.20259E-05 0.037363173 0.156461043 0.084808392 0.000822899
3.34526E-05 0.037150974 0.156562126 0.084963969 0.000888858
3.64802E-05 0.038366946 0.156522617 0.085222442 0.000921728
擬合模型：
r_(〖CH〗_4 )=(kp_(〖CO〗_2)^0.5 p_(H_2)^0.5)/(1+〖k_1 p〗_(〖CO〗_2)^0.5 p_(H_2)^0.5+k_2  (p_(〖CO〗_2)^0.5)⁄(p_(H_2)^0.5 )+k_3 p_CO )^2
感謝回答

:hand

S前幾天請教過版主，但我調(diào)了參數(shù)后還是很拉跨，麻煩了各位

回復此樓

» 猜你喜歡

西安電子科學大學杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲能器件方向博士已經(jīng)有13人回復
重慶交通大學光子學微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息已經(jīng)有1人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有112人回復
一志愿鄭大材料學碩298分，求調(diào)劑已經(jīng)有6人回復
尋合作：應(yīng)力腐蝕多尺度模擬已經(jīng)有3人回復
考研交流已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2021-11-17 17:59:34

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liufeisheng

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 134.3
帖子: 47
在線: 10.3小時
蟲號: 1379323

引用回帖:

2樓: Originally posted by xiaoqiuqiu20 at 2021-11-17 18:33:40

感謝大佬回復，還望具體指導啊

贊一下

回復此樓

3樓2021-11-17 18:48:45

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liufeisheng

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 134.3
帖子: 47
在線: 10.3小時
蟲號: 1379323

引用回帖:

4樓: Originally posted by yuekong at 2021-11-17 19:08:00

大佬回復太言簡意賅啊，敢問如何搞定？

贊一下

回復此樓

5樓2021-11-17 19:18:39

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liufeisheng

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 134.3
帖子: 47
在線: 10.3小時
蟲號: 1379323

引用回帖:

6樓: Originally posted by zzyxiang2012 at 2021-11-17 19:24:46
h

這惜字如金啊

發(fā)自小木蟲Android客戶端

回復此樓

7樓2021-11-17 19:53:28

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liufeisheng

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 134.3
帖子: 47
在線: 10.3小時
蟲號: 1379323

大神們，我理解的是，我這個問題就是解非線性超定方程組，要個最優(yōu)解，不知對不對

贊一下

回復此樓

9樓2021-11-17 21:30:28

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

★
liufeisheng(金幣+1): 謝謝參與

a(:,3)列數(shù)據(jù)擬合公式里沒用上吧，1stOpt穩(wěn)定唯一解：

Root of Mean Square Error (RMSE): 3.34612849534681E-6
Sum of Squared Residual: 2.68717821776925E-10
Correlation Coef. (R): 0.694037081019363
R-Square: 0.481687469829878

Parameter Best Estimate
---------- -------------
k1 0.000137780278403927
k2 0.490633000674274
k3 -0.805989362395384
k4 -53.9028500300834

贊一下

回復此樓

11樓2021-11-18 16:13:47

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

chenhuanlong

禁蟲 (文學泰斗)

★
liufeisheng(金幣+1): 謝謝參與

本帖內(nèi)容被屏蔽

12樓2021-11-20 17:07:42

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liufeisheng

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 134.3
帖子: 47
在線: 10.3小時
蟲號: 1379323

引用回帖:

12樓: Originally posted by chenhuanlong at 2021-11-20 17:07:42
gytik

這是什么？

發(fā)自小木蟲Android客戶端

回復此樓

13樓2021-11-20 23:25:59

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

簡單回復

xiaoqiuqiu202樓

2021-11-17 18:33 回復

liufeisheng(金幣+1): 謝謝參與

yuekong4樓

2021-11-17 19:08 回復

liufeisheng(金幣+1): 謝謝參與

發(fā)自小木蟲Android客戶端

zzyxiang20126樓

2021-11-17 19:24 回復

liufeisheng(金幣+1): 謝謝參與

h 發(fā)自小木蟲IOS客戶端

tzynew8樓

2021-11-17 20:06 回復

liufeisheng(金幣+1): 謝謝參與

nono200910樓

2021-11-17 23:53 回復

liufeisheng(金幣+1): 謝謝參與

。發(fā)自小木蟲Android客戶端

miaojiabing14樓

2021-11-23 17:09 回復

liufeisheng(金幣+1): 謝謝參與

jjijjj15樓

2021-11-24 10:04 回復

liufeisheng(金幣+1): 謝謝參與

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 化學 0703求調(diào)劑總分293 一志愿211 +8	土土小蟲 2026-03-03	8/400	2026-03-04 16:45 by xwxstudy
[考研] 085600材料與化工 298 調(diào)劑 +7	小西笑嘻嘻 2026-03-03	7/350	2026-03-04 16:31 by yuningshan
[考研] 301求調(diào)劑 +3	李LJR 2026-03-04	3/150	2026-03-04 16:15 by sslc1985
[考研] 材料328求調(diào)劑 +7	一個蘿卜02 2026-03-03	7/350	2026-03-04 16:02 by zhukairuo
[考研] 武理材料工程302調(diào)劑 +3	Doleres 2026-03-03	6/300	2026-03-04 15:54 by zhukairuo
[考研] 290求調(diào)劑 +3	Bananaiy 2026-03-04	3/150	2026-03-04 15:52 by sslc1985
[考研] 325求調(diào)劑 +5	學家科 2026-03-04	5/250	2026-03-04 14:04 by kakakapanpan
[考研] 291求調(diào)劑 +3	咕嚕咕嚕123123 2026-03-02	4/200	2026-03-04 08:39 by 每天只擺一小會
[基金申請] 請問大家，研究風險與應(yīng)對措施那里，大家都怎么寫呢？ +3	cauasen 2026-03-02	3/150	2026-03-03 19:43 by Tide man
[考研] 298求調(diào)劑 +3	人間唯你是清歡 2026-03-03	4/200	2026-03-03 18:09 by ms629
[考研] 290求調(diào)劑 +9	ErMiao1020 2026-03-02	9/450	2026-03-03 18:03 by linlonghao
[考研] 261求調(diào)劑 +3	陸lh 2026-03-01	3/150	2026-03-02 19:32 by zhukairuo
[考研] 一志愿華南理工大學材料與化工326分，求調(diào)劑 +3	wujinrui1 2026-02-28	3/150	2026-03-02 16:36 by chuocheng
[考研] 材料與化工328求調(diào)劑 +3	。，。，。，。i 2026-03-02	3/150	2026-03-02 13:09 by houyaoxu
[考研] 295求調(diào)劑 +8	19171856320 2026-02-28	8/400	2026-03-02 11:19 by yuchj
[考研] 材料學調(diào)劑 +10	提神豆沙包 2026-02-28	12/600	2026-03-02 09:26 by 李老師！
[考研] 化工299分求調(diào)劑一志愿985落榜 +5	嘻嘻(^ω^) 2026-03-01	5/250	2026-03-01 19:47 by 無際的草原
[考研] 尋找調(diào)劑 +4	LYidhsjabdj 2026-02-28	4/200	2026-03-01 10:56 by sunny81
[考研] 307求調(diào)劑 +4	73372112 2026-02-28	6/300	2026-03-01 00:04 by ll247
[考研] 304求調(diào)劑 +3	52hz~~ 2026-02-28	5/250	2026-03-01 00:00 by 52hz~~