基于Matlab实现多次最佳一致的函数逼近（类似求渐进函数）

%%%做系统识别很重要，方法上完全符合系统识别最基础的理论

function [sun]=main(n)
fplot('1/(x+2)',[-1,1],'r');
x=ones(n+2,1);
for j=0:n+1
    x(j+1)=cos(pi*(n+1-j)/(n+1));
end
first=ones(n+2,1); 
f=1./(x+2);  %原函数
last=first;
for j=2:n+2
   last(j)=(-1)*last(j-1);
end
A=ones(n+2,n+2);
A(:,1)=first;
A(:,n+2)=last;
for j=2:n+1
    for t=2:j
        A(:,j)=x.*A(:,j);
    end
end
e=(1e-15)*first; %精度控制条件
sun=A\f;
while (1)    
    at=''; 
    for i=1:n
        for t=1:i
            if (t==1)
                at=strcat('(',num2str(sun(t+1)),')');
            elseif (t>1)
                xt='t';j=2; 
                while j<t
                    xt=strcat('t*',xt);j=j+1;
                end 
                at=strcat(num2str(t),'*(',num2str(sun(t+1)),')*