Message Error Plot

Bonjour,

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

Voici l'erreur que j'obtiens :

Code:

---

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:

---

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;

---

Code:

---

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:

---

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

---

au lieu de:

Code:

---

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

---

?

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:

---

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:

---

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:

---

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:

---

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;

---

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

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.

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.