 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 980 | 回復(fù): 3 | ||
| 本帖產(chǎn)生 1 個 博學(xué)EPI ，點擊這里進行查看 | ||
閑人與懶貓木蟲 (小有名氣)
|
[求助]
matlab 非線性擬合 數(shù)據(jù)擬合
|
|
|
我在做一個動力學(xué)的方程擬合，用的是最小二乘法，下面是程序，不知道哪里有問題，結(jié)果偏差很大，求各位大神幫助！ 原方程 K=k0*exp(-E/T)*PH^a*PS^b*P，知K T PH PS P,求 k0 E a b function nihe global T0 PH0 PS0 P0 K0 T0=[5014.5891 5014.5891 5014.5891 5180.8691 5180.8691 5180.8691 5347.1491 5347.1491 5347.1491]; PH0=[788732.39 592592.59 396039.6 790123.46 394366.19 594059.41 792079.21 591549.29 395061.73]; PS0=[20.19 13.28 7.1 17.71 10.1 10.65 14.2 15.15 8.85]; P0=[24220.91 29411.10 20587.77 20587.77 29411.10 24220.91 29411.10 20587.77 24220.91]; K0=[1047 1122 1233 1188 1225 1148 1815 1185 1125]; y=log(K0); x1=T0 x2=log(PH0); x3=log(PS0); x4=log(P0); y=y'; X=[ones(size(y)) x1' x2' x3' x4']; [b bint]=regress(y,X,0.1); k=exp(b(1)); E=b(2) a=b(3); c=b(4); beta0=[k E a c]; lb=[];ub=[]; [beta,resnorm,resid,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,beta0,lb,ub); ci=nlparci(beta,resid,jacobian) % 殘差關(guān)于擬合值的殘差圖 K0c = Rate(beta,PH0,PS0,T0,P0); plot(K0c,resid,'*') xlabel('反應(yīng)速率擬合值, torr s^-^1') ylabel('殘差R, torr s^-^1') refline(0,0) % 參數(shù)辨識結(jié)果 fprintf('Estimated Parameters:\n') fprintf('\tk = %.4f ± %.4f\n',beta(1),ci(1,2)-beta(1)) fprintf('\tE = %.2f ± %.2f\n',beta(2),ci(2,2)-beta(2)) fprintf('\ta = %.2f ± %.2f\n',beta(3),ci(3,2)-beta(3)) fprintf('\tc = %.2f ± %.2f\n',beta(3),ci(3,2)-beta(3)) fprintf('\tThe sum of the squares is: %.1e\n\n',resnorm) % ---------------------------------------------------------------- function f=ObjFunc(beta) global T0 PH0 PS0 P0 K0 f=K0-Rate(beta,T0,PH0,PS0,P0); % ------------------------------------------------------------------ function K=Rate(beta,T,PH,PS,P) K=beta(1).*exp(-beta(2)./T).*PH.^beta(3).*PS.^beta(4).*P; 運行結(jié)果如下 ci = 1.0e+004 * -0.0000 0.0000 -2.4270 3.5058 -0.0000 0.0002 -0.0002 0.0001 Estimated Parameters: k = 0.0001 ± 0.0015 E = 5393.89 ± 29664.08 a = 0.64 ± 1.05 c = 0.64 ± 1.05 The sum of the squares is: 2.2e+005 |
回復(fù)此樓» 猜你喜歡
北科281學(xué)碩材料求調(diào)劑
已經(jīng)有13人回復(fù)
-
289求調(diào)劑
已經(jīng)有11人回復(fù)
-
網(wǎng)絡(luò)空間安全0839招調(diào)劑
已經(jīng)有5人回復(fù)
-
材料277求調(diào)劑
已經(jīng)有4人回復(fù)
-
274求調(diào)劑
已經(jīng)有7人回復(fù)
-
085600 材料與化工 329分求調(diào)劑
已經(jīng)有3人回復(fù)
-
302求調(diào)劑
已經(jīng)有3人回復(fù)
-
312求調(diào)劑
已經(jīng)有5人回復(fù)
-
329求調(diào)劑
已經(jīng)有3人回復(fù)
-
【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- matlab 非線性擬合
已經(jīng)有16人回復(fù)
- 用matlab，非線性擬合兩個參數(shù)之間的數(shù)量級差距比較大，參數(shù)怎么設(shè)置權(quán)重？
已經(jīng)有8人回復(fù)
- 如何用Matlab擬合大量散點的邊界曲線
已經(jīng)有5人回復(fù)
- 急！在用origin進行非線性擬合怎樣創(chuàng)建新的擬合函數(shù)？
已經(jīng)有14人回復(fù)
- matlab 多元函數(shù)擬合
已經(jīng)有9人回復(fù)
- 在用nlinfit對非線性函數(shù)擬合時遇到的一些問題請教
已經(jīng)有5人回復(fù)
- matlab如何在一個圖里擬合多條直線
已經(jīng)有5人回復(fù)
- 求大神幫忙擬合一個非線性方程，求出模型參數(shù)
已經(jīng)有15人回復(fù)
- 在進行多元非線性擬合的時候，用matlab擬合的值殘差都大于0
已經(jīng)有14人回復(fù)
- 求助用matlab 多元非線性擬合
已經(jīng)有4人回復(fù)
- matlab 擬合方法
已經(jīng)有4人回復(fù)
- 求Origin或Matlab高手非線性擬合公式
已經(jīng)有5人回復(fù)
- 關(guān)于matlab擬合
已經(jīng)有4人回復(fù)
- 用Matlab做多變量非線性擬合求參數(shù)值，求各位大俠鼎力相助！
已經(jīng)有3人回復(fù)
- matlab對數(shù)據(jù)進行最小二乘擬合，若果數(shù)據(jù)為負數(shù)怎么擬合？急急急！
已經(jīng)有3人回復(fù)
- 關(guān)于matlab微分、及非線性擬合問題
已經(jīng)有5人回復(fù)
- matlab數(shù)據(jù)擬合問題求助
已經(jīng)有7人回復(fù)
- 非線性數(shù)據(jù)擬合的數(shù)學(xué)原理�。。。。。�？？？
已經(jīng)有30人回復(fù)
- 用MATLAB或者ORIGIN怎么擬合復(fù)數(shù)數(shù)據(jù)
已經(jīng)有10人回復(fù)
- 求指點 langmuir線性擬合 與 非線性擬合 不一致
已經(jīng)有12人回復(fù)
- 1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？
已經(jīng)有13人回復(fù)
- MATLAB非線性優(yōu)化擬合怎么改才正確
已經(jīng)有3人回復(fù)
- matlab非線性參數(shù)擬合問題
已經(jīng)有7人回復(fù)
dingd
鐵桿木蟲 (職業(yè)作家)
- 博學(xué)EPI: 1
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
閑人與懶貓: 金幣+40, 博學(xué)EPI+1, ★★★很有幫助, 謝謝！ 2015-01-29 20:14:43
閑人與懶貓: 金幣+40, 博學(xué)EPI+1, ★★★很有幫助, 謝謝！ 2015-01-29 20:14:43
|
均方差(RMSE):145.616909851711 殘差平方和(SSE):190838.559912851 相關(guān)系數(shù)(R): 0.808436694770782 相關(guān)系數(shù)之平方(R^2): 0.653569889451906 確定系數(shù)(DC): 0.533303922309213 卡方系數(shù)(Chi-Square): 84.1949588985453 F統(tǒng)計(F-Statistic): 3.57120350124504 參數(shù) 最佳估算 ---------- ------------- a 0.000908846493918729 b 15363.6309483451 c 0.613132584682165 d -0.466403903584844 ====== 結(jié)果輸出 ====== No 實測值y 計算值y 1 1047 1044.65871777826 2 1122 1294.24697276882 3 1233 947.632511671439 4 1188 1042.59182404487 5 1225 1263.95133737708 6 1148 1305.46096405217 7 1815 1813.28448120857 8 1185 1029.71059671882 9 1125 1215.31618898618 |
2樓2015-01-29 20:07:38
閑人與懶貓
木蟲 (小有名氣)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 1515.2
- 帖子: 275
- 在線: 105.5小時
- 蟲號: 2111568
- 注冊: 2012-11-07
- 性別: GG
- 專業(yè): 能源化工
3樓2015-01-29 20:16:59
閑人與懶貓
木蟲 (小有名氣)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 1515.2
- 帖子: 275
- 在線: 105.5小時
- 蟲號: 2111568
- 注冊: 2012-11-07
- 性別: GG
- 專業(yè): 能源化工
4樓2015-01-29 22:38:09
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 閑人與懶貓 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +7 | 顧九笙要謙虛 2026-03-24 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
|
[考研] 359求調(diào)劑 +3 | 王了個楠 2026-03-25 | 3/150 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +17 | 陌の森林 2026-03-18 | 17/850 |
|
|
|
[考研]
材料學(xué)碩，求調(diào)劑
 6+4
|
糖葫蘆888ll 2026-03-22 | 9/450 |
|
|
[考研] 上海電力大學(xué)材料防護與新材料重點實驗室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +4 | 我愛學(xué)電池 2026-03-23 | 4/200 |
|
|
|
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +3 | gdlf9999 2026-03-24 | 3/150 |
|
|
|
[考研] 291求調(diào)劑 +8 | hhhhxn.. 2026-03-23 | 8/400 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達問津 2026-03-23 | 5/250 |
|
|
|
[考研] 291 求調(diào)劑 +4 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 5/250 |
|
|
|
[考研]
 一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
+7
|
2117205181 2026-03-21 | 8/400 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 考研調(diào)劑 +3 | 呼呼？~+123456 2026-03-21 | 3/150 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-20 | 7/350 |
|
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|