Discussion :

[Jeremy971]

Bonjour,
je suis en train de faire un TP sur les équations différentielles et j'ai du mal a voir par ou commencer et comment coder sous Matlab. J'essaie tout d'abord de rentrer les données du probleme mais je bloque a sur l'importation de ces données a partir d'un fichier .txt
Voici mon code :

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
%%Initialisation
h=0.005;
k=h;
xi=zeros(1,101);
yj=zeros(1,61);
uij=zeros(61,101);
xi(1)=0;
for i = 1:100
    xi(i+1) = i*h;
end
yj(1)=0;
for j=1:60
    yj(j+1)= j*k;
end
chargement = fopen('matlab.mat');
donnee = fscanf(chargement,'%g %g %g',[3 inf]);

En fait dans ce fichier .txt, il y a les valeur de f(x,y), j'essaie de les importer pour pouvoir l'utiliser avec mes équations aux différences. Pouvez-vous m'aider merci ?

[Jerome Briot]

Utilise plutôt les fonctions load ou dlmread pour lire le contenu du fichier.

[Jeremy971]

Ok donc j'obtiens une matrice 6161*5 :
i j x y f(x,y)
avec
i va de 1 a 101 (valeur)
j va de 1 a 61 (valeur)
Ce que j'aimerais c'est créer un maillage de ma plaque allant de i=1:101 et j=1:61, le probleme c'est que j'ai 6161 lignes pour i et j. Et je sais pas comment dire a Matlab quand j'ai i=1 et j=2, tu me donnes tel x, tel y, et tel f(x,y), vu que ma colonne i et j font 6161 lignes.

Voila ce que j'ai codé :

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
19
20
21
22
23
24
25
26
27
%%Initialisation
clear all
clc
M=dlmread('chargement.txt');
h=0.005;
k=h;
n = 5000
%%i=M(1:6161,1);
%%j=M(1:6161,2);
xi=M(1:6161,3);
yj=M(1:6161,4);
f=M(1:6161,5);
u=zeros(6161,6161);
E=10^-5;
N=5000;
%%Conditions limites
u(1:6161,6161)=0;
u(1:6161,6161)=0;
u(1,1:6161)=0;
u(6161,1:6161)=0;
for l = 1:n
    for i = 2:6161
        for j = 2:6161
            u(1,j) = 1/21 * ( 0.005 * f(1:61) - ( u(i+2,j) + u(i,j+2) + u(i,j-2) - 8*u(i+1,j) - 8*u(i-1,j) - 8*u(i,j+1) - 8*u(i,j-1) + 2*u(i+1,j+1) + 2*u(i+1,j+1) + 2*u(i-1,j) + 2*u(i-1,j-1)));
        end
    end
end

[Jeremy971]

Voici mon code final, j'ai reussi a resoudre la plupart des problèmes. La je bloque juste sur une partie qui consiste a dessiner un graphique avec la commande contour (a la fin de mon code), j'ai rentre mon code mais rien ne sort. Est-ce que vous pouvez m'aider ?

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
%%Construction du graphique
m = 101;
n = 61;
X = zeros(m,n);
Y = zeros(m,n);
for j = 1 : n
    X(:,j)= xi(:);
end
for j = 1 : m
    Y(j,:) = yj(:)';
end
[C,h] = contourf(X,Y,u');
clabel(C,h);