基于最小二乘法的线性回归分析

clear
clc
x=0:1/10:100;
y=x.^2+2*x+1;
y=y/max(y);
ye=y+0.2*rand(size(y));
ye=ye-0.1;
p=polyfit(x,ye,2);
yp=polyval(p,x);
plot(x,y);
hold on;
plot(x,ye,'r.');
plot(x,yp,'k');
figure
plot(x,y-yp);