URGENTE Algoritmo Diferencias Divididas (interpolacion Newton)
Necesito urgentemente que alguien me explique como puedo programar un algoritmo para diferencias divididas, es decir, interpolacion de newton :-(
Gracias
Gracias
Por favor enviame a mi correo si ya conseguiste algo sobre algoritmo diferencias divididas(interpolacion).Gracias
Julio
Julio
espero t ayude el programa de algoritmos de newton para dato
1.y=[0.7651977;0.6200860;0.4554022;0.2818186;0.1103623]'
x=[1.0;1.3;1.6;1.9;2.2]'
p=[1.7 5 2 1.5 0]
plot(x,y)
hold on
grid
fp=newton1(x,y,p)
2.function fp=newton1(x,y,p)
n1=length(p)
n=length(x);
a(1) = y(1);
for k = 1 : n - 1
d(k, 1) = (y(k+1) - y(k))/(x(k+1) - x(k));
end
for j = 2 : n - 1
for k = 1 : n - j
d(k, j) = (d(k+1, j - 1) - d(k, j - 1))/(x(k+j) - x(k));
end
end
d
for j=2:n
a(j) = d(1, j-1);
end
Df(1) = 1;
c(1) = a(1);
for i=1:n1
for j = 2 : n
Df(j)=(p(i) - x(j-1)) .* Df(j-1);
c(j) = a(j) .* Df(j);
end
fp(i)=sum(c);
end
1.y=[0.7651977;0.6200860;0.4554022;0.2818186;0.1103623]'
x=[1.0;1.3;1.6;1.9;2.2]'
p=[1.7 5 2 1.5 0]
plot(x,y)
hold on
grid
fp=newton1(x,y,p)
2.function fp=newton1(x,y,p)
n1=length(p)
n=length(x);
a(1) = y(1);
for k = 1 : n - 1
d(k, 1) = (y(k+1) - y(k))/(x(k+1) - x(k));
end
for j = 2 : n - 1
for k = 1 : n - j
d(k, j) = (d(k+1, j - 1) - d(k, j - 1))/(x(k+j) - x(k));
end
end
d
for j=2:n
a(j) = d(1, j-1);
end
Df(1) = 1;
c(1) = a(1);
for i=1:n1
for j = 2 : n
Df(j)=(p(i) - x(j-1)) .* Df(j-1);
c(j) = a(j) .* Df(j);
end
fp(i)=sum(c);
end
