UNIVERSIDAD POPULAR DEL CESAR

FACULTAD DE CIENCIAS BÁSICAS

DEPARTAMENTO DE MATEMÁTICAS Y ESTADÍSTICA
RED B

Programas: INGENIERIA DE ELECTRÓNICA E INGENIERÍA DE SISTEMAS

Asignatura: ANÁLISIS NUMÉRICO
Docente: JHONNY ANTONIO RIVERA V.
Semestre:
Grupos:
Fecha: ABRIL 2024
Ciudad: VALLEDUPAR

Tema: SOLUCIÓN NUMÉRICA DE ECUACIONES EN UNA VARIABLE Y SOLUCIÓN NUMÉRICA DE

SISTEMAS DE ECUACIONES LINEALES Y NO LINEALES.

1. Para los siguientes ejercicios:

1.1 Presente un estudio previo gráfico para analizar el comportamiento de la función
que define la ecuación dada y determin e un intervalo o un valor inicial cercano en
donde se en cuentre la solución deseada , según el método a emplear.
1.2 Seleccione un método apropiado y aplique los pasos manualmente del método
para determinar una aproximación de la solución requerida, con una tolerancia de 10−2 .

a. 𝑥 − 2−𝑥 = 0 ; única raíz.

b. 𝑒 𝑥 − 𝑥 2 + 3𝑥 − 2 = 0 ; raíz positiva más pequeña.
c. 2𝑥𝑐𝑜𝑠(2𝑥) − (𝑥 + 1)2 = 0 ; raíz negativa más cercana a cero.
d. (𝑥 − 2)2 − 𝑙𝑛𝑥 = 0 ; raíz más pequeña.
e. 𝑠𝑒𝑛𝑥 − 𝑒 −𝑥 = 0 ; raíz más pequeña.
f. 𝑒 𝑥 − 3𝑥 2 = 0 ; raíz más pequeña.
g. 𝑒 −0.5𝑥 = 5 − 5𝑥; raíz negativa más cercana a cero.
h. 𝑒 −0.5𝑥 (4 − 𝑥) − 2 = 0; raíz más pequeña.

2. Para cada uno de los siguientes problemas.

2.1 Haga un estudio previo y determine el método adecuado que le permita
encontrar la solución. (Gráfica en Matlab, análisis de la(s) condición(es) de
convergen cia, elección del intervalo o aproximación inicial según sea el caso, etc.).
2.2 Presente planteamien tos iniciales como: parámetros, adecuación de las ecuaciones,
sustituciones, cálculos, etc.
2.3 Ejecute el algoritmo del método que seleccionó implementado en Matlab , y
resuelva el problema con una tolerancia de 10−6.
2.4 Presente visualización completa del algoritmo implementado en Matlab y de los
resultados obtenidos (tabla). Haga su respectivo análisis y conclusión.

a. Como se ilustra en la figura 1. La velocidad del agua, v (m/s), en la descarga de un tanque

cilíndrico a través de un tubo largo se puede calcular como:

√2𝑔𝐻
𝑣 = √2𝑔𝐻𝑡𝑎𝑛ℎ⁡ ( 𝑡)
2𝐿

donde g=9.81m/s2, H=carga hidrostática inicial (m), L=longitud del tubo (m) y t=tiempo
transcurrido (s). Determine la carga hidrostática necesaria para obtener v=5m/s en 2.5s para un
tubo de 4 m de longitud. Utilice los valores iniciales de 𝑎 = 0⁡𝑦⁡⁡𝑏 = 2⁡𝑚, con una tolerancia de 10−2 .

1
 figura 1.

b. Una partícula inicia en reposo en un plano ligeramente inclinado (ver figura 2).,
𝑑𝜃
cuyo ángulo 𝜃 cambia a una velocidad con stante = 𝜔 < 0.
𝑑𝑡
Al final de t segundos, la posición del objeto está determinada por:

𝑔 𝑒 𝜔𝑡 − 𝑒 −𝜔𝑡
𝑥(𝑡) = − ( − 𝑠𝑒𝑛𝜔𝑡)
2𝜔 2 2

Suponga que la partícula se ha movido 1.7⁡𝑓𝑡 en 1 segundo. Encuentre, dentro de 10−5, la velocidad
𝜔 a la que 𝜃 cambia. Suponga que 𝑔 = 32.17⁡𝑓𝑡/𝑠 2 .

figura 2.

c. Un objeto que cae verticalmente a través del aire está sujeto a una resistencia viscosa, así
como a la fuerza de gravedad. Suponga que un objeto con masa m cae desde una altura 𝑠0 y que
la altura del objeto después de t segundos es:

𝑚𝑔 𝑚2 𝑔
𝑠(𝑡) = 𝑠0 − 𝑡 + 2 (1 − 𝑒 −𝑘𝑡/𝑚 )
𝑘 𝑘
donde ⁡𝑔 = 32.17𝑝𝑖𝑒𝑠/𝑠 2 y k representa el coeficiente de la resistencia del aire en 𝑙𝑏 −
𝑠/𝑝𝑖𝑒. Suponga 𝑠0 = 300⁡𝑝𝑖𝑒𝑠, 𝑚 = 0.25⁡𝑙𝑏⁡𝑦⁡𝑘 = 0.1𝑙𝑏 − 𝑠/𝑝𝑖𝑒. Encuentre, dentro de 0.01
segundos, el tiempo que tarda un cuarto de libra en golpear el piso.

d. Encuentre una aproximación para 𝜆, precisa dentro de 10−4, para la ecuación de crecimiento
poblacional
435⁡000
1⁡564⁡000 = 1⁡000⁡000⁡𝑒 𝜆 + ⁡(𝑒 𝜆 − 1),
𝜆

Use este valor para predecir la población al final del segundo año, al suponer que la tasa de
inmigración durante este año sigue siendo 435 000 individuos por año.

e. Los problemas relacionados con la cantidad de dinero requerida para pagar una hipoteca por un
período fijo incluye la formula

2
 𝑃
𝐴 = 𝑖 [1 − (1 + 𝑖)−𝑛 ],

Conocida como ecuación de anualidad ordinaria. En esta ecuación, A es la cantidad de la hipoteca,

P es la cantidad de cada pago e i es la tasa de interés por período para los n períodos de pago.
Suponga que se necesita una hipoteca por $135⁡000⁡𝑎⁡30 años para una vivienda y que el prestatario
puede efectuar pagos de máximo $1⁡000 por mes. ¿Cuál es la tasa de interés que el prestatario
puede pagar?

f. El volumen V de un líquido contenido en u tanque horizontal cilíndrico de radio r y longitud L está

relacionado con la profundidad del líquido h por

𝑟−ℎ
𝑉 = [𝑟 2 𝑐𝑜𝑠 −1 ( ) − (𝑟 − ℎ)√2𝑟ℎ − ℎ2 ] 𝐿
𝑟
Determine h dados r=2m, L=5m, y V=8m 3. Observe, que si usted está usando un lenguaje de
programación o una herramienta de software que no tenga abundancia de funciones
trigonométricas, el arco coseno se puede calcular con:

𝜋 𝑥
𝑐𝑜𝑠 −1 𝑥 = − 𝑡𝑎𝑛−1 ( )
2 √1 − 𝑥 2
g. Una carga total Q se encuentra distribuida en forma uniforme alrededor de un conductor en forma
de anillo con radio a. Una carga q se localiza a una distancia x del centro del anillo (véase la figura
3). La fuerza que el anillo ejerce sobre la carga está dada por la ecuación:

1 𝑞𝑄𝑥
𝐹=
4𝜋𝑒0 (𝑥 + 𝑎2 )3/2
2

donde 𝑒0 = 8.85⁡ × 10–12 ⁡𝐶 2 /(𝑁⁡𝑚2 ). Encuentre la distancia 𝑥 donde la fuerza es de 1⁡𝑁, si q y Q

son 2⁡ × ⁡ 10−5 𝐶 para un anillo con un radio de 0.9 m.

figura 3.

h. La resistividad ρ de un lubricante de sílice se basa en la carga q en un electrón, la densidad del

electrón n, y la movilidad del electrón 𝜇. La densidad del electrón está dada en términos de la
densidad del lubricante N, y la densidad intrínseca de acarreo ni. La movilidad del electrón está
descrita por la temperatura T, la temperatura de referencia T0, y la movilidad de referencia μ0. Las
ecuaciones que se requieren para calcular la resistividad son las siguientes:

1 1 𝑇 −2.42
𝜌 = 𝑞𝑛𝜇 , donde 𝑛 = 2 (𝑁 + √𝑁 2 + 4𝑛𝑖2 ) y 𝜇 = 𝜇0 (𝑇 )
0

3
Determine N, dado que T0 = 300 K, T = 1 000 K, μ0 = 1 300 cm2(V s)–1, q = 1.6 × 10–19 C,
ni = 6.21 × 109 cm–3, y un valor deseable de 𝜌 = 6 × 106 V s cm/C.

i. En la figura 4 se muestra un circuito con una resistencia, un inductor y un capacitor en paralelo.

Para expresar la impedancia del sistema se emplean las leyes de Kirchhoff, así:

1 1 1 2
= √ + (𝜔𝐶 − )
𝑧 𝑅2 𝜔𝐿

Donde z= impedancia (𝛺) y 𝜔=frecuencia angular. Encuentre la 𝜔 que da como resultado una
impedancia de 75⁡𝛺 , utilice un método cerrado con valores iniciales de 1 y 1000 y los parámetros
siguientes: 𝑅 = 225⁡𝛺, 𝐶 = 0.6 × 10−6 ⁡𝐹 y 𝐿 = 0.5𝐻.
Determine cuántas iteraciones son necesarias con la técnica a fin de encontrar la respuesta con una
𝑇𝑜𝑙 = 10−1. Utilice el enfoque gráfico para explicar cualesquiera dificultades que surjan.

figura 4.

j. La ecuación 𝑖 = 9𝑒 −𝑡 𝑠𝑒𝑛(2𝜋𝑡), describe una corriente oscilatoria en un circuito eléctrico, donde t

se expresa en segundos. Determine todos los valores de t de modo que i = 3.5

k. Considere el circuito eléctrico

La ecuación diferencial
𝑑2 𝑞(𝑡) 𝑑𝑞(𝑡) 1
𝐿 + + 𝑞=0
𝑑𝑡 2 𝑑𝑡 𝐶
Es el modelo matemático que describe el comportamiento de la variación de la carga q(t)
almacenada en el capacitor descrito, cuando se cierra el circuito y se genera un proceso de ajuste
hasta alcanzar un estado estacionario.
La solución de esta ecuación está dada por:

1 𝑅 2
𝑞(𝑡) = 𝑞0 𝑒 −𝑅𝑡/(2𝐿)
𝑐𝑜𝑠 ( √ − ( ) 𝑡)
𝐿𝐶 2𝐿

4
Si para t=0 segundos, 𝑞 = 𝑞0 = 𝑉0 𝐶 coulombios, 𝑉0 =voltaje de la batería en voltios. R=resistencia
dada en ohmios, L=inductancia, dada en henrios. =capacitancia dada en faradios.

k1. Se quiere determinar la resistencia R apropiada para disipar la carga a 1% de su valor original
𝑞
(𝑞 = 0.01) ⁡𝑒𝑛⁡𝑡 = 0.05𝑠⁡, 𝑐𝑜𝑛⁡𝐿 = 5𝐻⁡𝑦⁡𝐶 = 10−4 𝐹
0
k2. Determine el valor de L (inductancia) que se requiere en el circuito de manera que disipe el
1% del valor de la carga original en 𝑡 = 0.05𝑠, 𝑐𝑜𝑛⁡𝑅 = 280𝛺⁡𝑦⁡𝐶 = 10−4 𝐹 .

l. De acuerdo con el principio de Arquímedes la fuerza de flotación es igual al peso del fluido
desplazado por la porción sumergida de un objeto. Para el frústrum de un cono ilustrado en la figura
5, determine la altura h de la porción que queda encima del agua. Use los siguientes valores para
su cálculo 𝑟1 = 0.5⁡𝑚, 𝑟2 = 1⁡𝑚, ℎ = 1⁡𝑚, 𝜌𝑓 = 𝑑𝑒𝑛𝑠𝑖𝑑𝑎𝑑⁡𝑑𝑒𝑙⁡𝑓𝑟ú𝑠𝑡𝑟𝑢𝑚⁡ = 200⁡𝑘𝑔/𝑚3 ⁡𝑦⁡𝜌ℎ⁡ =
𝑑𝑒𝑛𝑠𝑖𝑑𝑎𝑑⁡𝑑𝑒𝑙⁡𝑎𝑔𝑢𝑎⁡ = 1000⁡𝑘𝑔/𝑚3 . Observe que el volumen del frústrum está dado por:

𝜋ℎ 2
𝑣= (𝑟 + 𝑟22 + 𝑟1 𝑟2 )
3 1

figura 5.

m. Suponga que se está diseñando un tanque esférico (ver figura 6)., para almacenar agua para
un poblado pequeño en un país en desarrollo. El volumen de líquido que puede contener se calcula
con
[3𝑅 − ℎ]
𝑉 = 𝜋ℎ2
3
donde V=volumen (𝑚3 ) , h=profundidad del agua en el tanque (m) y R=radio del tanque (m). Si R
= 3 m, ¿a qué profundidad debe llenarse el tanque de modo que contenga 30 𝑚3 ? Según el metodo
cerrado a emplear, utilice valores iniciales de 0 y R.

figura 6.

5
n. La ecuación 𝑖 = 9𝑒 −𝑡 𝑐𝑜𝑠(2𝜋𝑡), describe una corriente oscilatoria en un circuito eléctrico, donde t
se expresa en segundos. Determine todos los valores de t de modo que i = 3.

3. En cada uno de los siguientes ejercicios.

3.1 Realice un planteamiento inicial del sistema. (condicionamiento, coeficientes ceros,
etc.).
3.2 Haga un análisis de las condiciones de convergencia (comprobación de solución única, sistemas
equivalentes, etc.) y determine un vector inicial.
3.3 Realice paso a paso tres iteraciones manualmente, aplicando la formula del método de Jacobi
y la del método de Gauss-seidel y determin e una aproximación de la solución con una
tolerancia de 10−2 .
3.4 Ejecute el algoritmo del método de Jacobi y del método de Gauss-seidel implementado en
Matlab, y aproxime a la solución del sistema con una tolerancia de 10−6 en la norma 𝑙∞ .
3.5 Presente visualización completa del algoritmo implementado en Matlab y de los
resultados obtenidos (tabla). Haga su respectivo análisis y conclusión.

a. El sistema lineal
2𝑥1 − 𝑥2 + 𝑥3 = −1
2𝑥1 + 2𝑥2 + 2𝑥3 = 4
−𝑥1 − 𝑥2 + 2𝑥3 = −5

b. El sistema lineal
𝑥1 + 2𝑥2 − 2𝑥3 = 7
𝑥1 + 𝑥2 + 𝑥3 = 2
2𝑥1 + 2𝑥2 + 𝑥3 = 5

c. El sistema lineal
⁡⁡⁡⁡𝑥1⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡ ⁡⁡⁡ − 𝑥3 = 0.2
1 1
⁡⁡⁡⁡− 𝑥1 + ⁡ ⁡⁡𝑥2 − 𝑥3 = −1.425
2 4
1
⁡⁡⁡𝑥1 − 𝑥2 + ⁡⁡⁡ 𝑥3 = 2⁡⁡⁡⁡
2
d. Dado el sistema lineal 𝐴𝑥 = 𝑏

6
e. Dado el sistema lineal 𝐴𝑥 = 𝑏

f. Dado el sistema lineal 𝐴𝑥 = 𝑏

g. Dado el sistema lineal 𝐴𝑥 = 𝑏

h. Un cable coaxial está formado por un conductor interno de 0.1 pulgadas cuadradas y un
conductor externo de 0.5 pulgadas cuadradas. El potencial en un punto en la sección transversal
del cable se describe mediante la ecuación de Laplace.
Suponga que el conductor interno se mantiene en 0 volts y el conductor externo se mantiene en
110 volts. Aproximar el potencial entre los dos conductores requiere resolver el siguiente sistema
lineal.

7
i. Utilice las reglas de Kirchhoff para determinar el sistema de ecuaciones de la figura 7. Luego
calcule las corrientes del circuito aplicando lo indicado en el punto 3.

figura 7.

j. Utilice las reglas de Kirchhoff para determinar el sistema de ecuaciones de la figura 8. Luego
calcule las corrientes del circuito aplicando lo indicado en el punto 3.

figura 8.

k. Utilice las reglas de Kirchhoff para determinar el sistema de ecuaciones de la figura 9. Luego
calcule las corrientes del circuito aplicando lo indicado en el punto 3.

8
 figura 9.

l. Utilice las reglas de Kirchhoff para determinar el sistema de ecuaciones de la figura 10. Luego
calcule las corrientes y los voltajes del circuito aplicando lo indicado en el punto 3.

figura 10.

m. Utilice las reglas de Kirchhoff para determinar el sistema de ecuaciones de la figura 11. Luego
calcule las corrientes y los voltajes del circuito aplicando lo indicado en el punto 3.

figura 11.

n. Utilice las reglas de Kirchhoff para determinar el sistema de ecuaciones de la figura 12. Luego
calcule las corrientes para cada conductor del circuito aplicando lo indicado en el punto 3.

9
 figura 12.

o. Utilice las reglas de Kirchhoff para determinar el sistema de ecuaciones de la figura 13. Luego
calcule las corrientes para el circuito aplicando lo indicado en el punto 3.

figura 13.

4. En cada ejercicio.
4.1 Realice un planteamiento inicial (adecuación de las ecuaciones, variables o parámetros).
4.2 Haga un estudio previo gráfico (en matlab) y análisis de las condiciones de convergencia,
y determine un vector inicial apropiado que permita generar las aproximaciones a una de las
soluciones (intersecciones) del sistema.
4.3 Realice paso a paso dos iteraciones manualmente, aplicando la estrategia del método de
Newton para sistemas, y determine una aproximación de la solución con una tolerancia de
10−2 .
4.4 Ejecute el algoritmo del método de Newton para sistemas implementado en Matlab, y aproxime
a la solución del sistema con una tolerancia de 10−6 en la norma 𝑙∞ .
3.5 Presente visualización completa del algoritmo implementado en Matlab y de los
resultados obtenidos (tabla). Haga su respectivo análisis y conclusión.

a. 𝑒 𝑥+𝑦 + 𝑥 − 𝑦 = 3 b. 𝑠𝑒𝑛(𝑥) + 𝑒 𝑦 − 𝑥𝑦 = 5
𝑠𝑒𝑛(𝑥 + 𝑦) − 2𝑥 + 𝑦 = 1 𝑥 2 + 𝑦 3 − 3𝑥𝑦 = 7
1
c. 𝑐𝑜𝑠(𝑥 + 𝑦) + 𝑥𝑦 = 3 d. 4𝑥12 − 20𝑥1 + 𝑥22 + 8 = 0
4
1
3(𝑥 − 2)2 − 2(𝑦 − 3)2 = 5𝑥𝑦 𝑥 𝑥2
2 1 2
+ 2𝑥1 − 5𝑥2 + 8 = 0

10
e. 𝑠𝑒𝑛(4𝜋𝑥1 𝑥2 ) − 2𝑥2 − 𝑥1 = 0 f. 5𝑥12 − 𝑥22 = 0
(4𝜋−1)
(𝑒 2𝑥1 − 𝑒) + 4𝑒𝑥22 − 2𝑒𝑥1 = 0 𝑥2 − 0.25(𝑠𝑒𝑛𝑥1 + 𝑐𝑜𝑠𝑥2 ) = 0
4𝜋

g.⁡ 3𝑥12 − 𝑥22 = 0 h.⁡ 𝑙𝑛(𝑥12 + 𝑥22 ) − 𝑠𝑒𝑛(𝑥1 𝑥2 ) = 𝑙𝑛2 + 𝑙𝑛𝜋

3𝑥1 𝑥22 − 𝑥13 − 1 = 0 𝑒 𝑥1 −𝑥2 + 𝑐𝑜𝑠(𝑥1 𝑥2 ) = 0

11