Matrice et algorithme

Version imprimable

14/01/2008, 14h38
Boule de coco

Matrice et algorithme

Bonjour,

Je ne connais absolument pas matlab et nous avons un projet à rendre qui consiste à résoudre Ax=b pour deux matrices :
matrice 1 : n>=2, Aii=2, A(i-1,i)=A(i,i+1)=-1, Aij=0 si li-jl>=2
matrice 2 : n>=2, Aij=1/(i+j-1) , i,j=1,...,n

A est une matrice d'ordre n, b un vecteur de R^n, epsilon la précision, w un paramètre pour la matrice de préconditionnement C et K(e) le nb d'itérations pour que llr^kll soit inférieure à epsilon fois llr^0ll.

Il faut reporter dans un tableau pour chaque méthode le nombre d'itérations, le temps de calcul et le temps par nombre d'itérations.

Les méthodes utilisées sont l'algorithme SSOR; celui de Richardson et celui du gradient conjugué.

Si quelqu'un pouvait m'aider ca serait super, ou me dire comment écrire déjà tout ceci sur Matlab, je suis paumé.

Merci bcp
14/01/2008, 15h22
christophe_halgand
Bonjour,

Tout d'abord, nous ne sommes pas ici pour faire les projets des autres.

Un bon départ serait que tu nous montres ce que tu as fait.

Et voici une aide que je t'apporte :

Pour connaitre le temps qu'à pris ta machine pour calculer des lignes tu peux utiliser ceci :
Code:

1 2 3 tic bla bla bla toc
Dans la fenêtre command window matlab écrira le temps mis pour le calcul de "bla bla bla"...

Une deuxième choses, tu peux remplir une matrice comme ceci :
Code:

1 2 3 4 5 for i=1:4 for j=1:4 mat(i,j)=i*j; end end
Fais un copier/coller de ce bout de code dans matlab et tu comprendras vite...

Pur un vecteur, utilise qu'une seule boucle "for" mais n'oublie surtout pas l'indice de ton vecteur que tu crées sinon Matlab va croire que c'est un scalaire et là il remplacera sa valeur à chaque itération.

Si tu rencontres un problème particulier merci de mettre ton code + l'erreur renvoyé par matlab.

Merci

Un coup de pouce :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
n = 100;   % La taille
 
matrice1 = zeros(n,n);
matrice2 = zeros(n,n);
 
for i=1:n
  for j=1:n
    % Replire matrice1
    if j == i
      matrice1(i,j) = 2;
    elseif abs(i-j) == 1
      matrice1(i,j) = -1;
    end
 
    % Replire matrice2
    matrice2(i,j)= 1/(i+j-1);
  end
end

Et avec les conseils de christophe_halgand tu dois t'en sortir

++

15/01/2008, 18h54
Boule de coco

Bonsoir,

Merci mr_samurai et christophe_halgand pour vos conseils, c'est bon, j'ai réussi à débuter le programme.

bonne soirée

PS: le sujet peut etre fermé ou effacé.