線性回歸分析實驗報告

PAGEPAGE12012——2013學(xué)年第二學(xué)期實驗報告課程名稱：數(shù)學(xué)建模實驗項目：回歸分析實驗類別：綜合性□設(shè)計性eq\o\ac(□,√)驗證性□專業(yè)班級：姓名：學(xué)號：實驗地點：實驗時間：指導(dǎo)教師：成績：一.實驗?zāi)康氖炀氄莆斩嘣€性回歸模型的方法用matlab求解線性回歸問題。實驗內(nèi)容水泥凝固時放出的熱量y與水泥中的四種化學(xué)成分x1,x2,x3,x4有關(guān),今測得一組數(shù)據(jù)如下,試確定多元線性模型.分析：X1=[7,1,11,11,7,11,3,1,2,21,1,11,10];X2=[26,29,56,31,52,55,71,31,54,47,40,66,68];X3=[6,15,8,8,6,9,17,22,18,4,23,9,8];X4=[60,52,20,47,33,22,6,44,22,26,34,12,12];Y=[78.5,74.3,104.3,87.6,95.9,109.2,102.7,72.5,93.1,115.9,83.8,113.3,109.4];由式（9）可得X=[]（eT為單位列向量）Y=yT三．模型求解（含經(jīng)調(diào)試后正確的源程序）x1=[7,1,11,11,7,11,3,1,2,21,1,11,10];x2=[26,29,56,31,52,55,71,31,54,47,40,66,68];x3=[6,15,8,8,6,9,17,22,18,4,23,9,8];x4=[60,52,20,47,33,22,6,44,22,26,34,12,12];y=[78.5,74.3,104.3,87.6,95.9,109.2,102.7,72.5,93.1,115.9,83.8,113.3,109.4];X=[ones(length(y),1),x1',x2',x3',x4'];%把行向量轉(zhuǎn)軼為列向量Y=y';%把行向量轉(zhuǎn)軼為列向量[b,bint,r,rint,stats]=regress(Y,X);b,bint,statsrcoplot(r,rint)z=b(1)+b(2)*x1+b(3)*x2+b(4)*x3+b(5)*x4;plot(X,Y,'k+',X,z,'r')四．結(jié)果分析b=62.40541.55110.51020.1019-0.1441bint=-99.1786223.9893-0.16633.2685-1.15892.1792-1.63851.8423-1.77911.4910stats=0.9824111.47920.00005.8930因此我們可得=62.4054,的置信區(qū)間(-99.1786,223.9893);=1.5511,的置信區(qū)間(-0.1663,3.2685);=0.5102,的置信區(qū)間(-1.1589,2.1792);=0.1019,的置信區(qū)間(-1.6385,1.8423);=-1441,的置信區(qū)間(-1.7791,1.4910).r2=0.9824，F(xiàn)=111.4792，p=0.0000。p<0.05,可知回歸模型成立.五．實驗總結(jié)本次實驗主要是熟悉線性回歸的數(shù)學(xué)模型的建立與求解，通