Programme maple d'une representation d'état discrete

[invite00c17237]

Bonjour,

Le calcul de représentation d'état discrète à partir d'une représentation d'état continue se base sur le calcul d'une exponentielle de matrice.

Je possède A et je dois calcule F = exp(A*T) où T une constante d'échantillonnage.

Voici le programme :





J'ai quelques bugs et je ne comprends pas d'où çà provient. A est diagonalisable. Aussi, pourquoi est ce que je ne peux pas déterminer son inverse. Est-ce que mon souvenir comme quoi tout matrice diagonalisable est inversible est correct ? Et si oui, quel est le problème au niveau du programme.

Je vous remercie pour votre aide

Bien cordialement,

Benjamin BOUDON
-----

[invite551c2897]

Bonjour.
Plutôt :
evalm(M &* exp(DD*T)&*inverse(M));

[invite551c2897]

Bonjour.
Peux-être : :
DD:=diag(eigenvalues(A));de:=e xponential(DD);df:=evalm(de*T) ;
> dg:=evalm(M &* df);dh:=evalm( dg&* inverse(M));
with(exponential) n'a pas l'air d'ajouter quelque chose !