extraire une variable dans une fonction ODE

Bonjour,

Dans une fonction ode45, j'ai une variable que je veux afficher
une fois l'ode tourne.
Comment appeler cette variable qui se trouve dans un Mfile à partir
du workspace matlab?

Disons que ma fonction s'appelle "apple"

J'ai dans un Mfile:

Code:

---

1 2 3 4 5 6

function dx=apple(t,x) . . . G=x(1)*m*x(2) dx=[dx(1);dx(2)]

---

Ensuite quand je veux extraire la variable G, biensur elle n'est pas
connue. Je ne sais pas comment l'appeler car apple, je ne peux
l'appeler que pour faire tourner l'ode45.

Déjà pour une n'importe quelle fonction, comment peut on appeler
une variable?

Salut,

Que veux-tu faire exactement?
tu veux afficher G pendant que ode45 tourne? après que ode45 ait terminé?

Je ne comprends pas ton deuxième message...

Sinon en général pour récupérer la valeur d'une variable, on la passe en argument de sortie d'une fonction (function [sortie1,sortie2,sortie3] = mafunction()) mais j'imagine que pour la fonction ode45, ta fonction apple doit avoir une syntaxe particulière, non?

Exactement elle a une syntaxe particulière ce qui fait que si je l'appelle c'est juste pour la faire tourner (donc intégrer le x)

Code:

---

[T,x]=ode45(@apple,[0 100],x0)

---

je n'ai même pas accès à apple (qui est la fonction contenant G)

Or j'ai besoin de G en même temps que l'ode tourne, car elle est
fonction de x.

J'espère que je me suis bien expliquée.

Et même quand je fais plot (G) à l'intérieur de apple ça ne marche pas

Citation:

---

Envoyé par Sarra_

et même quand je fais plot (G) à l'intérieur de apple ça ne marche pas

---

Peux-tu être plus précis? pourquoi ça ne marche pas?

Quand je mets:

Code:

---

1 2 3

figure(1) plot(G)

---

à l'intérieur de la fonction apple, y'a rien, c'est à dire que la
figure ne se crée pas, et sans message d'erreur!

Bizarre autant qu'étrange!

Es-tu sûr de ne pas avoir plusieurs fonctions nommées apple? Tu n'as pas mis figure dans une partie du code qui ne serait pas éxécuté? (à l'intérieur d'une condition par exemple)

J'ai testé avec l'exemple de la doc à l'intérieur d'une fonction (qui s'appelle rigid dans la doc), et ça marche...

Pourrais-tu nous montrer ta fonction apple?

Bonjour,

Voici ma fonction:

Code:

---

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 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51

function dx=apple(t,x)   global J m g d k C B l0 q deltaa deltae epsilon gama R1 R2 r1 r2 k4   J=0.01; m=0.5; g=10; d=0.1; k=0.04; C=30; B=75; l0=0.2; q=0.1; deltaa=0.035; deltae=0.005; epsilon=-0.13; gama=0; r1=10; r2=10; k4=0.24;   R1=(r1*q*(-k*x(2))^epsilon)*(l0-k*x(1)-gama)*heaviside(l0-k*x(1)-gama)*heaviside(-k*x(2))+r1;   R2=(r2*q*(-k*x(2))^epsilon)*(l0+k*x(1)-gama)*heaviside(l0+k*x(1)-gama)*heaviside(k*x(2))+r2;   Ge=(R1-C*k*x(1))^2-B*k*x(2)*heaviside(k*x(2)); Gf=(R2+C*k*x(1))^2-B*k*x(2)*heaviside(-k*x(2));     if ((x(1)>0) & (x(2)>0))==1 Ge=k4*x(2); else Gf=-k4*x(2); end   if ((x(1)<0) & (x(2)>0))==1 Gf=0; end   if ((x(1)>0) & (x(2)<0))==1 Ge=0; end   dx(1)=x(2); dx(2)=(m*g/d)*sin(x(1))+(k/J)*(Ge-Gf);   figure(1) plot(Gf)   dx=[dx(1);dx(2)];

---

Et puis dans un autre Mfile je fais l'appel de fonction

Code:

---

1 2 3 4 5 6 7 8

x0=[1;1]; [T,x]=ode45(@apple,[0 30],x0); figure(2) plot(T,x(:,1),T,x(:,2))   figure(3) plot(x(:,1),x(:,2))

---

Voila, et la figure 1, ne se crée même pas!
j'ai essayé de la mettre avant les sections de conditions if, mais tjrs la meme chose.
Enfin, je crois que puisqu'elle dépend de la variable x, elle va pas se créer avant que l'ode tourne. Et dans ce cas là, comment faire pour visualiser cette variable? qq1 a une idée?

J'ai fais tourner le code que tu as donné, et pour ma part, la figure 1 se crée...

Que retourne :

Code:

---

which apple -all

---

?

Code:

---

Enfin, je crois que puisqu'elle dépend de la variable x, elle va pas se créer avant que l'ode tourne.

---

Quoi qu'il en soit, je ne comprends toujours pas bien ce que tu veux faire, et de plus je ne connais pas très bien la fonction ode45, ce qui rajoute à mon incompréhension...
néanmoins, c'est la fonction ode qui fait appelle à la fonction apple; elle fait d'ailleurs de multiples appels. Donc ce que je n'ai toujours pas compris :
- tu veux connaître les valeurs de Gf prises au cours de l'exécution de la fonction ode45?
- tu veux connaître la valeur de Gf à la fin de l'exécution de la fonction ode45?

La figure 1, contient elle une courbe ou c'est juste une figure vide?
Si oui, pourrais tu m'envoyer cette figure (et merci d'avance)
parceque pour ma part, rien ne s'affiche!! et c'est bizarre.

Concernant ta question. je veux visualiser Gf dans le même horizon
de temps sur lequel je fais tourner l'ode45, donc pendant que l'ode
tourne.
Je m'explique, je veux connaitre la valeur de Gf (c'est une force),
à chaque instant d'intégration, c'est à dire pendant que ode45 tourne
et appelle la fonction "apple". à chaque instant veut dire une courbe
d'où le plot (Gf).

Citation:

---

Envoyé par Sarra_

La figure 1, contient elle une courbe ou c'est juste une figure vide?
Si oui, pourrais tu m'envoyer cette figure (et merci d'avance)
parceque pour ma part, rien ne s'affiche!! et c'est bizarre.

---

Tu as quand même une figure qui est créé, non?
Parce que oui, ce n'est pas une courbe qui s'affiche, c'est juste un point car Gf est un scalaire, et à chaque fois que tu fais un plot, le nouveau point écrase les précédents...

Citation:

---

Envoyé par Sarra_

Concernant ta question. je veux visualiser Gf dans le même horizon
de temps sur lequel je fais tourner l'ode45, donc pendant que l'ode
tourne.
Je m'explique, je veux connaitre la valeur de Gf (c'est une force),
à chaque instant d'intégration, c'est à dire pendant que ode45 tourne
et appelle la fonction "apple". à chaque instant veut dire une courbe
d'où le plot (Gf).

---

Voir au-dessus, mais je ne vois pas comment tu veux tracer une courbe avec un seul point...

Voici une solution pour récupérer toutes les valeurs prises par Gff. Elle n'est pas très jolie car elle utilise des variables globales, et il faut éviter autant que possible d'utiliser des variables globales :
la fonctions apple devient :

Code:

---

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 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49

function dx=apple(t,x)   global J m g d k C B l0 q deltaa deltae epsilon gama R1 R2 r1 r2 k4 Gff compt   J=0.01; m=0.5; g=10; d=0.1; k=0.04; C=30; B=75; l0=0.2; q=0.1; deltaa=0.035; deltae=0.005; epsilon=-0.13; gama=0; r1=10; r2=10; k4=0.24;   R1=(r1*q*(-k*x(2))^epsilon)*(l0-k*x(1)-gama)*heaviside(l0-k*x(1)-gama)*heaviside(-k*x(2))+r1;   R2=(r2*q*(-k*x(2))^epsilon)*(l0+k*x(1)-gama)*heaviside(l0+k*x(1)-gama)*heaviside(k*x(2))+r2;   Ge=(R1-C*k*x(1))^2-B*k*x(2)*heaviside(k*x(2)); Gf=(R2+C*k*x(1))^2-B*k*x(2)*heaviside(-k*x(2));     if ((x(1)>0) & (x(2)>0))==1 Fe=k4*x(2); else Ff=-k4*x(2); end   if ((x(1)<0) & (x(2)>0))==1 Ff=0; end   if ((x(1)>0) & (x(2)<0))==1 Fe=0; end   dx(1)=x(2); dx(2)=(m*g/d)*sin(x(1))+(k/J)*(Ge-Gf);   Gff(compt) = Gf; compt = compt+1; dx=dx';

---

Et la l'appel de ode45 :

Code:

---

1 2 3 4

global Gff compt; compt = 1; x0=[1;1]; [T,x]=ode45(@apple,[0 30],x0);

---

Dans Gff tu auras alors toutes les valeurs de Gf prises au cours des apples de apple. Par contre, ce sont toutes des valeurs complexes, est-ce normal?

Citation:

---

Envoyé par magelan

Dans Gff tu auras alors toutes les valeurs de Gf prises au cours des apples de apple. Par contre, ce sont toutes des valeurs complexes, est-ce normal?

---

En fait, je ne peux pas répondre car jusqu'à maintenant je n'ai pas les valeurs
de Gff (le programme tourne encore). Il est busy depuis quelques minutes!!!

Citation:

---

Envoyé par magelan

Code:

---

1 2

dx=dx';

---

---

Enfin dx=[dx(1);dx(2)], pourquoi t'as mis dx=dx'?

Oui, il est un peu long à tourner...

Citation:

---

Enfin dx=[dx(1);dx(2)], pourquoi t'as mis dx=dx'?

---

Au début je l'avais enlevé car je me demandais pourquoi tu l'avais mis avant de me rendre compte que dx devait être une matrice colonne.
mais dx=[dx(1);dx(2)] est équivalent à dx=dx', l'opérateur ' permet de transposer une matrice.