Différences

Ci-dessous, les différences entre deux révisions de la page.

--- math:matrices:systeme_lineaire:directe [2019/03/17 22:54] – Ajout du code python pour LU root
+++ math:matrices:systeme_lineaire:directe [2020/04/28 23:00] (Version actuelle) – mhtml -> html root
@@ Ligne 174: / Ligne 174: @@
 ====Méthode de Gauss / Crout : LU (A non symétrique)====
-Valable pour une matrice carrée définie positive.
+A doit être inversible (pas de pivot nul).
 Étape 1 :
@@ Ligne 258: / Ligne 258: @@
 Méthode :
-[[https://en.wikipedia.org/wiki/Cholesky_decomposition#The_Cholesky_algorithm|The Cholesky algorithm]] {{ math:matrices:systeme_lineaire:directe:cholesky_decomposition_-_wikipedia.mhtml |Archive 17/10/2018}}
+[[https://en.wikipedia.org/wiki/Cholesky_decomposition#The_Cholesky_algorithm|The Cholesky algorithm]] {{ :math:matrices:systeme_lineaire:directe:cholesky_decomposition_-_wikipedia_2020-04-28_10_55_30_pm_.html |Archive du 13/04/2020 le 28/04/2020}}
 Soit
@@ Ligne 426: / Ligne 426: @@
 $L = \begin{pmatrix} \sqrt 4 = 2 & 0 & 0 \\ \frac{12}{2} = 6 & \sqrt{37-6^2} = 1 & 0 \\ \frac{-16}{2} = -8 & \frac{-43-(-8)*6}{1} = 5 & \sqrt{98-8^2-5^2} = 3 \end{pmatrix}$
+Implémentation en python :
+<code python>
+def cholesky(A):
+    n = len(A)
+    L = np.zeros((n,n), dtype=np.complex64)
+    L[0,0] = cmath.sqrt(A[0,0])
+    L[1:,0] = A[1:,0]/L[0,0]
+    for j in range(1,n):
+        if (not np.isreal(A[j, j])):
+# Matrice A non hermitienne
+            L.fill(np.nan)
+            return (L,-1)
+        carre = A[j,j] - np.dot(L[j,:j],L[j,:j].transpose().conjugate())
+        if (carre <= 0):
+# Matrice non définie strictement positive
+            L.fill(np.nan)
+            return (L,0)
+        L[j,j] = cmath.sqrt(carre)
+        for i in range(j+1,n):
+            if (A[i,j] != np.conjugate(A[j,i])):
+# Matrice A non hermitienne
+                L.fill(np.nan)
+                return (L,0)
+            L[i,j]=(A[i,j]-np.dot(L[i,:i],L[j,:i].transpose().conjugate()))/L[j,j]
+    return (L,1)
+A4 = np.array(
+        [[  4,  12, -16],
+         [ 12,  37, -43+2j],
+         [-16, -43-2j,  98]], dtype=np.complex64)
+(L, sol) = cholesky(A4)
+print ("L",L)
+print ("sol",sol)
+print (L.dot(L.transpose().conjugate()))
+</code>
 ====QR (A non symétrique) méthode Householder====
@@ Ligne 451: / Ligne 493: @@
 $$Q = \prod_{k=1}^{N-1}Q_i$$
-[[https://fr.wikipedia.org/wiki/D%C3%A9composition_QR|Décomposition QR]] {{ :math:matrices:systeme_lineaire:directe:decomposition_qr_wikipedia.mhtml |Archive du 06/02/2019}}
+[[https://fr.wikipedia.org/wiki/D%C3%A9composition_QR|Décomposition QR]] {{ :math:matrices:systeme_lineaire:directe:decomposition_qr_wikipedia_2020-04-28_10_55_29_pm_.html |Archive du 16/05/2019 le 28/04/2020}}
 [[http://metronu.ulb.ac.be/MATH-H-202/an_ch3.pdf|Factorisation QR et systèmes surdéterminés]] {{ :math:matrices:systeme_lineaire:directe:an_ch3.pdf |Archive du 06/02/2019}}