Cet article se veut de présenter un nouvel outil de calcul matriciel [1] développé grâce à Snap!. Nous y aborderons des notions d’algèbre linéaire essentielles pour manipuler des matrices.
Tout d’abord, nous présenterons quelques how to :
– Comment générer la matrice Identité de taille n ?
– Comment générer une matrice de nombres aléatoires ?
– Sélection d’une ligne ou d’une colonne dans une matrice A
– Comment extraire le coefficient a_i,j d’une matrice A ?
– Comment calculer la trace d’une matrice à l’aide de la puissance du bloc {map}
– Comment calculer le déterminant d’une matrice donnée à l’aide d’une formule récursive ?
– Comment obtenir la comatrice d’une matrice donnée ?
– Comment obtenir l’inverse d’une matrice ?
– Puissances d’une matrice
Nous avons essayé dans l’écriture des blocs de cet outil de toujours avoir une pensée fonctionnelle et non pas une pensée séquentielle. On pourra se référer à cet article : Tout est algorithme, tout est fonction, dans lequel on explique qu’amener les algorithmes aux données est l’essence de la programmation fonctionnelle… [2]
La mise au point de cet outil nous a permis ensuite de découvrir de manière empirique un résultat original d’algèbre linéaire. Nous avons d’abord émis une conjecture sur les matrices carrées d’ordre n d’entiers naturels de 1 à n².
– Comment générer une telle matrice carrée ?
– Calcul du déterminant de cette matrice carrée pour n appartenant à 2, 3, 4, 5, 6, 7
– Conjecture : Dès que n ≥ 3, le déterminant de la matrice carrée d’ordre n dont les coefficients sont les entiers de 1 à n² (dans l’ordre naturel) est égal à zéro.
– Preuve très courte de ce théorème.
Depuis la version 6.0.0 de Snap!, une librairie très puissante pour les manipulations de vecteurs lui a été ajoutée par l’un de ses développeurs, Brian Harvey. C’est la librarie APL adaptée par Brian Harvey du langage APL : initialement /A Programming Language/, officieusement Array Processing Langage. Ce langage très puissant a été créé dans les années soixante pour décrire commodément des opérations globales sur des tableaux (de booléens, de données numériques et de caractères). Et l’intégration de cette librairie à Snap! lui confère de nouvelles vastes possibilités pour traiter les vecteurs de toute sorte.
Les formules utilisées dans cet article découlent d’un magnifique livre d’algèbre linéaire Matrices et réduction de Cabane et Lebœuf édité chez Ellipses en 1990.
Déplier tous les blocs Replier tous les blocs
Découverte rapide de la librairie APL de Snap! 6 et ultérieur
Voici quelques exemples de fonctions qui nous seront utiles.
Quelques How To
Comment générer la matrice Identité de taille n ?
Le résultat :
Comment générer une matrice de nombres aléatoires ?
Ce qui donne ceci :
Sélection d’une ligne ou d’une colonne dans une matrice donnée
Prenons par exemple la matrice carré d’ordre 4, B = [ [6,8,9,4],[9,6,1,7],[8,8,2,3],[1,5,7,8] ].
La matrice transposée de B est : tB = [ [6,9,8,1],[8,6,8,5],[9,1,2,7],[4,7,3,8] ]
On extrait la ième ligne :
On extrait la jème colonne :
On obtient alors respectivement la 3e ligne et la 2e colonne de la matrice B :
Comment extraire le coefficient a_i,j d’une matrice donnée ?
Voici la fonction qui renvoie ce coefficient a_i,j :
Ainsi, a_3,1 de B est 8.
Comment calculer la trace d’une matrice à l’aide de la puissance du bloc {map}
Comme tout un chacun sait (au moins vous qui lisez cet article...), la trace d’une matrice est la somme des éléments diagonaux.
Nous avons donc deux manière de calculer la trace : 1. version « old school » 2. avec le bloc {map}
1. La version que j’appelle old school se fait classiquement avec une boucle for :
2. Avec le bloc {map}, c’est beaucoup plus élégant (et efficace ! ) :
Sur la matrice carré d’ordre 4, B = [ [6,8,9,4],[9,6,1,7],[8,8,2,3],[1,5,7,8] ],
que l’on peut importer en csv si on le souhaite
,
les coefficients diagonaux s’obtiennent en appliquant la fonction qui calcule les coefficients a_i,j de la matrice B à la matrice ρ(4) [3] grâce au bloc map.
|
|
Ainsi le bloc map permet d’appliquer la fonction a_i,j( )( )(B) aux éléments de la liste ρ(4)=[1,2,3,4] :
Il ne reste alors plus qu’à sommer les éléments diagonaux obtenus en appliquant l’opérateur + grâce au bloc combine.
Comment calculer le déterminant d’une matrice donnée à l’aide d’une formule récursive ?
On extrait d’abord la matrice Ɑ_i,j en enlevant les lignes i et j de la matrice de départ.
Étant donné une matrice carré d’ordre n, de terme général (a_i,j), son déterminant est donné par cette formule récursive :
[4]
faisant intervenir la notion de cofacteur.
D’où le bloc récursif qui calcule le déterminant d’une matrice carré.
Comment obtenir la comatrice d’une matrice donnée ?
La comatrice d’une matrice donnée est la matrice de ses cofacteurs.
Le cofacteur d’indice i, k d’une matrice A est :
[5]
D’où le bloc comatrice
La comatrice de la matrice B est donc :
Comment obtenir l’inverse d’une matrice ?
L’inverse d’une matrice est :
D’où le bloc qui calcule l’inverse d’une matrice donnée :
Comment obtenir la puissance d’une matrice ?
Nous proposons ce bloc récursif :
Exemple d’application :
|
|
Conjecture sur les matrices carrées d’ordre n d’entiers naturels de 1 à n²
Calcul du déterminant de cette matrice carrée pour n appartenant à 2, 3, 4, 5, 6, 7
Il suffit d’appliquer (mapper) la fonction det( ) à chacun des termes de la liste (2,3,4,5,6,7).
Il semblerait que l’on obtienne 0 dès que n ≥ 3.
Snap! 7 permet la création de nouvelles catégories de blocs : la catégorie Matrix tools
Depuis la version 7 de Snap!, sortie pour Noël 2021, l’utilisateur peut créer ses propres catégories de blocs personnalisés.
Cela permet de ranger ses blocs dans des catégories particulières, ayant des couleurs choisies par l’utilisateur.
Vous trouverez le projet complet de calcul matriciel que je propose à cette adresse.
Bien sûr, il est incomplet mais je vous incite vivement à poursuivre ce travail pour le rendre encore plus utile...
L’utilisation de certains blocs de ce projet nécessite d’activer les extensions Javascript de Snap! comme ceci :
Sinon vous aurez ce genre de retour :
Je vous souhaite une agréable ballade dans ce projet de calcul matriciel. Afin de vous aider à découvrir ce projet et afin de vous en faciliter son utilisation, j’ai créé un lutin par étape :
N’hésitez pas à l’utiliser et à le faire utiliser par vos élèves de terminale et post-bac !
Vous trouverez au bas de cet article en téléchargement le support de la présentation de cet article qui s’est déroulée lors d’un séminaire de l’IREM de la Réunion à l’Université de Saint-Denis le 6 avril 2022
Et ci-dessous le document XML qui contient les catégories Snap! de Calcul Matriciel à importer dans une instance de Snap! :
Livré sous licence CC BY-NC-SA 4.0
Commentaires