Discussion :

[carlosmoyas]

bonjour

mon pb consiste à résoudre une EDP avec la méthode de différences finies.
à la fin de la discrétisation j'ai trouvé une matrice creuse
A=
-4 1 0 1 0 . . . . . . . . . . . . . .
1 -4 1 0 1 0 . . . . . . . . . . . .
0 1 -4 1 0 1 0 . . . . . . . . . .
. . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . 1 -4 1
. . . . . . . . . . . . . . . . 0 1 -4



et je veux présenter cette matrice sur matlab avec des boucles sachant que sizes (A)=(50,50).

merci

[duf42]

Bonjour,

Pourrais-tu être plus précis sur ce que tu souhaites faire avec cette matrice?

As-tu déjà un code?

Duf

[carlosmoyas]

je veux seulement écrire cette matrice sur matlab .

[Jerome Briot]

Les termes "écrire" et "présenter" ne sont pas assez précis.

Cherches-tu à créer la matrice ou bien à en afficher le contenu dans le Command Window ou dans un fichier ?

[carlosmoyas]

je cherche à créer la matrice
merci

[Jerome Briot]

Voici une idée :

Code :

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

1
2
3
4
5
6
7
N = 50;
 
A = diag(ones(1,N)*-4,0)+ ...
    diag(ones(1,N-1),1)+ ...
    diag(ones(1,N-1),-1)+ ...
    diag(ones(1,N-3),-3)+ ...
    diag(ones(1,N-3),3)

[carlosmoyas]

la partie éléments finis me pose problème ,surtout comment choisir les fonctions tests :
cours : ils faut représenter les fonctions de Vh par des fonctions de base très simples (fonctions chapeau )