3元函数拟合及结果.docx

代码x1=[13.51,14.08,13.82,13.40,13.77,13.41,13.56,14.08,13.23,13.35,13.09,13.86,13.07,13.19,13.30,13.58,13.95,13.59,13.45,13.76,13.58,13.48,13.39,13.35,13.37,13.46,13.20,13.18,13.21,13.38];x2=[38.69,35.13,30.80,36.02,34.64,38.74,35.33,34.60,39.13,34.08,39.93,39.18,39.68,34.14,37.74,37.75,39.06,36.13,38.15,34.38,37.75,40.56,34.67,38.27,36.98,38.68,42.48,38.62,41.32,36.10];x3=[2.00,1.95,2.10,2.11,2.01,1.92,1.77,1.89,2.10,1.89,2.05,1.85,1.97,1.83,1.76,2.02,2.04,1.88,2.06,1.97,2.02,2.00,2.01,1.91,1.95,1.83,1.98,1.89,2.03,1.94];y=[20.30,21.76,20.49,20.24,20.84,20.02,20.10,21.58,19.84,19.26,19.50,19.83,19.17,19.62,19.76,20.76,21.66,20.40,20.43,20.90,20.78,20.33,19.85,19.81,19.62,19.59,19.58,19.36,19.82,19.71];%每3个x为一组对应一个y，请问怎样编写matlab程序，进行最小二乘法拟合，得到回归方程format short gY=yX11=[ones(1,length(y));x1;x2;x3]B1=regress(Y,X11)% 多元一次线性回归[m,n]=size(X11)X22=[];for i=2:n for j=2:n if i=j X22=([X22,X11(:,i).*X11(:,j)]); else continue end endendX=[X11,X22];B2=regress(Y,X)% 多元二次线性回归[Y X*B2 Y-X*B2]plot(Y,X11*B1,o,Y,X*B2,*)hold on,line([min(y),max(y)],[min(y),max(y)])axis([min(y) max(y) min(y) max(y)])legend(一次线性回归,二次线性回归)xlabel(实际值);ylabel(计算值)运行结果 qqpzY = 20.3 21.76 20.49 20.24 20.84 20.02 20.1 21.58 19.84 19.26 19.5 19.83 19.17 19.62 19.76 20.76 21.66 20.4 20.43 20.9 20.78 20.33 19.85 19.81 19.62 19.59 19.58 19.36 19.82 19.71X11 = 1 13.51 38.69 2 1 14.08 35.13 1.95 1 13.82 30.8 2.1 1 13.4 36.02 2.11 1 13.77 34.64 2.01 1 13.41 38.74 1.92 1 13.56 35.33 1.77 1 14.08 34.6 1.8