Discussion :

[Oistrakh]

Bonjour à tous,
Voilà j'avoue que ça fait un peu de temps déjà que j'ai plus touché à Matlab, j’ai retrouvé quelques bases mais là je suis totalement planté et donc j'espère que vous pourrez m'aider, sans doute que ça sera du gâteau pour vous.

J'ai fait un algorithme pour résoudre un système linéaire simple Ax=b par itération ou l'on donne en entrée, évidemment nos deux matrices A et b ainsi que epsilon (pour que l'itération s'arrête bien sûr)
Mon programme renvoie le vecteur x, le nombre d'itération.

Mais je voudrais tracer un graphe ou j'aurais en ordonné le nombre d'itération et en abscisse l'epsilon (à l'échelle logarithmique ça semble plus judicieux)
Ca sembla facile mais je sais plus du tout comment on fait (même avec l'aide j'ai pas retrouvé ) Je vois déjà pas comment je pourrais avoir plusieurs points alors que le programme renvoie à chaque fois une itération (donc un point sur le graphique )

[Oistrakh]

Personne pour m'aider?
je crois que c'est semilogs pour les tableaux avec une absicce d'echelle logarithmique si je me souviens bien non?

Merci

[Cheetor]

Tu as le choix entre

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

plot(epsilon, iteration)

et

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

semilog(epsilon, iteration)

Tu peux faire également des graphiques semi-logarithmiques avec
semilogy et semilogx.

[L'hermine]

En fait si j'ai bien compris, tu as fait une fonction qui te renvoi "epsilon" et le nombre total d'itérations, disons "e" et "n" ? C'est ça.
Si tu veut faire un graphe, dans ta fonction il faut que tu stocke chacune des valeurs d'e et de n.
Tu dois probablement avoir une boucle while dans ta fonction que tu stoppes quand "e<e0", non ?
Et bien dans cette while au lieu d'écrire e = ... et n = ..., tu écris simplement e(i) =... et n(i) = ...
Comme ça tu renverras deux vecteurs contenant toutes les infos que tu veux, sachant que le dernier élément des deux vecteurs(e(length(e) et n(length(n)) contient l'info que tu renvoyais précédemment.
Je ne sais pas si j'ai été très clair...

[tachmou]

Bonjour,

Les deux réponses précédentes te donnent les réponses à tes questions.

Par contre, pourquoi fais-tu ça ? Ton objectif est de développer un algorithme qui te permette de résoudre ce genre de problème, ou bien tu as juste besoin de résoudre cette équation pour faire autre chose après ? Si tu es dans le deuxième cas, alors tu as juste besoin de faire :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

x = inv(A)*b

Ou, mieux dans certains cas :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

x = pinv(A)*b

Cordialement,