IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

Pascal Discussion :

Matrice factorisation LU


Sujet :

Pascal

  1. #1
    Nouveau membre du Club
    Inscrit en
    Mars 2005
    Messages
    83
    Détails du profil
    Informations forums :
    Inscription : Mars 2005
    Messages : 83
    Points : 34
    Points
    34
    Par défaut Matrice factorisation LU
    Salut a tous;
    Je cherche un programme pascal qui fait l'élimination de Gauss pour donner la factorisation LU si elle existe, bien sur cette discutions est destinée a des mathématicien informaticien
    Merci d'avance.

  2. #2
    Membre averti Avatar de Tux++
    Étudiant
    Inscrit en
    Avril 2008
    Messages
    281
    Détails du profil
    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Avril 2008
    Messages : 281
    Points : 379
    Points
    379
    Par défaut
    Bonjour,

    si tu cherches un code je ne vois pas en quoi elle est destinée à des informaticiens mathématiciens.
    si par chercher tu entends qu'on te le fasse alors je comprends, mais non (bien qu'il ne faille pas être méthématicien pour utiliser une décomp LU)
    Sinon ce n'est pas très difficile à réaliser et ça procure un bon entrainement sur les matrices.
    Étant donné une matrice de dimension NxN
    A= (a indice n,n);
    on définit
    A^(0) := A;
    et on itère pour n = 1,...,N-1 comme ceci:
    sur nième colonne de A^(n-1), on élimine les éléments sous la diagonale en ajoutant à la ième ligne de cette matrice, la nième ligne multipliée par
    l(i,n) := -a(i,n)^(n-1)/a(n,n)^(n-1)
    pour i = n+1..N. Ceci peut être fait en multipliant par la gauche A^(n-1) avec la matrice triangulaire inférieure Ln
    A^(n) := Ln A^(n-1);
    Après N-1 itérations, on a éliminé tous les éléments sous la diagonale, par conséquent, tu as maintenant une matrice triangulaire supérieure A^(N-1).
    on a la décomposition

    A = L(1)^(-1)*L(1)A^(0) = L(1)^(-1) *A^(1) = L(1)^(-1) * L(2)^(-1)* L(2)*A^(1) = ... = L(1)^(-1)*...*L(N-1)^(-1)*A(N-1)

    Ensuite pour U beh, c'est la matrice triangulaire supérieure A(N-1) et L=L(1)^(-1)*...*L(N-1)^(-1).
    tu obtiens A = LU.
    en vérifiant bien que les pivots soient non-nuls, tu as ton algo

  3. #3
    Nouveau membre du Club
    Inscrit en
    Mars 2005
    Messages
    83
    Détails du profil
    Informations forums :
    Inscription : Mars 2005
    Messages : 83
    Points : 34
    Points
    34
    Par défaut Re
    Merci pour tout ces éclaircissement mais moi je cherche l'algorithme sinon je le ferai moi même mais j'ai pas le temps car je suis surmené et sa demande comme un temps de réflexion, en suite c’est a ça que serre les forum, alors si tu peut m’aider j’en est vraiment besoin car je doit en suite le faire en Delphi plus facile je te l’accorde

  4. #4
    Membre averti Avatar de Tux++
    Étudiant
    Inscrit en
    Avril 2008
    Messages
    281
    Détails du profil
    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Avril 2008
    Messages : 281
    Points : 379
    Points
    379
    Par défaut
    salut,

    la factorisation LU EST un algorithme...
    Si tu veux qu'il soit plus formel, propose-en un (bien qu'internet en regorge), et on te dira ce que tu peux éventuellement améliorer.

    On est pas là pour faire le boulot à ta place, mais pour t'aider si tu as un problème concret, voilà à quoi sert un forum

  5. #5
    Nouveau membre du Club
    Inscrit en
    Mars 2005
    Messages
    83
    Détails du profil
    Informations forums :
    Inscription : Mars 2005
    Messages : 83
    Points : 34
    Points
    34
    Par défaut Re
    J'ai trop cherché mais je trouve pas en pascal, je veut pas qu'on le face a ma place mais je suppose que c'est un programme qui existe déjà dans le net je veut juste qu'on m'oriente vers un lien ou bien le copier, mais le faire moi même je croix pas que je pourrais alors voila, si vous avez une idée ou je peut le trouver sa m'aidera bcp. en plus je doit le remettre demain donc vous comprendrai que je suis serré dans le temps.

  6. #6
    Expert confirmé
    Avatar de krachik
    Inscrit en
    Décembre 2004
    Messages
    1 964
    Détails du profil
    Informations forums :
    Inscription : Décembre 2004
    Messages : 1 964
    Points : 4 015
    Points
    4 015
    Par défaut
    Bonjour
    Citation Envoyé par amireve
    .... je veut pas qu'on le face a ma place ...
    Alors je pense que tu n'as pas très bien compris(Les règles ) et comme Tux++ te l'a dit précedemment on ne fera pas ton travail à ta place; non seulement ça ne te rendra pas service mais en plus ça sera ne pas faire d'effort.
    Tout ce que tu as à faire c'est de faire une petite recherche sur ce que tu cherches et vu que Tux++ t'a deja donné un pseudo algo il faut juste l'ecrire proprement sous forme de vrai algo je peux te dire que la traduction en Pascal ne devrait pas poser problème
    je veut juste qu'on m'oriente vers un lien ou bien le copier, mais le faire moi même je croix pas que je pourrais alors voila
    Moi je suis sûr que tu pourras et que tu peux aussi faire des efforts

    Bon courage
    Cordialement

  7. #7
    Nouveau membre du Club
    Inscrit en
    Mars 2005
    Messages
    83
    Détails du profil
    Informations forums :
    Inscription : Mars 2005
    Messages : 83
    Points : 34
    Points
    34
    Par défaut Re
    Décidément ;
    si c'est possible donner moi un utilitaire qui traduire c++/c# vers pascal
    sinon donner moi les mot que je dois taper pour la recherche car moi j'ai essayé sans succès en suite je suis débordé de travail et j'ai pas que ça j'ai fait le déterminant qui ma coûté un temps énorme et le produit matricielle ... Ax=B... il me reste que cette partie de factorisation que j'ai pas pu la faire alors si je demande un cou de pousse ça va m'éviter de donner mon compte rendu amputé d'une question alors comprenait moi

  8. #8
    Expert confirmé
    Avatar de krachik
    Inscrit en
    Décembre 2004
    Messages
    1 964
    Détails du profil
    Informations forums :
    Inscription : Décembre 2004
    Messages : 1 964
    Points : 4 015
    Points
    4 015
    Par défaut
    Bon je reprend : l'effort qu'on te demande c'est vu que Tux t'a deja donné un pseudo algo reécris le proprement et tu obtiendras un algo ensuite tu peux le traduire en Pascal.Une fois que tu as un algo sous la main le traduire est plus rapide (la bonne pratique c'est partir d'un algo pour en faire du code dans un langage et non le contraire )

    On est prêt à t'aider si tu essaies de traduire cet algo en code; et pour la recherche te fournira plusieurs liens vers des algo mais vers du code pascal ça serait trop facile

    Cordialement

  9. #9
    Nouveau membre du Club
    Inscrit en
    Mars 2005
    Messages
    83
    Détails du profil
    Informations forums :
    Inscription : Mars 2005
    Messages : 83
    Points : 34
    Points
    34
    Par défaut Re
    merci comeme de vos réponse, je dois me rendre a l'évidence les forum ne sont pas la pour faire le travail a notre place, tempi

  10. #10
    Membre régulier Avatar de Transgarp
    Profil pro
    Inscrit en
    Février 2008
    Messages
    99
    Détails du profil
    Informations personnelles :
    Localisation : Canada

    Informations forums :
    Inscription : Février 2008
    Messages : 99
    Points : 76
    Points
    76
    Par défaut
    J'ai réalisé ce petit programme codé BP7 en 1990.

    Il fonctionne pour des matrices de 9 x 9 et se compose d'à peine 158 lignes de code.

    Et je ne dois pas le montrer ici

  11. #11
    Nouveau membre du Club
    Inscrit en
    Mars 2005
    Messages
    83
    Détails du profil
    Informations forums :
    Inscription : Mars 2005
    Messages : 83
    Points : 34
    Points
    34
    Par défaut Re
    pitié juste pour cette foi

  12. #12
    Modérateur
    Avatar de ToTo13
    Homme Profil pro
    Chercheur en informatique
    Inscrit en
    Janvier 2006
    Messages
    5 793
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 45
    Localisation : Etats-Unis

    Informations professionnelles :
    Activité : Chercheur en informatique
    Secteur : Santé

    Informations forums :
    Inscription : Janvier 2006
    Messages : 5 793
    Points : 9 860
    Points
    9 860
    Par défaut
    Bonjour,

    si tu veux l'algorithme et une "idée" d'implémentation, tu peux regarder la bible du domaine qui est le "numerical recipse".
    La plupart des algorithmes de calculs matriciel sont dedans.

Discussions similaires

  1. Factorisation LU matrice tridiagonale par blocs
    Par djocin dans le forum Fortran
    Réponses: 3
    Dernier message: 26/10/2009, 22h27
  2. [CR] entête et pied sur page 1/B de matrice
    Par chloe.j3 dans le forum SAP Crystal Reports
    Réponses: 2
    Dernier message: 04/09/2002, 13h07
  3. Matrice de filtrage ?
    Par gimlithedwarf dans le forum Traitement d'images
    Réponses: 2
    Dernier message: 24/08/2002, 10h44
  4. Gestion de matrice
    Par bzd dans le forum C
    Réponses: 4
    Dernier message: 12/08/2002, 19h19
  5. Comment définir le type matrice ?
    Par charly dans le forum Langage
    Réponses: 7
    Dernier message: 15/06/2002, 22h01

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo