【有獎(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)問如何得到最優(yōu)初值→最優(yōu)結(jié)果？
若用其他軟件，有無類似教程（案例）推薦？希望案例詳細(xì)直接些，小白不勝感激

學(xué)校無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)教過版主，但我調(diào)了參數(shù)后還是很拉跨，麻煩了各位

回復(fù)此樓

» 猜你喜歡

幫我的英語口語老師找學(xué)生已經(jīng)有4人回復(fù)
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時(shí)參數(shù)如何選擇？臨近原子數(shù)怎么取值已經(jīng)有7人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有247人回復(fù)
2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有1人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有2人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)

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

查看全部散金貼

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

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

chenhuanlong

禁蟲 (文學(xué)泰斗)

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

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

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

已閱回復(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)意賅啊，敢問如何搞定？

贊一下

回復(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

引用回帖:

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的回帖

簡(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ā)表

[考研] 一志愿北京化工大學(xué)材料與化工296分求調(diào)劑 +15	稻妻小編 2026-03-09	17/850	2026-03-11 22:45 by 2020015
[考研] 311求調(diào)劑 +6	zchqwer 2026-03-10	6/300	2026-03-11 15:54 by 30660438
[考博] 26申博求助 +3	跳躍餅干 2026-03-10	4/200	2026-03-10 21:15 by Tntcnn
[考研] 0703化學(xué)調(diào)劑 +3	三dd. 2026-03-10	3/150	2026-03-10 15:45 by peike
[考研] 0860求調(diào)劑（272分） +3	lllllcsjsj 2026-03-05	4/200	2026-03-10 15:29 by circleffyy
[考研] 0856材料與化工309分求調(diào)劑 +4	ZyZy…… 2026-03-10	4/200	2026-03-10 15:20 by houyaoxu
[考研] 材料與化工求調(diào)劑 +11	與冬清寧 2026-03-07	14/700	2026-03-10 14:23 by Linda Hu
[考研] 研究生招生 +3	徐海濤11 2026-03-10	6/300	2026-03-10 12:43 by 徐海濤11
[考研] 一志愿武理314求調(diào)劑 +4	（￣～￣;） 2026-03-08	5/250	2026-03-08 23:12 by L135790
[考研] 269求調(diào)劑 +3	朔朔話 2026-03-08	4/200	2026-03-08 20:39 by 熱情沙漠
[考研] 334求調(diào)劑 +8	Trying] 2026-03-06	8/400	2026-03-08 14:51 by hmn_wj
[考研] 085701環(huán)境工程317分求調(diào)劑 +9	6汆尼9 2026-03-07	9/450	2026-03-08 06:41 by 劉兵
[考研] 085701環(huán)境工程專業(yè)，初試305，均過國(guó)家A區(qū)線 +7	卡卡來了@ 2026-03-07	8/400	2026-03-08 06:40 by 劉兵
[考研] 材料化工求調(diào)劑 +7	博斯特525 2026-03-06	8/400	2026-03-08 03:51 by lfhuang
[考研] 第一志愿上海大學(xué)，專業(yè)化學(xué)工程與技術(shù)，總分288，求調(diào)劑 +3	1829197082 2026-03-07	3/150	2026-03-07 19:14 by houyaoxu
[考研] 2026調(diào)劑】考試A區(qū)0703化學(xué)類323分誠(chéng)求接收 +3	卷柏卷柏 2026-03-05	4/200	2026-03-07 17:00 by yuandd_2001
[考研] 276求調(diào)劑本科天津工業(yè)大學(xué) 一志愿北京化工大學(xué)0855 +4	路lyh123 2026-03-06	6/300	2026-03-07 11:37 by shdliugang
[考研] 不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過來 +5	啊擺啊擺 2026-03-05	9/450	2026-03-06 12:06 by 啊擺啊擺
[考研] 287求調(diào)劑 +3	看看我. 2026-03-05	6/300	2026-03-06 10:40 by Iveryant
[考研] 材料調(diào)劑 +4	L9370 2026-03-05	4/200	2026-03-05 19:36 by caszguilin