function [M,C,de_error]=lspoly(X,Y,e)

% 調(diào)用形式：
%        [M,C]=lspoly(X,Y,e)
% 功能：
%       用最小二乘法多項(xiàng)式擬合非線性數(shù)據(jù)
% 輸入：
%       -- X   行向量，包含擬合數(shù)據(jù)點(diǎn)的橫坐標(biāo)
%       -- Y   行向量，包含擬合數(shù)據(jù)點(diǎn)的縱坐標(biāo)
%       -- e    標(biāo)量，誤差精度
% 輸出：
%     -- M  標(biāo)量，多項(xiàng)式的階數(shù)
%      --C  列向量，從低階到高階排列多項(xiàng)式的系數(shù)
%       其中調(diào)用子函數(shù)：[C,F]=sublspoly(X,Y,M)
% 袁怡圃 2003/04/28

n=length(X); 
Y=Y';
M=1;
de_error = 0;
while M<10
    [C,F]=sublspoly(X,Y,M);
    P=F*C;
    k=M;
    for i=1:n
        if abs((P(i)-Y(i))/P(i))>e
            M=M+1;
            break
        end
    end
    if k==M
        break
        de_error = 1;
    end
end
% min=min(X);
% max=max(X);
% u=min:0.01:max;
% n=length(u);
% f=zeros(n,M+1);
% for i=1:M+1
%     f(:,i)=u'.^(i-1);
% end
% p=f*C;
% plot(X,Y','*b',u,p,'r')