Discussion :

[UterP]

Bonjour à tous,

je dois utiliser un intégrateur donc une fonction handle.
Les deux fonctions à determiner sont de types : u(x,t) et v(x,t)

Or, j'ai deux problèmes pour cette fonction handle :

1) Sur tous les exemples que j'ai pu trouvé sur le Web, aucun ne mentionne la possibilité d'avoir directement la variable x dans le système d'équations.

2) Comment écrire dans MatLab une ocndition aux limites du type : 3.u(1,t) + u_x(1,t) = 3 ??

(avec u_x : dérivée de u par rapport à x).

Merci de vos réponses

[salseropom]

salut,

quelle méthode utilises-tu ? Différences finies, éléments finis, volumes finis ?

[UterP]

C'est la méthode des différences finies ...

[salseropom]

Salut, dans une EDP, tu dois te donner des conditions aux limite. Tu dois te donner u(1,t) et u(L,t), si par exemple ta variable spacialle vit dans l'intervalle [1, L] (je prends un exemple en 1D...)

2*u(1,t)+u_x(1,t) = 3 <=> 2*u(1,t) + (u(1+h,t)-u(1,t))/h = 3
<=> u(1+h,t) = 3 - (2-h)*u(1,t)

où u(1+h,t) est le second point de discrétisation de ton axe.

Ai-je été clair

[rostomus]

Bonjour,

est-ce que tu as essayé la fonction pdepe?