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;