Discussion :

[B_E_L]

Bonjour,

J' ai un petit soucis avec la realisation de l'affichage en graphique de mes codes.

Voici l'erreur que j'obtiens :

Code :

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

1
2
3
4
5
6
7
8
??? Error using ==> plot
Vectors must be the same lengths.
 
Error in ==> C:\Documents and Settings\ODEexp2.m
On line 10  ==> plot(x,y_time);
 
??? Error using ==> edit
Can't edit the built-in function 'plot'.

>> En fait mon but est de pouvoir voir si il existe un difference dans la resolution d' un ODE avec ou sans la fonction deval
et de le representer par graphe... Merci d'avance pour votre aide

je vous poste un cas simple d equation que j ai essayer de resoudre:

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
22
23

function ODE

tau = [1 3]
sol = ode45(@ODEexp1,tau ,4.2);


x = linspace(1,3,15);
y_time= deval(sol,tau,1);

plot(x,y_time);



plot(t,y);
xlabel('t');
ylabel('y');

function ydot = ODEexp1(t,y)

C_H2=y(1,1)

ydot=(t^3-2*y)/t;

[Uren]

Code :

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

1
2
3
4
5
6
 
??? Error using ==> plot
Vectors must be the same lengths.
 
Error in ==> C:\Documents and Settings\ODEexp2.m
On line 10 ==> plot(x,y_time);

Ce message d'erreur veut dire que les vecteurs x et y_time ne sont pas de même longueur.

En regardant ton code de plus près, si je fais pas d'erreur, je trouve que x est un vecteur de longueur 15 alors que y_time est un vecteur de la meme longueur que le vecteur "tau", donc 2.

Forcement tu ne peux pas dessiner y_time sur 15 points s'il n'en contient que 2, d'ou le message d'erreur.

Je ne suis pas sur de comprendre totalement ce que tu veux faire mais tu ne devrais pas plutot ecrire:

Code :

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

1
2
 
y_time= deval(sol,x,1);

au lieu de:

Code :

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

1
2
 
y_time= deval(sol,tau,1);

?

[B_E_L]

oki Merci pour l'erreur je l'ai rectifier.

En fait je voulais essayer de comprendre la difference qu il existe pour resoudre un ODE, en utilisant seulement la fonction ODE45 ou en utilisant un ODE +un deval. Je poste les deux exercices que j ai fait mais quand j essaye d obtenir 2 figure pour le 3eme code> ca ne fonctionne plus. L' erreur que j obtiens est pour le code 3 :

Code :

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

1
2
3
4
??? Undefined function or variable 't'.
 
Error in ==> C:\Documents and Settings\ODEexp2.m
On line 17  ==> plot(t,y)

CODE 1

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
22
23
 
 
function ODE
 
y_initial =4.2
tau = [1 3]
sol = ode45(@ODEexp1, tau , y_initial );
 
 
x = linspace(1,3,100);
y_time= deval(sol,x,1);
 
figure;
plot(x,y_time);
xlabel('t')
ylabel('y_time')
 
 
function ydot = ODEexp1(t,y)
 
C_H2=y(1,1);
 
ydot=(t^3-2*y)/t;

CODE 2

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
function ODE
 
[t,y] = ode45(@ODEexp1,[1 3],4.2);
 
figure;
plot(t,y)
xlabel('t')
ylabel('y')
 
function ydot = ODEexp1(t,y)
 
C_H2=y(1,1);
 
ydot=(t^3-2*y)/t;

CODE 3

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
22
23
24
25
26
27
function ODE
 
y_initial =4.2
tau = [1 3]
sol = ode45(@ODEexp1, tau , y_initial );
 
 
x = linspace(1,3,100);
y_time= deval(sol,x,1);
 
figure;
plot(x,y_time);
xlabel('t')
ylabel('y_time')
 
figure;
plot(t,y)
xlabel('t')
ylabel('y')
 
%figure ; bar3(y_time-y_initial),ylabel('i'),xlabel('j'),zlabel('\delta y^H^2,m.f.'), title (' Variation of mole fr. in anode chanels, (y^H^2-y^H^2_F))TSR model ');
 
function ydot = ODEexp1(t,y)
 
C_H2=y(1,1);
 
ydot=(t^3-2*y)/t;

[Jerome Briot]

Les variables t et y ne sont pas définies dans la fonction principale du 3ème code

[Merel]

Salut,

Quand tu définis une fonction, les variables sont internes à la fonction. Même si tu définis tes vecteurs t et y dans ton code 2, tu ne les as pas définis dans ton code 3. D'où le message d'erreur.

[B_E_L]

Mais en fait y est l inconnue mais alors comment est ce que je peux la définir?

dy/dt = (t^3 -2y)/t for t = [1,3] avec y= 4.2 a t =1

et donc je résous cette '' Ordinary Differential Eqaution '' avec ode 45 (code 2) et je veux comparer si ce que j obtiens avec la fonction deval (code 3) est pareil ou meilleure.

Je n utilise pas souvent Matlab merci pour l aide.