matlab插值算法实现.docxVIP

1.高斯插值 function f = Gauss(x,y,x0) if(length(x) == length(y)) n = length(x); else disp(x 和 y 的维数不相等！ ); return; end xx =linspace(x(1),x(n),(x(2)-x(1))); if(xx ~= x) disp(节点之间不是等距的！); return; end if( mod(n,2) ==1) if(nargin == 2) f = GStirling(x,y,n); else if(nargin == 3) f = GStirling(x,y,n,x0); end end else if(nargin == 2) f = GBessel(x,y,n); else if(nargin == 3) f = GBessel(x,y,n,x0); end end end function f = GStirling(x,y,n,x0) syms t; nn = (n+1)/2; f = y(nn); for(i=1:n-1) for(j=i+1:n) y1(j) = y(j)-y(j-1); end if(mod(i,2)==1) c(i) = (y1((i+n)/2)+y1((i+n+2)/2))/2; else c(i) = y1((i+n+1)/2)/2;