Discussion :

[antonb]

Bonjour,

je sais comment créer une matrice creuse diagonale; mais je ne parviens pas à créer une matrice creuse diagonale par blocs :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
% Création de la matrice B
 
U = ones(n,1);
B = spdiags([-U,4*U,-U],-1:1,n,n);
display(full(B));
 
% Création de la matrice A
 
V = eye(n,n);
display(-V);
A = spdiags([-V,B,-V],-1:1,n*n,n*n);
display(full(A));

Il y a une erreur à la ligne où je définis A,

Index exceeds matrix dimensions.

En lisant l'aide j'ai bien vu que spdiags ne permet pas (ou alors je ne vois pas comment) de créer une matrice diagonale à l'aide de blocs, alors je vous sollicite si vous avez une idée !
Ce que je cherche à avoir, c'est que la diagonale par bloc principale de A contienne les blocs B, et la sur et sous diagonale les blocs -V.

Merci d'avance !

[Invité]

Bonjour,

Je ne suis pas sûr de bien comprendre

Ce que je cherche à avoir, c'est que la diagonale par bloc principale de A contienne les blocs B, et la sur et sous diagonale les blocs -V.

et surtout le terme de sous-diagonale...
Cherches-tu à avoir ceci:

Code :

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

full([B,-V; -V B])

Si ce n'est pas ça, aurais-tu un petit exemple pour n=2?

[antonb]

Et bien soit A une matrice d'ordre n^2
Soit B une matrice d'ordre n
Soit In la matrice identité d'ordre n
Ce que je cherche à avoir c'est

A = 

B   In  0   0   0   0   0   0   0
In  B   In  0   0   0   0   0   0
0   In  B   In  0   0   0   0   0
0   0   In  B   In  0   0   0   0
0   0   0   In  B   In  0   0   0
0   0   0   0   In  B   In  0   0
0   0   0   0    0  In  B   In  0
0   0   0   0    0   0  In   B  In
0   0   0   0    0   0   0   In  B 

(En fait au post précédent j'avais dit -In à la place de In mais peu importe, c'est la forme générale qui compte)
Edit : il est assez important que la matrice obtenue soit creuse parce que je risque de travailler avec n assez grand et j'ai peur que les temps de calculs soient affreux avec une matrice non creuse...
J'ai modifié A, donc ici j'ai pris n=3 et A est donc de dimension 9

[FR119492]

Salut!

mais je ne parviens pas à créer une matrice creuse diagonale par blocs :

Pour le plaisir ou pour résoudre un problème particulier et, dans ce cas, lequel?
Jean-Marc Blanc

[antonb]

Bonjour !
J'implémente une méthode de résolution par différences finies et ma matrice pouvant atteindre de grandes tailles je préfère utiliser des matrices creuses, afin de limiter les temps de calculs

[FR119492]

Salut!
Alors, si tu veux qu'on puisse t'aider, tu dois poser ta question en commençant par le commencement:

Tu veux résoudre une équation de Laplace ou de Poisson sur un domaine bidimensionnel borné.

Alors, première question de ma part: quelle est la forme de ce domaine?
Jean-Marc Blanc