Scribd
Cargar
175 vistas2 páginas

Método de Broyden

El método de Broyden es un método iterativo para resolver sistemas de ecuaciones no lineales. Sustituye la matriz Jacobiana en cada iteración por otra matriz A(i) que aproxima el cambio en la función entre iteraciones. El método del descenso más rápido minimiza una función cuadrática moviéndose en la dirección negativa del gradiente. El método conjugado genera una sucesión ortogonal según el producto interno con A para resolver sistemas lineales en n pasos.

Cargado por

Andrea Rivera
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
175 vistas2 páginas

Método de Broyden

El método de Broyden es un método iterativo para resolver sistemas de ecuaciones no lineales. Sustituye la matriz Jacobiana en cada iteración por otra matriz A(i) que aproxima el cambio en la función entre iteraciones. El método del descenso más rápido minimiza una función cuadrática moviéndose en la dirección negativa del gradiente. El método conjugado genera una sucesión ortogonal según el producto interno con A para resolver sistemas lineales en n pasos.

Cargado por

Andrea Rivera
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd

Método de Broyden

La primera iteración del método se realiza como el método de Newton, pero a partir de la
segunda iteración, la matriz Jacobiana J f ( X (i ) ) es sustituida por otra matriz [ A( i) ] que se
[ ]
obtiene a partir de [ A( i−1 ) ].

Siendo [ A( 0) ] = J f ( X (0 ) ) y siguiendo el método iterativo:


[ ]
−1 (i )
X ( i+1 )= X ( i)-[ A( i) ] * f ( X ) ( i=0,1 , … )

Inspirada en el método de la secante donde el valor de f ´ ( X i ) se aproxima mediante la


expresión:

f ( X i )−f ( X i−1 )
f ´ ( X i ) ≈ f ´ i=
X i− X i−1

Se obtiene que:

f ´ i ( X i−X i−1 )= f ( X i )−f ( X i−1 )

Broyden propuso reemplazar en cada iteración del método de Newton la matriz tangente
[J f
( X (i ) ) ] por otra matriz [ A( i) ] que compruebe que:
[ A( i) ]*( X (i)− X (i−1) )= f ( X (i )) −f ( X (i−1) )
Siendo X ( i)y X ( i−1) dos vectores cercanos y x otro vector cercano también próximo a los
anteriores se tiene que:

f ( X ) ≈ f ( X (i ) ) + [ J f ( X (i )) ]∗( X−X (i) )

f ( X ) ≈ f ( X (i−1 )) + [ J f ( X (i−1) ) ]∗( X −X ( i−1 ) )

Restando las aproximaciones anteriores se obtiene que:

f ( X (i )) + [ J f ( X (i) ) ]∗( X −X ( i) )−( X (i−1 )) + [ J f ( X (i−1) ) ]∗( X −X (i−1 )) ≈ 0

En el método de Broyden las aproximaciones del vector solución se buscan de la forma:

X ( i+1 )=X (i )−[ A(i ) ]∗f ( X (i ) )


Método del descenso mas rápido

El método del descenso más rápido consiste en dar una sucesión de valores x k tal que
q ( x k +1 ) ≤ q ( x k ) y que la igualdad solo se produzca en algún caso especial. La idea genérica es
obtener x k+1 a partir de x k y de otro vector v k que puede ya estar dado o que se lo vaya
construyendo durante la evolución del método.

Lo que se busca es minimizar la cuadrática q ( x ) si nos movemos desde x en una dirección v, es


decir que se minimiza la función de una variable q ( x +tv )en la variable t; tenemos que el

minimo sobre este rayo se produce si t^ =


⟨ v ,b− Ax ⟩ y por lo tanto
⟨ v , Av ⟩
2
q ( x + t^ v )=q ( x )− ⟨ v , b− Ax ⟩ ⟨ v , Av ⟩
A partir de esto, se genera la sucesión

x k+1=x k +t k v k

Donde t k =⟨ v , b− Ax k ⟩ ⟨ v , Av k ⟩ y la sucesión de v k es una conveniente (esta elección es la que


da origen a los distintos métodos).

Descenso más rápido → El método del descenso más rápido se define tomando

v k =−∇ q ( x k )(O un múltiplo positivo de él). La practicidad del método reside en que el
rk
residuo r k =b−Ax k apunta en la dirección de ¿−∇ q ( x k ), por lo tanto tomamos v k =
‖r k‖
para generar la sucesión aproximante x k .

Conjugada La idea es considerar una base { u1 , … , un } ortogonal según el producto interno


⟨ u , v ⟩= ⟨ u , A v ⟩ entonces vale que la secuencia dada por
⟨ b− Axk −1 ,u k ⟩
x k =x k−1 + 1≤k≤n
⟨ uk , A u k ⟩
Arrancando desde un punto x 0 cualquiera. Una gran ventaja de este método es que, en

teoría, en n pasos alcanza la solución del sistema, pero en la práctica es muy sensible

a los errores.

También podría gustarte