整理Matlab多元回归regress实例.docVIP

下载本文档

8
0
约1.67千字
约 2页
2021-01-26 发布于广东
举报
版权申诉

整理Matlab多元回归regress实例.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

PAGE PAGE 1 问题： y=[7613.51 7850.91 8381.86 9142.81 10813.6 8631.43 8124.94 9429.79 10230.81 10163.61 9737.56 8561.06 7781.82 7110.97]; x1=[7666 7704 8148 8571 8679 7704 6471 5870 5289 3815 3335 2927 2758 2591]; x2=[16.22 16.85 17.93 17.28 17.23 17 19 18.22 16.3 13.37 11.62 10.36 9.83 9.25] 回归成y=ax1^2+bx1^2+cx1+dx2+ex1*x2+f(二次曲线） regress解法： y=[7613.51 7850.91 8381.86 9142.81 10813.6 8631.43 8124.94 9429.79 10230.81 10163.61 9737.56 8561.06 7781.82 7110.97]; x1=[7666 7704 8148 8571 8679 7704 6471 5870 5289 3815 3335 2927 2758 2591]; x2=[16.22 16.85 17.93 17.28 17.23 17 19 18.22 16.3 13.37 11.62 10.36 9.83 9.25]; X=[ones(size(y)) x1.^2 x2.^2 x1 x2 x1.*x2]; [b,bint,r,rint,stats] = regress(y,X); format long b%b为对应的参数 b(1)为f，b(2)为a，b(3)为b，b(4)为c，b(4)为d，b(5)为e %b = 1.0e+004 * % -1.353935450267797 % 0.000000089381408 % -0.005811190715468 % -0.000605427789545 % 0.479983626458520 % -0.000037869040292 bint%为b的95%置信区间 %bint = 1.0e+004 * % -2.621944842897243 -0.085926057638351 % 0.000000034253753 0.000000144509063 % -0.027588831662545 0.015966450231609 % -0.001309493882546 0.000098638303455 % 0.119564693553906 0.840402559363135 % -0.000105954336341 0.000030216255756 stats%stats的第三个参数为F检验的P值，p值很小P0.001,说明拟合模型有效 %stats = 1.0e+005 * %0.000008444011951 0.000086828553270 0.000000043344434 3.162249735298925 scatter3(x1,x2,y,filled)%以下绘图查看拟合效果 hold on x1fit = min(x1):100:max(x1); x2fit = min(x2):1:max(x2); [X1FIT,X2FIT] = meshgrid(x1fit,x2fit); YFIT=b(1)+b(2)*X1FIT.^2+b(3)*X2FIT.^2+b(4)*X1FIT+b(5)*X2FIT+b(6)*X1FIT.*X2FIT; mesh(X1FIT,X2FIT,YFIT) mesh(X1FIT,X2FIT,YFIT) view(10,10) xlabel(x1) ylabel(x2) zlabel(y) %希望能解决你的问题