Exponentiel de matrice: solution d'un systeme
Réponses à toutes vos questions après le Bac (Fac, Prépa, etc.)
-
fidi
- Messages: 3
- Enregistré le: 24 Juin 2008, 19:01
-
par fidi » 24 Juin 2008, 19:39
Bonsoir à tous,
Voilà, je me retrouve devant un "petit" problème que j'avoue avoir du mal à résoudre...
A partir de cette équation fonctionnelle:
M(x+y) = M(x).M(y)
avec M(x) la fonction de transfert d'un tuyau (matrice 2x2).
Si A, B et C désigne 3 sections du tuyau séparé par les distances x et y, on peut lier les vecteurs perturbation P (pression et débit) par:
Pb=M(x).Pa ; Pc=M(y).Pb et Pc = M(x+y).Pa
Comment peut on démontrer que les matrices M(x) solution de l'équation fonctionnelle sont de la forme:
M(x) = exp (x.K)
avec K une matrice (2.2) à coefficients complexes ne dépendant pas de x
Par avance merci
-
Aspx
- Membre Relatif
- Messages: 285
- Enregistré le: 25 Avr 2006, 13:12
-
par Aspx » 24 Juin 2008, 20:01
Dérivabilité automatique :La multiplication à gauche par
)
étant linéaire et continue elle commute avec l'intégration d'où en intégrant ta relation
[CENTER]
dy = M(x) \int_{0}^{t} M(y)dy)
[/CENTER]
i.e
[CENTER]
dy = M(x) \int_{0}^{t} M(y)dy)
[/CENTER]
Il reste plus qu'à trouver un

tel que
dy)
soit inversible (on aura alors le caractère

de

). Or (j'imagine que les conditions initiales doivent imposer
=I_2)
) vu

continue
 = I_2 + \epsilon(y))
puis en primitivant le DL
dy = t(I_2+\epsilon(t)))
qui est inversible pour

assez petit (
)
ouvert)
On peut dériver (d/dy) :[CENTER]
=M(x)M'(y))
[/CENTER]
Puis

[CENTER]
=M(x)M'(0))
[/CENTER]
Fonction auxiliaire :Posons
[CENTER]
 = M(x) \exp (-x M'(0)))
[/CENTER]
Sa dérivée est nulle, puis la conclusion vient ensuite.
-
fidi
- Messages: 3
- Enregistré le: 24 Juin 2008, 19:01
-
par fidi » 24 Juin 2008, 20:44
Bonsoir, merci de ta réponse !! Je suis vraiment impressionné par ta rapidité !!! merci. Cependant une chose m'échappe encore comment détermines tu la fonction auxiliaire et comment conclure?
Merci de ces éclairsicements
-
Aspx
- Membre Relatif
- Messages: 285
- Enregistré le: 25 Avr 2006, 13:12
-
par Aspx » 24 Juin 2008, 20:50
Eh bien
[CENTER]
=\exp (-x M'(0)) (M'(x)-M(x)M'(0))=0)
[/CENTER]
Donc

est constant donc
=Q(0)=M(0)=I_2)
Finalement en multipliant à droite par
))
il vient
[CENTER]
 = \exp (x M'(0)))
[/CENTER]
-
fidi
- Messages: 3
- Enregistré le: 24 Juin 2008, 19:01
-
par fidi » 24 Juin 2008, 21:56
Ok, désolé, je n'avais pas vu qu'il fallait dérivé Q...(mal lu :dodo: )
Je vais me permettre d'abuser un peu:
Je voudrais comprendre comment on arrive à l'expression suivante
Si on pose M'(0) = K

Encore merci pour la réactivité.
Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 17 invités