[go: up one dir, main page]
Scribd
Cargar
333 vistas4 páginas

Interpolación Lineal Matlab

Este documento presenta el código para realizar un ajuste polinomial de datos usando Matlab. Primero se grafican los datos y se realiza un ajuste polinomial de cuarto grado. Luego se hace un ajuste polinomial de segundo grado y se grafican ambas curvas. Finalmente, se calculan las raíces del polinomio de segundo grado y la integral de este polinomio.

Cargado por

Memo Leon
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
333 vistas4 páginas

Interpolación Lineal Matlab

Este documento presenta el código para realizar un ajuste polinomial de datos usando Matlab. Primero se grafican los datos y se realiza un ajuste polinomial de cuarto grado. Luego se hace un ajuste polinomial de segundo grado y se grafican ambas curvas. Finalmente, se calculan las raíces del polinomio de segundo grado y la integral de este polinomio.

Cargado por

Memo Leon
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
Está en la página 1/ 4

Guillermo León Nuñez

3EV4 METODOS NUMERICOS FLORES JAIME ARMANDO

INGENIERIA ELECTRICA
clc

clear

x=[0 1 2 3 4 5];

y=[-2 -3 8 67 234 593];

plot(x,y,'o')

grid

pause

hold on

P=polyfit(x,y,4,'r')

t=0:0.001:5;

Px=P(1)*t.^4+P(2)*t.^3+P(3)*t.^2+P(4)*t+P(5);

plot(t,Px)

P=polyfit(x,y,2)

Gx=P(1)*t.^2+P(2)*t+P(3);

plot(t,Gx)

hold off

xr=roots(P)

disp('Solucion Examen')

fprintf('xr1=%.3f\n', xr(1))

fprintf('xr2=%.3f\n', xr(2))

syms t

p=[x' y']

Gx=P(1)*t.^2+P(2)*t+P(3);

A=int(Gx,0,5);

fprintf('Gx=%.1f t^4+%.1f t^3+%.1f t^2+%.1f\n', P(1),P(2),P(3),P(4),P(5))

fprintf('A=%.3f u^2\n', A)
clc
clear
x=[0 1 2 3 4 5];
y=[-2 -3 8 67 234 593];
plot(x,y,'o')
grid
pause
hold on
P=polyfit(x,y,4)
t=0:0.001:5;
Px=P(1)*t.^4+P(2)*t.^3+P(3)*t.^2+P(4)*t+P(5);
plot(t,Px)
P=polyfit(x,y,2)
Gx=P(1)*t.^2+P(2)*t+P(3);
plot(t,Gx)
hold off
xr=roots(P)
disp('Solucion Examen')
fprintf('xr1=%.3f\n', xr(1))
fprintf('xr2=%.3f\n', xr(2))

syms t
p=[x' y']
Gx=P(1)*t.^2+P(2)*t+P(3);
A=int(Gx,0,5);
fprintf('Gx=%.1f t^4+%.1f t^3+%.1f t^2+%.1f\n', P(1),P(2),P(3),P(4),P(5))
fprintf('A=%.3f u^2\n', A)

También podría gustarte