Matlab多元回归regress实例.pdf

下载文档 降价啦

11
0
约1.93千字
约 2页
2021-01-31 发布于广东
举报
版权申诉
保障服务

Matlab多元回归regress实例.pdf

1、本文档共2页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

问题： 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; 1 mesh(X1FIT,X2FIT,YFIT) mesh(X1FIT,X2FIT,YFIT) view( 10,10) xlabel(x1) ylabel(x2) zlabel(y) %希望能解决你的问题 2