- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
院
专
系:
业:
数学学院
信息与计算科学
年
级:
2014 级
学生姓名: 王继禹
学 号: 201401050335
教师姓名: 徐霞
温度(℃)2025
温度(℃)
20
25
30
35
40
45
50
55
60
65
产量(kg)
13.2
15.1
16.4
17.1
17.9
18.7
19.6
21.2
22.5
24.3
求 y 关于 x 的线性回归方程,检验回归效果是否显著,并预测 x=42℃时产量的估值及预测区间(置信度 95%).
解:
输入数据:
x=[20 25 30 35 40 45 50 55 60 65];
X=[ones(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,bint,r,rint,stats]=regress(Y,X) 得结果:
b =
9.1212
0.2230
bint =
8.0211 10.2214
0.1985 0.2476
stats =
0.9821 439.8311 0.0000 0.2333
即 ??
0
? 9.1212, ??
1
? 0.223 , ??
0
的置信区间为[8.0211,10.2214], ??
1
的置信区间为[0.1985,
0.2476] , R2
? 0.9821, F ? 439.8311, p ? 0.0000 , p0.05, 可 知 回 归 模 型
y ? 9.1212 ? 0.223x 成立。 y 关于 x 的线性回归方程的回归效果是显著的。
残差分析,作残差图:
在(2)输入命令得出结果的基础上,再输入命令: rcoplot(r,rint)
得到残差图 1:
图 1
从残差图图 1 可以看出,所有数据的残差离零点均较近,且残差的置信区间均包含零点,这说明回归模型能较好地符合原始数据。
预测及作图
在(3)的命令基础上,再输入以下命令: z=b(1)+b(2)*x
再输入作图命令: plot(X,Y,k+,X,z,r)
得到各数据点及回归方程的图形如图 2.
图 2
结论:由图 2 可以看出回归直线很好的拟合了所有数据点。(5)计算当 x=42℃时,产量的估值及预测区间:
在(4)的命令基础上,输入以下程序: x=42;
z0=b(1)+b(2)*x 得结果:
z0 = 18.488
所以,当 x=42℃时,产量的估值为 18.488kg 及预测区间为[16.3581,20.6206] (置信度 95%)。
2、某零件上有一段曲线,为了在程序控制机床上加工这一零件,需要求这段曲线的解析表
达式,在曲线横坐标 x
i
处测得纵坐标 y
i
共 11 对数据如下:
x
x
i
0
2
4
6
8
10
12
14
16
18
20
y
i
0.6
2.0
4.4
7.5
11.8 17.1 23.3 31.2 39.6 49.7 61.7
求这段曲线的纵坐标 y 关于横坐标 x 的二次多项式回归方程. 解:
输入数据:
x=[0 2 4 6 8 10 12 14 16 18 20];
y=[0.6 2.0 4.4 7.5 11.8 17.1 23.3 31.2 39.6 49.7 61.7];
作二次多项式回归: [p,s]=polyfit(x,y,2) 得结果:
p =
0.1403 0.1971 1.0105
S =
R: [3x3 double] df: 8
normr: 1.1097
即这段曲线的纵坐标 y 关于横坐标 x 的二次多项式回归方程为
预测及作图
在 matlab 中输入的程序:
x=[0 2 4 6 8 10 12 14 16 18 20];
y=[0.6 2.0 4.4 7.5 11.8 17.1 23.3 31.2 39.6 49.7 61.7];
[p,s]=polyfit(x,y,2) 得出结果再输入: Y=polyconf(p,x,S); 得出结果再输入:
plot(x,y,’k+’,Y,’r’)
得到试验点与回归曲线的图形(图 3)。
图 3
3.某校 60 名学生的一次考试成绩如下:
93
93
75
83
93
91
85
84
82
77
76
77
95
94
89
91
88
86
83
96
81
79
97
78
75
67
69
68
84
83
81
75
66
85
70
94
84
83
82
80
78
74
73
76
70
86
76
90
89
71
66
86
73
80
94
79
78
77
63
53
55
计算均值,标准差,极差,偏度,峰度,画出直方图;
检验分布的正态性;
若检验符合正态分布,估计正态分布的
文档评论(0)