(完整word版)多元线性回归模型之数学建模实验报告

1、(完整word版)多元线性回归模型之数学建模实验报告(完整word版)多元线性回归模型之数学建模实验报告PAGE PAGE 4(完整word版)多元线性回归模型之数学建模实验报告%大学统计与数学模型分析实验中心经济数学实验与建模实验报告班级： 学号： 姓名： 实验时间 11.27 实验地点 2栋402 实验名称：多元线性回归模型 使用软件： matlab软件 实验目的1. 作出各个自变量与因变量的散点图，判断是否可以建立多元线性回归模型2. 建立多元线性回归模型3。 根据残差平方和与可决系数等统计量判断可否需要进行模型的改进4. 做残差及其置信区间图，剔除异常值5. 改进模型实验内容（具体题目

2、及程序）建立多元线性回归模型题目：下表为血压与年龄、体重指数、吸烟习惯， 已知体重指数 = 体重(kg）/身高（m）的平方，吸烟习惯: 0表示不吸烟，1表示吸烟。先作散点图，然后判断是否可以建立多元线性回归模型，如果可以，建立血压与年龄、体重指数、吸烟习惯之间的回归模型序号 血压年龄体重指数吸烟习惯 序号 血压年龄体重指数吸烟习惯11443924.2091585628。0122154731。11101545925.8031384522。60111623427.3141456724.01121504220.1051624224.01131404821。7061426719.70141104522

3、。2071705627.21151281827.4081246419。30 程序：y=144215138145162142170124158154 162150140110128130135114116124 136142120120160158144130125175;x1=39474547654667426756 64565934424845182019 36503921445363292569;x2=24。2 31.1 22.6 24.0 25。9 25。1 29.5 19。7 27.2 19。3 28.0 25.8 27。3 20.1 21。7 22.2 27。4 18.8 22。6

4、21.5 25.0 26。2 23.5 20.3 27.1 28.6 28。3 22。0 25。3 27.4； x3=0 1 0 1 1 0 1 0 1 0 1 0 0 0 0 1 0 0 0 0 0 1 0 0 1 1 0 1 0 1；plot(x1，y，)plot(x2，y，）x=ones（30,1）, x1，x2，x3b,bint，r，rint,s=regress(y,x)b,bint，s %y= 45。3636+0.3604*x1+3.0906*x2+11。8246x3rcoplot（r，rint）x1=x（a,:）y1=y（a)b1,bint1，r1，rint1,s1=regress

5、（y1，x1)b1,bint1,s1y= 58。5101+0。4303x1+2。3449x2+10.3065*x3rcoplot(r1，rint1）实验结果分析结果（分析): x = 1.0000 39。0000 24.2000 0 1。0000 47。0000 31。1000 1。0000 1。0000 45。0000 22。6000 0 1.0000 47.0000 24。0000 1。0000 1。0000 65。0000 25。9000 1.0000 1。0000 46。0000 25.1000 0 1.0000 67。0000 29.5000 1.0000 1。0000 42。000

6、0 19。7000 0 1.0000 67。0000 27。2000 1。0000 1。0000 56。0000 19。3000 0 1.0000 64。0000 28.0000 1。0000 1.0000 56。0000 25.8000 0 1.0000 59。0000 27。3000 0 1.0000 34.0000 20。1000 0 1.0000 42。0000 21。7000 0 1。0000 48。0000 22。2000 1。0000 1.0000 45。0000 27。4000 0 1。0000 18.0000 18.8000 0 1。0000 20.0000 22。6000

7、0 1。0000 19。0000 21。5000 0 1.0000 36。0000 25。0000 0 1。0000 50。0000 26。2000 1。0000 1。0000 39.0000 23。5000 0 1.0000 21.0000 20。3000 0 1。0000 44。0000 27.1000 1.0000 1。0000 53。0000 28.6000 1。0000 1.0000 63.0000 28.3000 0 1.0000 29.0000 22.0000 1。0000 1.0000 25。0000 25。3000 0 1.0000 69。0000 27.4000 1。000

8、0b = 45.3636 0。3604 3.0906 11.8246bint = 3。5537 87。1736 0。0758 0.7965 1.0530 5。1281 -0.1482 23.7973r = 9.7907 44.7583 6。5734 -3.2986 1。3429 2。4867 -2.5039 2。6172 7。3956 28.8084 -4。7870 4。7197 10。9972 -9.7362 0。4360 13。0960 11.2613 4.0473 -6。4176 5.3424 0.3993 -14。1789 12。0459 4。3303 3.2017 6。6774 11

9、。5292 -5.6311 7.5639 8.2656rint = -16。3703 35。9517 30.0184 59.4982 19。6030 32。7499 -28.9960 22.3987 -24。1742 26.8599 -23.8105 28.7839 -27.9825 22。9747 22。6411 27。8754 32。8481 18.0569 9。3635 48.2532 30。5838 21.0099 -20。9189 30.3583 35.7261 13。7317 35。3466 15。8742 25.8033 26.6754 37.3751 11.1832 36.26

10、15 13。7389 20.7883 28.8828 31。6163 18.7811 -19.9463 30.6311 -25.7575 26。5561 39。6230 11.2651 -38.1191 14。0273 21.1508 29。8115 -22。4081 28.8115 32.3024 18。9477 35.4770 12。4186 -29。5343 18。2722 32.2556 17。1277 16。9794 33.5106s = 0。6855 18。8906 0。0000 169。7917b = 45.3636 0。3604 3.0906 11。8246bint = 3。5

11、537 87。1736 -0。0758 0.7965 1。0530 5.1281 0.1482 23.7973s = 0.6855 18.8906 0。0000 169。7917x1 = 1。0000 39。0000 24。2000 0 1.0000 45.0000 22.6000 0 1。0000 47。0000 24.0000 1。0000 1。0000 65.0000 25.9000 1。0000 1.0000 46.0000 25。1000 0 1。0000 67.0000 29.5000 1。0000 1.0000 42。0000 19.7000 0 1.0000 67.0000 2

12、7.2000 1。0000 1。0000 64。0000 28.0000 1.0000 1。0000 56.0000 25。8000 0 1.0000 59。0000 27。3000 0 1.0000 34。0000 20。1000 0 1.0000 42.0000 21。7000 0 1.0000 48。0000 22。2000 1.0000 1。0000 45.0000 27.4000 0 1.0000 18。0000 18。8000 0 1.0000 20。0000 22.6000 0 1.0000 19。0000 21。5000 0 1.0000 36。0000 25.0000 0 1

13、.0000 50。0000 26.2000 1。0000 1。0000 39.0000 23.5000 0 1。0000 21。0000 20。3000 0 1.0000 44.0000 27。1000 1.0000 1.0000 53.0000 28。6000 1.0000 1.0000 63.0000 28。3000 0 1。0000 29.0000 22.0000 1。0000 1.0000 25.0000 25。3000 0 1.0000 69.0000 27.4000 1.0000y1 = Columns 1 through 26 144 138 145 162 142 170 12

14、4 158 162 150 140 110 128 130 135 114 116 124 136 142 120 120 160 158 144 130 Columns 27 through 28 125 175b1 = 58.5101 0。4303 2。3449 10。3065bint1 = 29。9064 87。1138 0.1273 0.7332 0.8509 3。8389 3.3878 17.2253r1 = 11.9628 7.1331 0.3168 4.4831 4.8406 3.1809 1.2241 -3.4258 0。0109 6.8965 -7。9116 10.2717

15、0。5343 11.5262 -7。1224 3。6609 -4。1103 6.8994 3。3777 -9。7664 -10。3957 4.8527 8。7048 -0。6850 7.9776 -2.8823 -3。5929 12。2447rint1 = 2。0764 26。0021 -7。1799 21。4462 -14。8406 14。2070 9。6962 18.6624 9。8711 19。5523 -11.0145 17。3763 -12.1944 14。6427 17。7814 10。9298 -14.5697 14.5479 -7.2335 21.0266 21。6675 5。8442 -23.9037 3.3603 -13。9980 15。0667 24.5035 1。4510 20。9898 6。7449 -10.2501 17。5718 18.1939 9。9734 -7.1030 20。9018 11.2647 18。0201 23.7990 4。2662 24.6930 3.9016 -9.4266 19。1321 -4。8324 22.2419 -14