Universidad Nacional Mayor De San Marcos

(Universidad Del Perú, DECÁNA DE AMÉRICA)

FACULTAD DE INGENIERÍA GEOLÓGICA MINERA METALÚRGICA Y GEOGRÁFICA

Escuela Profesional de Ingeniería Metalúrgica

ANÁLISIS NUMÉRICO

TRABAJO GRUPAL

Profesora:

Ing. Melgarejo Estremadoyro Gladys Giovanna

Alumnos:

• Collachagua Giron Jhonny 17160238

• Julca Vera Antony Brandi 17160264
• Obregón Aguilar Jose Luis 16160225
• Rimache Arotoma Jonathan 16160248

LIMA - PERÚ

2020 – I

[ANÁLISIS NUMÉRICO] UNMSM

2
 Método de punto fijo
Una chancadora secundaria produce 27% de partículas finas -3u , trabajando a una grado de
reducción de 5. Se ha muestreado al alimento. Determinando que obedece a una distribución
3
granulométrica dada por F(x) = x + cos2x – 14 ,con el valor inicial X0 = 0.1 y considerando la
función de iteración (14 − cos2x)1/3 e iterar hasta que el error estimado sea menor que
0.01%, trabajar con 6 dígitos en la parte decimal

3
Tenemos que f(x) = x + cos2x – 14 = 0 → x = (14 − cos2x)
→ g(x) = (14 − cos2x)

Tenemos : g(x) = (14 − cos2x)

X0 = g(X0) = 0.1

X1 = g(X0) = g(0.1) = ( 14 − cos2(0.1)) = 2.352536

X2 = g(X1) = g(2.352536) = (14 − cos2(2.352536)) = 2.410562

X3 = g(X2) = g(2.410562) = (14 − cos2(2.410562)) = 2.403899

X4 = g(X3) = g(2.403899) = (14 − cos2(2.403899)) = 2.404663

X5 = g(X4) = g(2.404663) = (14 − cos2(2.404663)) = 2.404575

i xi Es=|Vactual – Vanterior|
0 0.1
1 2.352536 2.252536
2 2.410562 0.058026
3 2.403899 0.006663
4 2.404663 0.000764
5 2.404575 0.000088 < 0.01%

Entonces el valor de la raíz X5 = 2.404575 u, ya el error estimado “Es” 0.000088 < 0.01%
 Método de punto fijo
En la empresa Antamina es una empresa minera metalúrgica polimetálico y extrae sus recursos
a tajo abierto y produce concentrados de cobre, zinc y molibdeno y observamos que la
recuperación del mineral valioso en la flotación de minerales sulfurados está en función del pH
del concentrado utilizando el método del punto fijo para poder localizar la raíz que esta en
2
función f(x) = x – eX + 10x – 4 con un valor inicial X0 = 1 e iterar hasta que el error estimado
sea menor o igual 0.1% y trabajar con 6 cifras en la parte decimal ,considerar que la función de
iteración es g(x) = (eX - x2 + 4)/10

Tenemos que f(x) = x2 – eX + 10x – 4 = 0 → x = (eX - x2 + 4)/10

2
→ g(x) = (eX - x + 4)/10

Tenemos : g(x) = (eX - 3) 1/3

X0 = g(X0) = 1.0

2
X1 = g(X0) = g(1.0) = (e1 - 1 + 4)/10 = 0.571828

0.571828
X2 = g(X1) = g(0.571828) = (e - 0.5718282 + 4)/10 = 0.544452

0.544452
X3 = g(X2) = g(0.544452) = (e - 0.5444522 + 4)/10 = 0.542724

0.542724
X4 = g(X3) = g(0.542724) = (e - 0.5427242 + 4)/10 = 0.542614

i Xi Es=|Vactual – Vanterior|

0 1

1 0.571828 0.428172

2 0.544452 0.027376

3 0.542724 0.001728

4 0.542614 0.000110

Entonces el valor de la raíz X4 = 0.542614 , ya el error estimado “Es” 0.000110 ≤ 0.001

 INTERPOLACIÓN NUMÉRICA

Problema 1:

La dureza de un acero 1040 (0.40%C – 0.30%Si – 0.70%Mn) disminuye en función de la

temperatura de revenido. Se dispone de los siguientes datos:
Puntos 0 1 2 3 4
Temperatura 200 300 400 500 600
(°C)
Dureza (HB) 410 338 260 231 220
brinell

TABLA DE DIFERENCIAS DIVIDIDAS FINITAS

i xi f (x i ) 1eras 2das 3ra 4ta 5ta
0 200 410
-0.720000
1 300 338 -0.000300
-0.780000 0.000009
0.000000
2 400 260 0.002450 0
-0.290000 -0.000005
3 500 231 0.000900
-0.110000
4 600 220

“DIFERENCIAS
DIVIDIDAS DE ORDEN CERO”

PRIMERAS DIFERENCIAS

f (x 0, x 1) = f (x 1) - f (x 0) = 338 – 410 = -0,720000

x 1- x 0 300 - 200

f (x 1, x 2) = f (x 2) - f (x 1) = 260 - 338 = -0,780000

x 2- x 1 400 - 300

f (x 2, x 3) = f (x 3) - f (x 2) = 231 - 260 = -0,290000

x 3- x 2 500 - 400
SEGUNDAS DIFERENCIAS

f (x 0, x 1, x 2) = f (x 1, x 2) - f (x 0, x 1) = -0,780000 – (-0,720000) = -0,000300

x 2- x 0 400 - 200

f (x 1, x 2, x 3) = f (x 2, x 3) - f (x 1, x 2) = -0,290000 – (-0,780000) = 0,002450

x 3- x 1 500 - 300

TERCERAS DIFERENCIAS

f (x 0, x 1, x 2, x 3) = f (x 1, x 2, x 3) - f (x 0, x 1, x 2) = 0,002450 - (-0,000300) = 0,000009

x 3- x 0 500 - 200

f (x 1, x 2, x 3, x 4) = f (x 2, x 3, x 4) - f (x 1, x 2, x 3) = 0,000900 – 0,002450 = -0,000005

x 4- x 1 600 - 300
Problema 2:

La resistencia a la tracción de un Acero 1040 disminuye en función de la temperatura de

revenido. Se tiene los siguientes datos obtenidos experimentalmente.
Puntos 0 1 2 3 4
Temperatura 300 400 500 600 700
(°C)
Ot 120 89 71 60 51
(resistencia a
la tracción)
Kg/ mm2

TABLA DE DIFERENCIAS DIVIDIDAS FINITAS

i xi f (x i ) 1eras 2das 3ra 4ta 5ta
0 300 120
-0,310000
1 444400444 89889 00,0006500
-0,180000 -0,000001
2 500 71 0,000350 0,000000 0
-0,110000 0,000000
3 600 60 0,000100
-0,090000
4 700 51

“DIFERENCIAS
DIVIDIDAS DE ORDEN CERO”

PRIMERAS DIFERENCIAS
f (x 0, x 1) = f (x 1) - f (x 0) = 89 – 120 = -0,310000
x 1 - x 0 400 - 300

f (x 1, x 2) = f (x 2) - f (x 1) = 71 – 89 = -0,180000
x 2 - x 1 500 - 400

f (x 2, x 3) = f (x 3) - f (x 2) = 60 – 71 = -0,110000
x 3 - x 2 600 - 500
SEGUNDAS DIFERENCIAS

f (x 0, x 1, x 2) = f (x 1, x 2) - f (x 0, x 1) = -0,180000 – (-0,310000) = 0,000650

x 2- x 0 500 - 300

f (x 1, x 2, x 3) = f (x 2, x 3) - f (x 1, x 2) = -0,110000 – (-0,180000) = 0,000350

x 3- x 1 600 - 400

TERCERAS DIFERENCIAS

f (x 0, x 1, x 2, x 3) = f (x 1, x 2, x 3) - f (x 0, x 1, x 2) = 0,000350 - 0,000650 = -0,000001

x 3- x 0 600 - 300

f (x 1, x 2, x 3, x 4) = f (x 2, x 3, x 4) - f (x 1, x 2, x 3) = 0,000100 – 0,000350 = 0,000000

x 4- x 1 700 - 400
 NEWTON HACIA ADELANTE
En la empresa minera southern que comprende un proceso amplio de extracción, molienda y
flotación se determinó que la liberación de partículas en la etapa de molienda estaba dada en
función al tiempo, sabemos que un factor muy importante es la liberación del mineral valioso
de la ganga que se llama “granulometría” para que pueda posteriormente pasar a la etapa de
concentración. a través de la flotación o también lixiviación.
Puntos 0 1 2 3 4
Tiempo 3 6 9 12 15
Liberación de 55 70 85 70 68
partículas(%)

Se considera optimo cuando la liberación se encuentra entre 74 a 86 % por lo que la planta

trabaja en un promedio de 7 horas si por error de manipulación de un operario trabaja el
mineral en tiempo de molienda 6.5 la pregunta del operario será satisfactoria en ese tiempo,
trabajar con 6 cifras en la parte decimal
i Xi f(xi) ∆f(x) ∆2 f(x) ∆3 f(x) ∆4 f(x)

0 3 55
15.000000
1 6 70 0.000000
15.000000 -30.000000
2 9 85 -30.000000 73.000000
-15.000000 43.000000
3 12 70 13.000000
-2.000000
4 15 68

∆f(x0) = f(x1) – f(x0) =70-55=15.000000

∆f(x1) = f(x2) – f(x1) =85-70=15.000000

∆f(x2) = f(x3) – f(x2) =70-85=-15.000000

∆2f(x0) = ∆f(x1)- ∆f(x0) = 15.000000-15.000000=0.000000

∆2f(x1) = ∆f(x2)- ∆f(x1) = -15.000000-15.000000=-30.000000

∆2f(x2) = ∆f(x3)- ∆f(x2) =-2.000000+15.000000=13.000000

∆3f(x0) = ∆2f(x1)- ∆2f(x0) =-30.00000-0.00000=-30.000000

∆3f(x1) = ∆2f(x2)- ∆2f(x1) =13.00000+30.00000=43.000000

∆4 f(x)= ∆3f(x1)- ∆3f(x0) = 43.00000+30.00000=73.000000

Veremos

x1 − x0 = h ⇒ h = 6 − 3 ⇒ h = 3
el valor que tenemos que interpolar es:6.5, entonces se obtiene el valor de “s” Es
decir; x = x0 +sh esto es: x x
𝑆=

6.
𝑆=

S=1.166667

Aproximaremos a un polinomio de 1er grado, entonces se toma los dos primeros

términos, es decir 0 y 1
p(x)=f(x0) +∆f(x0) s
p (6.5) =55+15.000000(1.166667)
p (6.5) =72.500005
vemos que estamos extrapolando, porque el valor de “x” queda fuera del intervalo de
los puntos que se usaron para formar el polinomio de aproximación
Que sucede si realizamos una aproximación mejor con los puntos 1 y 2, para esto
vamos a modificar: p(x) = f(x1) +Δf (x1) s
Los dos primeros términos dan la aproximación polinomial de 1er grado

xx
p(x) = f(x1) +Δf (x1) s 𝑆=

6. p(x)

= 70+15.000000(0.166667) 𝑆
==0.166667 p(x) =72.500005
Ahora se desea aproximar con un polinomio de 2do grado, se requiere 3 puntos, se
aconseja tomar los puntos 1,2 y 3 en lugar de 0,1 y 2, ya que el argumento por
interpolar está más al centro de los primeros.
 x x 6.
2
f(x1) s ,𝑆= = 𝑆 =
=0.166667 p(x) = f(x1) +Δf (x1) s+∆
! 3
p(x) =70+15.000000(0.166667) + (-30.000000) (0.166667(0.166667-1)/2)
p(x) = 74.583342

observamos que es eficiente ya que se encuentra en el intervalo de rendimiento

NEWTON HACIA ADELANTE

la empresa minera las bambas es uno de los mayores productores de concentrados
cobre en el Perú y en el mundo por eso tiene una visión muy importante en el
mejoramiento de los equipos y una mayor eficiencia y recuperación de minerales
valiosos con la finalidad de cumplir con las expectativas del mercado nacional y
mundial mejorando su proceso de recuperación en la fase de flotación de los bancos
de celdas rougher,scavenger y cleaner, se determinó que la recuperación estaba dada
en función de la temperatura.
Puntos 0 1 2 3 4
Tiempo 10 20 30 40 50
Recuperación% 45 70 95 98 98

Considerando que una buena recuperación se da cuando se obtiene entre 85 al 98, se

quiere saber por curiosidad del operario si en el tiempo 25 horas habrá una buena
recuperación trabajar con 6 cifras en la parte decimal
i Xi f(xi) ∆f(x) ∆2 f(x) ∆3 f(x) ∆4 f(x)

0 10 45
25.000000
1 20 70 0.000000
25.000000 -22.000000
2 30 95 -22.000000 41.000000
3.000000 19.000000
3 40 98 -3.000000
0.000000
4 50 98

∆f(x0) = f(x1) – f(x0) =70-45=25.000000

∆f(x1) = f(x2) – f(x1) =95-70=25.000000

∆f(x2) = f(x3) – f(x2) =98-95=3.000000

∆2f(x0) = ∆f(x1)- ∆f(x0) = 25.000000-25.000000=0.000000

∆2f(x1) = ∆f(x2)- ∆f(x1) = 3.000000-25.000000=-22.000000

∆2f(x2) = ∆f(x3)- ∆f(x2) =0.000000-3.000000=-3.000000

∆3f(x0) = ∆2f(x1)- ∆2f(x0) =-22.00000-0.00000=-22.000000

∆3f(x1) = ∆2f(x2)- ∆2f(x1) =-3.00000+22.00000=19.000000

∆4 f(x)= ∆3f(x1)- ∆3f(x0) = 19.00000+22.00000=41.000000

Veremos

x1 − x0 = h ⇒ h = 20 − 10 ⇒ h = 10
el valor que tenemos que interpolar es:25, entonces se obtiene el valor de “s” Es decir;
x = x0 +sh esto es:
𝑆=

𝑆=1.500000

Aproximaremos a un polinomio de 1er grado, entonces se toma los dos primeros

términos, es decir 0 y 1
p(x)=f(x0) +∆f(x0) s
p (25) =45+25.000000(1.500000)
p (6.5) =82.500000
vemos que estamos extrapolando, porque el valor de “x” queda fuera del intervalo de
los puntos que se usaron para formar el polinomio de aproximación
Que sucede si realizamos una aproximación mejor con los puntos 1 y 2, para esto
vamos a modificar: p(x) = f(x1) +Δf (x1) s
Los dos primeros términos dan la aproximación polinomial de 1er grado

x x
p(x) = f(x1) +Δf (x1) s 𝑆= p(x) =

70+25.000000(0.500000) 𝑆 = =0.500000 p(x)

=82.500000
Ahora se desea aproximar con un polinomio de 2do grado, se requiere 3 puntos, se
aconseja tomar los puntos 1,2 y 3 en lugar de 0,1 y 2, ya que el argumento por
interpolar está más al centro de los primeros.

x x
f(x1)
2
s , 𝑆= = 𝑆
==0.500000 p(x) = f(x1) +Δf (x1) s+∆
! 10
p(x) =70+25.000000(0.500000) + (-22.000000) (0.500000(0.500000-1) /2)

p(x) = 85.250000 observamos que es eficiente ya que se encuentra en el

intervalo de rendimiento