Discussion :

[Xx-Astro]

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

Out of memory

Voici le code:

De ma fonction

Code :

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

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 :

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

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 !

[Matersss]

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

[Jerome Briot]

Materss a sans doute raison

Que renvoi ceci avant le tracé ?

Code :

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

whos t Y

[Xx-Astro]

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é

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 :

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

function H=G(t,Y,lb)

Puis dans mon scripte

Code :

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

[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 ?

[Invité]

Bonsoir,

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

[Xx-Astro]

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

Voilà au final:

La fonction:

Code :

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

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 :

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
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 :

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

H=H'

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

Code :

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

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 ?

[Invité]

Parce que

Code :

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

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 :

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

1
2
3
4
5
6
>> H = [1 ; 2]
 
H =
 
     1
     2

La documentation stipule :

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).

[Xx-Astro]

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

Je te remercie !

A bientôt, je pense !

[Jerome Briot]

Ce qui revient simplement à

Code :

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

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);

[Xx-Astro]

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 ?