Comment puis-je trouver la réponse impulsionnelle d'un système à partir de sa représentation espace-état en utilisant la matrice de transition d'état?

15

Supposons que nous ayons un linéaire représenté dans la notation standard de l'espace d'état:

x˙(t)=Ax(t)+Bu(t)
y(t)=Cx(t)+Du(t)

Pour obtenir sa réponse impulsionnelle, il est possible de prendre sa transformée de Laplace pour obtenir

Y = C X + D U

sX=AX+BU
Y=CX+DU

puis résoudre la fonction de transfert qui est

YU=C(sIA)1B+D

De même, pour un système discret, la transformation de Z

x[n+1]=Ax[n]+Bu[n]
y[n]=Cx[n]+Du[n]

est

YU=C(zIA)1B+D

Ce processus semble un peu long et je me souviens qu'il existe un moyen de trouver la réponse impulsionnelle en utilisant la matrice de transition d'état qui est la solution pour des premières équations de chaque paire. Est-ce que quelqu'un sait comment faire ça?x

Phonon
la source

Réponses:

6

Vous pouvez aborder le problème en utilisant la matrice de transition d'état en résolvant l'ODE non homogène standard dans la première équation. La solution à estx˙(t)=Ax(t)+Bu(t)

x(t)=x0eAt+0teA(tt)Bu(t)dt

. La quantité e A t est appelée matrice de transition d'état (également la solution à l'OD homogène), que j'appellerai Ξ ( t ) (je ne me souviens pas de la notation standard pour cela). En prenant x 0 = 0 , l'équation pour y ( t ) devientx0=x(0)eAtΞ(t)x0=0y(t)

y(t)=C0tΞ(tt)Bu(t)dt+Du(t)

Ξ(t)=eAt(sIA)1

Y=C(sIA)1BU+DU

ce qui vous donne la même fonction de transfert que dans votre question.


En ce qui concerne votre commentaire sur l'approche entièrement transformée de Laplace étant long, je ne dirais pas nécessairement qu'il en est ainsi. Cependant, l'approche de la matrice de transition d'état peut être plus simple à mettre en œuvre , car plusieurs opérations la impliquant peuvent être calculées avec de simples multiplications matricielles et rien de plus.

Lorem Ipsum
la source
Très belle description.
Jason R