Durum geçiş matrisini kullanarak bir sistemin durum-uzay geri dönüşünden gelen dürtü yanıtını nasıl bulabilirim?


15

Standart durum alanı gösteriminde temsil edilen bir lineerimiz olduğunu varsayalım:

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

Dürtü tepkisini almak için Laplace dönüşümünü almak mümkündür.

sX=AX+BU
Y=CX+DU

ve sonra transfer fonksiyonu için çözmek

YU=C(sIA)1B+D

Benzer şekilde, ayrık bir sistem için, dönüşümü Z

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

dır-dir

YU=C(zIA)1B+D

Bu süreç biraz uzun görünüyor ve her bir çiftin ilk denklemlerinin çözümü olan durum geçiş matrisini kullanarak dürtü yanıtını bulmanın bir yolu olduğunu hatırlıyorum . Bunu nasıl yapacağını bilen var mı?x

Yanıtlar:


6

İlk denklemde standart homojen olmayan ODE'yi çözerek duruma geçiş matrisini kullanarak soruna yaklaşabilirsiniz. Çözüm olanx˙(t)=Ax(t)+Bu(t)

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

burada . miktarı , olarak anlatacağım durum geçiş matrisi (ayrıca homojen ODE'nin çözümü) olarak adlandırılır (bunun için standart gösterimi hatırlamıyorum). Alarak , denklemi olurx0=x(0)eAtΞ(t)x0=0y(t)

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

Yukarıdaki denklem, sistem impuls tepkisi ile kıvrık giriş olarak çıkışı verir ve gerçekten, doğrulamak için yukarıdaki denklemin Laplace dönüşümünü alabilirsiniz. Laplace transformu belirten olan s-etki zaman alanı haline ürünlerde ve kıvrım bulunur, elde ederizΞ(t)=eAt(sIA)1

Y=C(sIA)1BU+DU

bu da sorunuzdakiyle aynı aktarım işlevini sağlar.


Tam Laplace dönüşümü yaklaşımının uzun olması hakkındaki yorumunuzla ilgili olarak, bunun böyle olduğunu söylememeliyim. Bununla birlikte, durum geçiş matrisi yaklaşımının uygulanması daha basit olabilir , çünkü onu içeren çeşitli işlemler basit matris çarpımları ve daha fazlası ile hesaplanabilir.


Çok güzel bir açıklama.
Jason R
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.