| 2 | 1/1 | 返回列表 |
| 查看: 2610 | 回復(fù): 1 | |||
qlpeter銀蟲 (初入文壇)
|
[求助]
matlab多元非線性擬合 已有1人參與
|
|
希望蟲友能給與幫助 采用多元非線性回歸: data.txt 168.7 202.9 145.9 12.1 142.0 206.4 152.9 15.8 181.9 205.8 144.2 10.7 119.7 206.8 156.2 20.9 243.9 203.9 130.6 5.0 170.1 205.8 141.9 17.2 167.3 200.1 138.7 16.3 269.2 203.2 132.0 9.7 254.1 205.4 128.9 10.6 325.2 208.8 123.4 6.7 83.9 210.0 153.4 17.1 26.7 206.4 176.6 92.9 -1.3 207.1 190.8 149.8 79.3 209.2 157.2 21.4 68.2 210.6 167.9 44.7 41.1 205.2 162.9 61.3 8.4 198.8 176.1 115.7 50.1 209.8 167.3 64.9 189.7 216.8 160.2 20.9 50.7 200.4 158.2 46.4 35.9 207.1 167.9 64.0 function main clc clear data=load('data.txt'); y=data(:,1); x=data(:,2:4); beta=; % 回歸系數(shù)的初值,任意取 yy=fun(beta,x); =nlinfit(x,y,@fun,beta); beta_out mse betaci=nlparci(beta_out,r,'Jacobian',J); betaa= =nlpredci(@fun,x,beta_out,r,'Jacobian',J); alpha=0.05; nlintool(x,y,'fun',beta,alpha) r=corrcoef(y,yy) % 求相關(guān)系數(shù)r plotregression(y,yy) % 并且畫圖r regstats(y,yy) function yy=fun(beta,x) b1 = beta(1); b2 = beta(2); b3 = beta(3); b4 = beta(4); yy=b1*(x(:,1).^b2).*(x(:,2).^b3).*(x(:,3).^b4); 運行時顯示??? Error: File: D:\新建文件夾\main.m Line: 9 Column: 1 The expression to the left of the equals sign is not a valid target for an assignment. |
鐵桿木蟲 (職業(yè)作家)
|
均方差(RMSE):29.1845118412956 殘差平方和(SSE):17886.4503597092 相關(guān)系數(shù)(R): 0.950674501500641 相關(guān)系數(shù)之平方(R^2): 0.903782007803491 確定系數(shù)(DC): 0.899149439045897 卡方系數(shù)(Chi-Square): 116.489341933952 F統(tǒng)計(F-Statistic): 56.1887471230383 參數(shù) 最佳估算 ---------- ------------- b1 0.0690396368659353 b2 4.68095500291962 b3 -3.29925868731259 b4 -0.310074339479526 ====== 結(jié)果輸出 ====== No 實測值y 計算值y 1 168.7 145.834326238406 2 142 124.609811365089 3 181.9 168.291349599024 4 119.7 107.453774499923 5 243.9 282.878567758331 6 170.1 153.172181912177 7 167.3 147.231225521178 8 269.2 218.821981736768 9 254.1 242.150731537058 10 325.2 348.111406472481 11 83.9 130.430380974856 12 26.7 44.7214275898647 13 -1.3 30.3560270325071 14 79.3 110.242995026258 15 68.2 72.8394588615087 16 41.1 64.6176535703849 17 8.4 35.3797172907226 18 50.1 64.4980548296721 19 189.7 123.308050070583 20 50.7 69.4516535737008 21 35.9 60.2511608251563 |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0817化工學(xué)碩調(diào)劑 +7 | 努力上岸中! 2026-03-31 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-30 | 5/250 |
|
|
[考研] 085404 22408 315分 +3 | zhuangyan123 2026-03-31 | 3/150 |
|
|
[有機交流]
10+3
|
kaobao456 2026-03-29 | 4/200 |
|
|
[考研] 生物考研337分求調(diào)劑 +3 | cgxin 2026-03-30 | 4/200 |
|
|
[考研] 311求調(diào)劑 +7 | 藍月亮亮 2026-03-30 | 7/350 |
|
|
[考研] 085601材料工程找調(diào)劑 +17 | oatmealR 2026-03-29 | 18/900 |
|
|
[考研] 286求調(diào)劑 +5 | Faune 2026-03-30 | 5/250 |
|
|
[考研] 071010 323 分求調(diào)劑 +3 | Baekzhy 2026-03-27 | 3/150 |
|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分,英語已過四六級,六級521,求調(diào)劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|