 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 1190 | 回復(fù): 0 | |||
安之兩隅金蟲 (小有名氣)
|
[求助]
Matlab 非線性擬合問題
|
|
Matlab小白，參照黃華江matlab在化學(xué)工程中的應(yīng)用的教材來改寫模型，但是運(yùn)行的時(shí)候出現(xiàn)了，請(qǐng)大家?guī)兔匆幌�，麻煩了�?br />
function kinetics3 clear all clc %實(shí)驗(yàn)數(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]; % 用多變量線性回歸方法估計(jì)動(dòng)力學(xué)參數(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); % 殘差關(guān)于擬合值的殘差圖 rc = RateEqs(beta,B,C); plot(rc,resid,'*'); xlabel('反應(yīng)速率擬合值'); ylabel('實(shí)驗(yàn)速率'); refline(0,0); % 參數(shù)辨識(shí)結(jié)果 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); 錯(cuò)誤使用 snls (line 47) Objective function is returning undefined values at initial point. lsqnonlin cannot continue. 出錯(cuò) lsqncommon (line 167) snls(funfcn,xC,lb,ub,flags.verbosity,options,defaultopt,initVals.F,initVals.J,caller, ... 出錯(cuò) lsqnonlin (line 253) lsqncommon(funfcn,xCurrent,lb,ub,options,defaultopt,allDefaultOpts,caller,... 出錯(cuò) Kinetics3 (line 23) lsqnonlin(@ObjFun,beta0,lb,ub,[],B,C,r); |
回復(fù)此樓» 猜你喜歡
哈工大材料324求調(diào)劑
已經(jīng)有3人回復(fù)
-
304求調(diào)劑（085602一志愿985）
已經(jīng)有5人回復(fù)
-
環(huán)境調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工，291，求調(diào)劑
已經(jīng)有9人回復(fù)
-
085701環(huán)境工程295求調(diào)劑，一志愿A區(qū)211
已經(jīng)有9人回復(fù)
-
材料與化工 一志愿山大 321分 求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿鄭州大學(xué)，學(xué)碩，物理化學(xué)， 333求調(diào)劑
已經(jīng)有9人回復(fù)
-
一志愿北京化工大學(xué)材料與化工296分求調(diào)劑
已經(jīng)有9人回復(fù)
-
2026博士申請(qǐng)
已經(jīng)有6人回復(fù)
-
265求調(diào)劑
已經(jīng)有7人回復(fù)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 安之兩隅 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 環(huán)境調(diào)劑 +5 | 曉看天暮看云 2026-03-09 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 085701環(huán)境工程295求調(diào)劑，一志愿A區(qū)211 +8 | 等春來， 2026-03-06 | 9/450 |
|
|
|
[考研] 313分生物學(xué)求調(diào)劑 +4 | Yyt楊1 2026-03-09 | 6/300 |
|
|
[考研]
 295復(fù)試調(diào)劑
+4
|
簡(jiǎn)木ChuFront 2026-03-09 | 4/200 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 很老實(shí)人 2026-03-09 | 4/200 |
|
|
|
[考研] 297求調(diào)劑 +3 | 胡達(dá)靈 2026-03-05 | 5/250 |
|
|
|
[考研] 288求調(diào)劑085600材料與化工 +13 | Daunrin 2026-03-07 | 15/750 |
|
|
|
[考研] 材料科學(xué)（0805）338 求調(diào)劑 +7 | xiaokang3286 2026-03-07 | 7/350 |
|
|
|
[考研] 一志愿鄭大071000分?jǐn)?shù)282求調(diào)劑 +3 | 研研顏 2026-03-05 | 7/350 |
|
|
|
[考研] 0307化學(xué)求調(diào)劑 +6 | 0ok0k 2026-03-07 | 6/300 |
|
|
|
[考研] 誠(chéng)求調(diào)劑，323分 有專利/科研/實(shí)習(xí)經(jīng)歷 +6 | 橙子cyx 2026-03-06 | 8/400 |
|
|
|
[考研] 267化工調(diào)劑求助 +7 | 聰少OZ 2026-03-04 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +4 | 52hz~~ 2026-03-05 | 5/250 |
|
|
|
[考研] 275求調(diào)劑 +4 | 大爆炸難民 2026-03-06 | 5/250 |
|
|
|
[考研] 301求調(diào)劑 +4 | 李LJR 2026-03-04 | 4/200 |
|
|
|
[考研] 304分材料專碩求調(diào)劑 +11 | qiuzhigril 2026-03-03 | 14/700 |
|
|
|
[考研] 材料學(xué)碩080500復(fù)試調(diào)劑294 +3 | 四葉zjz 2026-03-04 | 3/150 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑，總分320分，一志愿中南民族大學(xué) +3 | 19286712440 2026-03-04 | 3/150 |
|
|
|
[考研] 一志愿314求調(diào)劑 +7 | 202111120625 2026-03-03 | 7/350 |
|
|
|
[考研] 能動(dòng)297求調(diào)劑，本科川大 +4 | 邵11 2026-03-04 | 4/200 |
|