实验十一 回归分析

1、辽宁工程技术大学上机实验报告成绩实验名称回归分析院系专业班级姓名学号日期实验 目的简述本次实验目的：1、了解回归分析基本内容2、掌握用matlab软件求解回归分析问题。实验准备你为本次实验做了哪些准备：复习书中相关的MATLAB函数知识，及相关的概率论的知识实验进度本次共有 6 个练习，完成6个。实验总结本次实验的收获、体会、经验、问题和教训：在现实生活中存在着很多的相关关系，我们熟悉也只是其中的关 系，这节课的练习，让我学会了利用matlab对数据进行处理，首先在 图像中画出散点图，从而确定相关的回归模型来，求出他们之间存在 的某种函数关系，来了解函数模型与回归模型之间的关系。以此利用数学知

2、识来说明相关的实际问题，我想这也是我们在日 后的生活有很大帮助的。教师评语1、考察温度x对产量j的影响，测得下列10组数据：温度（C）20253035404550556065产量（kg）13.215.116.417.117.918.719.621.222.524.3求j关于x的线性回归方程，检验回归效果是否显著，并预测x=42C时产量的估值及预测区间(置 信度95%). x=20 25 30 35 40 45 50 55 60 65;Y=13.2 15.1 16.4 17.1 17.9 18.7 19.6 21.2 22.5 24.3;s.：L-图2.豚rd时Id岫i蜘 rstool(x,Y,

3、linear)图1Variables have been created in the current workspace. beta, rmse beta =9.12120.2230rmse =0.4830结论：由图2知x=42C时产量的估值18.4885. j关于x的线性回归方程：y=9.1212+ 0.2230 x剩余标准差为0.4830，说明回归模型显著且显著性较好。2、某零件上有一段曲线，为了在程序控制机床上加工这一零件，需要求这段曲线的解析表达式，在 曲线横坐标xi处测得纵坐标ji共11对数据如下:x. i10 12 14 16 18 20;11.8 17.1 23.3 31.2

4、39.6 49.7 61.7;1.01050.62.04.47.511.817.123.331.239.649.761.7今测得一组数据如下表，试由此确定参数邕,.,p5,并给出置信区间., %的参考值为(1, 0.05, 0.02, 0.1, 2).序号反应速度j氢x1n戊烷x2异构戊烷x318.554703001023.79285801034.8247030012040.024708012052.754708010614.391001901072.54100806584.3547019065913.0010030054108.50100300120110.05100801201211.32

5、28530010133将拟合的非线性模型，建立m文件dongli.m如下function y=dongli(beta,x)y=(beta(1)*x(:,2)-x(:,3)./beta(5)./(1+beta(2)*x(:,1)+beta(3)*x(:,2)+beta(4 )*x(:,3);输入数据及求回归系数和置信区间(yy土delta)clearclcclose ally=8.55 3.79 4.82 0.02 2.75 14.39 2.54 4.35 13.00 8.50 0.05 11.32 3.13;x1=470 285 470 470 470 100 100

6、 470 100 100 100 285 285;x2=300 80 300 80 80 190 80 190 300 300 80 300 190;x3=10 10 120 120 10 10 65 65 54 120 120 10 120;x=x1 x2 x3;beta0=1, 0.05, 0.02, 0.1, 2;beta,r,J=nlinfit(x,y,dongli,beta0);betayy,delta=nlpredci(dongli,x,beta,r ,J);yydelta得出结果：beta =1.25260.06280.04000.11241.1914yy =8.41793.95

7、424.9109-0.01102.635814.34022.56624.038513.02928.3904-0.021611.47013.4326delta =0.28050.24740.17660.18750.15780.42360.24250.16380.34260.32810.36990.32370.1749 y=8.55 3.79 4.82 0.02 2.75 14.39 2.54 4.35 13.00 8.50 0.05 11.32 3.13;x1=470 285 470 470 470 100 100 470 100 100 100 285 285;x2=300 80 300 80

8、 80 190 80 190 300 300 80 300 190;x3=10 10 120 120 10 10 65 65 54 120 120 10 120;x=x1 x2 x3;beta0=1,0.05, 0.02, 0.1, 2;beta,r,J=nlinfit(x,y,donglixue,beta0); yy,delta=nlpredci(donglixue,x,beta,r ,J); beta,r,J=nlinfit(x,y,donglixue,beta0); betabeta =1.25260.06280.04000.11241.1914 yy,delta=nlpredci(do

9、nglixue,x,beta,r ,J); yydeltayy =8.41793.95424.9109-0.01102.635814.34022.56624.038513.02928.3904-0.021611.47013.4326delta =0.28050.24740.17660.18750.15780.42360.24250.16380.34260.32810.36990.32370.1749可以得出在显著性水平为1-0.05的时候，置信区间yydelta4、混凝土的抗压强度随养护时间的延长而增加，现将一批混凝土作成12个试块，记录了养护日期x (日)及抗压强度 y (kg/cm2)的数

10、据:养护时间x234579121417212856抗压强度y354247535965687376828699试求y = a + b ln x型回归方程.建立m文件：function yhat=(beta,x)yhat=beta(1)+beta(2)*log(x);x=234579121417212856;y=354247535965687376828699;beta0=20 22;beta,r,J=nlinfit(x,y,volum,beta0)beta =21.005819.5285所以有回归方程：y=21.0058+19.5285log(x)5、下表给出了某工厂产品的生产批量与单位成本(元

11、)的数据，从散点图，可以明显的发现，生产 批量在500以内时，单位成本对生产批量服从一种线性关系，生产批量超过500时服从另一种线性 关系，此时单位成本明显下降。希望你构造一个合适的回归模型全面地描述生产批量与单位成本的 关系。生产批量 650 340 400 800 300 600 720 480 440 540 750单位成本 2.48 4.45 4.52 1.38 4.65 2.96 2.18 4.04 4.20 3.10 1.50先画出其散点图： x=650340 400 800 300 600 720 480 440 540 750; y=2.484.45 4.52 1.38 4.6

12、5 2.96 2.18 4.04 4.20 3.10 1.50; plot(x,y,*)6、一矿脉有13个相邻样本点，人为设定一个原点，现测得各样本点对原点的距离x,与该样本点某种金 属含量y的一组数据如下,画出散点图观察二者的关系，试建立合适的回归模型,如二次曲线,双曲线，对 数曲线等.x23457810y106.42109.20109.58109.50110.00109.93110.49x111415151819y110.59110.60110.90110.76111.00111.20 x1=2 3 4 5 7 8 10 11 14 15 15 18 19;y=106.42 109.20

13、109.58 109.50 110.00 109.93 110.49110.59 110.60 110.90 110.76 111.00 111.20;plot(x1,y)线性模型y=b0+b1*x x=ones(13,1),x1;b bint r rint stats=regress(y,x); b,bint,statsb = 108.25810.1742bint = 107.2794 109.23670.08910.2593stats =0.648420.28660.0009即模型是 y=108.2581+0.1742*x二次曲线模型 y=b0+b1*x+b2*xA2x=ones(13,1

14、),x1,x1.”2;b bint r rint stats=regress(y,x);b,bint,statsb = 106.95220.5271-0.0170bint = 105.4769 108.42750.18960.8645-0.0329-0.0011stats =0.775917.31120.0006即模型为 y=106.9522+0.5271*x-0.0170*xA2双曲线模型y=b0+b1/xx=ones(13,1),1./x1;b bint r rint stats=regress(y,x);b,bint,statsb = 111.4405-9.0300bint = 111.

15、1068 111.7743-10.6711-7.3889stats =0.9302 146.67330.0000即y=111.4405-9.0300/x 由stats变量看，此种模型拟合效果较好。对数模型 y=b0+b1*log(x)x=ones(13,1),log(x1);b bint r rint stats=regress(y,x);b,bint,statsb = 106.71131.5663bint = 105.6382 107.78441.08282.0499stats =0.822150.82850.0000即模型是 y=106.7113+1.5663*log(x)实验分析：由各拟合效果的stats变量来看，双曲线模型y=111.4405-9.0300/拟合的效果最好求这段曲线的纵坐标j关于横坐标x的二次多项式回归方程. x=0 2468 y=0.62.0 4.4