数学建模回归分析实验报告.pdfVIP

实验目的 1、直观了解回归分析基本内容。 2 、掌握用数学软件求解回归分析问题。 实验内容 1、回归分析的基本理论。 2、用数学软件求解回归分析问题。 程序设计 1、考察温度x 对产量y 的影响，测得下列10 组数据： 温度（℃） 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 ℃时产量的 估值. 解： (1) 、输入数据： 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]; (2) 、回归分析及检验： [b,bint,r,rint,stats]=regress(Y,X) b,bint,stats (3) 、得结果： b = 9.1212 0.2230 bint = 8.0211 10.2214 0.1985 0.2476 stats = 0.9821 439.8311 0.0000 0.2333 ˆ ˆ ˆ ˆ 即 9.1212 ,  0.2230 ； 的置信区间为[8.0211，10.2214], 的置信区间 0 1  0  1 2 为[0.1985 ，0.2476]; r =0.9821, F=439.8311, p=0.0000 p0.05, 可知回归模型 y=9.1212+0.2230x 成立. （4 ）、残差分析，作残差图： rcoplot(r,rint) 从残差图可以看出，除第二个数据外，其余数据的残差离零点均较近，且残 差的置信区间均包含零点，这说明回归模型 y=9.1212+0.2230x 能较好的符合原 始数据，没有异常点. （5 ）、预测及作图： z=b(1)+b(2)*x plot(x,Y,k+,x,z,r) 预测x=42 ℃时产量的估值.y=18.4872 2 、某零件上有一段曲线，为了在程序控制机床上加工这一零件，需要求这段曲 线的解析表达式，在曲线横坐标xi 处测得纵坐标yi 共11 对数据如下： 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 求这段曲线的纵坐标y 关于横坐标x 的二次多项式回归方程. 解： (1) 、直接作二次多项式回归： t=[0 2 4 6 8 10 12 14 16 18 20] ; s=[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(t,s,2) (2) 、得到结果 p = 0.1403 0.1971 1.0105 S = R: [3x3 double] df: 8 normr: 1.1097 (3) 、得回归模型为 ： y= 0.1403x²+ 0.1971x+ 1.0105 (4) 、预测及作图 Y=polyconf(p,t,S) plot(t,s,k+,t,Y,r) 4 、混凝土的抗压强度随养护时间的延长而增加，现将一批混凝土作成12 个试块， 记录了养护