interpolacion de la grange
explicacion :
interpolacion es, apartir de una serie de puntos, obtener una ecuacion cuya curva
pase por todos ellos o lo mas cerca posible.
en el metodo de la grange, se uiliza una ecuacion que aunque se va alargado conforme
mas puntos se quieren unir, es siempre del mismo tamaño y de la misma forma por lo
que una de sus ventajas es que es mas claro y facil de hacer.
procedimiento:
primero se deben tener ya los puntos que se quieren interpolar :
x f(x) despues se empieza a rellenar la ecuacion de la grange
x0 f(x0) comenzando por el primer punto a tomarlo como el punto o
sea x0 de referencia.
x1 f(x1)
x2 f(x02) g(x)=((𝑥−𝑥1)(𝑥−𝑥2)
x3 f(x3) (𝑥−𝑥3))/((𝑥0−𝑥1)(𝑥0−𝑥2)
(𝑥0−𝑥3))f(x0)
el metodo de lagrange se compone por una serie de sumandos, cada uno se "arma"
tomando como referencia cada uno de los punto a interpolar empezando desde el
primero, si se analiza la ecuacion de arriba, se puede ver como es que se forma
cada sumando, en el caso de arriba con el primer punto como referencia.
asi el segundo sumando quedaria de la siguiente manera:
asi el segundo sumando quedaria de la siguiente manera:
g(x)=((𝑥−𝑥1)(𝑥−𝑥2)(𝑥−𝑥3))/((𝑥0−𝑥1)(𝑥0−𝑥2)
(𝑥0−𝑥3))f(x0)+ ((𝑥−𝑥𝑜)(𝑥−𝑥2)(𝑥−𝑥3))/((𝑥1−𝑥0)
(𝑥1−𝑥2)(𝑥1−𝑥3)) f(x1) el tomar un punto de referencia significa ignorarlo
como se puede observar,
en el numerador y utilizarlo en el denominador y utilizar su f(x) para
multiplicar al final.
Al final, si se tiene n numeros para unir, la ecuacion resultante tendra como
grado n-1 asi si se quieren unir 3 puntos, la ecuacion que resulta tiene grado 2.
Faltan 2 puntos, ahora se toma com referencia el x2 y la ecuacion va
quedando de la siguiente manera:
g(x)-((𝑥−𝑥1)(𝑥−𝑥2)(𝑥−𝑥3))/((𝑥0−𝑥1)(𝑥0−𝑥2)(𝑥0−𝑥3)) f(x0) + ((𝑥−𝑥0)(𝑥−𝑥2)
(𝑥−𝑥3))/((𝑥1−𝑥0)(𝑥1−𝑥2)(𝑥1−𝑥3)) f(x1) + ((𝑥−𝑥0)(𝑥−𝑥1)(𝑥−𝑥3))/((𝑥2−𝑥0)
(𝑥2−𝑥1)(𝑥2−𝑥3)) f(x2)
Por ultimo se toma como referencia el siguiente punto que es x3 y se termina de
escribir la ecuacion que queda de la siguiente manera:
interpolacion por el metodo de lagrange pdf
Interpolacion
AGUA
temp °c densidad relativa kg/m3 viscocidad cinematica m2/seg
5 1 1.52
10 1 1.308
15 0.999 1.142
20 0.998 1.007
25 0.997 0.897
30 0.995 0.804
35 0.993 0.727
40 0.991 0.661
50 0.99 0.556
65 0.98 0.442
se desea conocer la densidad relativa y la viscocidad
cinematica cuando la temperatura sea de 27.5°c.
colocacion del valor a
interpolar
valores interpolados
correspondientes
𝑌𝑛=𝑌𝑜+ (𝑋−𝑋𝑜)/(𝑋−𝑋
(Y1 - Yo)
Formula para
interpolacion en
excel
tabla para valores necesarios para realizar la interpolacion
acion
informacion o datos
iniciales
viscocidad
e 27.5°c.
temp densidad viscosidad
27.5
+ (𝑋−𝑋𝑜)/(𝑋−𝑋𝑜)
(Y1 - Yo)
temp densidad viscocidad
Lim inferior
Lim superior
valor
en las celdas escribir"=BUSCAR (valor_buscad
resultados
temp densidad viscocidad
27.5 0.996 0.8505
en esta celda escribir"= MULTIPLO. INFERIOR(numero,cifra_significaiva)"
cifra_significativa=es la
diferencia que hay entre x y y
cada dato de la temperatura temp densidad viscocidad
lim inferior 25 0.997 0.897
lim superior 30 0.995 0.804
valor 27.5 0.996 0.8505
en esta celda escribir "=F16+5"
introducir formula
𝑌𝑛=𝑌𝑜+
(𝑋−𝑋𝑜)/(𝑋1−𝑋𝑜) (
ir"=BUSCAR (valor_buscado, valor_de_comparacion,[vector_resultado])"
valor_buscando=F16 y/o F17
tempdensidad
°c relativa kg/m3
5 1
valor_de_comparacion = datos de la columna de temp
10 1
15 0.999
vector_resultando = datos de la columna de la densidad 20 0.998
25 0.997
30 0.995
Realiza en esta s celdas, las mismas operaciones 35 0.993
realizadas en la columna de la densidad; pero con 40 0.991
la columna de la viscocidad. 50 0.99
65 0.98
ntroducir formula
𝑌𝑛=𝑌𝑜+
𝑜)/(𝑋1−𝑋𝑜) (Y1 - Yo)
viscocidad cinematica m2/seg
1.52
1.308
1.142
1.007
0.897
0.804
0.727
0.661
0.556
0.442
geopresiones
profundidad (ft) presion de formacion lb/in2 presion de sobrecarga lb/in 2
9500 4273.1 3835
9750 7605 9106.5
10000 6999.2 9380
10250 8453.4 9655.5
10500 9134.58 9933
10750 9670.7 10212.5
11000 10336 10494
12000 11519.64 11446
determinar las tres presiones geologicas cuando se esta a una
profundidad de:
a) 10400 ft
b) 11500 ft
c) 10600 ft
lim inferior
lim superior
resultados valor
profundidad [Link] psc p fract
10400 4959.61432 4556.04210526316 8335.69157894737
presion de fractura lb/in2
8004.7
8898.2
8999.3
9427.7
𝑌𝑛=𝑌𝑜+
9819.6 (𝑋−𝑋𝑜)/(𝑋1−𝑋𝑜) (Y1 - Yo)
10141.5
10474.8
11498.5
esta a una
x y y y
profundidad presion de formacion presiom de sobrecarga presion de fractura
9500 4273.1 3835 8004.7
19000 11519.64 11446 11498.5
10400 4959.61431578947 4556.04210526316 8335.69157894737
volumen especifico m3/kg energia interna kj/kg
temp °c presion de saturada kpa liquido saturado vapor saturado liquido saturado evaporacion
0.01 0.6117 0.001 206 0 2374.9
5 0.8725 0.001 147.03 21.019 2360.8
10 1.2281 0.001 106.32 42.02 2346.6
15 1.7057 0.001001 77.885 62.98 2332.5
20 2.3392 0.001002 57.762 83.913 2318.4
25 3.1698 0.001003 43.34 104.83 2304.3
30 4.2469 0.001004 32.879 125.73 2290.2
35 6.6291 0.001006 25.205 146.63 2276
40 7.3851 0.001008 19.515 167.53 2261.9
se desea conocer las propiedades termodinamicas de volumen especifico,
energia interna y entalpia del estado del agua a 38°c
38 °c lim inferior
lim superior
valor
𝑌𝑛=𝑌𝑜+
(𝑋−𝑋𝑜)/(𝑋1−𝑋𝑜) (Y1 - Yo)
a interna kj/kg entalpia kj/kg
vapor saturado evaporacion vapor saturado
2374.9 2500.9 2500.9
2381.8 2489.1 2510.1
2388.7 2477.2 2519.2
2395.5 2465.4 2528.3
2402.3 2453.5 2537.4
104.83 2441.7 2546.5
125.74 2429.8 2555.6
146.64 2417.9 2564.6
167.53 2406 2573.5
temp p. saturada kpa liquido saturado vapor saturado liquido saturado evaporacion
35 6.6291 0.001006 25.205 146.63 2276
40 7.3851 0.001008 19.515 167.53 2261.9
38 7.0827
vapor saturado evaporacion vapor saturado
146.64 2417.9 2564.6
167.53 2406 2573.5
en la siguiente tabla se indican los limites maximos de solidos recomendados por la
asiciacion de fabricantes de calderas de los estados unidos
presion salida de la caldera (psi) solidos totales alcalinidad total (PPM)
225.5 3500 700
375.5 3000 600
525.5 2500 500
675.5 2000 400
825.5 1500 300
950.5 1250 250
1250.5 1000 200
determinar los limites maximos de solidos cuando la presion de salida de la caldera
es de:
a) 300 psi
b) 380.5 psi
c) 790 psi
d) 1000 psi
mos de solidos recomendados por la
dos unidos
solidos suspendidos PPM silice PPM
300 125 𝑌𝑛=𝑌𝑜+
250
150
90
50
(𝑋−𝑋𝑜)/(𝑋1−𝑋𝑜) (Y1 - Yo)
100 35
60 20
40 8 150
20 2.5
de la caldera
x y
presion salida de la caldera (psi) solidos totales
lim inferior 300 3500
lim superior 450 3000
valor 300 3500
resultados
presion salida de la caldera (psi) solidos totales
300 3500
=𝑌𝑜+
−𝑋𝑜) (Y1 - Yo)
y y y
alcalinidad total (PPM) solidos suspendidos PPM silice PPM
700 300 125
600 250 90
700 300 125
resultados
alcalinidad total (PPM) solidos suspendidos PPM silice PPM
700 300 125
valor a b c d
22 2 50 3.5 0.01
28.5 6 45 4.8 0.5
35 10 42.5 5.6 1.7
41.5 14 30 5.9 1.95
48 20 28 6.7 2.34 6.5
54.5 26 21 7.1 2.9 6.5
61 32 16 7.7 4.8
67.5 40 10 9.2 5.85
a) 30
b)
c)
45.6
58.7
𝑌𝑛=𝑌𝑜+
d) 64.35 (𝑋−𝑋𝑜)/(𝑋1−𝑋𝑜) (Y1
x y y y y
valor a b c d
lim inferior 26 2 50 3.5 0.01
lim superior 32.5 6 45 4.8 0.5
valor 45.6 14.0615385 34.9230769 7.42 1.48753846
𝑌𝑛=𝑌𝑜+
𝑜)/(𝑋1−𝑋𝑜) (Y1 - Yo)
el siguiente conjunto de mediociones de temperatura tomadas de la
cabeza de cilindro de un motor nuevo que se esta probando con miras a
usase en un auto de carreras
°C
Tiempo tem1 tem2 tem3
0 0 0 0
1 20 25 52
2 60 62 90
3 68 67 91
4 77 82 93
5 110 103 96
resultado
valor a b c
566 4572 2882 654
Determinar las temperaturas cuando an transcurrido:
a) 2.3 seg
b) 3.8 seg
c) 4.25 seg
x y y y
valor a b c
lim inferior 2 60 62 90
lim superior 3 68 67 91
valor 566 4572 2882 654
