最小二乘法拟合

x = [28.8 42.8 57.4 71.5 85.6 100 115 144  173];
y = [20   30   40   50   60   70  80  100  120];
plot(x,y,'o');
grid on
p = polyfit(x,y,1)  %1为你要拟合的阶数，一阶直线拟合，二阶抛物线拟合
hold on             %新图不会覆盖旧图

x1 = 0:0.01:180;
y1 = polyval(p,x1);
plot(x1,y1,'k')
 
xlabel('方波频率f(Hz)');
ylabel('车速V(km/h)');
title('黑线为最小二乘法拟合')