(完整word版)管理決策實(shí)驗(yàn)二(答案)

1、魯東大學(xué)信息與電氣工程學(xué)院 管理決策模型與方法實(shí)驗(yàn)報(bào)告（二） 實(shí)驗(yàn)名稱(chēng)：MATLAB歸分析命令及其應(yīng)用 專(zhuān) 業(yè)： _ 班 級(jí)： _ 姓 名： _ 學(xué) 號(hào)： _ 一、 實(shí)驗(yàn)?zāi)康模赫莆誱atlab統(tǒng)計(jì)工具箱中的回歸分析命令及其應(yīng)用。 二、 實(shí)驗(yàn)要求：將相關(guān)程序（M文件）和運(yùn)行結(jié)果粘貼到 word文檔中，說(shuō)明 操作過(guò)程中出現(xiàn)的問(wèn)題及改正方法，以及后續(xù)學(xué)習(xí)應(yīng)注意的問(wèn)題。 三、 實(shí)驗(yàn)內(nèi)容： 1 考察溫度x對(duì)產(chǎn)量y的影響，測(cè)得下列10組數(shù)據(jù)： 溫度C） 20 25 30 35 40 45 50 55 60 65 產(chǎn)量（kg） 13.2 15.1 16.4 17.1 17.9 18.7 19.6 21.2

2、22.5 24.3 求y關(guān)于x的線性回歸方程，檢驗(yàn)回歸效果是否顯著，并預(yù)測(cè)x=42C時(shí)產(chǎn)量的估 值及預(yù)測(cè)區(qū)間(置信度95%). 主程序： x=20 25 30 35 40 45 50 55 60 65; X=o nes(10,1) x; y=13.2 15.1 16.4 17.1 17.9 18.7 19.6 21.2 22.5 24.3; b,b in t,r,ri nt,stats=regress(y,X); b, stats rstool(x,y, li near ) P,S=polyfit(x,y,1); Y,DELTA=polyco nf(P,x,S); plot(x,y, k+ ,

3、x,Y, r) Y2,DELTA2=polyco nf(P,42,S) 運(yùn)行結(jié)果： b = 9.1212 0.2230 stats =0.9821 439.8311 0.0000 0.2333 Y2 =18.4885 DELTA2 = 1.1681 425 X1 L L i i 24 - - 22 - - 20 - - 18 - - 16 - - 14 - - f r r r r r r t Predicted 13.& 0.45609 25 30 35 40 45 50 55 60 線性預(yù)測(cè) 預(yù)測(cè)值與真實(shí)值比較26 1220 1 1 i L 十 4- + + - 4- + _ + 1

4、! f r f r r r 65 24 22 20 18 16 14 25 30 35 40 45 50 55 60 (1) y關(guān)于x線性回歸方程為y=9.1212+0.223x,由于相關(guān)系數(shù)R2=0.9821，接近 1,并且與F對(duì)應(yīng)的概率p=0.00000.05所以回歸方程成立，回歸效果顯著。 (2) 預(yù)測(cè) x =42E 時(shí)產(chǎn)量為 18.4885 ( kg) (3) 預(yù)測(cè)x =42E時(shí)預(yù)測(cè)區(qū)間(置信度 95%)為17.3204, 19.6566. 2某零件上有一段曲線，為了在程序控制機(jī)床上加工這一零件，需要求這段曲 線的解析表達(dá)式， 在曲線橫坐標(biāo) Xi處測(cè)得縱坐標(biāo)yi 共11對(duì)數(shù)據(jù)如下： x

5、i 0 2 4 6 8 10 12 14 16 18 20 yi 0.6 2.0 4.4 7.5 11.8 17.1 23.3 31.2 39.6 49.7 61.7 求這段曲線的縱坐標(biāo)y關(guān)于橫坐標(biāo)x的二次多項(xiàng)式回歸方程. 方法一，利用n次多項(xiàng)式函數(shù)polyfit(xi,yi,n)，得到多項(xiàng)式系數(shù) xi=0 2 4 6 8 10 12 14 16 18 20; yi=0.6 2.0 4.4 7.5 11.8 17.1 23.3 31.2 39.6 49.7 61.7; plot(xi,yi) p=polyfit(xi,yi,2) %根據(jù)plot圖，猜想是y關(guān)于x的二次幕方程 運(yùn)行結(jié)果： p =

6、 0.1403 0.1971 1.0105 即曲線擬合方程為 y=0.1403xZ0.1971x+1.0105 方法二，利用多元線性回歸命令regress(y,X)得到多項(xiàng)式系數(shù) xi=0 2 4 6 8 10 12 14 16 18 20; yi=0.6 2.0 4.4 7.5 11.8 17.1 23.3 31.2 39.6 49.7 61.7; x=xi,y=yi, n=len gth(xi); X=x.A2 x on es( n,1); b=regress(y,X) 0.1403 0.1971 1.0105 即曲線擬合方程為 y=0.1403*x0.1971*x+1.0105 3.在研

7、究化學(xué)動(dòng)力學(xué)反應(yīng)過(guò)程中，建立了一個(gè)反應(yīng)速度和反應(yīng)物含量的數(shù)學(xué)模 型，形式為 X3 y 臥飛 y 1 心2冷；3X2 -4X3 其中是未知參數(shù)，X1，X2，X3是三種反應(yīng)物（氫，n戊烷，異構(gòu)戊烷）的含 量，y是反應(yīng)速度.今測(cè)得一組數(shù)據(jù)如下表，試由此確定參數(shù) -,，一：5，并給出置 信區(qū)間.亠：5的參考值為（1,0.05, 0.02, 0.1,2 序號(hào) 反應(yīng)速度y 氫X1 n戊烷X2 異構(gòu)戊烷X3 1 8.55 470 300 10 2 3.79 285 80 10 3 4.82 470 300 120 4 0.02 470 80 120 5 2.75 470 80 10 6 14.39 100

8、190 10 7 2.54 100 80 65 8 4.35 470 190 65 9 13.00 100 300 54 10 8.50 100 300 120 11 0.05 100 80 120 12 11.32 285 300 10 13 3.13 285 190 120 (1) 對(duì)回歸模型建立M文件model.m如下： fun cti on y= model( betaO,x ) b1=betaO(1); b2=beta 0( 2); b3=beta0(3); b4=beta 0； b5=beta0(5); y=(b1.*x(:,2)-(1/b5).*x(:,3).*(1+b2.*x(:,1)+b3.*x(:,2)+b4.*x(:,3).A(-1); end (2) 主程序： x1=470 285 470 470 470 100 100 470 100 100 100 285 285; x2=300 80 300 80 80 190 80 190 300 300 80 300 190; x3=10 10 120 120 10 10 65 65 54 120 120 10 120; x=x1 x2 x3; y=8.55 3.79 4.82 0.02 2.75 14.39 2.54 4.35 13.00 8.50 0.05 11.32