 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 1722 | 回復(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ù)此樓» 猜你喜歡
0856調(diào)劑，是學(xué)校就去
已經(jīng)有7人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有14人回復(fù)
-
申博26年
已經(jīng)有3人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有10人回復(fù)
-
0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷
已經(jīng)有12人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有13人回復(fù)
-
288求調(diào)劑
已經(jīng)有14人回復(fù)
-
307求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑
已經(jīng)有5人回復(fù)
-
288求調(diào)劑，一志愿華南理工大學(xué)071005
已經(jīng)有5人回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +7 | Ymlll 2026-03-18 | 10/500 |
|
|---|---|---|---|---|
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +8 | zlingli 2026-03-13 | 8/400 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-18 | 3/150 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 一志愿福大288有機化學(xué)，求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|