Discussion :

[The_FD]

bonjour à tous

je suis sur une boucle qui contient deux équa diff assez simples à résoudre, et à chaque fois la condition initiale de l'une est la valeure finale de l'autre

Jusque là tout va bien sauf qu'avant même de faire le premier tour il me sort

??? Error using ==> ode45
Solving eq2 requires an initial condition vector of length 0

je vois pas ce que je peux faire pour sortir de là (un prob avec les global? un clear all ne fait pas avancer le schmilblik...)

le code de la boucle:

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
t=1;
while t<10
 
   [Lcpp,P]=ode45('eq1',[0 Lcp],[Pv]);
 
   figure(1)
   plot(Lcpp,P)
   Pa=P(end);
 
   [Lcmm,P]=ode45('eq2',[0 Lcm],[Pa]);
 
   figure(2)
   plot(Lcmm,P)
   Pv=P(end);
 
end

Auriez vous des idées pour me sortir d'affaire?

[salseropom]

Que vallent Pa et P ?

[The_FD]

Avant d'entrer dans la boucle je fixe Pv arbitrairement à 100

Après la première équa diff Pa=100 (ce qui est tout à fait ce que j'attends avec Pv=100) par contre à partir de la deuxième boucle ça doit changer...

[salseropom]

Peux-tu envoyer ton code en entier ?

[The_FD]

depuis j'ai résolu le problème du "initial condition vector".
C'est un des parametres de l'eq2 (Ptm qu'il s'appelle ce coquin ) qu'il fallait declarer en global. Ce parametre évolue quand la boucle tourne et que eq1 et eq2 se résolvent l'une après l'autre. Le seul hic c'est qu'en declarant ce parametre en tant que global le modèle devient insensible à ce même parametre et ça ne m'arrange pas du tout mais alors pas du tout!

Je n'ai pas encore de stratégie pour sortir de ce cercle vicieux, peut être faut-il que je crée 2 Ptm: un pour permettre au modèle de tourner et l'autre pour la sensibilité du modèle... qu'est ce que vous en pensez?

merci de l'aide en tout cas