[求助] matlab曲线拟合 basic fitting 与polyfit

function [th,err,yi]=polyfits(x,y,N,xi,r)
%x,y  :数据点系列
%N    :多项式拟合的系统
%r    :加权系数的逆矩阵

M=length(x);
x=x(:);
y=y(:);

%判断调用函数的格式，考虑权阵
if nargin==4
%当用户调用函数格式为(x,y,N,r)
    if length(xi)==M
        r=xi;
        xi=x;
%当用户调用函数格式为(x,y,N, xi)，权阵为1
    else r=1;
    end
%当用户调用函数格式为(x,y,N)，不考虑权阵
elseif nargin==3
    xi=x;
    r=1;
end
%求解系数矩阵
A(:,N+1)=ones(M,1);
for n=N:-1:1
    A(:,n)=A(:,n+1).*x;
end

if length(r)==M
    for m=1:M
        A(m,:)=A(m,:)/r(m);
        y(m)=y(m)/r(m);
    end
end

%计算拟合系数
th=(A\y)';
ye=polyval(th,x);
err=norm(y-ye)/norm(y);
yi=polyval(th,xi);

clear;
clf;
x=input('请输入x值:');
y=input('请输入y值:');
% x=[ 10.5  13  13.8  16  18  19 20.8];
% y=[1190  1040 1000 800 600 430 200];
% eb=[0.2*ones(5,1);ones(5,1)];
[x,i]=sort(x);
y=y(i);
xi=min(x)+[0:100]/100*(max(x)-min(x));
for i=1:4
N=2*i-1;
[th,err1,yi]=polyfits(x,y,N,xi);
  subplot(220+i);
plot(x,y,'k*')
hold on
plot(xi,yi,'g:','LineWidth',1.5)
title(['The ',num2str(N),'th Polynomial Curve Fitting'])
grid on
end