minpolyfit9.m

实现最小二乘的仿真函数

clc
clear
n=20;
l=9;
x=-1:1/n:1;
y=exp(x);
disp('多项式系数，x降幂排列到0');
P=polyfit(x,y,l)
%xx=linspace(log(cond(P)),n,60);
xx=-1:0.01:1;
yy=polyval(P,xx);
plot(xx,yy,'-b',x,y,'.r','MarkerSize',20)
legend('拟合曲线','原始数据','Location','SouthEast')
yf=polyval(P,x);
ol=0;
for i=1:n
ol=ol+(yf(i)-y(i))^2;
end
disp('偏差的最小值');ol
title('6次曲线拟合')
xlabel('x') 
ylabel('y')