 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 1720 | 回復(fù): 7 | ||
o452794544新蟲 (初入文壇)
|
[求助]
二元非線性擬合，請大神幫忙看下程序哪有問題，謝謝！ 已有1人參與
|
|
兩個自變量x1、x2，一個因變量y。通過繪圖可以看出各個自變量和因變量的關(guān)系分布類似于正態(tài)分布，因此選用二維正太分布模型擬合。 下面是二維正態(tài)分布模型：  http://baike.baidu.com/view/9306579.htm 可詳見 http://baike.baidu.com/view/9306579.htm 百度百科 下面是數(shù)據(jù)： y x1 x2 0 1000 2 0 1200 5 1 1400 15 1 1570 20 2 1720 30 2 1790 40 3 1860 60 4 1880 90 6 1980 172 7 2250 200 7 2530 250 8 2850 350 9 3100 450 10 3300 550 11 3500 650 12 3600 800 13 3800 1100 14 4000 1440 14 4200 1501 14 4500 1563 14 4800 1637 6 6800 2863 下面是函數(shù)代碼： function yy=fun(beta,x) b1 = beta(1); b2 = beta(2); b3 = beta(3); b4 = beta(4); b5 = beta(5); yy=(2*pi*b3*b4*(1-b5^2)^(1/2))^(-1)*exp(-1/2*(1-b5^2)^(-1)*((x(:,1)-b1)^2*b3^(-2)-2*b5*(x(:,1)-b1)*(x(:,2)-b2)*(b3*b4)^(-1)+(x(:,2)-b2)^2*b4^(-2))); 下面是程序運行代碼： clc clear data=load('data.txt'); y=data(:,1); x=data(:,2:3); beta=[1,0.1,3,1,0.5]; [beta_out,r,J,COVB,mse]=nlinfit(x,y,@fun,beta); beta_out; mse; betaci=nlparci(beta_out,r,'Jacobian',J); betaa=[beta_out',betaci]; [yy,delta]=nlpredci(@fun,x,beta_out,r,'Jacobian',J); alpha=0.05; nlintool(x,y,'fun',beta,alpha) r=corrcoef(y,yy); plotregression(y,yy) regstats(y,yy) 運行結(jié)果： Warning: Rank deficient, rank = 0, tol = 0.000000e+00. > In nlinfit>LMfit at 351 In nlinfit at 184 In nihe at 7 Warning: Some columns of the Jacobian are effectively zero at the solution, indicating that the model is insensitive to some of its parameters. That may be because those parameters are not present in the model, or otherwise do not affect the predicted values. It may also be due to numerical underflow in the model function, which can sometimes be avoided by choosing better initial parameter values, or by rescaling or recentering. Parameter estimates may be unreliable. > In nlinfit at 253 In nihe at 7 Warning: Matrix is singular to working precision. > In nlinfit at 274 In nihe at 7 Warning: Matrix is singular to working precision. > In nlparci at 104 In nihe at 10 Warning: Matrix is singular to working precision. > In nlpredci at 181 In nihe at 12 Warning: Matrix is singular to working precision. > In nlintool>PredictionsPlusError at 409 In nlintool at 257 In nihe at 14 Warning: Matrix is singular to working precision. > In nlintool>PredictionsPlusError at 409 In nlintool at 257 In nihe at 14 Warning: Matrix is singular to working precision. > In nlintool>PredictionsPlusError at 409 In nlintool at 271 In nihe at 14 Warning: Matrix is singular to working precision. > In regstats at 193 In regstats at 155 In nihe at 17 Warning: Matrix is singular to working precision. > In regstats at 208 In regstats at 155 In nihe at 17 出現(xiàn)很多錯誤，擬合度很低或者出現(xiàn)NaN。 求大神看是什么問題。 |
回復(fù)此樓» 猜你喜歡
材料專碩英一數(shù)二306
已經(jīng)有5人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有9人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有12人回復(fù)
-
0703化學(xué) 305求調(diào)劑
已經(jīng)有4人回復(fù)
-
0703化學(xué)調(diào)劑，求各位老師收留
已經(jīng)有10人回復(fù)
-
271材料工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有16人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復(fù)
dingd
鐵桿木蟲 (職業(yè)作家)
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
o452794544: 金幣+6, ★★★很有幫助 2016-03-21 16:39:36
感謝參與，應(yīng)助指數(shù) +1
o452794544: 金幣+6, ★★★很有幫助 2016-03-21 16:39:36
|
換這個公式要好的多： y=a*exp(-0.5*(((x1-b)/c)^2))+d*exp(-0.5*(((x2-f)/g)^2)); 均方差(RMSE):0.401981621559676 殘差平方和(SSE):3.55496292957842 相關(guān)系數(shù)(R): 0.996661569977283 相關(guān)系數(shù)之平方(R^2): 0.993334285069583 確定系數(shù)(DC): 0.993333686971469 卡方系數(shù)(Chi-Square): 0.838652138061828 F統(tǒng)計(F-Statistic): 480.02546329649 參數(shù) 最佳估算 ---------- ------------- a -5.90625065613002 b 1204.03546366578 c -526.613682770912 d 14.0704287034134 f 1448.75058061013 g -1080.76690762784 |
2樓2016-03-15 14:50:58
dingd
鐵桿木蟲 (職業(yè)作家)
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
3樓2016-03-15 14:51:47
o452794544
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 14.6
- 帖子: 28
- 在線: 44.5小時
- 蟲號: 2100017
- 注冊: 2012-11-01
- 專業(yè): 安全科學(xué)與工程
4樓2016-03-15 15:26:48
o452794544
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 14.6
- 帖子: 28
- 在線: 44.5小時
- 蟲號: 2100017
- 注冊: 2012-11-01
- 專業(yè): 安全科學(xué)與工程
5樓2016-03-15 15:34:13
dingd
鐵桿木蟲 (職業(yè)作家)
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
6樓2016-03-15 17:09:23
o452794544
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 14.6
- 帖子: 28
- 在線: 44.5小時
- 蟲號: 2100017
- 注冊: 2012-11-01
- 專業(yè): 安全科學(xué)與工程
7樓2016-03-21 16:38:05
科院內(nèi)馬爾
銀蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 958.6
- 散金: 60
- 帖子: 32
- 在線: 35.4小時
- 蟲號: 3485137
- 注冊: 2014-10-19
- 性別: GG
- 專業(yè): 影像醫(yī)學(xué)與生物醫(yī)學(xué)工程其
8樓2016-08-31 21:35:23
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 o452794544 的主題更新
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 271材料工程求調(diào)劑
+5
|
.6lL 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考研] 321求調(diào)劑 +3 | 何潤采123 2026-03-18 | 3/150 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機用戶 2026-03-17 | 7/350 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +6 | 慕寒mio 2026-03-16 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | pupcoco 2026-03-17 | 7/350 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|