Bonne soirée à tous! Je cherche à résoudre avec matlab une équation différentielle sous la forme
x"+w^2 x=a*cos(gt)+u(t)
le terme u se calcule avec la fonction
function U=comm()
syms alp a gam T t s w
A=[0 1;-w^2 0]
B=[0;1]
G=[0;1]
a1=expm((T-s)*A)
f=a*exp(-alp*s)*cos(gam*s);
a2=a1*B*B'*a1.';
M=int(a2,'s',0,T);
N1=B'*expm((T-t)*A);
N2=int(expm((T-s)*A)*G*f,'s',0,T);
k=inv(M)
U=-N1*k*N2
end
mais j'ai aucune idée de la façon de faire ça.mercii de m'aider car c'est important pour moi
-----