pade_mine.m

对函数进行PADE近似程序可以对多种函数进行近似

function ya1=pade_mine(nn,mm,ya,xx,t)
    
       ya0=ya;
       p=ya(t);  
       yy=xx(t);

for i=1:nn
    for j=1:nn
        if j<=(nn-1)/2
            A(i,j)=p(i)*(yy(i))^(j);
        else
            A(i,j)=-(yy(i))^(j-(nn+1)/2);
        end
    end
end
 C=inv(A)*(-p');%这就是系数矩阵
 C1=C(1:(nn-1)/2);
 C2=C((nn+1)/2:nn);
   xx0=[];
for i=1:(length(xx)-1)
    ww=linspace(xx(i),xx(i+1),mm+2);%在第i到第i+1间插入mm个点的横坐标
    xx0=[xx0 ww(1:(mm+1))];
   end

   xx0=[xx0 xx(length(xx))];%补上最后一个点的值。
  
for i=1:length(xx0)
           g=www(xx0(i),(nn-1)/2);
               p0=1+C1'*g';
             g2=[1 g];
             q0=C2'*g2';
            ya1(i)=q0/p0;
           end

   cccc=1; 
 
  function ww1=www(ww,Q)
        ww1=[];
for ii=1:Q
    ww1=[ww1 ww^ii];
end