matlab代碼--函数逼近.docVIP

函数逼近 Chebyshev用切比雪夫多项式逼近已知函数 function f = Chebyshev(y,k,x0) 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 Legendre用勒让德多项式逼近已知函数 function f = Legendre(y,k,x0) syms t; P(1:k+1) = t; P(1) = 1; P(2) = t; c(1:k+1) = 0.0; c(1)=int(subs(y,findsym(sym(y)),sym(t))*P(1),t,-1,1)/2; c(2)=int(subs(y,findsym(sym(y)),sym(t))*P(2),t,-1,1)/2; f = c(1)+c(2)*t; for i=3:k+1 P(i) = ((2*i-3)*P(i-1)*t-(i-2)*P(i-2))/(i-1); c(i) = int(subs(y,findsym(sym(y)),t)*P(i),t,-1,1)/2; f = f + c(i)*P(i); if(i==k+1) if(nargin == 3) f = subs(f,t,x0); else f = vpa(f,6); end end end Pade用帕德形式的有理分式逼近已知函数 function f = Pade(y,n,x0) syms t; A = zeros(n,n); q = zeros(n,1); p = zeros(n+1,1); b = zeros(n,1); yy = 0; a(1:2*n) = 0.0; for(i=1:2*n) yy = diff(sym(y),findsym(sym(y)),n); a(i) = subs(sym(yy), findsym(sym(yy)), 0.0)/factorial(i); end; for(i=1:n) for(j=1:n) A(i,j)=a(i+j-1); end; b(i,1) = -a(n+i); end; q = A\b; p(1) = subs(sym(y),findsym(sym(y)),0.0); for(i=1:n) p(i+1) = a(n)+q(i)*subs(sym(y),findsym(sym(y)),0.0); for(j=2:i-1) p(i+1)=p(i+1)+q(j)*a(i-j); end end f_1 = 0; f_2 = 1; for(i=1:n+1) f_1 = f_1 + p(i)*(t^(i-1)); end for(i=1:n) f_2 = f_2 + q(i)*(t^i); end if(nargin == 3) f = f_1/f_2; f = subs(f,t,x0); else f = f_1/f_2; f = vpa(f,6); end lmz用列梅兹算法确定函数的最佳一致逼近多项式 function [coff,err]= lmz(func,m,a,b,eps) if(nargin == 4) eps=1.0e-6; end syms v; maxv = 0.0; max_x = a;