Discussion :

[fgallais]

Bonjour,

Nous essayons de résoudre un système d'EDP d'ordre 1 sur MATLAB. Le voici :




Pour l'instant on a fait les fichiers suivants :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
function [c,f,s] = pdex2pde(x,t,u)
gamma=0.2;
delta=0.05;
d1=-1;
 
d2 = -1;
 
c = [1; 1]; 
f = [d1; d2];
f1 = -gamma.*u(1);
f2 =-delta*u(2) - f0*u(2)/2;
s = [f1; f2];

Code :

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

1
2
3
4
5
6
7
8
9
10
 
m = 0;
 
 
x=[0:1:100];
t=[0:0.1:10];
 
sol = pdepe(m,@pdex2pde,@pdex4ic,@pdex4bc,x,t);
u1 = sol(:,:,1);
u2 = sol(:,:,2);

On ne sait pas si cela est correct et on ne sais pas comment ecrire les conditions aux bords car celle ci dépendent de nos variables p et r.


Merci