数学建模社区-数学中国

标题: 用切比雪夫多项式逼近已知函数 [打印本页]

作者: 鬼谷子天明 时间: 2011-10-28 18:16
标题: 用切比雪夫多项式逼近已知函数
function f = Chebyshev(y,k,x0)
%用切比雪夫多项式逼近已知函数
%已知函数y
%逼近已知函数所需项数：k
%逼近点的x坐标：x0
%求得的切比雪夫逼近多项式或是在x0处的逼近值：f
syms t;
T(1:k+1) = t;
T(1) = 1;
T(2) = t;
c(1:k+1) = 0.0;
c(1)=int(subs(y,findsym(sym(y)),sym('t'))*T(1)/sqrt(1-t^2),t,-1,1)/pi;
c(2)=2*int(subs(y,findsym(sym(y)),sym('t'))*T(2)/sqrt(1-t^2),t,-1,1)/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(y,findsym(sym(y)),sym('t'))*T(i)/sqrt(1-t^2),t,-1,1)/2;
f = f + c(i)*T(i);
f = vpa(f,6);

if(i==k+1)
      if(nargin == 3)
         f = subs(f,'t',x0);
      else
         f = vpa(f,6);
      end
end
end
调用时：
f = Chebyshev('1/(2-x)',6)
Warning: Explicit integral could not be found.
??? The following error occurred converting from sym to double:
Error using ==> mupadmex
Error in MuPAD command: DOUBLE cannot convert the input
expression into a double array.

If the input expression contains a symbolic variable, use the
VPA function instead.

Error in ==> Chebyshev at 7
c(1)=int(subs(y,findsym(sym(y)),sym('t'))*T(1)/sqrt(1-t^2),t,-1,1)/pi;

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)