 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1124 | 回復: 4 | ||
l1003785517新蟲 (小有名氣)
|
[求助]
Matlab 用多變量線性回歸方法估計動力學參數 已有1人參與
|
|
各位蟲友好！ 我在用多變量線性回歸方法估計動力學參數的時候寫了如下代碼（根據書上例題進行修改的），但是程序運行不了，說是目標函數返回了未定義參數。 有哪位幫忙看一下是哪里寫錯了嗎？我自己一直沒檢查出來哪里錯了，謝謝大家！ function KineticsEst % 動力學參數估計 clear all clc C_DC=[6.31317 6.29412 6.27507 6.26237 6.24713 6.23316 6.221095 6.210935]; C_TT=[0.029302 0.04459 0.059878 0.069433 0.080899 0.091091 0.099372 0.107016]; C_NP=[0.007644 0.011466 0.015288 0.018473 0.022295 0.026117 0.029939 0.032487]; r1=[0.0574 0.0407 0.0241 0.0278 0.0315 0.0268 0.022 0.0172]; % 用多變量線性回歸方法估計動力學參數 R = C_DC/r1; y = R; X = [ones(size(y)) C_DC C_TT C_NP]; b=X\y; % 或b = X\y k1= 1/b(1); k2 = b(2)*k1; k3 = b(3)*k1; k4 = b(4)*k1; % 用lsqnonlin()--求解非線性最小二乘法(非線性數據擬合)問題 beta0 = [k1 k2 k3 k4]; lb = [0 0 0 0]; ub = [+inf +inf +inf +inf]; [beta,resnorm,resid,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFun,beta0,lb,ub,[], C_DC,C_TT,C_NP,r1); ci = nlparci(beta,resid,jacobian); % 模型適定性判別 Ne = length(r1); Np = length(beta); [rho2, F] = rho2_F(k1, r1, resnorm, Ne, Np); % 殘差關于擬合值的殘差圖 rc = RateEqs(beta,C_DC,C_TT,C_NP); plot(rc,resid,'*') xlabel('反應速率擬合值, mol/(min g催化劑)') ylabel('殘差R, mol/(min g催化劑)') refline(0,0) % 參數辨識結果 fprintf('\n\nEstimated Parameters:\n') fprintf('\tk1 = %.4f ± %.4f\n',beta(1),ci(1,2)-beta(1)) fprintf('\tk3 = %.4f ± %.4f\n',beta(2),ci(2,2)-beta(2)) fprintf('\tk4 = %.4f ± %.4f\n',beta(3),ci(3,2)-beta(3)) fprintf('\tk5 = %.4f ± %.4f\n',beta(4),ci(4,2)-beta(4)) fprintf(' 決定性指標ρ^2: %.3f\n',rho2) fprintf(' F比: %.3f\n\n',F) % ------------------------------------------------------------------ function f = ObjFun(beta,C_DC,C_TT,C_NP,r1) rc =RateEqs(beta,C_DC,C_TT,C_NP); f = r1 - rc; % ------------------------------------------------------------------ function rc =RateEqs(beta,C_DC,C_TT,C_NP) % Rate equation rc = beta(1)*C_DC./(1+beta(2)*C_DC+beta(3)*C_TT+beta(4)*C_NP); |
回復此樓» 猜你喜歡
調劑
已經有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經有0人回復
-
物理學I論文潤色/翻譯怎么收費?
已經有65人回復
-
中德博士后2026項目5月30日截止- 國內博士后入站、德國兩年全職
已經有0人回復
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內博后入站德國兩年全職
已經有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經有2人回復
-
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結構與器件課題組。
已經有0人回復
-
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！
已經有1人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人
已經有0人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調劑，10人
已經有0人回復
2樓2018-07-11 19:52:37
3樓2018-07-13 06:46:24
l1003785517
新蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 11.4
- 散金: 30
- 帖子: 59
- 在線: 22.1小時
- 蟲號: 4326948
- 注冊: 2016-01-01
- 專業(yè): 能源化工
4樓2018-07-15 14:41:43
l1003785517
新蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 11.4
- 散金: 30
- 帖子: 59
- 在線: 22.1小時
- 蟲號: 4326948
- 注冊: 2016-01-01
- 專業(yè): 能源化工
5樓2018-07-15 14:42:05
相關版塊跳轉
我要訂閱樓主 l1003785517 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華理，數一英一285求A區(qū)調劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|---|---|---|---|---|
|
|
[考研] 本科新能源科學與工程，一志愿華理能動285求調劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
|
[考研] 283求調劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 化學0703 調劑 306分 一志愿211 +4 | 26要上岸 2026-03-28 | 4/200 |
|
|
|
[考研] 復試調劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研] 346求調劑 一志愿070303有機化學 +3 | 蘿卜燉青菜 2026-03-28 | 3/150 |
|
|
|
[考研] 086000調劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
|
[有機交流]
高溫高壓反應求助
 10+4
|
chibby 2026-03-25 | 4/200 |
|
|
|
[考研] 287求調劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
|
[考研] 085600，材料與化工321分，求調劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
|
[考研] 327求調劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研] 材料與化工304求B區(qū)調劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
|
[考研]
 085602 289分求調劑
+8
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
|
[考研] 347求調劑 +4 | L when 2026-03-25 | 4/200 |
|
|
|
[考研] 086003食品工程求調劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
|
[考研] 生物學學碩求調劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
|
[考研] 調劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
|
[考研] 328求調劑 +4 | LHHL66 2026-03-23 | 4/200 |
|
|
|
[考研]
 一志愿070300浙大化學358分，求調劑！
+4
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
|
[考研] 315分，誠求調劑，材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|