资源描述:
《《matlab实验》word版》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、成绩辽宁工程技术大学上机实验报告实验名称实验十一回归分析院系理学院专业信科班级姓名学号日期20131121实验目的简述本次实验目的:1、了解回归分析基本内容2、掌握用matlab软件求解回归分析问题。实验准备你为本次实验做了哪些准备:经过看书和看ppt完成了此次实验。实验进度本次共有6个练习,完成6个。实验总结日本次实验的收获、体会、经验、问题和教训:通过本次实验学会了有关回归分析的有关知识,怎么求线性与非线性回归模型,知道了怎么找β系数和置信区间,知道了有关回归分析有关的matlab命令,在实验的过程中遇
2、到了很多的难题,比如怎么观察模型是不是线性的,比如怎么找到最合适的模型,这时通过看书和认真思考学会了怎么做。本次实验收获了很多,知道了自己的不足,以后还需要勤加训练!教师评语1、考察温度x对产量y的影响,测得下列10组数据:温度(℃)20253035404550556065产量(kg)13.215.116.417.117.918.719.621.222.524.3求y关于x的线性回归方程,检验回归效果是否显著,并预测x=42℃时产量的估值及预测区间(置信度95%).>>x=[2025303540455055
3、6065]';>>X=[ones(10,1)x];>>y=[13.215.116.417.117.918.719.621.222.524.3]’;r=-0.38180.40300.58790.1727-0.1424-0.4576-0.6727-0.1879-0.00300.6818rint=-1.28580.5221-0.56751.3736-0.36391.5397-0.92931.2748-1.26320.9783-1.51230.5972-1.61790.2725-1.25630.8806-1.035
4、21.0291-0.07631.4399>>corrcoef(x,y)ans=1.00000.99100.99101.0000>>[b,bint,r,rint,stats]=regress(y,X)b=9.12120.2230bint=8.021110.22140.19850.2476stats=0.9821439.83110.00000.2333置信区间[8.0211,10.2214]置信区间[0.1985,0.2476]r2=0.9821F=439.8311p=0.0000p<0.05回归模型:y=9.
5、1212+0.2230x成立r2=0.9821接近1回归方程显著X=42时y=18.48722、某零件上有一段曲线,为了在程序控制机床上加工这一零件,需要求这段曲线的解析表达式,在曲线横坐标xi处测得纵坐标yi共11对数据如下:xi02468101214161820yi0.62.04.47.511.817.123.331.239.649.761.7求这段曲线的纵坐标y关于横坐标x的二次多项式回归方程.(并画出图形)>>x=[02468101214161820];>>y=[0.62.04.47.511.817
6、.123.331.239.649.761.7];X=100124141616361864110100112144114196116256118324120400>>[p,S]=polyfit(x,y,2)p=0.14030.19711.0105S=R:[3x3double]df:8normr:1.1097>>Y=polyconf(x,y,S)回归模型:y=0.1403x2+0.1971x+1.0105>>X=[ones(11,1)x'(x.^2)'][b,bint,r,rint,stats]=regress
7、(y',X);>>Y=polyconf(p,x,S)>>plot(x,y,'k+',x,Y,'r')Y=1.01051.96604.04417.244911.568317.014223.582831.274040.087850.024261.08323、在研究化学动力学反应过程中,建立了一个反应速度和反应物含量的数学模型,形式为其中是未知参数,是三种反应物(氢,n戊烷,异构戊烷)的含量,y是反应速度.今测得一组数据如下表,试由此确定参数,并给出置信区间.的参考值为(1,0.05,0.02,0.1,2).序号
8、反应速度y氢x1n戊烷x2异构戊烷x318.554703001023.79285801034.8247030012040.024708012052.754708010614.391001901072.54100806584.3547019065913.0010030054108.50100300120110.05100801201211.3228530010133.13285190120解:>>x1=[47028