Macierz przejścia (automatyka)

Z testwiki
Przejdź do nawigacji Przejdź do wyszukiwania

Macierz przejścia stanu (lub krótko macierz przejścia), macierz tranzycji, macierz transformacji, macierz fundamentalna, macierz podstawowa (Szablon:Ang.) – macierz, której iloczyn z wektorem stanu x z chwili początkowej t0 daje stan x w późniejszej chwili t. Macierz przejścia stanu może być wykorzystana do uzyskania ogólnego rozwiązania dla liniowych układów dynamicznych. Macierz ta znana jest też jako eksponenta macierzy.

Rozwiązanie równań stanu

Niech dany będzie ogólny liniowy model przestrzeni stanów w postaci równań stanu:

𝐱˙(t)=𝐀(t)𝐱(t)+𝐁(t)𝐮(t),
𝐲(t)=𝐂(t)𝐱(t)+𝐃(t)𝐮(t).

Rozwiązanie ogólne dane jest wówczas równaniem (jest to tak zwany wzór Cauchy-Bellmana):

𝐱(t)=Φ(t,t0)𝐱(t0)+t0tΦ(t,τ)𝐁(τ)𝐮(τ)dτ,

gdzie Φ(t,τ) jest macierzą przejścia określoną poniżej.

Innymi słowy: stan układu przedstawiany jest zwykle jako wektor 𝐱=[x1,x2,,xn]Rn i przedstawia pamięć układu. Znając stan układu oraz sterowanie jesteśmy w stanie określić stan, który osiągnie układ po zadanym czasie.

Dla układu regulacji opisanego układem równań różniczkowych przyjmuje on postać:

x(t)=etAx(0)+0te(tr)ABu(r)dr,

gdzie etAx(0) nazywana jest składową swobodną (zależną od warunków początkowych), a 0te(tr)ABu(r)dr składową wymuszoną (która jest splotem odpowiedzi impulsowej i wejścia). W przypadku układu swobodnego postać rozwiązania sprowadza się do składowej swobodnej (tzw. rozwiązanie swobodne).

Wyprowadzenie wzoru dla układu jednowymiarowego

Wzór na stan x(t) układu jednowymiarowego, opisanego równaniami stanu:

dxdt=ax+bu(t),
y=cx,

gdzie u(t) to zadane sterowanie.

Wyznacza się go w dwóch krokach:

  1. Obliczane jest rozwiązanie bez części sterującej
    dxdt=ax.
    Przekształca się powyższy wzór tak, aby po jednej stronie znalazło się dx oraz x, a po drugiej stronie dt
    dxx=adt.
    Uzyskany wzór całkuje się obustronnie uzyskując:
    lnxS=at, gdzie S to stała całkowania.
    Na koniec następuje pozbycie się logarytmu naturalnego używając eksponenty dla obydwu stron równania:
    x(t)=Seat.
  2. Uzyskany x(t) podstawia się do równań podanych na wstępie i oblicza pochodną x po czasie.
    dxdt=aeatS+eatdSdt=ax+bu(t),
    dSdt=bu(t)eat.
    Przenosi się dt na prawą stronę i całkuje obustronnie:
    S(t)=S0+0terabu(r)dr,
    x(0)=S0.
    Na koniec wstawia się uzyskane S(t) do wzoru x(t)=Seat.
    x(t)=etax(0)+0te(tr)abu(r)dr.

Macierz przejścia

Macierz przejścia Φ(t,τ) określona jest jako:

Φ(t,τ)𝐔(t)𝐔1(τ),

gdzie 𝐔(t) jest podstawową macierzą rozwiązania, która spełnia zależność:

𝐔˙(t)=𝐀(t)𝐔(t)

jest macierzą o wymiarach n×n, która stanowi liniowe mapowanie na siebie samą, na przykład z 𝐮(t)=0, przy danym stanie 𝐱(τ) w dowolnej chwili czasu τ, stan w dowolnej innej chwili t określony jest przez mapowanie:

𝐱(t)=Φ(t,τ)𝐱(τ).

Podczas gdy macierz przejścia stanu ϕ nie jest całkowicie nieznana, to zawsze musi spełniać następujący związek:

ϕ(t,t0)t=A(t)ϕ(t,t0) i
ϕ(τ,τ)=I dla każdego τ i gdzie I jest macierzą jednostkową.

Ponadto ϕ musi posiadać następujące właściwości:

1. ϕ(t2,t1)ϕ(t1,t0)=ϕ(t2,t0)
2. ϕ1(t,τ)=ϕ(τ,t)
3. ϕ1(t,τ)ϕ(t,τ)=I
4. dϕ(t,t0)dt=A(t)ϕ(t,t0)

Jeśli układ jest niestacjonarny, można zdefiniować ϕ jako:

ϕ(t,t0)=eA(tt0).

W przypadku niestacjonarnym, istnieje wiele różnych funkcji, które spełniają te wymagania, a rozwiązanie uzależnione jest od struktury układu. Macierz przejścia stanu musi zostać określona przed dalszą analizą rozwiązania dla układu niestacjonarnego.

Zobacz też