Discussion :

[Mathieu37000]

Bonjour,

J'ai un projet de fin d'année, et ce dernier me demande d'utiliser Matlab afin de réaliser un programme.

Je n'ai pas beaucoup de connaissances et suis un peu bloqué, je ne sais pas trop par ou commencer. Le problème est le suivant :

Je dois calculer le système différentiel suivant :


dr/dt = 2*r-a*r*f

df/dt = -f+a*r*f

a = 0,01 sur l'intervalle [0,25] ; r0 = 100 ; f0 = 400

On me demande d'utiliser la fonction ode23.


Avez-vous une idée pour m'aider ?

Merci d'avance,

Bien cordialement,

Mathieu

[le fab]

bonjour

tu as regardé l'aide de ode23 ?
qu'as tu commencé à coder ?

[Mathieu37000]

Bonjour,

Pour l'instant j'ai réalisé ce code, qui a l'air de fonctionner ? :

Code :

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

1
2
3
4
5
6
7
8
9
t0 =0 ;
tfinal=25 ;
y0=[100,400] ;
[t,y]=ode23(@springb,[t0,tfinal],y0)
plot(t,y);
title('Lapin/Renard Population Over Time')
xlabel('t')*
ylabel('Population');
legend('lapin','renard','Location','North');

Fichier lokta *:

Code :

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

1
2
3
4
5
6
7
function yp = springb(t,y);
alpha = 0.01;
yp = zeros(size(y));
*
yp(1) =(2-alpha*y(2))*y(1);
yp(2) =(-1+alpha*y(1))*y(2);
end

Cependant, j'ai un deuxième code ou j'ai un peu plus de mal :

[Mathieu37000]

x_z_vac=dsolve('D2x=0', 'D2z=-g', 'x(0)=0', 'Dx(0)=v0*cos(theta*pi/180)', 'z(0)=h', 'Dz(0‡=v0*sin(theta*pi/180)', 't');
xt=x_z_vac.x
zt=x_z_vac.z
g=9.8;
h=1;
theta=15;
v0=25;
xt=vectorioze(xt);
zt=(vectorize(zt));
t=linspace(0,1.5,200);
xt=eval(xt);
zt=eval(zt);
z_i=find(zt==min(abs(zt)));
touch_gr=xt(z_i);
t_touch=t(z_i);
plot(xt, zt), grid on xlim([0, touch, gr])†;


Pour le deuxième , mais je suis un peu bloqué