【有獎(jiǎng)交流】積極回復(fù)本帖子，參與交流，就有機(jī)會(huì)分得作者 liufeisheng 的 18 個(gè)金幣，回帖就立即獲得 1 個(gè)金幣，每人有 1 次機(jī)會(huì)

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

liufeisheng

新蟲 (初入文壇)

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

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

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

學(xué)校無(wú)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前幾天請(qǐng)教過(guò)版主，但我調(diào)了參數(shù)后還是很拉跨，麻煩了各位

回復(fù)此樓

» 猜你喜歡

考研交流已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有177人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
幫我的英語(yǔ)口語(yǔ)老師找學(xué)生已經(jīng)有3人回復(fù)
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時(shí)參數(shù)如何選擇？臨近原子數(shù)怎么取值已經(jīng)有5人回復(fù)
2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)

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

查看全部散金貼

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

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

liufeisheng

新蟲 (初入文壇)

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

引用回帖:

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

這惜字如金啊

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

回復(fù)此樓

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

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

查看全部 15 個(gè)回答

liufeisheng

新蟲 (初入文壇)

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

引用回帖:

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

感謝大佬回復(fù)，還望具體指導(dǎo)啊

贊一下

回復(fù)此樓

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

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

liufeisheng

新蟲 (初入文壇)

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

引用回帖:

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

大佬回復(fù)太言簡(jiǎn)意賅啊，敢問(wèn)如何搞定？

贊一下

回復(fù)此樓

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

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

liufeisheng

新蟲 (初入文壇)

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

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

贊一下

回復(fù)此樓

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

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

簡(jiǎn)單回復(fù)

yuekong4樓

2021-11-17 19:08 回復(fù)

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

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

nono200910樓

2021-11-17 23:53 回復(fù)

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

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

zzyxiang20126樓

2021-11-17 19:24 回復(fù)

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

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

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

[考研] 材料專碩323求調(diào)劑 +3	李白26 2026-03-07	3/150	2026-03-07 14:24 by 楊楊楊紫
[教師之家] 售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】 +3	vs9ky3szrj 2026-03-06	3/150	2026-03-07 12:24 by 68vlfw6g5q
[考研] 售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】 +3	ipys00zm58 2026-03-06	4/200	2026-03-07 03:57 by vwchyxxzm4
[考研] 266材料化工求調(diào)劑 +5	哇塞王帥 2026-03-03	5/250	2026-03-06 16:44 by 學(xué)員FnSWZj
[考研] 復(fù)試調(diào)劑 +7	呼呼？~+123456 2026-03-05	10/500	2026-03-06 16:07 by @颯颯颯颯
[考研] 085600材料調(diào)劑總分330 +6	池池丶 2026-03-03	6/300	2026-03-06 15:50 by @颯颯颯颯
[考研] 化學(xué)專業(yè)調(diào)劑 +5	好好好1233 2026-03-04	6/300	2026-03-06 15:38 by @颯颯颯颯
[考研] 289求調(diào)劑 +3	馳光絕景 2026-03-04	4/200	2026-03-06 12:34 by 馳光絕景
[考研] 304求調(diào)劑 +3	52hz~~ 2026-03-05	4/200	2026-03-06 00:27 by wutongshun
[考研] 一志愿西交化工專碩288專業(yè)課93求調(diào)劑 +6	好運(yùn)好運(yùn)接接 2026-03-04	6/300	2026-03-05 09:36 by houyaoxu
[考研] 281電子信息求調(diào)劑 +5	jhtfeybgj 2026-03-02	9/450	2026-03-05 09:36 by 085400調(diào)劑黨
[考研] 材料328求調(diào)劑 +10	一個(gè)蘿卜02 2026-03-03	10/500	2026-03-05 09:31 by oxidpl
[考研] 320材料與化工，求調(diào)劑 +6	鶴遨予卿 2026-03-04	8/400	2026-03-04 20:47 by wutongshun
[考研] 本科太原理工采礦工程，求調(diào)劑 +3	onlx 2026-03-01	3/150	2026-03-04 15:57 by Stephen_ym
[考研] 281求調(diào)劑 +3	Y?l?h 2026-03-04	3/150	2026-03-04 10:51 by 汪��？！
[考研] 化工專碩調(diào)劑 +4	利好利好. 2026-03-03	7/350	2026-03-03 21:30 by L135790
[考研] 278求調(diào)劑 +3	滿天星11_22 2026-03-02	3/150	2026-03-03 13:51 by Iveryant
[考研] 299求調(diào)劑 +5	kkcoco25 2026-03-02	9/450	2026-03-03 12:55 by 公瑾逍遙
[考研] 調(diào)劑 +3	13853210211 2026-03-02	4/200	2026-03-02 10:16 by 13853210211
[考研] 291分工科求調(diào)劑 +9	science餓餓 2026-03-01	10/500	2026-03-01 18:55 by 18137688336