0

he estado tratando de resolver un ejercicio de interpolación de Hermite con matlab pero no logro ejecutarlo, pueden ayudarme con las correcciones? no entiendo por que el código no se ejecuta y me muestra demasiados errores. para el segundo Código, he intentado calcular la cota de error sin ningún éxito. es el siguiente:

% Definir la función f(x)
f = @(x) 3*x*exp(x) - exp(2*x);

% Definir los puntos conocidos
x0 = 1;
x1 = 1.05;

% Calcular las derivadas de f en los puntos x0 y x1
df_x0 = (3*exp(x0) + 3*x0*exp(x0)) - 2*exp(2*x0);
df_x1 = (3*exp(x1) + 3*x1*exp(x1)) - 2*exp(2*x1);

% Calcular el polinomio interpolante de Hermite de grado 3
syms x;
p = ((1 - 2*(x - x0)*subs(diff((x - x0)^2*(x - x1)^2), x0)) * f(x0) + (x - x0)*subs(diff((x - x0)*(x - x1)^2), x0)*df_x0 + (1 - 2*(x - x1)*subs(diff((x - x0)*(x - x1)^2), x1)) * f(x1) + (x - x1)*subs(diff((x - x0)^2*(x - x1)), x1)*df_x1);

% Aproximación de f(1.03) usando el polinomio interpolante
aproximacion = double(subs(p, 1.03));

% Calcular el error real
valor_exacto = f(1.03);
error_real = abs(valor_exacto - aproximacion);

% Calcular la cota de error
cuarta_derivada = diff(f, 4);
xi = 1.03; % Punto en el que se evalúa la cota de error
cota_error = abs(subs(cuarta_derivada, xi) / (4*factorial(4)) * (xi - 1)*(xi - 1)*(xi - x0)*(xi - x1));

% Mostrar resultados
disp(['Aproximación del polinomio a 1.03: ', num2str(aproximacion)]);
disp(['Error real: ', num2str(error_real)]);
disp(['Cota de error: ', num2str(cota_error)]);

% Graficar el polinomio interpolante
fplot(p, [1, 1.05]);
title('Polinomio interpolante de Hermite de grado 3');
xlabel('x');
ylabel('p(x)');

tambien he intentado con este otro código, el detalle es que este no me muestra la cota de error, no logro entender.

%--------------------------------------------------
function P = HPoly(x,y,yp)
    n = length(x) ;
    for k=1:(2*n)
        Bx(:,k)=x'.^(k-1) ;
    end
    Bp(:,1)=zeros(n,1) ;
    for k=2:(2*n)
        Bp(:,k)=(k-1)*x'.^(k-2) ;
    end
    P = flip([Bx;Bp]\[y,yp]')';
end
%--------------------------------------------------
function y = f(x)
    y = 3*x.*exp(x)-exp(2*x) ;
end
function y = fp(x)
    y = 3*exp(x).*(x + 1) - 2*exp(2*x);
end
%--------------------------------------------------
% GRADO 3
%
x  = [1, 1.05] ;
y  = f(x) ;
yp = fp(x);

xtest = 1.03 ;
%--------------------------------------------------

P = HPoly(x,y,yp) ;
Err3 = polyval(P,xtest) - f(xtest)
function Lagrange_error_bound = max(abs((xtest-x(1))*(xtest-x(2)))*(abs(fp(t))/(2*factorial(3))), abs((xtest-x(1))*(xtest-x(2))*(xtest-x(3)))*(abs(fp(t))/(6*factorial(4))));
  end
xb = [1:.001:1.05];
yb = f(xb) ;
H = polyval(P,xb) ;

subplot(2,2,1);plot(xb,yb,"r",xb,H,"g",x,y,"r*")
subplot(2,2,2);plot(xb,H-yb,xtest,Err3,"r*")

%--------------------------------------------------

% GRADO 5
%
x  = [1, 1.05, 1.07] ;
y  = f(x) ;
yp = fp(x);

xtest = 1.03 ;
%--------------------------------------------------

P = HPoly(x,y,yp) ;
Err5 = polyval(P,xtest) - f(xtest)

xb = [1:.001:1.07];
yb = f(xb) ;
H = polyval(P,xb) ;

subplot(2,2,3);plot(xb,yb,"r",xb,H,"g",x,y,"r*")
subplot(2,2,4);plot(xb,H-yb,xtest,Err5,"r*")
%--------------------------------------------------

necesito calcular la cota de error tanto para el polinomio de grado tres, como para el de grado 5. Espero encontrar ayuda.

0

Tu Respuesta

By clicking “Publica tu respuesta”, you agree to our terms of service and acknowledge you have read our privacy policy.

Examina otras preguntas con la etiqueta o formula tu propia pregunta.