Discussion :

[barns_89]

Bonjour,

j'aimerai résoudre une équation différentiel avec le solver ode45, mais il faut que tout mon code soit dans un seul script, je sais que c'est possible mais je n'ai pas trouvé comment

Pour l'instant mon code est séparé en deux script :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
t_init=0;
t_final=10;
x1_0=1;
x2_0=1;
 
options_struct = odeset('RelTol', 1E-9, 'AbsTol', 1E-6);
[t, x] = ode45(@func,[t_init t_final],[x1_0 x2_0], options_struct); 
 
plot(t,x(:,1),'-',t,x(:,2),'.'); 
 
inits='x1(0)=1,x2(0)=1';
[x1,x2]=dsolve('Dx1=-x1+x2','Dx2=-x2',inits)

et la fonction :

Code :

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

1
2
3
4
5
function x_dot = func(t, x)
x_dot(1) = -x(1) + x(2);  
x_dot(2) = -x(2);
 
x_dot = x_dot';

mon problème se trouve surtout dans la syntaxe :s

Quelqu'un pourrait m'aider svp? c'est un peu pressant

[le fab]

techniquement ton problème est défini dans deux fichiers, pas dans deux scripts : il y a un script et une fonction
si ton pb est de n'avoir qu'un fichier, alors il te suffit de transformer ton script en fonction et de mettre l'autre fonction à la fin du même fichier

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
function main
t_init=0;
t_final=10;
x1_0=1;
x2_0=1;
 
options_struct = odeset('RelTol', 1E-9, 'AbsTol', 1E-6);
[t, x] = ode45(@func,[t_init t_final],[x1_0 x2_0], options_struct); 
 
plot(t,x(:,1),'-',t,x(:,2),'.'); 
 
inits='x1(0)=1,x2(0)=1';
[x1,x2]=dsolve('Dx1=-x1+x2','Dx2=-x2',inits)
 
function x_dot = func(t, x)
x_dot(1) = -x(1) + x(2);  
x_dot(2) = -x(2);
 
x_dot = x_dot';

et tu peux aussi défini ta fonction comme fonction anonyme dans ton premier script :

Code :

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

func = @(t,x) [-x(1)+x(2);-x(2)];

[barns_89]

Merci bcp pour la réponse!

Et encore une petite question en passant. Est-ce que vous savez si on peut choisir l'intervalle de temps entre lesquels notre système évalué?

[tachmou]

Salut !

"Est-ce que vous savez si on peut choisir l'intervalle de temps entre lesquels notre système évalué?"

==> Euh ?? Sauf erreur de ma part, tu les changes avec ton "t_init" et ton "t_final". Ou bien j'ai loupé un truc ? Si tu veux spécifier un pas de temps spécifique "dt", tu peux utiliser [t_init:dt:t_final] à la place de [t_init t_final].

Par ailleurs, je te conseille d'utiliser des fonctions imbriquées : plus de clarté (avis personnel)...

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
function main
    t_init=0;
    t_final=10;
    x1_0=1;
    x2_0=1;
 
    options_struct = odeset('RelTol', 1E-9, 'AbsTol', 1E-6);
    [t, x] = ode45(@func,[t_init t_final],[x1_0 x2_0], options_struct); 
 
    plot(t,x(:,1),'-',t,x(:,2),'.'); 
 
    inits='x1(0)=1,x2(0)=1';
    [x1,x2]=dsolve('Dx1=-x1+x2','Dx2=-x2',inits)
 
    function x_dot = func(t, x)
        x_dot(1) = -x(1) + x(2);  
        x_dot(2) = -x(2);
        x_dot = x_dot';
    end
end

Petite question : pourquoi as-tu besoin de le mettre dans un seul fichier ?

Cordialement,