- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
3.残差分析,作残差图: rcoplot(r,rint) 从残差图可以看出,除第二个数据外,其余数据的残差离零点均较近,且残差的置信区间均包含零点,这说明回归模型 y=-16.073+0.7194x能较好的符合原始数据,而第二个数据可视为异常点. 4.预测及作图: z=b(1)+b(2)* plot(x,Y,k+,x,z,r) 返回 To MATLAB(liti12) 第三十页,编辑于星期五:十三点 二十分。 多 项 式 回 归 (一)一元多项式回归 (1)确定多项式系数的命令:[p,S]=polyfit(x,y,m) (2)一元多项式回归命令:polytool(x,y,m) 1.回归: y=a1xm+a2xm-1+…+amx+am+1 2.预测和预测误差估计: (1)Y=polyval(p,x)求polyfit所得的回归多项式 在x处 的预测值Y; (2)[Y,DELTA]=polyconf(p,x,S,alpha) 求polyfit所得的回归多项式在x处的预测值Y及预测值的显 著性为1-alpha的置信区间Y DELTA;alpha缺省时为0.5. ± 第三十一页,编辑于星期五:十三点 二十分。 法一 直接作二次多项式回归: t=1/30:1/30:14/30; s=[11.86 15.67 20.60 26.69 33.71 41.93 51.13 61.49 72.90 85.44 99.08 113.77 129.54 146.48]; [p,S]=polyfit(t,s,2) To MATLAB(liti21) 得回归模型为 : 第三十二页,编辑于星期五:十三点 二十分。 法二 化为多元线性回归: t=1/30:1/30:14/30; s=[11.86 15.67 20.60 26.69 33.71 41.93 51.13 61.49 72.90 85.44 99.08 113.77 129.54 146.48]; T=[ones(14,1) t (t.^2)]; [b,bint,r,rint,stats]=regress(s,T); b,stats To MATLAB(liti22) 得回归模型为 : Y=polyconf(p,t,S) plot(t,s,k+,t,Y,r) 预测及作图 To MATLAB(liti23) 第三十三页,编辑于星期五:十三点 二十分。 (二)多元二项式回归 命令:rstool(x,y,’model’, alpha) n?m矩阵 显著性水平 (缺省时为0.05) n维列向量 第三十四页,编辑于星期五:十三点 二十分。 例3 设某商品的需求量与消费者的平均收入、商品价格的统计数 据如下,建立回归模型,预测平均收入为1000、价格为6时 的商品需求量. 法一 直接用多元二项式回归: x1=[1000 600 1200 500 300 400 1300 1100 1300 300]; x2=[5 7 6 6 8 7 5 4 3 9]; y=[100 75 80 70 50 65 90 100 110 60]; x=[x1 x2]; rstool(x,y,purequadratic) 第三十五页,编辑于星期五:十三点 二十分。 在画面左下方的下拉式菜单中选”all”, 则beta.rmse和residuals都传送到MATLAB工作区中. 将左边图形下方方框中的“800”改成1000,右边图形下方的方框中仍输入6.则画面左边的“Predicted Y”下方的数据由原来的“86.3791”变为88.4791,即预测出平均收入为1000.价格为6时的商品需求量为88.4791. 第三十六页,编辑于星期五:十三点 二十分。 在MATLAB工作区中输入命令: beta, rmse To MATLAB(liti31) 第三十七页,编辑于星期五:十三点 二十分。 第一页,编辑于星期五:十三点 二十分。 优选matlab多元回归分析 第二页,编辑于星期五:十三点 二十分。 实验目的 实验内容 2.掌握用数学软件求解回归分析问题. 1.直观了解回归分析基本内容. 1.回归分析的基本理论. 3.实验作业. 2.用数学软件求解回归分析问题. 第三页,编辑于星期五:十三点 二十分。 一元线性回归 多元线性回归 回归分析 数学模型及定义 *模型参数估计 *检验、预测与控制 可线性化的一元非线 性回归(曲线回归) 数学模型及定义 *模型参数估计 逐步回归分析 * 多元线性回归中的 检验与预测 第四页,编辑于
文档评论(0)