| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1069 | 回復(fù): 3 | |||||
| 【懸賞金幣】回答本帖問(wèn)題,作者kk1424將贈(zèng)送您 50 個(gè)金幣 | |||||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | |||||
kk1424主管區(qū)長(zhǎng) (著名寫(xiě)手)
中藥研究者
|
[求助]
求高手幫忙為小弟做一個(gè)四元回歸曲線(xiàn),高分懸賞! 已有1人參與
|
||||
|
遇到棘手的問(wèn)題,我有MATLAB軟件,但不是很會(huì)。求高手達(dá)人,幫忙給個(gè)編程。。。。感激不盡. 這回歸方程: y1=0.227220-0.002859x1+0.018825x2-0.003369x3-0.021391x4 拜托了! |

榮譽(yù)版主 (文壇精英)
![]() |
專(zhuān)家經(jīng)驗(yàn): +518 |
|
Matlab多元線(xiàn)性回歸 [ b , bint , r , rint , stats ]=regress ( y , x ) , 其中b 是回歸方程中的參數(shù)估計(jì)值,bint 是b 的置信區(qū)間,r 和rint 分別表示殘差及殘差對(duì)應(yīng)的置信區(qū)間。StatS 數(shù)組包含三個(gè)數(shù)字,分別是相關(guān)系數(shù),F(xiàn) 統(tǒng)計(jì)量及對(duì)應(yīng)的概率p 值。 擬合結(jié)果: Y=b(1)x(1)+b(2)x(2)+b(3)x(3)+…+b(n)x(n) b(1)是系數(shù),x(1)為全1的一個(gè)列向量。 注意:不是插值。 x=[1097 1284 1502 1394 1303 1555 1917 2051 2111 2286 2311 2003 2435 2625 2948 3155 3372];%因變量時(shí)間序列數(shù)據(jù) y=[698 872 988 807 738 1025 1316 1539 1561 1765 1762 1960 1902 2013 2446 2736 2825];%自變量時(shí)間序列數(shù)據(jù) X=[ones(size(x')),x']; [b,bint,r,rint,stats]=regress(y',X,0.05);%調(diào)用一元回歸分析函數(shù) rcoplot(r,rint)%畫(huà)出在置信度區(qū)間下誤差分布. 舉例: x = 1 2 4 9 1 4 3 7 1 5 9 0 1 9 1 8 >> y=[10 3 90 48]'; >> [ b , bint , r , rint , stats ]=regress ( y , x ) 得到的結(jié)果 b = -186.8333 16.0238 21.8571 8.5952 bint = NaN NaN NaN NaN NaN NaN NaN NaN r = 1.0e-013 * -0.5684 -0.9948 -0.1421 0.1421 rint = NaN NaN NaN NaN NaN NaN NaN NaN stats = 1 NaN NaN NaN 另外,用b=inv(x)*y 得到的結(jié)果和上面用regression得到的一樣?梢(jiàn),求逆的問(wèn)題也是用了最小二乘的原理。 但是,regress更優(yōu)于inv,體現(xiàn)在當(dāng)x的元素存在缺陷時(shí)。舉例: x = 1 0 0 4 1 0 0 3 1 0 0 5 1 0 0 8 y = 3 4 3 9 [ b , bint , r , rint , stats ]=regress ( y , x ) Warning: X is rank deficient to within machine precision. b = -0.9643 0 0 1.1429 bint = -12.0010 10.0724 0 0 0 0 -0.9245 3.2102 r = -0.6071 1.5357 -1.7500 0.8214 rint = -9.2316 8.0173 -1.9141 4.9855 -7.4997 3.9997 0.2465 1.3964 stats = 0.7388 5.6575 0.1405 3.2321 |
榮譽(yù)版主 (文壇精英)
![]() |
專(zhuān)家經(jīng)驗(yàn): +518 |
榮譽(yù)版主 (文壇精英)
![]() |
專(zhuān)家經(jīng)驗(yàn): +518 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|---|---|---|---|---|
|
[文學(xué)芳草園] 伙伴們,祝我生日快樂(lè)吧 +17 | myrtle 2026-03-10 | 26/1300 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +7 | Ncdx123456 2026-03-13 | 8/400 |
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
[考研] 0856專(zhuān)碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
[考研] 279求調(diào)劑 +3 | 抓著星星的女孩 2026-03-10 | 3/150 |
|
|
[考研] 材料專(zhuān)碩288分求調(diào)劑 一志愿211 +4 | 在家想你 2026-03-11 | 4/200 |
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7 | mysdl 2026-03-11 | 9/450 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[考研] 【考研調(diào)劑求收留】 +3 | Ceciilia 2026-03-11 | 3/150 |
|
|
[考研] 材料專(zhuān)碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[考研] 270求調(diào)劑 085600材料與化工專(zhuān)碩 +3 | YXCT 2026-03-11 | 3/150 |
|