FACULTAD DE CIENCIAS NATURALES

Y MATEMÁTICAS 
UNIVERSIDAD NACIONAL DEL CALLAO
CIENCIA Y TENCNOLOGÍA RUMBO AL TERCER MILENIO 

 ASIGNATURA
 ASIGNATURA
Métodos Numéricos de la Física
TEMA
Método de Cholesky
DOCENTE
Mg. Jorge Luis Godier Amburgo
 ALUMNO
Yemerson Encarnación Rosales
ESCUELA
Profesional de Física

CALLAO, LIMA 2020

 
1
 

ÍNDICE

INTRODUCCIÓN
INTRODUC CIÓN ..................................
.....................................................................
.......................................................................
.............................................
.........3
CAPÍTULO
CAPÍTULO 1 ....................................
.......................................................................
.......................................................................
...................................................
...............4
MÉTODO CHOL
CHOLESKY
ESKY .............................
..................................................................
.........................................................................
.......................................
...4
ANDRE-LOU
ANDRE-LOUIS
IS CHOL
CHOLESKY
ESKY.............................
..................................................................
...............................................................
..........................4
DESCOMPOSICIÓN DE CHOLESKY ............................................................................. 4
PROPIEDA
PROPIEDADES
DES ....................................
.......................................................................
......................................................................
.......................................
.... 4
CÁLCULO
CÁLCULO ...............................
....................................................................
........................................................................
..................................................
...............5
APLICAC
APLICACIÓN
IÓN REAL
REALES...........
ES..............................................
......................................................................
...................................................
................6
VENTAJAS
VENTAJAS ..............................
...................................................................
........................................................................
..................................................
...............6
DESVENTA
DESVENTAJAS
JAS....................................
.......................................................................
......................................................................
.......................................
.... 6
PASOS PARA LA FACTORIZACION DE CHOLESKY............................................. 6
CAPITULO
CAPITULO 2 ....................................
.......................................................................
.......................................................................
...................................................
...............7
APLICACIÓN EN EL SISTEMA DE CIRCUITOS ............................................................. 7
MÉTODO DE RESOLUCIÓN EN FORTRAN ..................................................................... 9
SOLUCION SEGÚN EL MÉTODO DE FORTRAN ...................................................... 11
CONCLUCI
CONCLUCIÓN.........
ÓN.............................................
......................................................................
......................................................................
........................................ 12
BIBLIOGRAF
BIBLIOGRAFÍA
ÍA .............................
..................................................................
........................................................................
................................................
............. 13

 
2
 

INTRODUCCIÓN

El método de factorización o descomposición de cholesky lleva por nombre gracias al

matemático
LU, haciendo Andre-louis cholesky. Francés
que la descomposición que logro
matricial tengarealizar una variante
una pequeña en el método
variación de LU
,siendo así que al momento de considerarse aplicab
aplicable
le este es dos veces más
más efectiva que
la anteriormente menci
mencionada.
onada.

El método de cholesky considerado con más similitud al lower-upper (LU), logra

desprender de este último una matriz U en una matriz transpuesta Lt  , por lo que este
desarrollo a realizar va mostrar el proceso que se realizara para tener la precisión de
lograr este procedimiento, que este gran matemá
matemático
tico francés logro desarrollar, es así que
estamos ante el análisis del trabajo histórico de un intelectual por lo que veremos la
 precisión
 precisión que nos brinda este desarroll
desarrolloo en cuanto a una solución
solución de una forma más
efectiva brindando respuesta a problemas de cada ámbito de estudio que cuente con las
características necesarias
necesarias para utilizar el
e l método.

Donde presentaremos la efectividad y diferencias notorias que puede tener con respecto
a los otros métodos, en especial a la descomposición más similar de lower-upper; como
también mostrar las carencias que se logra obtener en el desarrollo.
desarrollo.

 
3
 

CAPÍTULO 1

MÉTODO CHOLESKY 

ANDRE-LOUIS CHOLESKY
Matemático francés nacido en Montguyon el 15 de octubre en 1875, Francia. Estudió en
la École polytechnique y trabajó en geodesia y cartografía.
Sirvió en el ejército francés como oficial de ingeniería y murió en una batalla a pocos
meses del final de la primera guerra mundial, siendo su trabajo de la descomposición o
método de Choles
Choleskyky publicado póstumamente.

DESCOMPOSICIÓN DE CHOLESKY

La factorización o descomposición de Cholesky toma nombre del matemático André-

Louis Cholesky
Cholesky,, quien encontró que una matriz simétrica definida positiva puede
ser descompuesta como el producto de una matriz triangular inferior y la traspuesta de
esta matriz triangular.

Forzando un poco más la factorización LU de una matriz, en casos especiales

encontramos que pueden elegirse L = Ut, es lo que se conoce como la factorización de
Cholesky para una matriz cuadrada simétrica: A = LLt.
Sólo para matrices simétricas definidas positivas.

Cualquier matriz cuadrada A con pivotes no nulos puede ser escrita como el producto de
una matriz triangular inferior L y una matriz triangular superior U; esto recibe el
nombre de factorización LU sin embargo, si A es simétrica y definida positiva, se puede
escoger los factores tales que U es la transpuesta
t ranspuesta de L, y esto se llama descomposición
descomposición o
factorización
factorización de cholesky.

Tanto la descomposición LU como la descomposición de Cholesky son usadas para

resolver sistemas de ecuaciones lineales. Cuando es aplicable, la descomposición de
Cholesky es dos veces más eficiente que la descomposición LU.

Este método sirve para resolver ecuaciones de la forma     .

PROPIEDADES
a)  Matriz simétrica.
 b)  Si la matriz  A siempre es definida positiva At  también lo es.
c)  La inversa de una matriz definida positiva siempre existe.
d)  Todas las submatrices principales de una matriz definida
def inida positiva son
 positivas.
 
4
 

e)  Una matriz simétrica es definida positiva si y solo si todos los menores
 principales son
son positivos (teorem
(teoremaa de Sylvester).

CÁLCULO

Para el cálculo efectivo de la factorización

factorización de Cholesky
Cho lesky,, se parte de la igualdad:

              
                    

(     ) (    ) (     )

Con    .

Con este desarrollo, podemos identificar los elementos de la primera fila de  A respecto a
 LLt , obtenemos:

 

  
 
             

Identificamos los elementos de la segunda fila de  A respecto a  LLt , obtenemos:

           

Entonces tenemos como resultado:


 

{             

Generalizando la propiedad, suponemos las conocidas k-1 columnas L, entonce
t 
entonces,
s,
identificamos los k-ésima fila de A correspondientes a  LL , obtenemos:

               

En consecuencia:

 
5
 


     ∑   
 
     ∑          
La factorización se calcularía con la consecuente relación:


 
 
  {  
  Utilizada para los el
elementos
ementos de llaa diagonal principal.
Utilizada para los demás elementos de la matriz  L.

APLICACIÓN REALES
•   Circuitos eléctricos
•   Ley de Hooke
•   Conservación de la energía
•   Sistemas de potencia
•
  Etc. 

VENTAJAS
•   Proporciona una resolución práctica y muy sencilla de realizar para sistemas no
muy comunes.

DESVENTAJAS
•   Tiendes a ser un método selectivo solo se usa para sistemas poco comunes y con
 propiedades específicas.
específicas.

PASOS PARA LA FACTORIZACIÓN DE CHOLESKY

1.  Tomamos los coeficientes del sistema de ecuaciones    

  y lo formamos
como la forma de matriz  A y b.
2.  Posteriormente partimos de la igualdad .    
3.  Se encuentra la matriz triangular inferior L y la matriz triangular superior  Lt ,
 paso seguido
seguido se calcula llos
os coeficientes de la m

matriz
atriz triangular inferi
inferior
4.  Teniendo los coeficientes de L automáticamente se tiene los valores de .
5.  Por la igualdad , hallamos el valor de y.
or  L.

 
6. relación
Teniendoa la
la matriz  es posible
. hallar el valor de nuestras incógnitas  x   en
 y 
igualdad
  
 6
 
 

CAPITULO 2

APLICACIÓN EN EL SISTEMA DE CIRCUITOS

Del siguiente circuito determinar los valores de las intensidades de la corriente I 1, I2, I3 
aplicando el método de cholesky en el sistema de ecuaciones.

Sea el sistema de ecuaciones:


   

        

 

Paso 1.-

Definido de la forma:      

  

Paso 2.-

Realizamos la igualdad:
[   ][ ]  []
 

[  ]  [   ] [  ]

          

Paso 3.-

Desarrollamos la matriz L:

•
 
 
  √ √   √ √  7
 
 

•        √ √ √ √ 

 
      
•

•
 
 
   
 
 

•    

Paso 4.-
        
Con la solución de L automá
automáticamente
ticamente tenem os Lt .
t enemos

Paso 5.-

Por la siguiente igualdad hallamos los valores de y:


        

•    
[ 
 
 
 
 
  
  
  ] [
       ] []  
•

•
 
           
 
 
 

Paso 6.-

Hallamos los valores de X con la siguiente relación:

     

 
[  ] [ ]  [] ]
•          
  2.314 
•

•
 
          
 
 
 

Solución:

Por lo tanto
t anto las soluciones para las intensidades son:
 
 

 
 
8
 
 

MÉTODO DE RESOLUCIÓN EN FORTRAN

PROGRAM METODO_DE_CHOLESKY
!DEFINICIÓN DE VARIABLES
REAL(4):: A(100,100),L(100,100)
A(100,100),L(100,100),LT(100,100)
,LT(100,100)
REAL(4):: B(100),Y(100),X(
B(100),Y(100),X(100)
100)
REAL(4):: SUM2,SUM3
INTEGER::T,I,J,K,N,COLUMNA
!INGRESO DE DATOS
WRITE(*,*)" METODO DE CHOLEKSKY PARA MATRICES SIMETRICAS POSITIVAS"
WRITE(*,*)" ================================"
WRITE(*,*)""
WRITE(*,*)" INGRESO DE DATOS"
WRITE(*,*)" ============================================="
WRITE(*,*)""
WRITE(*,*)" INGRESE LA CANTIDAD DE VARIABLES DEL SISTEMA"
READ(*,*)
READ(* ,*) T
WRITE(*,*)" INGRESE LOS COEFICIENTES DE LA MATRIZ A"
DO I=1,T
READ(*,*)(A(I,J),J=1,T)
END DO
WRITE
DO (*,*)' INGRESE MATRIZ DE TERMINOS
I=1,T T ERMINOS INDEPEN
INDEPENDIENTES
DIENTES b'
READ(*,*)B(I)
END DO

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!PROCESO PARA EL CALCULO CALCULO DE LA MATRIZ "L"
L(1,1)=SQRT(A(1,1))!
L(1,1)=SQRT(A(1,1) )! CALCULO DE DE L(1,1)

!PRIMERA COLUMNA
DO I = 1,T
L(I,1) = A(I,1)/L(1,1)
END DO

!DEMAS COLUMNAS
DO WHILE (COLUMNA <= T)
DO I=COLUMNA, T
E=0
DO J=1, COLUMNA-1
E=E+L(I,J)*L(COLUMNA,J)
IF(I == COLUMNA
COLUMNA)THEN
)THEN
L(I,I)=SQRT(A(I,I)-E)
END IF
L(I,COLUMNA) = (A(I,COLUMNA) -
E)/L(COLUMNA,COLUMNA)
END DO
END DO
COLUMNA = COLUMNA + 1
END DO

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

 
9
 

LT=TRANSPOSE(L) ! HALLAR LA MATRIZ "L" TRANSPUESTA

!PROCESO PARA HALLAR EL VETOR "Y"

Y(1)=B(1)/L(1,1)
 NK=2
26 SUM2=0

DO J=1,NK-1
SUM2=SUM2+L(NK,J)*Y(J)
END DO
Y(NK)=(B(NK)-SUM2)/L(NK,NK)
 NK=NK+1
IF (NK.GT.T) THEN
 NK=T
GOTO 27
ELSE
GOTO 26
END IF

!CALCULO DEL VECTOR X

27 X(T)=Y(T)/LT(T,T)
JK=NK-1
29
DOSUM3=0
I=JK+1,T
SUM3=SUM3+LT(JK,I)*X(I)
END DO
X(JK)=(Y(JK)-SUM3)/LT(JK,JK)
JK=JK-1
IF (JK.LE.0) THEN
GOTO 28
ELSE
GOTO 29
END IF

!SALIDA DE DATOS

28 WRITE(*,*)" MATRIZ L"

WRITE(*,*)" ========"
DO I=1,T
WRITE(*,*)(L(I,J),J=1,T)
END DO
WRITE(*,*)" MATRIZ LT"
WRITE(*,*)" ========"
DO K=1,T
WRITE(*,*)(LT(K,J),J=1,T)
END DO
WRITE(*,*)" MATRIZ Y"
WRITE(*,*)" ========"
DO N=1,T
WRITE(*,*)Y(N)
END DO
WRITE(*,*)" ======================================"
WRITE(*,*)" LOS VALORES DE X PARA EL SISTEMA SON:"

 
10
 

WRITE(*,*)" ======================================"
DO N=1,T
WRITE(*,*)N,X(N)
END DO

END PROGRAM

SOLUCIÓN SEGÚN EL MÉTODO DE FORTRAN

 
11
 

CONCLUCIÓN
Concluimos afirmando la cercana relación que existe entre el método de lower-upper
Concluimos
con el método de cholesky, pero aclarando las variaciones y las facilidades que nos

 brinda
mucho este
más método,
métod o, esy decir
sencillos la aaccesibilidad
prácticos, ccesibilidad ppara
pero objetandoara que
la solución
es muy por los procedimientos
pren
limitante ocedimientos
cuando a
distintos sistemas; ya que son  seleccionados  los sistemas que cumplen con este
“  ”

método para su solución.

 
12
 

BIBLIOGRAFÍA

Dpto. Ecuaciones Diferenciales y Análisis Numérico, U. d. (s.f.). Cálculo Numérico II. 

Recuperado el 10 de Agosto de 2020, de Cálculo Numérico II:
https://personal.us.es/pmr/images
https://personal.us.es/pmr/images/pdfs/1314gm-cnii-practica1bis.pdf
/pdfs/1314gm-cnii-practica1bis.pdf
LEON, F. (2020). StuDocu. Recuperado el 10 de AGOSTO de 2020, de StuDocu:
https://www.studocu.com/co/docum
https://www.studocu.c om/co/document/universidad-industrial-de
ent/universidad-industrial-de--
santander/dispositivos-electronicos/practica/elizguer-metodo-de-
cholesky/8855227/view?utm_c
cholesky/885 5227/view?utm_campaign=sha
ampaign=shared-
red-
document&utm_source=studocu-
document&utm_medium=social
documen t&utm_medium=social_sharing&utm_con
_sharing&utm_content=elizguer-met
tent=elizguer-met
Trevisan., P. M. (s.f.). CURSO Análisis Numérico 2. Recuperado el 10 de AGOSTO de
2020, de CURSO Análisis Numérico 2:
http://www.ciens.ula.ve/matematica/publicaci
http://www.ciens.ula.ve/matematica/publicaciones/contenidoPub
ones/contenidoPub.html
.html

 
13