 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1127 | 回復: 4 | |||
l1003785517新蟲 (小有名氣)
|
[求助]
Matlab 用多變量線性回歸方法估計動力學參數(shù) 已有1人參與
|
|
各位蟲友好！ 我在用多變量線性回歸方法估計動力學參數(shù)的時候?qū)懥巳缦麓a（根據(jù)書上例題進行修改的），但是程序運行不了，說是目標函數(shù)返回了未定義參數(shù)。 有哪位幫忙看一下是哪里寫錯了嗎？我自己一直沒檢查出來哪里錯了，謝謝大家！ function KineticsEst % 動力學參數(shù)估計 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]; % 用多變量線性回歸方法估計動力學參數(shù) 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()--求解非線性最小二乘法(非線性數(shù)據(jù)擬合)問題 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) % 參數(shù)辨識結果 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); |
回復此樓» 猜你喜歡
調(diào)劑
已經(jīng)有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有239人回復
-
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復
-
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結構與器件課題組。
已經(jīng)有0人回復
-
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！
已經(jīng)有1人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有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
相關版塊跳轉(zhuǎn)
我要訂閱樓主 l1003785517 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 295求調(diào)劑 +4 | wei-5 2026-03-26 | 4/200 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學學碩304分，有大創(chuàng)項目，四級已過 +7 | 幸運哩哩 2026-03-22 | 11/550 |
|
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4 | 小張zxy張 2026-03-26 | 8/400 |
|
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
|
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +10 | 幸運的醬醬 2026-03-22 | 12/600 |
|
|
|
[考研] 333求調(diào)劑 +3 | question挽風 2026-03-23 | 3/150 |
|
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實習經(jīng)歷 +10 | 材料學257求調(diào)劑 2026-03-23 | 11/550 |
|
|
|
[考研]
材料學碩，求調(diào)劑
 6+5
|
糖葫蘆888ll 2026-03-22 | 10/500 |
|
|
|
[考研] 342求調(diào)劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
|
[考研] 349求調(diào)劑 +4 | 李木子啊哈哈 2026-03-25 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
|
[考研] 303求調(diào)劑 +6 | 藍山月 2026-03-25 | 6/300 |
|
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
|
[考研] 環(huán)境學碩288求調(diào)劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
|
[考研] 315分，誠求調(diào)劑，材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|