Résoudre un ODE Matlab

Question

je l'ODE suivante:

b'(t) + k16*b(t) = k15*a(t) 

où K15 et K16 sont des constantes.

Une idée sur la façon de le résoudre?

Merci! Amit

Answer 1

C'est une ODE de premier ordre. Il y a une solution analytique pour cela (il suffit d'utiliser un facteur d'intégration). Aucune intégration requise. http://www.math.hmc.edu/calculus/tutorials/odes/

Cependant, si vous voulez résoudre dans Matlab:

>> k15 = 0.2; k16 = 0.3; % type your constants here 
>> a = @(t) t^2; % type your expression for a here 
>> dbdt = @(t,b) -k16*b + k15*a(t); 
>> tf = 10; % final time of integration 
>> b0 = 1; % initial value of b 
>> [t,y] = ode45(@dbdt,[0 tf],b0) 
>> plot(t,y) % display solution.