【有獎(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í)
蟲(chóng)號(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ú)類(lèi)似教程（案例）推薦？希望案例詳細(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ù)此樓

» 猜你喜歡

重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息已經(jīng)有2人回復(fù)
一志愿鄭大材料學(xué)碩298分，求調(diào)劑已經(jīng)有6人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有164人回復(fù)
尋合作：應(yīng)力腐蝕多尺度模擬已經(jīng)有5人回復(fù)
考研交流已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
幫我的英語(yǔ)口語(yǔ)老師找學(xué)生已經(jīng)有1人回復(fù)

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

查看全部散金貼

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

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

liufeisheng

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 134.3
帖子: 47
在線: 10.3小時(shí)
蟲(chóng)號(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的回帖

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

liufeisheng

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 134.3
帖子: 47
在線: 10.3小時(shí)
蟲(chóng)號(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

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 134.3
帖子: 47
在線: 10.3小時(shí)
蟲(chóng)號(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

新蟲(chóng) (初入文壇)

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

引用回帖:

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

這惜字如金啊

發(fā)自小木蟲(chóng)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ā)自小木蟲(chóng)Android客戶端

nono200910樓

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

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

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

zzyxiang20126樓

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

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

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

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

[考研] 302材料工程求調(diào)劑 +6	Doleres 2026-03-01	7/350	2026-03-04 20:09 by caszguilin
[考研] 環(huán)境工程專碩307 一志愿211 四六級(jí)已過(guò) 求調(diào)劑 +5	ccc! 2026-03-03	5/250	2026-03-04 17:39 by ZYNG
[考研] 一志愿中科大080500總分324求調(diào)劑 +3	jorna 2026-03-03	6/300	2026-03-04 16:20 by 花開(kāi)富貴幸福人?/a>
[考研] 武理材料工程302調(diào)劑 +3	Doleres 2026-03-03	6/300	2026-03-04 15:54 by zhukairuo
[考研] 歡迎采礦、地質(zhì)、巖土、計(jì)算機(jī)、人工智能等專業(yè)的同學(xué)報(bào)考 +9	pin8023 2026-02-28	12/600	2026-03-04 15:27 by rorr
[考研] 281求調(diào)劑 +3	我是小小蔥蔥 2026-03-03	5/250	2026-03-04 14:23 by kakakapanpan
[考研] 347求調(diào)劑 +5	啊歐歐歐 2026-03-03	5/250	2026-03-04 14:10 by kakakapanpan
[考研] 0703化學(xué)學(xué)碩 288分 +4	旅行中的紫葡萄 2026-03-03	4/200	2026-03-03 23:41 by ch1979
[考研] 0703 總分319求調(diào)劑 +4	Xinyuu 2026-03-02	4/200	2026-03-03 22:58 by 讓我們一起快樂(lè)
[考研] 環(huán)境調(diào)劑 +8	chenhanheng 2026-03-02	8/400	2026-03-03 22:13 by 劉兵
[考研] 298求調(diào)劑一志愿中海洋 +3	lour. 2026-03-03	3/150	2026-03-03 20:41 by gxg2025
[考研] 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
[碩博家園] 2025屆雙非化工碩士畢業(yè)，申博 +4	更多的是 2026-02-27	5/250	2026-03-03 17:44 by wenium
[考研] 清華大學(xué) 材料與化工 353分求調(diào)劑 +5	awaystay 2026-03-02	6/300	2026-03-03 09:03 by houyaoxu
[考研] 化工京區(qū)271求調(diào)劑 +7	11ing 2026-03-02	7/350	2026-03-03 07:30 by 利好利好.
[考研] 306分材料調(diào)劑 +5	chuanzhu川燭 2026-03-01	6/300	2026-03-02 19:51 by 張曉芳0105
[考研] 一志愿中石油（華東）本科齊魯工業(yè)大學(xué) +3	石能偉 2026-03-02	3/150	2026-03-02 18:54 by caszguilin
[考研] 303求調(diào)劑 +5	今夏不夏 2026-03-01	5/250	2026-03-02 15:01 by 向上的胖東
[考研] 尋找調(diào)劑 +4	LYidhsjabdj 2026-02-28	4/200	2026-03-01 10:56 by sunny81