Universidad Nacional de Ingeniería

Facultad de Ingeniería Civil 21-04-2013

Sección Post Grado

MÉTODOS NUMÉRICOS EN INGENIERÍA

Primer Trabajo 2013-1

1. Prepare una hoja de cálculo en Excel para resolver un sistema de 3 ecuaciones lineales, con coeficientes
reales, por eliminación directa. Suponga que no se requerirán intercambios de filas y/o columnas.
Verifique el correcto funcionamiento empleando Matlab, Mathcad, o con un ejemplo resuelto en algún
texto.

2. Utilice la hoja de cálculo para resolver el sistema de ecuaciones Ax b por eliminación directa:
1 1 1 x1 1
2 3

1 1 1 x2 0
2 3 4

1 1 1 x3 0
3 4 5

3. Luego repita nuevamente todas las operaciones, pero trabajando con sólo tres cifras significativas. Esto
significa que en lugar de 1 3 debe escribir 0.333. Al hacer operaciones como

debe escribir:

El resultado exacto sería 112 0.08333 ... Hay un error como resultado de la aritmética limitada a tres
cifras significativas. No agregue más cifras “ayudando” a la hipotética computadora, ni en este ni en los
pasos siguientes. Compare y comente los resultados

4. Obtenga la inversa de la matriz de coeficientes, A 1 . Para ello puede utilizarse la misma hoja de
cálculo, considerando en el segundo miembro del sistema de ecuaciones las columnas de la matriz
identidad. Alternativamente, puede emplear Matlab o Mathcad.
5. Determine el número de condicionamiento (A ) . Comente sobre las soluciones antes obtenidas.

6. Haga un pequeño cambio en un coeficiente de la matriz A y observe como cambia la inversa. Comente
los resultados. Llas operaciones pueden ser realizadas empleando nuevamente la hoja de cálculo o
utilizando la fórmula de Shermann-Morrison:
(A u v T ) 1
A 1
A 1u v T A 1

siendo (1 v T A 1u) 1 . Puede verificar fácilmente la validez de la fórmula efectuando el producto:

(A u v T ) (A u v T ) 1
I

El trabajo no deberá exceder de 10 páginas. Se espera un documento parcialmente manuscrito, con

algunos resultados listados a partir de los programas de cómputo cuando eso resulte más simple. Será
entregado el sábado 27 de abril de 2013.
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Civil 27-04-2013
Sección Post Grado

MÉTODOS NUMÉRICOS EN INGENIERÍA

Segundo Trabajo 2013-1

1. Escriba una hoja de cálculo para resolver un sistema de cuatro ecuaciones lineales con un procedimiento
de relajación. Utilice la hoja de cálculo para resolver un sistema de ecuaciones lineales en el que la
diagonal principal de la matriz de coeficientes sea fuertemente dominante. Utilice la misma hoja de
cálculo para resolver otro sistema de ecuaciones, en el que los coeficientes de la diagonal principal sean
menores que la suma de los valores absolutos de los otros coeficientes de la correspondiente fila. Sólo
presente los cómputos iniciales y los finales. Comente sobre la convergencia.

2. Se han medido los tres ángulos internos de un triángulo plano, 1 2 3, teniéndose los valores
(aproximados) A, B, C . Suponga que A B C . Sin embargo, 1 2 3 . Esta expresión
es exacta, por lo que hay sólo 2 valores independientes y se puede escribir:
1 A
2 B
1 2 C
Determine los valores de 1 2 3 con un criterio de mínimos cuadrados. ¿Deben los errores en cada
ángulo corregirse en proporción con los valores medidos? o ¿Las correcciones en cada ángulo deben
ser iguales?
Suponga luego que se tiene un triángulo rectángulo, cuyos catetos son L1 , L2 y cuya hipotenusa es H .
Al medir los tres lados se obtienen valores aproximados A, B, C respectivamente, pero A 2 B 2 C 2 .
Para corregir las mediciones, considere x1 L12 y x 2 L22 y por tanto:
x1 A2
x2 B2
x1 x2 C2

Determine x1 y x 2 con un criterio de mínimos cuadrados y a partir de éstas obtenga valores corregidos
para L1 , L2 , H .

3. (Adaptado de S. M. Stiegler. “Gauss and the invention of least squares”. Ann. Statist., 9:465–474, 1981)
En 1793 se decidió establecer el nuevo sistema métrico sobre la base de una unidad, el metro, igual a la
diez millonésima parte de la distancia del polo Norte al Ecuador medida a lo largo del meridiano que
pasa por París. En 1795 se publicaron resultados de mediciones en cuatro sectores, desde Dunkerque
hasta Barcelona. Para cada sector se midió la longitud del arco (con referencia a un cierto patrón) S , la
diferencia de latitudes en grados d , y la latitud al centro L (determinada con observaciones
astronómicas):

Sector Longitud de arco S Diferencia d Latitud L

Dunkerque a Pantheon 62472.59 2.18910º 49º 56’ 30”
Pantheon a Évaux 76145.74 2.66868º 47º 30’ 46”
Évaux a Carcasona 84424.55 2.96336º 44º 41’ 48”
Carcasona a Barcelona 52749.48 1.85266º 42º 17’ 20”

Suponiendo que el meridiano es un elipsoide, podría plantearse la aproximación:

 z y sen 2 L S/d
en la que z e y son parámetros desconocidos, expresados en unidad de longitud (con el patrón de
referencia) por grado. La longitud de un cuadrante del meridiano resulta entonces: M 90 ( z y / 2) y
1
la excentricidad e se obtiene mediante: e 3 (z / y 0.5) . Determine z e y por mínimos cuadrados
y luego M y e . Note que la constante 90 está en grados.

4. (Adaptado de Dalquist y Björck, Numerical Methods. 2008) Se tienen m pares de valores ( xi , yi ) para
los que se postula la relación lineal y x . Si x e y son los valores medios de los x i e y i ,
respectivamente:
m
xi y i mxy
i 1
m
xi2 mx 2
i 1

luego de lo cual: y x
Un procedimiento más preciso consiste en restar primero los valores medios de los datos, escribiendo:
(y y) (x x)
m
( xi x ) ( yi y)
i 1
m
( xi x)2
i 1

Suponga que se tienen los valores:

t 1001 1003 1004 1006 1007
f (t ) -3 -1 -0.5 0.5 1

Para los que se postula: f t o bien ( f f ) (t t ) . Obtenga trabajando primero con

no menos de 10 cifras significativas. Luego emplee las dos alternativas antes mencionadas, pero
trabajando en ambos casos con sólo 4 cifras significativas. Comente los resultados.

Este trabajo no deberá exceder de 10 páginas. Deberá ser entregado a más tardar el sábado 11 de
mayo de 2013.
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Civil 11-05-2013
Sección Post Grado

MÉTODOS NUMÉRICOS EN INGENIERÍA

Tercer Trabajo 2013-1

1. Prepare hojas de cálculo para resolver el problema de valores y vectores característicos A B , por
iteración directa y por iteración inversa con traslación. Suponga que las matrices son de orden 4 y que
B es diagonal. Verifique las hojas de cálculo, obteniendo los vectores característicos correspondientes a
los valores característicos máximo y mínimo de un problema de la forma:

a b 0 0
c a b 0
Φ Φ
0 c a 
0 0  

Para el cual:
1
c 2

sen s
b
1 c
c 2 sen 2 s s
s a 2b cos s y s
b siendo s
b 3 n 1
c 2

sen 3 s
b


donde n es el orden de la matriz y s un entero en el rango entre 1 y n .

2. Para la misma matriz tridiagonal considerada en el ítem 1, determine otro vector característico,
correspondiente a un valor característico intermedio.

3. Determine los valores característicos máximo y mínimo de la matriz:

4 1 1 1
1 4 1 1
A
1 1 5 1
1 1 1 3

4. Considerando la misma matriz A , determine los valores característicos máximo y mínimo de la matriz:
1 1
H D 2 A D 2 donde D es una matriz diagonal que contiene a los elementos de la diagonal principal
de A . Note que hij aij / aii a jj

5. Suponga que se tiene un sistema de ecuaciones de la forma: A x b con la misma matriz de

coeficientes considerada en el ítem 3. Note que la matriz A es simétrica. Suponga valores para las x y
calcule los correspondientes valores b , o alternativamente suponga valores b arbitrarios y determine
las correspondientes x empleando una hoja de cálculo en Excel, Matlab o Mathcad. Luego considere la
aproximación inicial x 0 0 y efectúe 4 pasos del método de gradiente conjugada. Compare los
resultados aproximados obtenidos en el tercer paso con los exactos.
6. Repita lo hecho en el ítem 4, pero con un procedimiento de gradiente conjugada pre condicionado, con
pre condicionador diagonal, que sería equivalente a resolver D 1 Ax D 1b o alternativamente
~ ~
H~ b con ~
1 1
x x D 2x y b D 2
b:
Con la aproximación inicial x 0 0 y por lo tanto r0 b A x0 b determine s 0 p0 D 1r0 e itere
para k 0,1, 2 n 1:
qk A sk
p Tk rk
k
s Tk q k
xk 1 xk k sk
rk 1 rk k qk
pk 1 D 1 rk 1

p Tk 1rk 1
k
p Tk rk
sk 1 pk 1 k sk

Compare los resultados aproximados obtenidos en el tercer paso con los exactos.

7. A partir de los resultados obtenidos en los ítems 3 y 4, determine los números de condicionamiento
2 ( A) y 2 ( H) . Luego, considerando que en un proceso de gradiente conjugada se observa que:
k
(k ) 1
x x 2 x ( 0) x
1

Comente sobre los resultados obtenidos en los ítems 4 y 5.

Este trabajo será entregado el sábado 25 de mayo de 2013. No deberá exceder de 20 páginas. Deberá
presentarse por escrito e/o impreso; no se aceparán trabajos presentados en disco o enviados por
correo, salvo circunstancias excepcionales debidamente justificadas.
Universidad Nacional de Ingeniería
Facultad de Ingeniería Civil 25-05-2013
Sección Post Grado

MÉTODOS NUMÉRICOS EN INGENIERÍA

Cuarto Trabajo 2013-1

1. Proponga una ecuación no lineal, preferentemente una que corresponda a su especialidad, y

determine una de sus raíces empleando el método de Isaac Newton (1642-1727).

2. Resuelva el sistema de ecuaciones no lineales:

x2 10 x y2 8 0

x y2 x 10 y 8 0
Se sugiere iterar con ecuaciones escritas en la forma x f ( x, y) e y g ( x, y)

3. Prepare una hoja de cálculo que permita obtener las raíces de un polinomio de cuarto grado por el
método de Edmond Laguerre (1834-1886) y utilícela para determinar las raíces de un polinomio.
Experimente con distintas aproximaciones iniciales y presente sus conclusiones.

4. En 1225 Leonardo Fibonacci (1175-1240) propuso la serie de números definida por la ecuación de
yk 1
diferencias: yn yn 1 yn 2 con y0 0 e y1 1 . Determine Lim (la “proporción áurea”).
k yk

5. Una deuda P0 A se amortiza en m cuotas mensuales iguales, R . La cantidad adeudada después

de k 1 meses es Pk 1 (1 i) Pk R , donde i es el interés mensual. Esta es una ecuación de
diferencias (no homogénea). Determine el monto de las cuotas para que la deuda se pague en
exactamente m meses Pm 0 . Proponga luego valores numéricos para las distintas variables y
determine el porcentaje que representa el monto total de los intereses en proporción a la deuda
inicial.
4 dx
6. Evalúe la integral: por el método de Werner Romberg (1909-2003), a partir de los
o 1 x2
resultados obtenidos con la regla de los trapecios con intervalos h iguales a 4, 2, 1, 0.5 y 0.25.

7. Para la misma integral, utilice el método de Gauss Legendre con distintos números de puntos de
integración y compare los resultados.

8. Se pretende desarrollar una fórmula de integración numérica del tipo:

1 1
f ( x, y) dx dy w0 f (0,0) w1{ f ( a, a) f (a, a) f ( a, a) f (a, a) }
1 1

Determine las relaciones entre a, w0 , w1 para que la expresión integre exactamente funciones
polinómicas del mayor grado posible.

El trabajo no deberá tener más de 15 páginas. Será entregado el día sábado 08 de junio de 2013.
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Civil 01-06-2013
Sección Post Grado

MÉTODOS NUMÉRICOS EN INGENIERÍA (C 502)

Examen Parcial 2013-1

1 1
0
1 1 1 2 2
1 1 1 1
1. Dada la matriz: A 1 y su inversa: A
2 4 4
1 1 1 1
2 4 4
Determine el número de condicionamiento (A) suponiendo que 0 1.

2. Dadas las mediciones:

x 2 3 4 5 6 8
y 4.46 7.82 10.76 15.59 18.37 23.89
Se proponen las expresiones:
y ax
y b cx
Determine a, b, c con el criterio de mínimos cuadrados y en cada uno de los dos casos calcule el error
total E r2
1 2 1 0
3. Determine los valores y vectores característicos de φ φ
2 8 0 4
usando el método de Jacobi (puede ser convirtiendo primero a la forma clásica o alternativamente con el
método generalizado).

4. Dados:
x 0.5 1.0 1.5
f (x) 2.625 0.750 -0.625
Determine por interpolación una raíz de f (x) 0 con la mayor precisión posible.

5. Utilice el método de Newton para determinar una de las raíces de la ecuación x cosx 0 . Haga
primero una gráfica de las funciones y cos x e y x para determinar una aproximación inicial.

3
6. Determine una raíz del polinomio x 2 x 5 0 en el intervalo [1,5] utilizando el método de Laguerre.
Luego determine un polinomio de grado 2 cuyas raíces sean las dos raíces restantes del polinomio
original.

7. Considere la ecuación de diferencias: yk 2 2 yk 1 yk 1. Determine la solución de la

correspondiente ecuación homogénea. Verifique que la solución particular es: y k c k 2 . Obtenga c y
luego determine la solución general sabiendo que y 0 y100 0.

2 sen x
8. Evalúe la integral: dx por el método de Gauss Legendre, con 3 puntos de integración.
0 1 x
Los ceros del polinomio de Legendre de grado 3 son 0.6 , 0 y 0.6 , a los que corresponden
  
pesos 0.5 , 0.8 y 0.5 respectivamente.

Tiempo: 2 horas. Con libros y/o apuntes. No se permite el préstamo de cualquier clase de material.
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Civil 15-06-2013
Sección Post Grado

MÉTODOS NUMÉRICOS EN INGENIERÍA (C 502)

Examen Parcial (R) 2013-1

17 18 1 448 504
1. Verifique que la inversa de la matriz simétrica A es A .
18 19 504 576
Determine los valores característicos de A y luego los números de condicionamiento 2 ( A) y (A) .

100 1 x1 3
2. Resuelva: utilizando un procedimiento iterativo, con la aproximación inicial
1 100 x2 2
x1 0
. Obtenga x1 y x 2 con un error relativo menor que 0.1%.
x2 0

1 1 1
x1
3. Resuelva el sistema sobre determinado de ecuaciones lineales: 0 0 en el que 1
x2
0 0
m
Si 10 ( m es un entero) ¿Cuál sería el mínimo número de cifras significativas con el que debería operarse
para que sea factible obtener una solución?

4 1 x1 1 0 x1
4. Determine los valores y vectores característicos del sistema:
1 3 x2 0 3 x2
utilizando el método que considere más conveniente.

4 1 2
5. Dada A 1 1 1 determine una matriz de la forma Hessemberg con los mismos valores propios.
2 1 4

6. Considere la aproximación inicial: x 0 1.5 y determine una mejor aproximación a una raíz de:
4
x x 8 efectuando 2 pasos del método de Laguerre.
n
7. La suma de n términos de una progresión geométrica puede expresarse como: y n aq k
k 0
n
Por lo tanto: yn yn 1 aq . Esta es una ecuación de diferencias. Suponiendo que la razón q es
distinta de 1 (en cuyo caso la solución sería obvia) obtenga la solución particular de la ecuación. Luego,
sabiendo que y 0 a , obtenga la solución general para y n .
1 x2
8. Evalúe: dx por el método de Gauss – Tchebycheff, con 2 puntos de integración.
1
1 x2

Tiempo: 2 horas
Con libros y/o apuntes. No se permite el préstamo de cualquier clase de material.
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Civil 08-06-2013
Sección Post Grado

MÉTODOS NUMÉRICOS EN INGENIERÍA

Quinto Trabajo 2013-1

1. Al resolver una ecuación diferencial de la forma y' f ( x, y) , un método de Runge-Kutta de cuarto orden
(global) requiere cuatro evaluaciones de la función f en cada paso, mientras que el método de Euler
requiere sólo una. Por tanto, si el método de Runge-Kutta de cuarto orden es mejor, debe producir con
intervalo h resultados más precisos que el método de Euler con intervalo h / 4 . Similarmente, debe
producir resultados mejores que el método de Heun (u otro Runge-Kutta de segundo orden) con intervalo
h / 2 . Compare con la solución exacta las soluciones obtenidas para la ecuación diferencial que le
corresponda según la tabla de variantes, empleando los métodos de Euler, Runge-Kutta de segundo
orden y Runge-Kutta de cuarto orden, con intervalos h / 4 , h / 2 y h , respectivamente.

En la tabla siguiente, n es el último dígito de su código UNI. Si no tiene un código UNI, considere n 0.

n Ecuación diferencial Condición inicial Intervalo Solución exacta

0 y' xy2 y(0) 2 0 x 4 y 2 /(1 x 2 )

x
1,2,3 y' y(0) 1
y
0 x 4 y 1 x2

4,5,6 y
y' 1 y(1) 2 1 x 5 y 2x x ln( x)
x
7,8,9 y' 2 xy y(0) 1 x2
0 x 4 y e

2. Considere la ecuación diferencial: y ' 1y ) 2 con condición inicial y(2) 1 . Suponga que se ha
(x
obtenido y(2.1) 1.19090909 Determine y (3) empleando un proceso Predictor – Corrector, con
x 0.1

3. Un perro, en el campo, observa que su amo corre a lo largo de un camino recto con velocidad constante
v =2m/s. El hombre está inicialmente en coordenadas (0,0) y se desplaza hacia (0,100) . La posición
inicial del perro es (200,0) . El perro corre hacia su amo siempre en la dirección en la que lo observa,
con velocidad u , de donde resulta su trayectoria dada por la ecuación:

xy c 1 ( y )2
donde c v / u 0.4 . Proponga un procedimiento numérico para determinar la trayectoria y prepare una
hoja de cálculo que permita determinar el punto de encuentro y el tiempo transcurrido.

Este trabajo será entregado el día sábado 22 de junio de 2013.

 Universidad Nacional de Ingeniería
Facultad de Ingeniería Civil 22-06-2013
Sección Post Grado

MÉTODOS NUMÉRICOS EN INGENIERÍA

Sexto Trabajo 2013-1

u 3 1 u 1
1. Para resolver las ecuaciones diferenciales
v 1 3 v 2
¿Puede usarse el método de Euler? ¿Y el método de Euler inverso? ¿Y el método de Crank Nicholson? En cada
caso afirmativo, ¿Cuál sería el máximo x para que el procedimiento sea estable?
Prepare una hoja de cálculo para resolver las ecuaciones con uno de los procedimientos antes mencionados y
compare las soluciones con distintos x.

u 1
5 3 u 0
2. La solución de las ecuaciones 2
puede basarse en la descomposición modal:
v 3 5 v 0
u 1 1
f ( x) g ( x)
v 1 1
f (x) y g (x) . Suponga que u (0) 1 ; v(0) 0 y
Escriba ecuaciones diferenciales que permitan determinar
u (0) v (0) 0 . Luego, recordando que la solución de y a2 y 0 es de la forma
y A cos ax B sen ax , determine las funciones u (x) y v(x) .

3. Determine la ecuación de Euler – Lagrange del funcional:

L
P (v) ( 12 EI (v ) 2 1
2
kv2 qv) dx mínimo
0
con las restricciones v(0) v( L) 0 y v (0) v ( L) 0 . Suponga que EI , k , q son constantes. La expresión
precedente es la energía potencial de una viga sobre cimentación elástica; la ecuación de Euler – Lagrange es la
correspondiente ecuación de equilibrio. Luego considere la aproximación:
x 3 x
v( x) a1 sen a 2 sen
L L
y obtenga los parámetros a1 , a 2 (en función de EI , k , q ) empleando el método de Rayleigh Ritz. Suponiendo
valores numéricos para EI , k , q grafique a la misma escala los dos modos componentes de la solución.

2 2
v v
4. La ecuación:
2
0 , con condiciones iniciales v( x,0) A sen ( x / L) y v(x,0) 0 y condiciones de
T t x2
borde v(0, t ) v( L, t ) 0 , describe la vibración libre de una cuerda. Suponga la aproximación:
v( x, t ) c(t ) sen ( x / L) . Esta expresión cumple las condiciones de borde. Obtenga una ecuación diferencial
ordinaria para c (t ) con el método de Galerkin.

b
5. Determine la ecuación de Euler-Lagrange del funcional: I ( y) [1 p( x)( y ) 2 1
q ( x) y 2 r ( x) y] dx
a 2 2

donde las funciones p(x) , q (x) , r (x) son conocidas. Si y e y no fueran conocidas en ambos extremos, ¿Qué
condiciones de borde debería cumplir la función que hace estacionario el funcional?
Suponga que el intervalo entre a y b se divide en subregiones o elementos (finitos), para cada uno de los cuales
s s
se hacen las aproximaciones: u (1 ) ui uj
l L
Con lo cual se obtiene un sistema de ecuaciones lineales: Au b , donde u es una matriz columna que agrupa los
valores nodales u1 , u 2, u 3, A y b se obtienen ensamblando las contribuciones de los elementos. Obtenga
expresiones para las contribuciones de un elemento típico a las matrices A y b .

Este trabajo será entregado el día sábado 13 de julio de 2013.