bonjour, voici le programme pour :
Écrire une fonction permettant de résoudre une équation différentielle ordinaire du second ordre : a¨ x + b˙ x + cx = cos(t), x(0) = X0, ˙ x(0) = V0. 3.1. Déterminer les racines de l’équation caractéristique. 3.2. Déterminer la solution de l’équation homogène. 3.3. Déterminer la solution particulière. 3.4. Tracer l’évolution de la solution générale en fonction du temps.
(a,b,c) sont remplacé par m,c,k
je ne comprends pas a quoi sert la matrice B
function [x] = TP3_EX3(m,c,k,X0,V0,pas,tf)
t=0as:tf;
Y=[X0;V0];
F=[1;0];
if m==0
disp("Ordre de l'équation <2");
else
[r1,r2] = TP3_EX1(m,b,c);
if r1~=r2
A=[1,1;r1,r2];
cste=Y\A;
xh=cste(1)*exp(r1*t)+cste(2)*e xp(r2*t);
else
A=[1,0;r1,1];
cste=Y\A;
xh=cste(1)*exp(r1*t)+cste(2)*t *exp(r1*t);
end
end
figure(1)
plot(t,xh)
B=[k-m,c;-c,k-m];
ALPHA=F\B;
figure(2)
xp=ALPHA(1)*cos(t)+ALPHA(2)*si n(t);
x=xp+xh;
plot(t,xp)
figure(3)
plot(t,x)
end
-----