Résoudre une équation différentiel avec Matlab

[invited1c96ae7]

bonjour,

est-ce que vous pouvez m'aider de résoudre une équation différentiel?


t=0:0.1:10

x est un fonction de t
diff(x,2)=diff(x)-5x

je veux avoir x(t) et x(0)=1

Merci d'avance
-----

[invite551c2897]

Bonjour.
PP :

clear
[t,x]=ode45('exeqd',[0:0.1:10],[1 0]);
subplot(2,1,1)
plot(t,x(:,1))
grid
subplot(2,1,2)
plot(t,x(:,2))
grid

SP :

function dx=exeqd(t,x)
dx(1)=x(2);
dx(2)=x(2)-5*x(1);
dx=[dx(1);dx(2)];

[invited1c96ae7]

Bonjour.
PP :

clear
[t,x]=ode45('exeqd',[0:0.1:10],[1 0]);
subplot(2,1,1)
plot(t,x(:,1))
grid
subplot(2,1,2)
plot(t,x(:,2))
grid

SP :

function dx=exeqd(t,x)
dx(1)=x(2);
dx(2)=x(2)-5*x(1);
dx=[dx(1);dx(2)];

Merci beaucoup
Est-ce que c'est possible qu'il renvoie les valeur de x1(t) et x2(t)?

[invite551c2897]

x1(t)=x(:,1)
et x2(t)=x(:,2)

[A voir en vidéo sur Futura]

[invited1c96ae7]

Merci beaucoup
Est-ce que c'est possible qu'il renvoie les valeur de x1(t) et x2(t)?

J'ai trouvé merci

[invited1c96ae7]

x1(t)=x(:,1)
et x2(t)=x(:,2)

Comment je peux afficher x(:,1)^2 ? (comme ça ne marche pas :s)

[invite551c2897]

Il faut mettre un point :
x(:,1).^2

[invited1c96ae7]

Merci beaucoup

Est-ce que quelqu'un peut dire comment je peux faire un animation avec 2 vecteur (afficher les point (x,y) un après l'autre ) ?

[invite551c2897]

Bonjour.
Un exemple d'animation :

clear
N=10;pas=0.5;
h=plot(0,0,'.r');
grid
set(gcf,'DoubleBuffer','on')
set(h,'Markersize',23);
axis([0 20 0 20])
nbf=0;
 for t=0:pas:N
    nbf=nbf+1;
     y=2*t;
    set(h,'XData',t,'YData',y)
   M(nbf)=getframe;
   end
movie(M,1)
save M

[invited1c96ae7]

il affiche 2 fois là meme chose
je peux le fafficher juste une fois?

[invite551c2897]

movie(M,0)