 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 1184 | 回復: 0 | ||
安之兩隅金蟲 (小有名氣)
|
[求助]
Matlab 非線性擬合問題
|
|
Matlab小白，參照黃華江matlab在化學工程中的應用的教材來改寫模型，但是運行的時候出現(xiàn)了，請大家?guī)兔匆幌�，麻煩了�?br />
function kinetics3 clear all clc %實驗數(shù)據(jù) B = [1.319 1.3161 1.3112 1.2892 1.2695 1.2601 1.2485 1.2389 1.2245 1.2052 1.181 1.1456 1.1402 1.1134 1.0827 1.036 0.9634]; C= [0.0813 0.0957 0.1188 0.1532 0.1877 0.2089 0.2333 0.2621 0.2891 0.3323 0.3738 0.4275 0.4464 0.4916 0.5393 0.6141 0.7114]; r = [0.02286 0.02244 0.02202 0.0216 0.02118 0.02076 0.02034 0.01992 0.0195 0.01908 0.01866 0.01796 0.01726 0.01656 0.01586 0.01376 0.01166]; % 用多變量線性回歸方法估計動力學參數(shù) R =B./r; y = R; X = [ones(size(y)) B C]; b = X\y % 或[b,bint] = regress(y,X,0.05); KZ = 1/b(1); KB = b(2)*KZ; KC = b(3)*KZ; % 用lsqnonlin()--求解非線性最小二乘法(非線性數(shù)據(jù)擬合)問題 beta0 = [KZ KB KC]; lb = [1 1 1]; ub = [+inf +inf +inf]; [beta,resnorm,resid,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFun,beta0,lb,ub,[],B,C,r); ci = nlparci(beta,resid,jacobian); % 模型適定性判別 Ne = length(r); Np = length(beta); [rho2, F] = rho2_F(KZ, r, resnorm, Ne, Np); % 殘差關于擬合值的殘差圖 rc = RateEqs(beta,B,C); plot(rc,resid,'*'); xlabel('反應速率擬合值'); ylabel('實驗速率'); refline(0,0); % 參數(shù)辨識結果 fprintf('\n\nEstimated Parameters:\n') fprintf('\tKZ= %.4f ± %.4f\n',beta(1),ci(1,2)-beta(1)) fprintf('\tKB = %.2f ± %.2f\n',beta(2),ci(2,2)-beta(2)) fprintf('\tKC = %.2f ± %.2f\n',beta(3),ci(3,2)-beta(3)) % ------------------------------------------------------------------ function f = ObjFun(beta,B,C,r) rc = RateEqs(beta,B,C); f = r - rc; % ------------------------------------------------------------------ function rc = RateEqs(beta,B,C) % Rate equation rc = beta(1)*B./(1+beta(2)*B+beta(3)*C); 錯誤使用 snls (line 47) Objective function is returning undefined values at initial point. lsqnonlin cannot continue. 出錯 lsqncommon (line 167) snls(funfcn,xC,lb,ub,flags.verbosity,options,defaultopt,initVals.F,initVals.J,caller, ... 出錯 lsqnonlin (line 253) lsqncommon(funfcn,xCurrent,lb,ub,options,defaultopt,allDefaultOpts,caller,... 出錯 Kinetics3 (line 23) lsqnonlin(@ObjFun,beta0,lb,ub,[],B,C,r); |
回復此樓» 猜你喜歡
085600 材料與化工 298
已經有8人回復
-
266求調劑
已經有5人回復
-
0854總分272
已經有6人回復
-
EST拒稿重投
已經有3人回復
-
歡迎采礦、地質、巖土、計算機、人工智能等專業(yè)的同學報考
已經有11人回復
-
0703化學學碩 288分
已經有4人回復
-
接收調劑
已經有23人回復
-
化學 0703求調劑 總分293 一志愿211
已經有3人回復
-
304分材料專碩求調劑
已經有5人回復
-
0703化學 學碩 理工科均可 不區(qū)分研究方向 總分279求調劑
已經有5人回復
相關版塊跳轉
我要訂閱樓主 安之兩隅 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 266求調劑 +5 | 哇塞王帥 2026-03-03 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 304分材料專碩求調劑 +3 | qiuzhigril 2026-03-03 | 5/250 |
|
|
[考研] 理學，工學，農學調劑，少走彎路，這里歡迎您！ +8 | likeihood 2026-03-02 | 11/550 |
|
|
|
[考研] 085700資環(huán)求調劑，初始279，六級已過，英語能力強 +3 | 085700資環(huán)調劑 2026-03-03 | 4/200 |
|
|
|
[考研] 264求調劑 +6 | 26調劑 2026-03-03 | 6/300 |
|
|
|
[考研] 085600材料與化工 298 調劑 +3 | 小西笑嘻嘻 2026-03-03 | 3/150 |
|
|
|
[考研] 289求調劑 +7 | BrightLL 2026-03-02 | 9/450 |
|
|
|
[考研] 0805總分292，求調劑 +12 | 幻想之殤 2026-03-01 | 12/600 |
|
|
|
[考研] 278求調劑 +3 | 滿天星11_22 2026-03-02 | 3/150 |
|
|
|
[考研] 26考研報考西工大材料308分求調劑 +4 | weizhong123 2026-03-01 | 5/250 |
|
|
|
[考研] 材料工程求調劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
|
[論文投稿]
通訊作者寫誰，問題是你意想不到的問題
 15+3
|
阿爾法啊 2026-03-01 | 3/150 |
|
|
|
[考研] 282求調劑 +4 | 2103240126 2026-03-02 | 7/350 |
|
|
|
[考研] 298求調劑 +10 | 人間唯你是清歡 2026-02-28 | 14/700 |
|
|
|
[考研] 材料化工調劑 +13 | 今夏不夏 2026-03-01 | 16/800 |
|
|
|
[考研] 289求調劑 +8 | yang婷 2026-03-02 | 9/450 |
|
|
|
[考研] 275求調劑 +3 | L-xin? 2026-03-01 | 6/300 |
|
|
|
[基金申請]
 本子寫完了，給DS兄弟看了，得了92分
+3
|
Doma 2026-03-01 | 7/350 |
|
|
|
[考研] 299求調劑 +3 | Y墨明棋妙Y 2026-02-28 | 5/250 |
|
|
|
[考研] 化工299分求調劑 一志愿985落榜 +5 | 嘻嘻(*^ω^*) 2026-03-01 | 5/250 |
|