Matlab多元回归regress实例.pptx

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