 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2029 | 回復(fù): 14 | |||
| 【有獎交流】積極回復(fù)本帖子，參與交流，就有機(jī)會分得作者 liufeisheng 的 18 個金幣 ，回帖就立即獲得 1 個金幣，每人有 1 次機(jī)會 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
[交流]
【求助】參數(shù)擬合；用lsqnonlin擬合氣固催化動力學(xué)模型如何得到最優(yōu)初值、結(jié)果？
|
|||
 【求助】參數(shù)擬合；用lsqnonlin擬合氣固催化動力學(xué)模型時不同參數(shù)初值得到不同的結(jié)果，請問如何得到最優(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前幾天請教過版主，但我調(diào)了參數(shù)后還是很拉跨，麻煩了各位 |
回復(fù)此樓» 猜你喜歡
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有2人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有290人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有6人回復(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ù)
-
幫我的英語口語老師找學(xué)生
已經(jīng)有3人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
江西師范大學(xué)陶端健教授團(tuán)隊(duì)2026年碩士研究生招生公告
+1/231
-
中科院強(qiáng)磁場中心高薪誠聘博士后（有機(jī)化學(xué)）
+3/95
-
2026年藥物化學(xué)博士招生
+1/83
-
北京信息科技大學(xué)智能機(jī)器人技術(shù)方向，聯(lián)合清華培養(yǎng)碩生、博生，并長期招博士后
+1/81
-
青島科技大學(xué)山東省流程模擬與優(yōu)化重點(diǎn)實(shí)驗(yàn)室/系統(tǒng)工程團(tuán)隊(duì)招收碩士研究生，誠邀加入
+1/36
-
西安工程大學(xué)能源化工新材料研發(fā)項(xiàng)目組接收2026年研究生
+1/35
-
26屆計算機(jī)、電子信息類、電科、控制、通信考研T j信息pp駿：74+08+76+6+20/
+1/14
-
電子科技大學(xué)材料學(xué)院SFT創(chuàng)新中心招收博士生及科研助理 長期有效
+1/12
-
武漢科技大學(xué)周志輝教授團(tuán)隊(duì)招師資博后（膜分離方向）
+1/8
-
深圳大學(xué)何傳新課題組誠聘博士后、（副）研究員
+1/6
-
齊魯工業(yè)大學(xué)輕工學(xué)部、綠色造紙與資源循環(huán)全國重點(diǎn)實(shí)驗(yàn)室招收一志愿或調(diào)劑研究生
+1/6
-
溫州大學(xué)王娟課題組招收2026年博士生
+1/6
-
歡迎報考中山大學(xué)課題組，提供2026級碩士研究生名額
+1/6
-
香港科技大學(xué) 招生 2026 Fall全獎博士 -- 機(jī)械/電子/材料
+1/5
-
26屆計算機(jī)、電子信息類、電科、控制、通信考研T j信息pp駿：74+08+76+6+20/
+1/5
-
中南林業(yè)科技大學(xué)生物質(zhì)綠色轉(zhuǎn)化與功能材料課題組現(xiàn)招收2026年秋季入學(xué)博士研究生
+1/5
-
歡迎報考中山大學(xué)課題組，提供2026級碩士研究生名額
+1/5
-
徐工-環(huán)境工程學(xué)院-招收調(diào)劑碩士
+1/2
-
【招聘科研助理】四川大學(xué)華西醫(yī)院Rajeev Kumar助理研究員課題組
+1/1
-
深圳大學(xué)何傳新課題組誠聘博士后、（副）研究員
+1/1
3樓2021-11-17 18:48:45
5樓2021-11-17 19:18:39
這惜字如金啊7樓2021-11-17 19:53:28
9樓2021-11-17 21:30:28
簡單回復(fù)
yuekong4樓
2021-11-17 19:08
回復(fù)
nono200910樓
2021-11-17 23:53
回復(fù)
liufeisheng(金幣+1): 謝謝參與
。 發(fā)自小木蟲Android客戶端
2021-11-17 19:24
回復(fù)
liufeisheng(金幣+1): 謝謝參與
h 發(fā)自小木蟲IOS客戶端
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085602高分子方向求調(diào)劑 +3 | tlgudy 2026-03-04 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 0856材料與化工求調(diào)劑！ +5 | 化工考生111 2026-03-04 | 9/450 |
|
|
[考研] 334求調(diào)劑 +4 | Trying] 2026-03-05 | 5/250 |
|
|
|
[考研] 材料328求調(diào)劑 +10 | 一個蘿卜02 2026-03-03 | 10/500 |
|
|
|
[考研] 材料與化工，291分，求調(diào)劑 +9 | 咕嚕咕嚕123123 2026-03-04 | 11/550 |
|
|
|
[考研] 070300化學(xué) 280 一志愿太原理工 求調(diào)劑 +8 | 拾玖壹 2026-03-04 | 8/400 |
|
|
|
[考研] 0703化學(xué) 一志愿東北大學(xué) 279分求調(diào)劑 +4 | 聶樂平 2026-03-03 | 4/200 |
|
|
|
[考研] 0857求調(diào)劑 +4 | CDDD1013 2026-03-04 | 4/200 |
|
|
|
[考博] 26申博 +5 | north， 2026-02-28 | 5/250 |
|
|
|
[考研] 293求調(diào)劑 +4 | 是樂渝哇 2026-03-03 | 4/200 |
|
|
|
[考研] 主題 +3 | realstar2006 2026-02-27 | 3/150 |
|
|
|
[考研] 0703化學(xué)306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|
|
|
[考研] 266材料化工求調(diào)劑 +3 | 哇塞王帥 2026-03-03 | 3/150 |
|
|
|
[考研]
材料270求調(diào)劑
6+6
|
Eiiiio 2026-03-01 | 11/550 |
|
|
|
[考研] 一志愿中科大能動297求調(diào)劑，本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
|
[考研] 268求調(diào)劑 +6 | 好運(yùn)連綿不絕 2026-03-02 | 6/300 |
|
|
|
[考研] 【2026 碩士調(diào)劑】課題組 招收調(diào)劑生 +3 | 考研版棒棒 2026-03-02 | 5/250 |
|
|
|
[考研] 306分材料調(diào)劑 +5 | chuanzhu川燭 2026-03-01 | 6/300 |
|
|
|
[考博] 博士自薦 +4 | kkluvs 2026-02-28 | 5/250 |
|
|
|
[考研] 化工299分求調(diào)劑 一志愿985落榜 +5 | 嘻嘻(*^ω^*) 2026-03-01 | 5/250 |
|