古道西風11

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 1245
帖子: 6
在線: 8.6小時
蟲號: 1699256
注冊: 2012-03-18
專業(yè): 生物化工與食品化工

[求助] 求助用matlab擬合三元線性回歸方程及分析

用matlab擬合三元線性回歸方程，我自己弄了一個如下
x1=[23.00 23.00 23.00 23.00 23.00 28.00 28.00 28.00 28.00 28.00]';
x2=[0.30 0.30 0.30 0.30 0.30 0.10 0.10 0.10 0.10 0.10]';
x3=[0.00 0.69 1.09 1.38 1.61 0.00 0.69 1.09 1.38 1.61]';
y=[-0.91 -0.35 -0.05 0.15 0.29 -0.58 -0.04 0.23 0.41 0.54];
x=[ones(10,1) x1 x2 x3];
>> [b,bint,r,rint,stats]=regress(y,x);
??? Error using ==> or
Matrix dimensions must agree.

Error in ==> regress at 70
wasnan = (isnan(y) | any(isnan(X),2));

>> b,bint,stats,rcoplot(r,rint)

b =

      0
-0.0115
-1.7878
0.6751

bint =

      0       0
-0.0168 -0.0062
-2.2782 -1.2973
0.5701 0.7800

stats =

0.9747  134.6812 0.0000 0.0059
問題：
1.我編寫的程序對嗎？
2.輸入程序中，總是會出現(xiàn)warning的那段文字，不知道怎么修改？

回復此樓

» 猜你喜歡

一志愿北京化工大學材料與化工（085600）296求調(diào)劑已經(jīng)有22人回復
322求調(diào)劑：一志愿湖南大學材料與化工（085600），已過六級。已經(jīng)有7人回復
332求調(diào)劑已經(jīng)有5人回復
310求調(diào)劑已經(jīng)有10人回復
26考研-291分-廈門大學（085601）-柔性電子學院材料工程專業(yè)求調(diào)劑已經(jīng)有5人回復
085600 材料與化工 329分求調(diào)劑已經(jīng)有15人回復
化學0703 調(diào)劑 306分一志愿211 已經(jīng)有8人回復
342求調(diào)劑已經(jīng)有4人回復
334分一志愿武理材料求調(diào)劑已經(jīng)有14人回復
材料科學與工程 317求調(diào)劑已經(jīng)有5人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

MATLAB用于三元回歸參數(shù) 已經(jīng)有3人回復
matlab擬合方程參數(shù)時初值的選擇已經(jīng)有15人回復
matlab 求指點動力學方程擬合過程中導數(shù)的獲取已經(jīng)有10人回復
求助SAS 多元非線性回歸擬合編程已經(jīng)有3人回復
求助EXCEL 用LINEST函數(shù)求解三元線性回歸方程已經(jīng)有3人回復
菜鳥求助：關于Origin多元線性回歸擬合y=x1*（a*x2+b*x3+c) 已經(jīng)有10人回復
【求助】使用MATLAB擬合出按照公式的系數(shù) 已經(jīng)有19人回復
【求助】matlab曲面擬合表達式已經(jīng)有10人回復
【求助】使用MATLAB怎么實現(xiàn)擬合力場參數(shù)的程序？【已解決】已經(jīng)有15人回復
【求助】使用Matlab預估動力學方程問題已經(jīng)有13人回復

1樓 2012-07-11 09:18:59

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
應助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

【答案】應助回帖

★
感謝參與，應助指數(shù) +1
臭水溝: 金幣+1, 謝謝交流~~ 2012-07-11 12:11:55

x,y的行數(shù)應該一樣，所以y應該是列向量

贊一下

回復此樓

showmethemoney

2樓2012-07-11 09:50:45

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

古道西風11

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 1245
帖子: 6
在線: 8.6小時
蟲號: 1699256
注冊: 2012-03-18
專業(yè): 生物化工與食品化工

請問y寫成列向量后又出現(xiàn)問題了，這是怎么回事��？？
>> x1=[23.00 23.00 23.00 23.00 23.00 28.00 28.00 28.00 28.00 28.00]';
x2=[0.30 0.30 0.30 0.30 0.30 0.10 0.10 0.10 0.10 0.10]';
x3=[0.00 0.69 1.09 1.38 1.61 0.00 0.69 1.09 1.38 1.61]';
>> y=[-0.91
   -0.35
   -0.05
   0.15
   0.29
   -0.58
   -0.04
   0.23
   0.41
   0.54];
>> x=[ones(10,1) x1 x2 x3];
>> [b,bint,r,rint,stats]=regress(y,x);
Warning: X is rank deficient to within machine precision.
> In regress at 82
>> b,bint,stats,rcoplot(r,rint)

贊一下

回復此樓

3樓2012-07-11 10:17:45

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
應助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

【答案】應助回帖

★ ★
臭水溝: 金幣+2, 謝謝交流~~ 2012-07-11 12:12:08

引用回帖:

3樓: Originally posted by 古道西風11 at 2012-07-11 10:17:45
請問y寫成列向量后又出現(xiàn)問題了，這是怎么回事�。�？
>> x1=';
x2=';
x3=';
>> y=;
>> x=;
>> =regress(y,x);
Warning: X is rank deficient to within machine precision.
> I ...

Most likely you have too many predictor variables, and too few observations. As an alaogy, imagine trying to fit a cubic polynomial regression with only two observations. REGRESS will make a choise about which coefficients to set to zero, but it can't possibly know what you really want.

變量太多，你可以試一下，分別使用x(:,1:2), x(:,1:3), x(:,2:3)算一下，使用3列時總會出現(xiàn)一個b=0。

贊一下(1人)

回復此樓

showmethemoney

4樓2012-07-11 10:38:00

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

古道西風11

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 1245
帖子: 6
在線: 8.6小時
蟲號: 1699256
注冊: 2012-03-18
專業(yè): 生物化工與食品化工

我上面編寫的程序是在別人的二元線性回歸方程的基礎上自己改的，我沒學過matlab，所以，對您的所說的"分別使用x(:,1:2), x(:,1:3), x(:,2:3)算一下，使用3列時總會出現(xiàn)一個b=0"不是十分明白什么意思，我主要就是想依據(jù)上面的數(shù)據(jù)擬合出來一個三元線性回歸方程，如果您愿意的話，請您幫寫一個程序行嗎？萬分感謝！

贊一下

回復此樓

5樓2012-07-11 11:29:51

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
應助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
臭水溝: 金幣+1, 謝謝交流~~ 2012-07-11 12:12:17
古道西風11: 金幣+10, ★★★很有幫助, 解決了問題 2012-07-11 12:33:22

引用回帖:

5樓: Originally posted by 古道西風11 at 2012-07-11 11:29:51
我上面編寫的程序是在別人的二元線性回歸方程的基礎上自己改的，我沒學過matlab，所以，對您的所說的"分別使用x(:,1:2), x(:,1:3), x(:,2:3)算一下，使用3列時總會出現(xiàn)一個b=0"不是十分明白什么意思，我主 ...

上面那樣的就可以了。問題是你的x的第2列和第3列，也就是x1和x2在回歸計算時有問題，所以matlab默認把常數(shù)項b(1)為零。
你可以把x1改成另外的一系列數(shù)試試，或者把x1刪掉，或者再做多組(x1,x2,x3,y)的數(shù)據(jù)用作回歸。

贊一下(1人)

回復此樓

showmethemoney

6樓2012-07-11 11:51:47

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

古道西風11

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 1245
帖子: 6
在線: 8.6小時
蟲號: 1699256
注冊: 2012-03-18
專業(yè): 生物化工與食品化工

謝謝，問題解決！

回復此樓

7樓2012-07-11 12:33:46

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

08小貝

金蟲 (正式寫手)

應助: 4 (幼兒園)
金幣: 162.3
散金: 1377
紅花: 35
帖子: 463
在線: 92.6小時
蟲號: 1202162
注冊: 2011-02-12
性別: MM
專業(yè): 環(huán)境工程

本人菜鳥一枚，請問樓主你這個的回歸剛才怎么寫呀？我在做一個題目，跟你的有點類似。。但是不知道怎么寫回歸方程。。請樓主賜教~~

贊一下

回復此樓

8樓2012-09-28 20:31:01

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhaowl764

禁蟲 (初入文壇)

本帖內(nèi)容被屏蔽

9樓2013-05-14 13:15:22

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

singingp

金蟲 (小有名氣)

應助: 2 (幼兒園)
金幣: 1043.6
散金: 100
紅花: 1
帖子: 124
在線: 58.1小時
蟲號: 2135793
注冊: 2012-11-19
性別: GG
專業(yè): 流體力學

迭代數(shù): 33
計算用時(時:分:秒:微秒): 00:00:00:241
優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法
計算結束原因: 達到收斂判斷標準
均方差(RMSE): 0.0274286872091944
殘差平方和(SSE): 0.00752332882019822
相關系數(shù)(R): 0.998018933040592
相關系數(shù)之平方(R^2): 0.996041790707481
決定系數(shù)(DC): 0.996041790707481
卡方系數(shù)(Chi-Square): -0.0155244265746169
F統(tǒng)計(F-Statistic): 2013.11596653563

參數(shù) 最佳估算
---------- -------------
a       -37.3038021226221
b       -934.025053065556
c       0.724150995543315
d       1137.33012469023

====== 結果輸出 =====

No 實測值y 計算值y
1 -0.91 -0.8648400
2 -0.35 -0.3651759
3 -0.05 -0.0755155
4 0.15 0.1344883
5 0.29 0.3010431
6 -0.58 -0.5788400
7 -0.04 -0.0791759
8 0.23 0.2104845
9 0.41 0.4204883
10 0.54 0.5870431

贊一下

回復此樓

10樓2015-04-28 12:12:09

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 334分一志愿武理材料求調(diào)劑 +14	李李不服輸 2026-03-26	14/700	2026-03-30 16:39 by oooqiao
[考研] 291求調(diào)劑 +7	HanBeiNingZC 2026-03-24	7/350	2026-03-30 16:36 by 梨子～
[考研] 求調(diào)劑323材料與化工 +10	1124361 2026-03-24	10/500	2026-03-30 16:26 by 690616278
[考研] 329求調(diào)劑 +8	星野? 2026-03-26	8/400	2026-03-30 13:41 by chemdavid
[考研] 085602 化學工程專碩 340分求調(diào)劑 +4	qianbai11 2026-03-29	4/200	2026-03-30 11:34 by 唐沐兒
[考研] 332求調(diào)劑 +6	Lyy930824@ 2026-03-29	6/300	2026-03-30 10:53 by 1172367218
[考研] 285求調(diào)劑 +5	AZMK 2026-03-29	8/400	2026-03-30 10:41 by AZMK
[考研] 282求調(diào)劑 +4	wcq131415 2026-03-24	4/200	2026-03-30 10:39 by Delta2012
[考研] 337求調(diào)劑 +6	《樹》 2026-03-29	6/300	2026-03-30 10:15 by herarysara
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導師 +7	傾____蕭 2026-03-27	8/400	2026-03-30 09:37 by longlotian
[考研] 299求調(diào)劑 +10	15188958825 2026-03-25	10/500	2026-03-29 17:51 by 王亮_大連醫(yī)科大
[考研] 289求調(diào)劑 +5	BrightLL 2026-03-29	5/250	2026-03-29 17:24 by zhyzzh
[考研] 349求調(diào)劑 +6	李木子啊哈哈 2026-03-25	6/300	2026-03-29 12:47 by 無際的草原
[考研] 315求調(diào)劑 +4	akie... 2026-03-28	5/250	2026-03-28 21:05 by zhq0425
[考研] 340求調(diào)劑 +5	jhx777 2026-03-27	5/250	2026-03-28 04:18 by fmesaito
[考研] 化學調(diào)劑 +4	愛吃番茄的旭 2026-03-24	5/250	2026-03-27 17:50 by kiokin
[考研] 281求調(diào)劑 +6	Koxui 2026-03-24	7/350	2026-03-26 15:37 by 無際的草原
[考研] 一志愿河工大 081700 276求調(diào)劑 +4	地球繞著太陽轉 2026-03-23	4/200	2026-03-26 14:27 by zzll406
[考研] 材料專碩 335 分求調(diào)劑 +4	拒絕冷暴力 2026-03-25	4/200	2026-03-25 18:45 by haxia
[考研] 一志愿武理085500機械專業(yè)總分300求調(diào)劑 +3	an10101 2026-03-24	7/350	2026-03-25 00:00 by 山鬼0-