实验十三回归分析.pdfVIP

实验十三 回归分析 材料系 材93 2009011976 邓陟 一． 实验目的 1. 了解回归分析的基本原理。 2. 根据问题的要求提出模型。 3. 对已确定的模型，使用MATLAB 确定参数。 二． 实验内容 1. 用切削机床加工时，为实时地调整机床需测定道具的磨损速度，每隔一小时测量刀 具的厚度得到以下数据，建立刀具厚度对于切削时间的回归模型，对模型和回归系 数进行检验，并预测7.5h 和15h 后的刀具厚度，用两种方法计算预测区间，解释计 算结果。 时间/h 0 1 2 3 4 5 6 7 8 9 10 刀具厚度/cm 30.6 29.1 28.4 28.1 28.0 27.7 27.5 27.2 27.0 26.8 26.5 记刀具厚度（因变量）为 ，时间（自变量）为 。 y x 编写程序录入数据并作图观察 y=[30.6 29.1 28.4 28.1 28.0 27.7 27.5 27.2 27.0 26.8 26.5]; x=0:10; plot(x,y,+); 时间与刀具厚度散点图 x y   x 从图形直观地看， 与 大致呈线性关系，建立一元线性回归模型 。 y 0 1 检验模型的有效性 编写程序 n=11; X=[ones(n,1),x]; [b,bint,r,rint,s]=regress(y,X); b,bint,s,rcoplot(r,rint) F0=finv(0.95,1,n-2) 得到如下结果 b = 29.5455 -0.3291 bint = 28.9769 30.1140 -0.4252 -0.2330 s = 0.8696 60.0018 0.0000 0.1985 F0 = 5.1174 列表显示结果 回归系数 估计值 置信区间 0 29.5455 [28.9769 30.1140] 1 -0.3291 [-0.4252 -0.2330] 2 2 R 0.8696 F 60.0018 p 0.0001 s 0.1985 检验结果显示： 的置信区间不含零点； ； 。说明模 1 p  F 1,9 ，0.95 5.1174 F   型有效。 残差及其置信区间图 图中第一个点残差的置信区间不包含零点，为异常点，予以剔除。 剔除异常点后重新检验模型的有效性 编写程序 x=x(1,2:11); y=y(1,2:11); n=10; X=[ones(n,1),x];