[ode45] out of memory !

Bonjour à tous !

Je dois faire un code avec ode45, mais il est absolument trop long. Il fonctionne car si mon intervalle de temps est assez petit ça fonctionne, du moins j'ai un bout de courbe, mais lorsque je prend l'intervalle de temps complet, il m'indique

Citation:

---

Out of memory

---

Voici le code:

De ma fonction

Code:

---

1 2 3 4 5 6 7 8 9

function H=G(t,Y) T=5700; lb=1000; l=log(2)/T;   H(1)=-l*Y(1); H(2)=l*Y(1)-lb*Y(2);   H=H';

---

Et mon script

Code:

---

1 2 3 4 5 6 7 8 9 10 11

clear all close all clc   No=100; T=5700; lb=1000;   [t,Y]=ode45(@G,[0 5*T],[No 0]);   plot(t,Y(:,1),t,Y(:,2))

---

Comment résoudre ce soucis ?

En vous remerciant ! ;)

Le problème est que le pas de temps est trop fin (5*10^-4 s) et que comme ta simulation est longue, tu as beaucoup trop de résultats.
Tu peux toujours essayer de changer les paramètres du solveur (voire changer de solveur) en indiquant un pas de temps maximal ou prendre un solveur discret où tu pourras régler le pas de temps comme tu le souhaites.

Maters

Materss a sans doute raison

Que renvoi ceci avant le tracé ?

Code:

---

whos t Y

---

Ou dois je le mettre le "whos t Y" ?

Dans mon script ? Car ça ne me renvois rien, quelle que soit sa position ...

J'ai regardé dans les options de ode45 pour choisir le pas, mais je n'ai rien trouvé :roll:

Ok l'erreur est humaine ... ma constant lambdaB, lb dans ma fonction, devait être beaucoup plus petite que 1, alors que le j'ai mise beaucoup plus grande ......

Donc ça fonctionne !

Mais j'en profite pour vous poser une autre question:

Dans ma fonction, j'ai lb qui est fixe, or je doit regarder l'influence de ce paramètre, donc le faire varier.

Si je fait une fonction

Code:

---

function H=G(t,Y,lb)

---

Puis dans mon scripte

Code:

---

[t,Y]=ode45(@G,[0 10*T],[No 0],lb);

---

Afin de faire varier lb aisément, ode45 ne veut rien savoir ...

Suis je clair ? :?

Bonsoir,

Voir la FAQ Que représente le symbole @ ? 3ème point.

Cool ! Ca marche ! J'ai cherché 30min mais ça fonctionne :mouarf:

Voilà au final:

La fonction:

Code:

---

1 2 3 4 5 6 7 8

function H=G(t,Y,lb) T=5700; l=log(2)/T;   H(1)=-l*Y(1); H(2)=l*Y(1)-lb*Y(2);   H=H';

---

Le script:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21

clear all close all clc   No=1000; T=5700;   for i=1:10; lb=(i/2)*10^(-4);   [t,Y]=ode45(@(t,Y) G(t,Y,lb),[0 10*T],[No 0]); S=Y(:,1)+Y(:,2);   [t1,y]=ode45(@f,[0 10*T],No);   plot(t,Y(:,1),t,Y(:,2),t,S,t1,y,'+') title('Désintégration multiple de noyaux radioactifs') legend('Element A','Element B (A->B)','somme de la désintégration multiple','Désintégration simple de A')   pause(0.7) end

---

Encore une question, qui me tracasse, pourquoi à la fin de la fonction je dois mettre

Code:

---

H=H'

---

alors que l'an passé, avec un autre enseignant, on aurait mis ça:

Code:

---

H=[H(1);H(2)]

---

?

Pourtant les deux fonctionne dans le cas présent. Je comprend pas ... Dans le premier cas je transpose, alors que dans le deuxième non ?

Parce que

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14

H(1) = 1; H(2) = 2; >> H   H =   1 2   >> H.'   ans =   1 2

---

et

Code:

---

1 2 3 4 5 6

>> H = [1 ; 2]   H =   1 2

---

La documentation stipule :

Citation:

---

The function, f = odefun(t,y), for a scalar t and a column vector y, must return a column vector f corresponding to f(t,y).

---

Ok ! C'est aussi simple que ça ...

Je te remercie !

A bientôt, je pense !

Ce qui revient simplement à

Code:

---

1 2 3 4 5 6

function H = G(t,Y,lb) T = 5700; l = log(2)/T;   H(1,1) = -l*Y(1); H(2,1) = l*Y(1)-lb*Y(2);

---

Ok !

Et une question rapide pour la route:

Hier j'ai vu une proposition de stage, je vais surement y candidater, mais le langage est en Python. Est ce très différent de matlab ? Devrais je tout réapprendre de zéros ?