新冠插值法matlab_MATLAB插值源代码(19个)

)*(-((abs(x1))^3)/6+(x1)^2-2*abs(x1)+4/3);

x2 = (a+index*h-x0)/h;

m2 = c(index+2)*((abs(x2))^3/2-(x2)^2+2/3);

x3 = (a+index*h+h-x0)/h;

m3 = c(index+3)*((abs(x3))^3/2-(x3)^2+2/3);

x4 = (a+index*h+2*h-x0)/h;

m4 = c(index+4)*(-((abs(x4))^3)/6+(x4)^2-2*abs(x4)+4/3);

f0 = m1+m2+m3+m4; %求出插值

function f = DCS(x,y,x0)

syms t;

if(length(x) == length(y))

n = length(x);

c(1:n) = 0.0;

else

disp('x和y的维数不相等！');

return;

end

c(1) = y(1);

for(i=1:n-1)

for(j=i+1:n)

y1(j) = (x(j)-x(i))/(y(j)-y(i));

end

c(i+1) = y1(i+1);

y = y1;

end

f = c(n);

for(i=1:n-1)

f