Matlab : repousser les limites de linsolve ???

[invite6f730007]

Bonjour,

dans le cadre d'un exercice, je suis amené à résoudre avec matlab des systèmes linéaires de taille très importante (en fait je pourrais me contenter de systèmes moins importants, mais autant optimiser).

Je précise que ces matrices sont méchamment creuses (en fait des matrices d'adjacence avec des nœuds d'un degré maximum de quatre)... linsolve n'accepte plus de résoudre pour des valeurs de l'ordre de 15 000 équations à 15 000 inconnues...

Mes questions sont les suivantes :

- Existe-t-il une option de linsolve ou une autre fonction tout aussi efficace pour spécifier que la matrice est creuse et accélérer la convergence (temps correct, de l'ordre de 5 minutes pour la résolution de 6000 équations, mais le reste du temps utilisé est-négligeable par rapport à cette résolution.)

-Est-il possible de traiter des systèmes de la taille mentionnée ci-dessus avec Matlab ? Si oui, comment procéder ?



Merci d'avance
-----

[invite6f730007]

J'ai trouvé...

Il faut abandonner 'linsolve' qui n'est pas compatible avec 'sparse' et utiliser \...

Il faut envoyer les matrices du système Ax=b dans la fonction sparse puis on applique une élimination gaussienne classique...

Désolé d'ennuyer le monde à cause d'un truc aussi débile...

[lou_ibmix_xi]

Désolé d'ennuyer le monde à cause d'un truc aussi débile...

Ce n'est pas débile puisque tu t'ai posé la question... Et si quelqu'un d'autre se la pose il trouvera ta réponse donc tout le monde y gagne: c'est la loi du forum

[predator06]

Salut

Je fais un petit déterrage du Topic

J'ai écrit un programme dans lequel j'utilise linsolve.
T = linsolve (A,B)
A est une matrice 800*800
B est une matrice 800*1
et donc T est une matrice 800*1

Ma matrice A est essentiellement composée de zéro.

J'utilise ce programme pour une simulation, c'est à dire qu'à chaque pas de temps, j'appelle ce programme.
Donc j'aimerais accélérer si possible le calcul.

J'ai essayé la proposition de FinchMulateur :
- J'ai réussi à utiliser "sparse" sur ma matrice A.
- Par contre je ne comprends la fin de sa 2ème phrase "et utiliser \...". Ca veut dire quoi ?
- J'ai fais quelques recherches sur l'éminination gaussiènne proposé par Finchmulateur. J'ai trouvé quelques lignes de code que j'ai essayé dans mon programme mais mon temps de calcul a explosé.

Un peu d'aide serait la bienvenue
Merci

[A voir en vidéo sur Futura]

[lou_ibmix_xi]

Par contre je ne comprends la fin de sa 2ème phrase "et utiliser \...". Ca veut dire quoi ?

c'est la division matricielle dans l'autre sens, si tu veux résoudre le système "Ax = y", alors la solution est "x = A\y"

[predator06]

Ah oui d'accord
J'ai réussi à l'utiliser sous matlab.
Mais bon, simulink ne comprends pas les matrice sparse. Finalement, ca ne me servira pas plus que ca.

merci

[Dlzlogic]

Bonjour,
Je voudrait aborder cette question tout à fait par un autre bout : pourriez vous expliquer, si c'est pas indiscret, pourquoi vous voulez résoudre un système linéaire aussi grand. C'est à dire quel est le problème que vous devez résoudre ?