newton interpolacion
INGENIERÌA INDUSTRIAL Y DE SISTEMAS
METODOS NUMÈRICOS
“MÈTODO DE INTERPOLACIÓN DE NEWTON EN MATLAB
PRACTICA Nº4”
PRESENTA: LUIS DANIEL ESPINOZA DÀVILA
Cuatrimestre: 5°Grupo:901-5CV
PROFESOR: DANIEL FLORES RODRÍGUEZ
INTRODUCCIÓN
LA INTERPOLACION DE NEWTON SE BASA EN LA OBTENCIÓN DE UN POLINOMIO A PARTIR DE UN CONJUNTO DE PUNTOS DADO, APROXIMANDOSE LO MAS POSIBLEAL VALOR BUSCADO
EN MATLAB ESTA INTERPOLACION SE HARA DE MANERA RAPIDA Y EFICAZ
LA ECUACIÓN GENERAL PARA LA OBTENCION DE LA FUNCIÓN POR ESTE MÉTODO ES:Ejercicio 2.-
1.- se transcribe el código que resolverá el problema de interpolación de newton en matlab
clear;clc;
disp('metodos numericos');
disp('interpolacion de newton');
n=input('ingrese elgrado del polinomio, n=');
fprintf('Se necesitan %.0f puntos\n',n+1);
disp('ingrese los puntos');
for i=1:n+1
fprintf('x%.0f=',i-1);
X(i)=input(' ');
fprintf('y%.0f=',i-1);Y(i)=input(' ');
end
DD=zeros(n+1);
DD(:,1)=Y;
for k=2:n+1
for J=k:n+1
DD(J,k)=[DD(J,k-1)-DD(J-1,k-1)]/[X(J)-X(J-k+1)];
end
end
disp('La matriz de diferencias divididas es:');disp(DD);
disp('El polinomio de newton es');
syms x;
polnew=DD(1,1);
P=1;
for i=1:n
P=P*(x-X(i));
polnew=polnew+P*DD(i+1,i+1);
end
polnew=expand(polnew);
pretty(polnew);x=input('ingrese el valor de x a interpolar,x=');
vi=eval(polnew);
fprintf('el valor interpolado es %.2f\n',vi);
hold on;
ezplot(polnew,[X(1) X(n+1)]);
plot(x,vi,'r+');
2.-dado el código nos pide ingresarel grado del polinomio esto se traduce como el nivel de datos que tenemos en la tabla de nuestro ejercicio en este caso 5
3.-ingresar los puntos x, y
4.- acontinuación el programa nos resuelve la interpolación con la matriz dada en forma de bandera donde se van originando los valores de las divisiones de los datos interpolados y nos pide evaluar un...
Regístrate para leer el documento completo.