Version imprimable
Bonjour, je sais je deviens chiant mais bon!
Si quelqu'un sait où je pourai trouver une fonction réalisant une pseudo-inversion, qui plus est sur des nombres complexes! J'ai réalisé une fonction qui fait cela, avec l'aide de monsieur vip mais c'est très lent. Donc merci!
:lol:
T'entends quoi par "pseudo-inversion"?
Une pseudo inverse est une inversion de matrice non carrée d'après la formule S+=(Sh S)^-1*Sh, S+ est la pseudo inverse de S, Sh l'opérateur hermitien sur S, l'équivalent d'une transposée de conjuguée.
http://perso.wanadoo.fr/marco.butte/images/shadok1.gifCitation:
Envoyé par petdelascar
http://gigi46.canalblog.com/images/drole078.jpg
:lol:
Bon!
salut, je n'ai pas d'algo à te proposer, mais je sais qu'il est déconseillé d'inverser une matrice. Pour résoudre un système linéaire (réel ou complexe), tu peux essayer la méthode de jacobi, gauss seidel , relaxation, LU, etc... Tout dépend des propriétés de ta matrice transposée(A) * A (qui est donc carrée, symétrique donc diagonalisable et si je ne m'abuse, toutes les valeurs propres sont réelles).
Je ne crois pas qu'il y ait énormément de solutions à ton problème. Le calcul de l'inverse a été suffisament étudié pour dire que la méthode à utiliser est généralement un pivot de Gauss...Citation:
Bonjour, je sais je deviens chiant mais bon!
Si quelqu'un sait où je pourai trouver une fonction réalisant une pseudo-inversion, qui plus est sur des nombres complexes! J'ai réalisé une fonction qui fait cela, avec l'aide de monsieur vip mais c'est très lent. Donc merci!
Si tu veux tenter de gagner du temps et si tes matrices sont des matrices 3*3 ou 4*4, fais travailler ta carte graphique, son processeur est conçu pour des calculs matriciels... Ca pourrait te faire gagner du temps... Par contre si tu joues avec des nombres complexes, alors il faudra le faire à la main...
Y'a une autre technique pour l'inversion matricielle : http://www.sciencedirect.com/science?_ob=GatewayURL&_method=citationSearch&_uoikey=B6TYH-43WV09W-J&_origin=SDEMFRHTML&_version=1&md5=2ec2327386691e7d860ddf73c3c531a9
Elle est en gros en N³
ou alors tu utilises lapack la librairie disponible sur le net qui est utilisé par Matlab pour resoudre tout les calculs matriciels.
Elle est en fortran mais tu trouveras difficillement plus efficace.
XXiemeciel