Universidad de Santiago de Chile Facultad de Ciencia Departamento de Matemática y Ciencias de la Computación EJERCICIOS RESUELTOS EJERCICIOS RESUELTOS EJERCICIOS RESUELTOS EJERCICIOS RESUELTOS DE DE DE DE ECUACI ECUACI ECUACI ECUACIONES NO LINEALES NES NO LINEALES NES NO LINEALES NES NO LINEALES Profesor: Jaime Álvarez Maldonado Ayudante: Rodrigo Torres Aguirre Ejercicios: Ejercicios: Ejercicios: Ejercicios: 1) Sea la ecuación ), donde g satisface | )|1 ,,. a) Probar que también se cumple | ) )|| | , ,. b) Demostrar que si se cumple la condición en a), la ecuación ), tiene a lo mas, una solución en el intervalo , . Sol: Se tiene que ) ), y que la función gx) satisface la condición de | )|1 ,,. a) Del problema de punto fijo ) se desprende que: ) ) Si se restan queda: | | | ) )| Ocuparemos el T.V.M. Teorema de Valor Medio), el cual está dado por: ) ) ) | || )| | ) )| Se sabe por el enunciado que | )| , por lo que podemos relacionar la resta entre 2 problemas de punto fijo con el T.V.M., entonces se tendría que: | ) )|| )|| || | Por lo tanto | ) )|| | , ,. Queda probado. b) Se desea demostrar que ) tiene 1 sola solución en , . Suponiendo que existen en gx) 2 puntos fijos o soluciones), entonces: ) ) Al ser restados queda: | ) )|| | , ahora aplicamos TVM | || )| | ) )|). | || )|| | , y por el enunciado principal | )| ), se tiene que: | || )|| | | | | || | Siendo que L está entre 0,1, la ecuación anterior resulta ser una contradicción Entonces, se concluye que para gx) existe un único punto fijo en el intervalo , , siempre y cuando se cumpla la condición de que | )|1 ,,.
18
Embed
EJERCICIOS RESUELTOS Profesor: Jaime Álvarez …rorro.ublog.cl/archivos/1445/ejercicios_resueltos_de_ecuaciones_no... · EJERCICIOS RESUELTOS EJERCICIOS RESUELTOS DE DE DE ECUACIECUACIECUACIOOOONES
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Universidad de Santiago de Chile
Facultad de Ciencia
Departamento de Matemática y Ciencias de la Computación EJERCICIOS RESUELTOSEJERCICIOS RESUELTOSEJERCICIOS RESUELTOSEJERCICIOS RESUELTOS DE DE DE DE ECUACIECUACIECUACIECUACIOOOONES NO LINEALESNES NO LINEALESNES NO LINEALESNES NO LINEALES Profesor: Jaime Álvarez Maldonado Ayudante: Rodrigo Torres Aguirre Ejercicios:Ejercicios:Ejercicios:Ejercicios: 1) Sea la ecuación � � !�), donde g satisface | (!�)| ) * ) 1 +�, � , -., /0. a) Probar que también se cumple | !�8) 9 !�:)| ; *|�8 9 �:| +�8, �: , -., /0. b) Demostrar que si se cumple la condición en !a), la ecuación � � !�), tiene a lo mas, una solución en el intervalo -�8, �:0. Sol: Se tiene que � � !�) !=>?/@AB. CA =DEF? GHI?), y que la función g!x) satisface la condición de | (!�)| ) * ) 1 +�, � , -., /0. a) Del problema de punto fijo � � !�) se desprende que: �8 � !�8) �: � !�:) Si se restan queda: |�8 9 �:| � | !�8) 9 !�:)| Ocuparemos el T.V.M. !Teorema de Valor Medio), el cual está dado por: (!�) � O!PQ)RO!PS)PQRPS � |�8 9 �:|| (!�)| � | !�8) 9 !�:)| Se sabe por el enunciado que | (!�)| ) *, por lo que podemos relacionar la resta entre 2 problemas de punto fijo con el T.V.M., entonces se tendría que: | !�8) 9 !�:)| � | (!�)||�8 9 �:| ) *|�8 9 �:| Por lo tanto | !�8) 9 !�:)| ) *|�8 9 �:| +�8, �: , -., /0. Queda probado. b) Se desea demostrar que � � !�) tiene 1 sola solución en -�8, �:0. Suponiendo que existen en g!x) 2 puntos fijos !o soluciones), entonces: !�8) � �8 !�:) � �: Al ser restados queda: | !�8) 9 !�:)| � |�8 9 �:| , ahora aplicamos TVM !|�8 9 �:|| (!�)| � | !�8) 9 !�:)|). �|�8 9 �:|| (!�)| � |�8 9 �:| , y por el enunciado principal !| (!�)| ) *), se tiene que: |�8 9 �:| � | (!�)||�8 9 �:| ) *|�8 9 �:| �|�8 9 �:| ) *|�8 9 �:| Siendo que L está entre -0,10, la ecuación anterior resulta ser una contradicción � Entonces, se concluye que para g!x) existe un único punto fijo en el intervalo -�8, �:0, siempre y cuando se cumpla la condición de que | (!�)| ) * ) 1 +�, � , -., /0.
2) La ecuación 2�Z [ 24�] [ 61�: 9 16� [ 1 � 0 tiene dos raíces cerca de 0.1 (0.1213203436; 0.1231056256), encuéntrelas mediante el método de Newton-Raphson.
Sol: El método de N-R, es el método iterativo que requiere de la función, su derivada y un punto de inicio, la formula está dada por: �de8 � �d 9 f(Pg)
8�d] [ 72�d: [ 122�d 9 16 Entonces las iteraciones son: �₁ � 0.1111328125 �₂ � 0.11664780053787 �₃ � 0.11936143263559 �₄ � 0.12064808476922 �₅ � 0.12117604663885 �₆ � 0.12131031004941 �₇ � 0.12132028783201 �₈ � 0.12132034355771 � Aproximación a la raíz buscada (0.1213203436) Ahora buscaremos la otra raíz, tomando como punto de inicio �m � 0.13 �₁ � 0.12616290927433 �₂ � 0.1242900624559 �₃ � 0.12344358909839 �₄ � 0.12315206375549 �₅ � 0.123106774549 �₆ � 0.12310562635671 �₇ � 0.1231056256179 � Aproximación a la raíz buscada (0.1231056256)
3) Demuestre que al usar el método de Newton-Raphson, para aproximar el reciproco de un numero S, S>0 se obtiene la formula iterativa �ye8 � �y!2 9 z { �y), | � 0,1 …
Calcular 8
8} usando el algoritmo.
Sol:
El algoritmo de N-R, esta dado por �ye8 � �y 9f!P~)
fh!P~)
El reciproco de un número S, es: � �8
� � z �
8
P � G!�) � z 9
8
P� 0
Gj!�) �1
�:
Ahora reemplazamos los datos en la ecuación de N-R
�ye8 � �y 9�R
Q
�~Q
�~S
� �ye8 � �y 9!�{P~R8)
P~{
P~S
8
� �y 9 !z { �y: 9 �y)
� 2�y 9 z { �y:
��y { !2 9 z { �y)
Entonces: �ye8��y { !2 9 z { �y)
Para calcular 1/17, por el algoritmo encontrado, se considera S�17, por lo tanto, la ecuación quedaría:
�ye8��y { !2 9 17 { �y)
Se sabe que 1/17 está entre 0 y 0.1.
Por lo que consideraremos como punto de inicio �m � 0.05
Por lo tanto el intervalo de convergencia es [9 �� ��� : , ∞]
5) Encuentre el punto positivo donde la función G!�) � ��� PPS alcanza su valor mínimo calculando los ceros de Gj!�) con el método de Newton-Raphson. Calcule dicho valor mínimo. Sol: Se necesita saber la función a la cual se le aplicará el método de N-R, su derivada y el punto de inicio del algoritmo !dentro de un intervalo, para disminuir la cantidad de iteraciones). Entonces: G!�) � ��� PPS G(!�) � 8PS !��� P)S 9 : ��� PP� � Función buscada para aplicar el algoritmo G((!�) � � ��� PP� 9 ZP�!��� P)S [ : ��� PPS!��� P)� �Derivada de la función a analizar Ahora sustituimos los datos en la formula de N-R, la cual en nuestro casi es así: �de8 � �d 9 G(!�)Gjj!�) �de8 � �d 9 � 1
�d: (cos �d): 9 2 tan �d�d] ��6 tan �d�dZ 9 4�d]!cos �d): [ 2 sin �d�d:!cos �d)]�
El intervalo !para obtener el punto de inicio del algoritmo) lo podemos obtener por medio de un barrido entre 0 y 2 !aplicando el TVI), de sugerencia � � [0.9; 1]. Al tener el intervalo, elegimos algún punto que esté en su interior, como �m � 0.9. Las iteraciones son: �₁ � 0.94775114536635 �₂ = 0.94774713352959 �₃ = 0.94774713351695 �₄ = 0.94774713351695 El error absoluto de �₃ es: ��(�₃) = |�₄ 9 �₃| = |0.94774713351695 9 0.94774713351695| � 0 �� (�₃) � 0 Después de 4 iteraciones se llego a �₄ = 0.94774713351695 que es una aproximación a la raíz de G((�), con un �� (�₃) � 0. Esta aproximación es el punto mínimo de la función G(�) = ��� PPS . G((�₄ = 0.94774713351695) = -1.7250028381586*10R8] � 0 G(�₄ = 0.94774713351695) = 1.5494400344836 El valor mínimo de la función es 1.5494400344836.
6) Determinar algoritmos de punto fijo para obtener una solución aproximada de la ecuación !1 [ �)�AE!�) � 1.
Sol: Ordenamos la ecuación G!�) � �AE!�) [ ��AE!�) 9 1 � 0 Ahora le sumamos �: en ambos lados de la ecuación, quedando: �AE!�) [ ��AE!�) 9 1 [ �: � �: Despejando una de las x, para formar el algoritmo de punto fijo !�de8 � !�d)) �de8 � �AE!�d) [ �d�AE!�d) 9 1 [ �d: � Algoritmo de punto fijo Evaluamos la función en algunos puntos, para obtener un intervalo que contenga la raíz que buscamos. Por lo que nos queda un intervalo !después del “barrido”), � � -2.85; 2.90. Elegimos un punto de inicio, �m � 2.875 Las iteraciones del algoritmo son: �₁ � 2.8786243631782 �₁₁ � 2.8809861083846 �₂ � 2.8800553193257 �₁₂ � 2.8809862372817 �₃ � 2.8806194966361 �₁₃ � 2.8809862880556 �₄ � 2.8808418107431 �₁₄ � 2.8809863080559 �₅ � 2.8809293947499 �₁₅ � 2.8809863159342 �₆ � 2.8809638968662 �₁₆ � 2.8809863190376 �₇ � 2.880977487889 �₁₇ � 2.88098632026 �₈ � 2.8809828415741 �₁₈ � 2.8809863207415 �₉ � 2.8809849504512 �₁₉ � 2.8809863209312 �₁₀ � 2.88098578116 �₂₀ � 2.8809863210059 El error absoluto de �₁₉ es: ��! �₁₉) � |�₂₀ 9 �₁₉| � |2.8809863210059 9 2.8809863209312|�7.47*10R88 �� ! �₁₉) � 7.47*10R88< ℰ�10R8m Después de 20 iteraciones se llego a �₂₀ � 2.8809863210059 que es una aproximación a la raíz de G!�) � �AE!�) [ ��AE!�) 9 1, con un �� ! �₁₉) < ℰ�10R8m f !�₂₀ � 2.8809863210059)�1.697*10R8m � 0 � f !�₂₀)� 0
7) Obtener un algoritmo de punto fijo para calcular log8/] 5 . Sol: La ecuación es � � logQ
�5, que es igual a !8
])P � 5,entonces G!�) � !8])P 9 5.
Trabajamos un poco la ecuación, para pasarla a un problema de punto fijo.
!13)P � 5 /{ �
1 = 3 { √5� /[x � � 3 √5� [ � 9 1 �Problema de punto fijo �de8 � 3 √5�g [ �d 9 1 � !�d) � Algoritmo de punto fijo Para aplicar el algoritmo, necesitamos un punto de inicio, este punto se obtiene de un intervalo, el cual responde al TVI. Este intervalo es � � [91.5; 91.4] !entre más pequeño mejor), y el punto de inicio será �m � 91.45. Las iteraciones quedan: �₁ � -1.4612807817018 �₆ = -1.464969925092 �₂ = -1.4640531867328 �₇ = -1.4649726215252 �₃ = -1.4647435586268 �₈ = -1.4649732958476 �₄ = -1.4649160238965 �₉ = -1.4649734644823 �₅ = -1.4649591426499 �₁₀ = -1.4649735066545 El error absoluto de �₉ es: ��(�₉) = |�₁₀ 9 �₉| = | 9 1.4649735066545 9 91.4649734644823|=4.21722*10R� �� (�₉) = 4.21722*10R�< ℰ=10R} Después de 10 iteraciones se llego a �₁₀ = -1.4649735066545 que es una aproximación a la raíz de la función, con un �� (�₉) < ℰ=10R}
f (�₁₀ = -1.4649735066545)=-7.72512*10R� � 0 � f (�₁₀)� 0 (La raíz exacta es 9 �� §�� ])
8) La ecuación 2�] [ 4�: 9 � 9 5 = 0 tiene una raíz cercana a x�1. Obtener tres algoritmos iterativos de la forma � � !�), siendo !�) un radical, tales que converjan a la raíz, comenzando con �m � 1.Encontrar la solución indicando cual es el algoritmo que más rápidamente converge. Sol: La función es G!�) � 2�] [ 4�: 9 � 9 5�0 Tres problemas de punto fijo podrían ser: a) despejando �: � � � √:Pe§R:P�
: � (�)
b) despejando �] � � � §e:PRZPS:
� � (�) c) Aplicamos N-R (es un tipo de problema de punto fijo) � � � 9 (2�] [ 4�: 9 2 9 5)
(6�: [ 8� 9 2) = (�) Las iteraciones para (a), son: �₁ = 1.1180339887499
�₂ = 1.0536819972868
….
�₃₀ = 1.0781625824823
Las iteraciones para (b), son:
�₁ = 1.1447142425533
�₂ = 1.0079279312584
�₃ = 1.1385995170201
�₄ = 1.015033461394
�₅ = 1.1330072625887
…Converge muy lentamente
Las iteraciones para (c), son:
�₁ = 1.08333333333333
�₂ = 1.0781830462682
�₃ = 1.0781625876515
�₄ = 1.0781625873293
� Este es el algoritmo que converge más rápidamente a la raíz de f(x).
9) Resolver 656112√P�� 6P usando un algoritmo de punto fijo.
Sol:
Para poder usar un algoritmo de punto fijo, necesitamos pasar la ecuación a algo más trabajable, entonces:
656112√P�
� 6P /{ ln !)
ln 656112√PZ � ln 6P
√�4
{ ln 656112 � � { ln 6
G!�) � √PZ
{ ln 656112 9 � { ln 6 � 0
Ahora hacemos un barrido para ver en que intervalo esta la raíz de la función f!x).
G!0) � 0 �x�0 es raíz de la función
G!3) � 0.424531300494
G!3.4) � 0.08238673361
G!3.5) � 90.006647672497
G!4) � 90.469994484382
El intervalo donde se encuentra la raíz de f!x) es � � -3.4; 3.50, y el punto de inicio será �m � 3.45.
El problema de punto fijo con que se trabajara será: �de8 � !�d) � Pg
Después de 12 iteraciones se llego a �₁₂ � 3.4925731992805 que es una aproximación a la raíz de de la función, con un �� !�₁₁) < ℰ�10RZ
f !�₁₂ � 3.4925731992805)�9.3711503*10R� � 0 � f !��® ¯)� 0
10) Resolver . 9 √. [ � � � usando un algoritmo de punto fijo. Sol: Al ordenar la ecuación !igualándola a 0), queda expresado de la siguiente forma: G!�) � . 9 √. [ � 9 ��0 !la idea es obtener la raíz de f!x)) Para poder obtener esa raíz, tomaremos la ecuación inicial, como nuestro problema de punto fijo. El algoritmo de punto fijo queda: �de8 � . 9 . [ �d Y el punto de inicio será: �m � √. Las iteraciones quedan: �8 � . 9 . [ √. !Cantidad de a: 3) �: � ±. 9 ². [ . 9 . [ √. !Cantidad de a: 5)
�] � ³́́́́´́́́́µ. 9 ¶. [ ±. 9 ². [ . 9 . [ √. !Cantidad de a: 7)
…… �d � ³́́́́
´́́́́µ. 9 ³́́́́´́́́́µ. [ · … . ³́́́́
´́́́́µ. 9 ¶. [ ±. 9 ². [ . 9 . [ √. !Cantidad de a: -2n[10) El error absoluto de �d es: ��!�d) � |�d 9 �dR8| � 0 Después de n iteraciones se llego a que �d es la raíz de de la función G!�) � . 9 √. [ � 9 �, con un �� !�d) � 0 G !�d)� 0 Recomendación: compruebe con a�3 !la raíz es x�1), y con a�7 !la raíz es x�2).
11) Cuando a>1 y 0<b<1/2 la ecuación 8
¸�/¹R89
º
¸�R89 !. 9 1)/ � 0 tiene una
raíz positiva. Demostrarlo y calcularla, con cuatro cifras decimales exactas, para a�5 y b�1/4.
Sol:
La ecuación al multiplicarla por !AP/º 9 1)! AP 9 1) queda:
Después de 3 iteraciones se llego a �₃ � 2.8304817995381 que es una aproximación a la raíz de de la función, con un �� !�₂ ) < ℰ�10R�, y con 7 cifras significativas.
f !�₃ � 2.8304817995381)�-3*10R8: � 0 � f !�₃)� 0
12) Resolver la ecuación A√PSR8 � 2√� [ 1.24 , usando un algoritmo de punto fijo, partiendo con �m � 1.5 , operar con cinco decimales. Sol: Al ordenar la ecuación queda: √�: 9 1 � ln!2√� [ 1.24) /*!):
�: 9 1 � ln:!2√� [ 1.24) G!�) � �: 9 1 9 ln:!2√� [ 1.24) � 0 De esta función se desprende el problema de punto fijo. � � ¼1 [ ln:!2√� [ 1.24) El algoritmo de punto fijo queda: �de8 � ¼1 [ ln:!2 �d [ 1.24) Las iteraciones son: �m � 1.5 �₁ � 1.6444754042634 �₆ = 1.6737725207851
Después de 10 iteraciones se llego a �₁₀ = 1.6737758077194 que es una aproximación a
la raíz de de la función, con un �� (�₉ ) < ℰ=10R}
f (�₁₀= 1.6737758077194)=-6.3463*10R½ � 0 � f (�₁₀)� 0
13) ¿En qué valor de � tiene la curva ¿ � ARP ln � un punto de inflexión? Sol: Para calcular el punto de inflexión de una curva se necesita la segunda derivada de esta, entonces: ¿ � ARP ln � ¿( � ARP� 9 ARP ln � ¿(( � ARP ln � 9 2ARP
� 9 ARP�:
La idea es obtener la raíz de ¿(( , ya que esta es el punto de inflexión de la función ¿ . Entonces la segunda derivada se iguala a cero y se aplica algún método iterativo. ¿(( = ARP ln � 9 :¸Á�
P 9 ¸Á�PS = 0 /*�:
¿(( � G(�) � �: ln � 9 2� 9 1 � 0 El intervalo donde se encuentra la raíz es � � -2.5; 2.60, el cual fue obtenido por medio del T.V.I. Usando un algoritmo de punto fijo: �de8 = ²2�d [ 1ln �d Con el punto de inicio �m � 2.55 Las iteraciones son: �₁ = 2.5527326631178
Después de 10 iteraciones se llego a �8m=3.1925823986334 que es una aproximación a la solución de H con un �� ; Ë = 10R}
16) Encontrar las Raíces de P(x) = 1+ 3x +5x²+6x³
Sol:
Para encontrar las raíces del polinomio debemos primero ver en que intervalo se encuentran. El intervalo se puede encontrar por medio de un “anillo”, dado por la siguiente fórmula:
Entonces el intervalo que contiene a la raíces del polinomio es
I = [-2; 1/7] U [1/7; 2] = I₁ U I₂
Ahora hacemos un “barrido” en los 2 sub intervalos.
-Sub intervalo
P (-2) = 33 P (-1.2) = 5.768 P (-0.6) = 0.296 P (-0.4) = 0.216 P (-1/7) 0.6559767
Según el análisis, hay una raíz real y está ubicada en
I= [-0.6;-0.4]. Ya que P (-0.6)*P (
Según la derivada P’(x) =3+10x+18x²
P’(x)>0, ∀ x ∊ I. La función es creciente en el
� Ahora aplicamos el método de Horner
Considerado el punto inicial de
P(x) = 1+ 3x +5x²+6x³ mediante el método de Horner.
Para encontrar las raíces del polinomio debemos primero ver en que intervalo se encuentran. El intervalo se puede encontrar por medio de un “anillo”, dado por la
≤|x|≤ (|an|+am)/|an|
� 1/7≤ |x| ≤ 2 � 0.142 ≤ |x| ≤ 2
Entonces el intervalo que contiene a la raíces del polinomio es
₁ U I₂
Ahora hacemos un “barrido” en los 2 sub intervalos.
-Sub intervalo
P (1/7) 1.548105 P (0.4) =3.384 P (0.6) =5.896 P (1.2) =22.168 P (2) =75
Grafico del polinomio P(x) = 1+ 3x+5x²+6x³. Se puede ver que la raíz del polinomio es pero la idea del ejercicio es calcular esa raíz por medio de un método iterativo.
Según el análisis, hay una raíz real y está ubicada en el intervalo
0.6)*P (-0.4) <0
derivada P’(x) =3+10x+18x²
. La función es creciente en el intervalo I. En el intervalo existe una raíz
Ahora aplicamos el método de Horner
Considerado el punto inicial de iteración como x₀=-0.45
mediante el método de Horner.
Para encontrar las raíces del polinomio debemos primero ver en que intervalo se encuentran. El intervalo se puede encontrar por medio de un “anillo”, dado por la
1.548105 P (0.4) =3.384
0.6) =5.896 P (1.2) =22.168
P(x) = 1+ 3x+5x²+6x³. Se puede ver que la raíz del polinomio es -0.5, pero la idea del ejercicio es calcular esa raíz por medio de un método iterativo.