Discussion :

[mehdi1991]

Bonjour tout le monde,

Alors voilà, je dois résoudre cet exercice sur matlab :

http://spiral.univ-lyon1.fr/ensccf/f...ts_2011-12.pdf

Projet numéro 26, c'est à dire le tout dernier.

Mon soucis est que je n'arrive pas à bien entrer les equa diff dans matlab, je reçois toujours une erreur. Voici mon code:

Ma Fonction crée:

Code :

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

1
2
3
4
5
function z = ff(t,y)
 
z(t,y) = -k*y(1)*y(2)*Ca0 - y(1)/tau
z(t,y) = -k*y(1)*y(2)*Ca0 + y(1) + (Cb0-y(2))/tau
z(t,y) = k*y(1)*y(2)*Ca0 - y(3)/tau

Mon programme:

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
%% Données de départ
 
global Ca0 Ca0p Cb0 k tau Cc y0 
Ca0 = 0.1; % Concentration initiale en A [mole/L]
Ca0p = 1; % Concentration ??
Cb0 = 0.4;
k = 2;
tau = 100;
Cc=0; % juste pour fixer valeur
 
y0 = [Ca0 Cb0 Cc]; % on définit le vecteur avec les concentrations initiales en A et B, et la conc finale en C
 
 
% graphique
 
[t,y]=ode45(@ff,tspan,y0)

Si quelqu'un saurai m'aider, car je suis complètement perdu !
Merci

[Invité]

Bonsoir,

Le lien du pdf est incorrect.

Je peux déjà dire que ta fonction ff est fausse :

• utilisant des variables Ca0, Cb0 et tau globales, il faut les redéclarer globales à l'intérieur (cf. ce message).
• ce n'est pas z(t,y), mais z(1,1), z(2,1) et z(3,1) qu'il faut assigner.

[mehdi1991]

Ah désole, je pensais que le lien fonctionnerait :

http://spiral.univ-lyon1.fr/ensccf/f...ts_2011-12.pdf

C'est le tout dernier projet .

Sinon merci de ta réponse rapide, je vais faire les changements .

[mehdi1991]

Voilà, j'ai apporté les changements.
Mais pour ma fonction, il me dit maintenant que :

Undefined function or variable 't' .

Apparemment il ne reconnait pas la variable t...

[Invité]

Pourrais-tu re-fournir le code corrigé ? t étant en paramètre d'entrée, je ne vois à priori aucune raison à cette erreur.

[mehdi1991]

Voila le code :

Code :

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

1
2
3
4
5
6
function z = ff(t,y)
 
global Ca0 tau Cb0 
z(1,1) = -k*y(1)*y(2)*Ca0 - y(1)/tau
z(2,1) = -k*y(1)*y(2)*Ca0 + y(1) + (Cb0-y(2))/tau
z(3,1) = k*y(1)*y(2)*Ca0 - y(3)/tau

et je remet le code de mon script :

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
%% Données de départ
 
global Ca0 Ca0p Cb0 k tau Cc y0 
Ca0 = 0.1;  % Concentration initiale en A [mole/L]
Ca0p = 1;   % Concentration ??
Cb0 = 0.4;
k = 2;
tau = 100;
Cc=0; % juste pour fixer valeur
 
y0 = [Ca0 Cb0 Cc]; % on définit le vecteur avec les concentrations initiales en A et B, et la conc finale en C
 
 
% graphique
 
[t,y]=ode45(@ff,tspan,y0)

[Invité]

À moins que tu ne l'aies renommée, c'est de la variable tspan de cette ligne [t,y]=ode45(@ff,tspan,y0) dont il est question, variable que tu n'as effectivement pas définie.

[mehdi1991]

Justement, je ne vois pas trop comment définir la fonction tspan ?
Je ne l'ai jamais utilisé..

[Invité]

Ce n'est pas une fonction, mais une variable, qui désigne le domaine d'arrivée souhaité.

[mehdi1991]

Ah oui désolé je voulais dire variable...
Faut t-il la définir avant la fonction ode45?

tspan = [t0 tf]

[Invité]

Avant son utilisation oui

[mehdi1991]

Ok ok merci .

Mais dans mon cas, que je mette ça ou par exemple tspan[1 50], ça ne fonctionne toujours pas :s.
Enfin tu saurai me montre ce qui ne va pas dans mon code? car j'ai toujours un problème :s.
Merci

[Invité]

Problème que je dois encore deviner ?
Si c'est un autre undefined sur la variable k, j'espère que tu a quand même su le corriger tout seul ?

[tachmou]

Bonjour,

tspan correspond à ton "temps".

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
% Resolution des equations
 
t0 = 0; % instant initial en secondes
tf = 10; % instant final en secondes
dt = 0.01; % intervalle de temps
 
tspan = [t0:dt:tf]; % definition de ta base temporelle
 
[t,y]=ode45(@ff,tspan,y0); % resolution
 
% Graphique
 
plot(t,y);

Cordialement,