\"Accueil\"





Sauvons le net Européen












Worst EU Lobbying Awards 2008
Votez dès le 15/10/2008








Plus de 200 000 signatures pour l'abandon d'Edvige recueillies depuis le 10 juillet 2008


Communiqué de presse du 29 octobre 2008
Le Conseil d'Etat rejette pour défaut d'urgence la demande de suspension du décret portant création du fichier « Edvige »..







Calcul matriciel

Description

Les calculs utiliseront généralement des matrices carrées de mêmes dimensions (recherche du polynôme caractéristique, du déterminant ...) et même parfois uniquement certaines de ces matrices (recherche de l'inverse). Mais on peut effectuer la somme et la différence de matrices rectangulaires de mêmes dimensions ou encore le produit d'une matrice (m,n) par une matrice (n,p).

L'application ci-dessous permet
  • De construire une matrice carrée particulière (nulle, unité, scalaire) ou encore au hasard (les coefficients sont entiers et choisis très proches de 0).
  • De construire une matrice carrée stochastique (matrice de transition) ou symétrique ou triangulaire supérieure ou orthogonale (Gram-Schmidt) ou encore circulante.
    De calculer par exemple le produit d'une matrice orthogonale et de sa transposée, vérifiez que l'on a obtenu (en arrondissant un peu les valeurs) une matrice unité.
  • De transposer une matrice quelconque, de permuter ses lignes ou ses colonnes,
  • D'obtenir le déterminant d'une matrice carrée et, s'il n'est pas nul de calculer l'inverse de la matrice,
  • De calculer une puissance entière de la matrice lorsque c'est possible (matrice carrée et inversible lorsque l'exposant est négatif),
  • De multiplier une matrice par un réel,
  • D'additionner, de soustraire, de multiplier deux matrices lorsque les dimensions des matrices sont compatibles,
  • De calculer le polynôme caractéristique d'une matrice carrée (utilise l'algorithme de Jean-Marie Souriau).
    En évaluant P(X) on retrouve les valeurs propres réelles 'lambda' lorsqu'on réussit à obtenir P(lambda)=0.
    (Essayer d'obtenir une racine par essais successifs et dichotomie. Essayez les trois intervalles délimités par -3, -2, 2, 10. À l'aide de GP/PARI vous pouvez obtenir les racines en calculant   polroots(X^3 - 6*X^2 - 21*X - 1), toutes les trois sont réelles).
  • Les valeurs propres (réelles ou complexes) peuvent être localisées à l'aide d'un théorème de Gerschgörin-Hadamard dans un disque de rayon r=min(L,C) où L (resp. C) sont les maximums des sommes des valeurs absolues des éléments de chacune des lignes (resp. des colonnes). Temporaire, en attendant une implémentation des complexes : En choisissant un pas suffisamment fin et en balayant l'intervalle réel [-r, r] on obtient par la méthode de dichotomie les valeurs propres réelles lorsqu'elles sont simples ou d'ordre de multiplicité impair. Les autres v.p. sont obtenues par détection des extremum, elles sont placées entre parenthèses comme (ceci). Toutes les valeurs obtenues sont des résultats approchés et incertains à vérifier soigneusement avant de conclure.
    Lorsque la matrice est symétrique, toutes les valeurs propres sont réelles.
    Essayer aussi une matrice triangulaire supérieure
  • De calculer la valeur de P(M) où M est la matrice affichée et P est donné par ses coefficients.
    Comme cas particulier on peut prendre pour P le polynôme caractéristique et vérifier que P(M) est la matrice nulle (th. de Cayley-Hamilton).
L'historique permet de retrouver aisément toutes les matrices utilisées précédemment et de les intégrer dans de nouveaux calculs.

Calculs

Transformations d'une matrice

Nouvelle matrice
Matrice carrée n =
Éléments réels ou entiers dans l'intervalle a b :
Matrice
Historique      

Opérations   

, ou les deux matrices no

ou par k = la matrice.

Permute les      
Autres      
     X= P(X)=
Racine de P(X) par en prenant x0 x1 = pas =
Majoration des modules des val. propres |lambda|
coeffs

Les résultats ci-dessus sont donnés sans garantie d'exactitude (en raison des très nombreux arrondis effectués et non maitrisés, des bugs de conception ou de programmation toujours possibles...). N'utilisez cette application qu'à vos risques et périls, en sachant ce que vous faites ! Prévenez-moi des erreurs relevées ou des améliorations souhaitables, j'en tiendrai le plus grand compte dans une prochaine version de cette page.


Autres pages

Documents, références, liens















Pour un premier contact, écrivez-moi en utilisant ce formulaire.
Les correspondances suivantes pourront se faire par messagerie électronique.
Important : Si votre question a un quelconque rapport avec un travail personnel (Devoir TIPE Master...) , vous devez absolument me le préciser dès maintenant et m'indiquer très précisément les limites des informations demandées. Vous devez aussi avertir la personne qui dirige votre travail ou le corrige de cette communication et lui montrer les documents fournis.

© (Copyright) Jean-Paul Davalan 2002-2008




J-P. Liens Th. des Jeux liens Location maison vacances Île Balanec Bretagne Jeux de Nim et autres