Practica 3.2..

PRACTICA 3 (parte2) METODOS NUMERICOS
METODO INTERPOLACIÓN NEWTON

function pol_newton
X=input('Ingrese los valores de X= ');
Y=input('Ingrese los valores de Y= ');
n=length (X);
D=zeros (n);
D(:,1)=Y
for j=2:n
for i=1:n-j+1
D(i,j)=(D(i+1,j-1)-D(i,j-1))/(X(i+j-1)-X(i));
end
end
syms x
pol=D(1,1);
p=1;
for i=1:n-1
p=p*(x-X(i));
pol=pol+D(1,i+1)*p;
end
disp(pol)
g=expand(pol)
disp(vpa(g,6))
x1=1:0.1:10;
y1=subs(pol,x,x1);
plot(x1,y1,'r')
hold on
plot (X,Y,'*'),grid
1.- Estime el logaritmo natural de 8 por medio de la interpolación lineal dado por el polinomio de
Newton a partir de:
a) log 6 = 1.7918 𝑦 log 10= 2.3026
>> pol_newton
Ingrese los valores de X= [6,10]
Ingrese los valores de Y= [1.7918,2.3026]
D=
1.7918 0
2.3026 0
(1277*x)/10000 + 641/625
g=
(1277*x)/10000 + 641/625
0.1277*x + 1.0256
b) calcule el error relativo porcentual con base al valor verdadero

Reemplazamos 8 en el polinomio obtenido:
A= 0.1277*x + 1.0256
A=2.0472
C) Error porcentual:
A1=Ln(8)=2.0794
A 1− A
%E= | A1 |
∗100
%E= |2.0794−2.0472
2.0794 |∗100=1.55 %
2.- Estime el logaritmo natural de 8 por medio de la interpolación cuadrática dado por el
polinomio de Newton a partir de:
a) uso de 𝑥 = 4, 6 𝑦 9 en la función logaritmo
b) calcule el error relativo porcentual con base al valor verdadero
 ln(4)=1.3863
 ln(6)=1.7917
 ln(9)= 2.1972
>> pol_newton
Ingrese los valores de X= [4,6,9]
Ingrese los valores de Y= [1.3863,1.7917,2.1972]
D=
1.3863 0 0
1.7917 0 0
2.1972 0 0
(2027*x)/10000 - (1013*(x - 4)*(x - 6))/75000 + 1151/2000
g=
(10133*x)/30000 - (1013*x^2)/75000 + 12567/50000
- 0.0135067*x^2 + 0.337767*x + 0.25134
Ahora reemplazamos 10 en el polinomio:
A= - 0.0135067*8^2 + 0.337767*8 + 0.25134

A=2.0890
A) Error porcentual:
A1=Ln(8)=2.0794
A 1− A
%E= | A1 |
∗100
%E= |2.0794−2.0890
2.0794 |∗100=0.46 %
3.- Dados los datos
𝑥 1.4 2.2 2.6 3 4 4.5
𝑓(𝑥) 2 8 14 15 8 2
a. Calcule 𝑓(2.8)𝑐𝑜𝑛 𝑒𝑙 𝑢𝑠𝑜 𝑑𝑒 𝑝𝑜𝑙𝑖𝑛𝑜𝑚𝑖𝑜 𝑑𝑒 𝑖𝑛𝑡𝑒𝑟𝑝𝑜𝑙𝑎𝑐𝑖ó𝑛 𝑑𝑒 Newton 𝑑𝑒 𝑜𝑟𝑑𝑒𝑛 1,
elija los valores para la construcción del polinomio entorno al valor buscado (un
valor inferior y otro superior)
b. Calcule 𝑓(2.8)𝑐𝑜𝑛 𝑒𝑙 𝑢𝑠𝑜 𝑑𝑒 𝑝𝑜𝑙𝑖𝑛𝑜𝑚𝑖𝑜 𝑑𝑒 𝑖𝑛𝑡𝑒𝑟𝑝𝑜𝑙𝑎𝑐𝑖ó𝑛 𝑑𝑒 Newton 𝑑𝑒 𝑜𝑟𝑑𝑒𝑛 3,
elija los valores para la construcción del polinomio entorno al valor buscado (dos
valores inferiores y dos valores superiores)
c. Determine los coeficientes de la parábola que pasa por los últimos tres puntos,
dados por el polinomio de Newton
 A. Polinomio de orden 1 (un valor inferior y un valor superior)

X=(2.5,3.2) f(x) = (14,3.2)
>> pol_newton
Ingrese los valores de X= [2.2,3]
Ingrese los valores de Y= [8,15]
D=
8 0
15 0
(35*x)/4 - 45/4
g=
(35*x)/4 - 45/4
8.75*x - 11.25
Ye=8.75*2.8-11.25
Ye=13.25
 B. Polinomio de orden 3,(dos valores superiores y dos valores inferiores)
>> pol_newton
Ingrese los valores de X= [2.2,2.6,4,4.5]
Ingrese los valores de Y= [8,14,8,2]
D=
8 0 0 0
14 0 0 0
8 0 0 0
2 0 0 0
15*x - (75*(x - 11/5)*(x - 13/5))/7 + (8850*(x - 4)*(x - 11/5)*(x - 13/5))/3059 - 25
g=
(8850*x^3)/3059 - (110655*x^2)/3059 + (423747*x)/3059 - 466436/3059
2.8931*x^3 - 36.1736*x^2 + 138.525*x - 152.48
Ye= 2.8931*2.8^3 - 36.1736*2.8^2 + 138.525*2.8 - 152.48
Ye=15.30
4.- Dados los datos

𝑥 1 2 3 5 7 8
𝑓(𝑥) 2.5 5 17 99 291 444
a.- Calcule (6.1)𝑠𝑜 𝑑𝑒 𝑝𝑜𝑙𝑖𝑛𝑜𝑚𝑖𝑜 𝑑𝑒 𝑖𝑛𝑡𝑒𝑟𝑝𝑜𝑙𝑎𝑐𝑖ó𝑛 𝑑𝑒 𝐿𝑎𝑔𝑟𝑎𝑛𝑔𝑒 𝑑𝑒 𝑜𝑟𝑑𝑒𝑛 1, elija los valores
para la construcción entorno al valor buscado
b.- Calcule (6.1)𝑠𝑜 𝑑𝑒 𝑝𝑜𝑙𝑖𝑛𝑜𝑚𝑖𝑜 𝑑𝑒 𝑖𝑛𝑡𝑒𝑟𝑝𝑜𝑙𝑎𝑐𝑖ó𝑛 𝑑𝑒 𝐿𝑎𝑔𝑟𝑎𝑛𝑔𝑒 𝑑𝑒 𝑜𝑟𝑑𝑒𝑛 3, elija los valores
para la construcción entorno al valor buscado
A) X=(1,5) f(x)=(2.5,99)
>> pol_newton
Ingrese los valores de X= [1,5]
Ingrese los valores de Y= [2.5,99]
D=
2.5000 0
99.0000 0
(193*x)/8 - 173/8
g=
(193*x)/8 - 173/8
24.125*x - 21.625
Reemplazando 6.1 en x
A=24.125*x – 21.625
A=24.125*6.1-21.625
ye =125.54
El valor estimado es: 627/5= 125.54
B) X=(1,3,5,8) f(x)=3,9,99,444)
>> pol_newton
Ingrese los valores de X= [1,3,5,8]
Ingrese los valores de Y= [2.5,17,99,444]
D=
2.5000 0 0 0
17.0000 0 0 0
99.0000 0 0 0
444.0000 0 0 0
(29*x)/4 + (135*(x - 1)*(x - 3))/16 + (509*(x - 1)*(x - 3)*(x - 5))/560 - 19/4
g=
(509*x^3)/560 + (9*x^2)/35 - (3133*x)/560 + 97/14
0.908929*x^3 + 0.257143*x^2 - 5.59464*x + 6.92857
Reemplazando 6.1 en x
A= 0.908929*6.1^3 + 0.257143*6.1^2 - 5.59464*6.1 + 6.92857
El valor estimado es: 188.68
5.- a) Determine un polinomio dado por interpolación Lagrange a partir de los valores
𝑥 2 2 4
𝑓(𝑥) 6 25 62
b) Observe que los valores de la tabla se generaron con la función 𝑓(𝑥) = 𝑥3 − 2, determine el
valor 𝑓(2.4), representado el valor correcto (verdadero o exacto)
c) Use el polinomio de interpolación de Lagrange para determinar 𝑓(2.4)
d) Calcule el error relativo porcentual verdadero
A) X=(1,2,3) f(x)= (-1,6,25)
>> pol_newton
Ingrese los valores de X= [2,2,4]
Ingrese los valores de Y= [6,25,62]
D=
6 0 0
25 0 0
62 0 0
Inf*(x - 2) - Inf*(x - 2)^2 + 6
g=
Inf*(x - 2) - Inf*(x - 2)^2 + 6
Inf*x - Inf - Inf*(x - 2.0)^2
B) f(x)=x^3-2
A1= f(2.6)=2.6^3-2 = 15.576
C) A= F(2.6)= Inf*2.6 - Inf - Inf*(x - 2.0)^2
El valor estimado es
304/25= 12.16
A 1− A
D) |
%E=
A1| ∗100
11.8240−12.16
%E= | 11.8240 |
∗100=2.8417 %
6.- a) Determine un polinomio lineal dado por interpolación Lagrange a partir de los valores
𝑥 1 ℮
𝑓(𝑥) 0 1
b) observe que los valores de la tabla se generaron con la función 𝑓(𝑥) = ln 𝑥, determine el valor
𝑓(3) = ln(3), representado el valor correcto (verdadero o exacto)
c) Use el polinomio de interpolación de Lagrange para determinar 𝑓(3)
d) Calcule el error relativo porcentual verdadero
A) X=(1,2.7183) F(x)=(0,1)
>> pol_newton
Ingrese los valores de X= [1,2.7183]
Ingrese los valores de Y= [0,1]
Ingrese el valor a estimar xe= log(3)
D=
0 0
1 0
(10000*x)/17183 - 10000/17183
la respuesta es
379/6604 ==0.05738
B) A1=F(x)=ln(x)
A1=F(3)=ln(3) = 1.09
C) A=F(3)=(10000*3)/17183 - 10000/17183
Ye=
El valor estimado es: 1.16
D) | A 1−
%E=
A
A1 |
∗100
%E= |0.693147−0.581970
0.693147 |∗100=16.0439%
7.- Emplee la porción de la tabla de vapor que se da para el 𝐻2𝑂 supercalentada a 200 MPa, para
a) encontrar la entropía correspondiente 𝑠 para un volumen específico 𝑣 de 0.108 𝑚3/𝑘𝑔 con
interpolación lineal Lagrange usando los dos primeros valores,
b) encontrar la misma entropía correspondiente con el uso de interpolación cuadrática dado por
Lagrange
𝑣[𝑚3/𝑘𝑔] 0.10377 0.11144 0.1254

𝑠[𝑘𝑙/𝑘𝑔 ∙ 𝐾] 6.4147 6.5453 6.7664
A) X=(0.10377,011144) f(x)=(6.4147,6.5453)
>> pol_newton
Ingrese los valores de X= [0.10377,0.11144]
Ingrese los valores de Y= [6.4147,6.5453]
Ingrese el valor a estimar xe= 0.108
D=
2444/381 0
3829/585 0
(13060*x)/767 + 35648387/7670000
la respuesta es
733/113 =5.51
B)
>> pol_newton
Ingrese los valores de X= [0.10377,0.11144,0.1254]
Ingrese los valores de Y= [6.4147,6.5453,6.7664]
Ingrese el valor a estimar xe= 0.108
D=
2444/381 0 0
3829/585 0 0
1651/244 0 0
(13060*x)/767 - (3869046362498253*(x - 1393/12500)*(x - 10377/100000))/70368744177664 +

35648387/7670000
la respuesta es
3380/521 ==6.48

Practica 3.2..

Cargado por

Practica 3.2..

Cargado por

PRACTICA 3 (parte2) METODOS NUMERICOS

METODO INTERPOLACIÓN NEWTON

b) calcule el error relativo porcentual con base al valor verdadero

Ahora reemplazamos 10 en el polinomio:

A= - 0.01350678^2 + 0.3377678 + 0.25134

 A. Polinomio de orden 1 (un valor inferior y un valor superior)

4.- Dados los datos

0.908929x^3 + 0.257143x^2 - 5.59464*x + 6.92857

𝑣[𝑚3/𝑘𝑔] 0.10377 0.11144 0.1254

(13060x)/767 - (3869046362498253(x - 1393/12500)*(x - 10377/100000))/70368744177664 +

También podría gustarte