 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 975 | 回復(fù): 3 | |||
| 本帖產(chǎn)生 1 個(gè) 博學(xué)EPI ，點(diǎn)擊這里進(jìn)行查看 | |||
閑人與懶貓木蟲 (小有名氣)
|
[求助]
matlab 非線性擬合 數(shù)據(jù)擬合
|
||
|
我在做一個(gè)動(dòng)力學(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ù)辨識(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; 運(yùn)行結(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ù)此樓» 猜你喜歡
317求調(diào)劑
已經(jīng)有12人回復(fù)
-
一志愿武漢理工材料工程專碩調(diào)劑
已經(jīng)有9人回復(fù)
-
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑
已經(jīng)有7人回復(fù)
-
324求調(diào)劑
已經(jīng)有5人回復(fù)
-
288求調(diào)劑
已經(jīng)有16人回復(fù)
-
317求調(diào)劑
已經(jīng)有9人回復(fù)
-
287求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。
已經(jīng)有4人回復(fù)
-
求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有4人回復(fù)
-
085600材料與化工
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- matlab 非線性擬合
已經(jīng)有16人回復(fù)
- 用matlab，非線性擬合兩個(gè)參數(shù)之間的數(shù)量級(jí)差距比較大，參數(shù)怎么設(shè)置權(quán)重？
已經(jīng)有8人回復(fù)
- 如何用Matlab擬合大量散點(diǎn)的邊界曲線
已經(jīng)有5人回復(fù)
- 急！在用origin進(jìn)行非線性擬合怎樣創(chuàng)建新的擬合函數(shù)？
已經(jīng)有14人回復(fù)
- matlab 多元函數(shù)擬合
已經(jīng)有9人回復(fù)
- 在用nlinfit對(duì)非線性函數(shù)擬合時(shí)遇到的一些問題請(qǐng)教
已經(jīng)有5人回復(fù)
- matlab如何在一個(gè)圖里擬合多條直線
已經(jīng)有5人回復(fù)
- 求大神幫忙擬合一個(gè)非線性方程，求出模型參數(shù)
已經(jīng)有15人回復(fù)
- 在進(jìn)行多元非線性擬合的時(shí)候，用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對(duì)數(shù)據(jù)進(jìn)行最小二乘擬合，若果數(shù)據(jù)為負(fù)數(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ù)
- 求指點(diǎn) 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小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 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)計(jì)(F-Statistic): 3.57120350124504 參數(shù) 最佳估算 ---------- ------------- a 0.000908846493918729 b 15363.6309483451 c 0.613132584682165 d -0.466403903584844 ====== 結(jié)果輸出 ====== No 實(shí)測(cè)值y 計(jì)算值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小時(shí)
- 蟲號(hào): 2111568
- 注冊(cè): 2012-11-07
- 性別: GG
- 專業(yè): 能源化工
3樓2015-01-29 20:16:59
閑人與懶貓
木蟲 (小有名氣)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 1515.2
- 帖子: 275
- 在線: 105.5小時(shí)
- 蟲號(hào): 2111568
- 注冊(cè): 2012-11-07
- 性別: GG
- 專業(yè): 能源化工
4樓2015-01-29 22:38:09
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 閑人與懶貓 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 材料學(xué)碩297已過四六級(jí)求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 13/650 |
|
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +19 | 不知道叫什么！ 2026-03-15 | 21/1050 |
|
|
|
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +6 | 懷瑾握瑜l 2026-03-20 | 6/300 |
|
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！ +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
|
[考研] 321求調(diào)劑 +8 | 何潤采123 2026-03-18 | 10/500 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|