假如你有一组实测数据,例如
x=[53 56 60 67.5 75 90 100];
y=[109 120.5 130 141.1 157.5 180 185];
求其回归直线,画回归直线图形并计算最小误差平方和,给出你的源程序
你可以这样子做:
x=[53 56 60 67.5 75 90 100];
y=[109 120.5 130 141.1 157.5 180 185];
p = polyfit(x,y,1);
yp = polyval(p,x);
plot(x,y,'ro',x,yp,'b-'); %画回归直线图形
legend('实测数据','回归直线')
R2 = sum((yp-y).^2) % 计算最小误差平方和
结果
误差平方和
R2 =
1.605360042163911e+02