@@ -22,6 +22,8 @@ Inversion de matrices
22
22
23
23
#. Calculer `M^{-1} ` par pivot de Gauß.
24
24
25
+ #. Refaire le même calcul avec une matrice par blocs `M=\begin{pmatrix}A&B\\C&D\end{pmatrix} `!
26
+
25
27
.. TOPIC :: Solution
26
28
27
29
La matrice inverse::
@@ -35,6 +37,7 @@ Inversion de matrices
35
37
[ d/(-b*c + a*d) (-b)/(-b*c + a*d)]
36
38
[(-c)/(-b*c + a*d) a/(-b*c + a*d)]
37
39
40
+
38
41
Par pivot de Gauß::
39
42
40
43
sage: I2 = matrix(2,2,1); I2
@@ -43,9 +46,10 @@ Inversion de matrices
43
46
sage: M = M.augment(I2, subdivide=True); M
44
47
[a b|1 0]
45
48
[c d|0 1]
46
- sage: M[1] = a*M[1] - c *M[0]; M
47
- [ a b| 1 0]
48
- [ 0 -b*c + a*d| -c a]
49
+ sage: M[1] = M[1] - c/a *M[0]; M
50
+ [ a b| 1 0]
51
+ [ 0 (-b*c + a*d)/a| (-c)/a 1]
52
+
49
53
sage: M[1] = M[1]/M[1,1]; M
50
54
[ a b| 1 0]
51
55
[ 0 1|(-c)/(-b*c + a*d) a/(-b*c + a*d)]
@@ -56,7 +60,6 @@ Inversion de matrices
56
60
[ 1 0| d/(-b*c + a*d) (-b)/(-b*c + a*d)]
57
61
[ 0 1|(-c)/(-b*c + a*d) a/(-b*c + a*d)]
58
62
59
-
60
63
.. TOPIC :: Théorème: formule d'inversion de matrice par blocs
61
64
62
65
Soit `M=\begin{pmatrix}A&B\\C&D\end{pmatrix} ` une matrice par
@@ -382,7 +385,7 @@ de degré `2^r` en appliquant récursivement l'étape précédente.
382
385
383
386
.. TOPIC :: En pratique: implantation
384
387
385
- L'algorithme de Karatsuba, étant plus compliqué en particulier à
388
+ L'algorithme de Karatsuba étant plus compliqué, en particulier à
386
389
cause de la récursion, est moins performant en petit degré que
387
390
l'algorithme naïf. Aussi les implantations utilisent l'étape de
388
391
récurrence en haut degré, et basculent sur un produit naïf en deçà
@@ -484,7 +487,7 @@ Transformée de Fourier Discrète
484
487
#. Donner la matrice inverse.
485
488
486
489
487
- Indication: `\sum_{k=0}^{n-1} \omega^{ik} = \begin{cases}n&\text{si $ i\equiv 0[n]$} \0円&\text{sinon}\end{cases} `
490
+ Indication: `\sum_{k=0}^{n-1} \omega^{ik} = \begin{cases}n&\text{si } i\equiv 0[n]\0円&\text{sinon}\end{cases} `
488
491
489
492
.. TOPIC :: Proposition
490
493
0 commit comments