Discussion :

[feynman]

Bonjour,
J ai un truc que j arrives pas a le resoudre , merci pour m aider svp
Je veux ecrire une matrice tridiagonale par bloc (sur scilab) mais de dimension N c ad la matrice est de la forme suivante

A D 0......0
D A D 0...0
.
.
0....0 D A D
0.........D A


Merci

[Jerome Briot]

Une solution sous MATLAB est donnée dans le tutoriel Gestion des matrices : Génération de matrices usuelles

Une autre solution consiste à utiliser la fonction DIAG :

Code :

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

1
2
3
4
5
A = 2
D = 3
 
m = 7;
M = diag(ones(1,m)*A)+diag(ones(1,m-1)*D,-1)+diag(ones(1,m-1)*D,1)

[feynman]

Atention A et D ne sont pas des nombres, mais des matrices aussi!

[FR119492]

Salut !
Si tu écris toute ta matrice dans un seul tableau à 2 indices, tu gaspilles une place énorme à stocker des termes nuls. Il est bien préférable d'utiliser 2 tableaux à 3 indices, l'un pour les A et l'autre pour les D. Un indice correspond à la matrice A ou D que tu es en train de traiter, et les deux autres, la position d'un terme à l'intérieur d'une de ces matrices.
Jean-Marc Blanc

[feynman]

Salut Jean Marc,
tu peut me donner un petit exemple ou le code svp?
merci

[Jerome Briot]

Atention A et D ne sont pas des nombres, mais des matrices aussi!

Cette information aurait dû être précisée dans le premier message.

Un exemple pour des matrices A et D carrées :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
N = 4;
 
A = ones(N);
D = ones(N)*2;
 
for n=1:3
 
    r = (n-1)*N+1;
    c = (n-1)*N+1;
    M(r:r+2*N-1,c:c+2*N-1) = [A D ; D A];
 
end
 
M

Ce code devrait être facilement modifiable pour accepter des matrices A et D de taille quelconque...

Salut Jean Marc,
tu peut me donner un petit exemple ou le code svp?
merci

Je pense que Jean-Marc voulait parler des matrices creuses (voir la fonction SPARSE)

[FR119492]

Salut!

tu peut me donner un petit exemple ou le code svp?

Malheureusement, je ne l'ai jamais fait dans MatLab, parce que, lorsque les problèmes dépassent une certaine taille, je reviens toujours à mon bon vieux Fortran.
Jean-Marc Blanc