用Matlab做曲线拟合的最小二乘法

来源:百度知道 编辑:UC知道 时间:2024/06/21 04:54:43
xl=[350 350 350 350 350 350 350 350 350 500 500 500 500 500 500 500 500 500 650 650 650 650 650 650 650 650 650];
x2=[2 2 2 3.5 3.5 3.5 5 5 5 2 2 2 3.5 3.5 3.5 5 5 5 2 2 2 3.5 3.5 3.5 5 5 5];
x3=[2.33 1 0.428 2.33 1 0.428 2.33 1 0.428 2.33 1 0.428 2.33 1 0.428 2.33 1 0.428 2.33 1 0.428 2.33 1 0.428 2.33 1 0.428];
y=[90.27 93.68 92.65 87.4 91.71 93.94 86.49 93.9 93.6 82.4 92.27 93.98 87.58 93.01 93.24 82.96 92.61 94.83 90.32 92.6 92.77 89.97 90.31 93.47 90.51 92.61 94.04];
每3个x为一组对应一个y,请问怎样编写matlab程序,进行最小二乘法拟合,得到回归方程……很迷茫,看书都不知道从哪里看起,还没有入门……

x1=[350 350 350 350 350 350 350 350 350 500 500 500 500 500 500 500 500 500 650 650 650 650 650 650 650 650 650];
x2=[2 2 2 3.5 3.5 3.5 5 5 5 2 2 2 3.5 3.5 3.5 5 5 5 2 2 2 3.5 3.5 3.5 5 5 5];
x3=[2.33 1 0.428 2.33 1 0.428 2.33 1 0.428 2.33 1 0.428 2.33 1 0.428 2.33 1 0.428 2.33 1 0.428 2.33 1 0.428 2.33 1 0.428];
y=[90.27 93.68 92.65 87.4 91.71 93.94 86.49 93.9 93.6 82.4 92.27 93.98 87.58 93.01 93.24 82.96 92.61 94.83 90.32 92.6 92.77 89.97 90.31 93.47 90.51 92.61 94.04];
%每3个x为一组对应一个y,请问怎样编写matlab程序,进行最小二乘法拟合,得到回归方程……很迷茫,看书都不知道从哪里看起,还没有入门……
format short g
Y=y'
X11=[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
end
end
X=[X11,X22];
B2=regress(Y,X)% 多元二次线性回归
[Y X*B2 Y-X*B2]

plot(Y,X11*B1,'