多元线性回归模型之数学建模实验报告.docx
- 文档编号:4837134
- 上传时间:2022-12-10
- 格式:DOCX
- 页数:13
- 大小:42.42KB
多元线性回归模型之数学建模实验报告.docx
《多元线性回归模型之数学建模实验报告.docx》由会员分享,可在线阅读,更多相关《多元线性回归模型之数学建模实验报告.docx(13页珍藏版)》请在冰豆网上搜索。
多元线性回归模型之数学建模实验报告
%%大学统计与数学模型分析实验中心
《经济数学实验与建模》实验报告
班级:
学号:
姓名:
实验时间11.27实验地点2栋402
实验名称:
多元线性回归模型使用软件:
matlab软件
实
验
目
的
1.作出各个自变量与因变量的散点图,判断是否可以建立多元线性回归模型
2.建立多元线性回归模型
3.根据残差平方和与可决系数等统计量判断可否需要进行模型的改进
4.做残差及其置信区间图,剔除异常值
5.改进模型
实
验
内
容(具体题目及程序)
建立多元线性回归模型
题目:
下表为血压与年龄、体重指数、吸烟习惯,已知体重指数=体重(kg)/身高(m)的平方,
吸烟习惯:
0表示不吸烟,1表示吸烟。
先作散点图,然后判断是否可以建立多元线性回归模型,如果可以,建立血压与年龄、体重指数、吸烟习惯之间的回归模型
序
号
血
压
年
龄
体重
指数
吸烟
习惯
序
号
血
压
年
龄
体重
指数
吸烟
习惯
1
144
39
24.2
0
9
158
56
28.0
1
2
215
47
31.1
1
10
154
59
25.8
0
3
138
45
22.6
0
11
162
34
27.3
1
4
145
67
24.0
1
12
150
42
20.1
0
5
162
42
24.0
1
13
140
48
21.7
0
6
142
67
19.7
0
14
110
45
22.2
0
7
170
56
27.2
1
15
128
18
27.4
0
8
124
64
19.3
0
…
…
…
…
…
程序:
y=[144215138145162142170124158154162150140110128130135114116124136142120120160158144130125175];
x1=[394745476546674267566456593442484518201936503921445363292569];
x2=[24.231.122.624.025.925.129.519.727.219.328.025.827.320.121.722.227.418.822.621.525.026.223.520.327.128.628.322.025.327.4];
x3=[010110101010000100000100110101];
plot(x1,y,’>’)
plot(x2,y,’*’)
x=[ones(30,1),x1',x2',x3']
[b,bint,r,rint,s]=regress(y',x)
b,bint,s
%y=45.3636+0.3604*x1+3.0906*x2+11.8246*x3
rcoplot(r,rint)
x1=x(a,:
)
y1=y(a)
[b1,bint1,r1,rint1,s1]=regress(y1',x1)
b1,bint1,s1
%y=58.5101+0.4303*x1+2.3449*x2+10.3065*x3
rcoplot(r1,rint1)
实
验
结
果
分
析
结果(分析):
x=
1.000039.000024.20000
1.000047.000031.10001.0000
1.000045.000022.60000
1.000047.000024.00001.0000
1.000065.000025.90001.0000
1.000046.000025.10000
1.000067.000029.50001.0000
1.000042.000019.70000
1.000067.000027.20001.0000
1.000056.000019.30000
1.000064.000028.00001.0000
1.000056.000025.80000
1.000059.000027.30000
1.000034.000020.10000
1.000042.000021.70000
1.000048.000022.20001.0000
1.000045.000027.40000
1.000018.000018.80000
1.000020.000022.60000
1.000019.000021.50000
1.000036.000025.00000
1.000050.000026.20001.0000
1.000039.000023.50000
1.000021.000020.30000
1.000044.000027.10001.0000
1.000053.000028.60001.0000
1.000063.000028.30000
1.000029.000022.00001.0000
1.000025.000025.30000
1.000069.000027.40001.0000
b=
45.3636
0.3604
3.0906
11.8246
bint=
3.553787.1736
-0.07580.7965
1.05305.1281
-0.148223.7973
r=
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.5292
-5.6311
-7.5639
8.2656
rint=
-16.370335.9517
30.018459.4982
-19.603032.7499
-28.996022.3987
-24.174226.8599
-23.810528.7839
-27.982522.9747
-22.641127.8754
-32.848118.0569
9.363548.2532
-30.583821.0099
-20.918930.3583
-35.726113.7317
-35.346615.8742
-25.803326.6754
-37.375111.1832
-36.261513.7389
-20.788328.8828
-31.616318.7811
-19.946330.6311
-25.757526.5561
-39.623011.2651
-38.119114.0273
-21.150829.8115
-22.408128.8115
-32.302418.9477
-35.477012.4186
-29.534318.2722
-32.255617.1277
-16.979433.5106
s=
0.685518.89060.0000169.7917
b=
45.3636
0.3604
3.0906
11.8246
bint=
3.553787.1736
-0.07580.7965
1.05305.1281
-0.148223.7973
s=
0.685518.89060.0000169.7917
x1=
1.000039.000024.20000
1.000045.000022.60000
1.000047.000024.00001.0000
1.000065.000025.90001.0000
1.000046.000025.10000
1.000067.000029.50001.0000
1.000042.000019.70000
1.000067.000027.20001.0000
1.000064.000028.00001.0000
1.000056.000025.80000
1.000059.000027.30000
1.000034.000020.10000
1.000042.000021.70000
1.000048.000022.20001.0000
1.000045.000027.40000
1.000018.000018.80000
1.000020.000022.60000
1.000019.000021.50000
1.000036.000025.00000
1.000050.000026.20001.0000
1.000039.000023.50000
1.000021.000020.30000
1.000044.000027.10001.0000
1.000053.000028.60001.0000
1.000063.000028.30000
1.000029.000022.00001.0000
1.000025.000025.30000
1.000069.000027.40001.0000
y1=
Columns1through26
144138145162142170124158162150140110128130135114116124136142120120160158144130
Columns27through28
125175
b1=
58.5101
0.4303
2.3449
10.3065
bint1=
29.906487.1138
0.12730.7332
0.85093.8389
3.387817.2253
r1=
11.9628
7.1331
-0.3168
4.4831
4.8406
3.1809
1.2241
-3.4258
-0.0109
6.8965
-7.9116
-10.2717
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.2447
rint1=
-2.076426.0021
-7.179921.4462
-14.840614.2070
-9.696218.6624
-9.871119.5523
-11.014517.3763
-12.194414.6427
-17.781410.9298
-14.569714.5479
-7.233521.0266
-21.66755.8442
-23.90373.3603
-13.998015.0667
-24.50351.4510
-20.98986.7449
-10.250117.5718
-18.19399.9734
-7.103020.9018
-11.264718.0201
-23.79904.2662
-24.69303.9016
-9.426619.1321
-4.832422.2419
-14.829013.4591
-21.30265.3474
-16.337910.5734
-17.08959.9038
-1.080825.5702
s1=
0.846244.00870.000053.6604
b1=
58.5101
0.4303
2.3449
10.3065
bint1=
29.906487.1138
0.12730.7332
0.85093.8389
3.387817.2253
s1=
0.846244.00870.000053.6604
成绩
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 多元 线性 回归 模型 数学 建模 实验 报告
![提示](https://static.bdocx.com/images/bang_tan.gif)