Discussion :

[Boule de coco]

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

[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 :

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

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 :

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

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

[mr_samurai]

Un coup de pouce :

Code :

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

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

++

[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é.