切比雪夫（Chebyshev）多项式展开式的matlab程序

切比雪夫(Chebyshev)多项式展开式的matlab程序，输入一个函数 可以很方便直接得到切比雪夫的展开式

function f=Chebyshev（ykx0）
%用切比雪夫多项式逼近已知函数?%已知函数：y?
%逼近已知函数所需项数：k?%逼近点的x坐标：x0?
%求得的切比雪夫逼近多项式或在x0处的逼近值：f? 

syms t;
T（1:k+1）=t;
T（1）=sym（‘1‘）;
T（2）=t;
c（1:k+1）=sym（‘0‘）;
c（1）=int（subs（yfindsym（sym（y））sym（‘t‘））*T（1）/sqrt（1-t^2）t-11）/pi;
c（2）=2*int（subs（yfindsym（sym（y））sym（‘t‘））*T（2）/sqrt（1-t^2）t-11）/pi;f=c（1）+c（2）*t;
    for i=3:k+1
     T（i） = 2*t*T（i-1）-T（i-2）; 
      c（i） = 2*int（subs（yfindsym（sym（y））sym（‘t‘））*T（i）/sqrt（1-t^2）t-11）/2;    
      f = f + c（i）*T（i）;     
      f = vpa（f6）;      
    if（i==k+1） 
        if（nargin == 3） 
            f = subs（f‘t‘x0）;         
           else 
            f = vpa（f6）;         
              end     
                 end 
                    end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        800  2014-03-28 19:57  Chebyshev.m

     文件        886  2014-04-02 15:24  ChebyshevPoly.m

     文件        131  2014-04-14 21:41  redeme.txt

----------- ---------  ---------- -----  ----

                 1817                    3