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;