Podręcznik
1. Metoda równań różniczkowych w analizie stanów nieustalonych w obwodach
1.8. Wartości własne i wektory własne macierzy kwadratowej
Załóżmy, że A jest macierzą kwadratową stopnia n. Macierz (s1-A) nazywana jest macierzą charakterystyczną A, przy czym 1 oznacza macierz jednostkową stopnia n, to jest macierz diagonalną 1=diag(1, 1,..., 1). Wyznacznik macierzy charakterystycznej det(s1-A) nazywamy wielomianem charakterystycznym macierzy, a równanie
|
det(s1-A) = 0 |
(1.18) |
nazywamy równaniem charakterystycznym macierzy A. Równanie to po rozwinięciu wyrażenia wyznacznika przyjmuje postać wielomianu n-tego stopnia
| (1.19) |
Pierwiastki tego równania s1, s2, ..., sn nazywamy wartościami własnymi macierzy A. Mogą one przyjmować wartości rzeczywiste lub zespolone, pojedyncze lub wielokrotne. Z każdą wartością własną si skojarzony jest wektor własny xi o niezerowej wartości i wymiarze n, spełniający równanie
| (1.20) |
Jeśli wszystkie wartości własne są różne to na podstawie równania (1.20) można napisać n równań liniowych o postaci
|
................... |
(1.21) |
z rozwiązania których można wyznaczyć wszystkie wektory własne xi. W Matlabie wartości i wektory własne sa wyznaczane przy użyciu funkcji eig.
Dla macierzy stanu
wyznaczyć wartości i wektory własne
Rozwiązanie
Równanie charakterystyczne
Pierwiastki tego równania będące wartościami własnymi A są równe s1=-4 oraz s2=-1. Wektory własne spełniają relację (1.25), która w naszym przypadku przyjmie postać
Powyższym równaniom odpowiadają cztery równania skalarne o postaci
Biorąc pod uwagę, że dwa spośród nich są zależne, dwie zmienne można przyjąć dowolnie, na przykład x11=1 oraz x22=-1. Z rozwiązania pozostałych 2 równań otrzymuje się wektory własne o postaci
Napisać układ równań stanu dla obwodu elektrycznego przedstawionego na rys. 1.1
Rys. 1.1. Schemat obwodu do przykładu 1.2
Rozwiązanie
Z praw Kirchhoffa napisanych dla obwodu z rys. 1.1 wynikają następujące równania
Biorąc pod uwagę, że
oraz
równania Kirchhoffa można przekształcić do równoważnej postaci równań różniczkowych
które przyjmują uporządkowaną formę odpowiadającą postaci (1.5)
Równania powyższe można zapisać w postaci zależności macierzowej równania stanu, w której zmiennymi stanu są prąd cewki i napięcie kondensatora.
Wektor stanu x jest równy
a wektor wymuszeń
Obwód liniowy zawierający dwa elementy reaktancyjne (cewka i kondensator) opisuje się więc macierzowym równaniem stanu drugiego rzędu. Macierz stanu A jest macierzą również drugiego rzędu o współczynnikach uzależnionych od wartości rezystancji, pojemności oraz indukcyjności. Macierz B zawiera dwa wiersze (liczba zmiennych stanu) oraz dwie kolumny (liczba wymuszeń w obwodzie). Przyjmując w analizie wartości liczbowe obwodu: R=2W, L=1H, C=1F otrzymuje się macierz stanu A o postaci
Równanie charakterystyczne tej macierzy jest równe
Wartości własne (pierwiastki równania charakterystycznego) są w tym przypadku sobie równe i wynoszą s1 = s2 = -1. Dla rozważanego obwodu RLC są one położone w lewej półpłaszczyźnie zmiennej zespolonej s na ujemnej osi rzeczywistej.





![\mathbf{A}=\left[\begin{matrix}-2&-2\\-1&-3\\\end{matrix}\right] \mathbf{A}=\left[\begin{matrix}-2&-2\\-1&-3\\\end{matrix}\right]](https://esezam.okno.pw.edu.pl/filter/tex/pix.php/d46f91e8c8b205838fa0c14c85c4c657.gif)
![det{(}s\mathbf{1}-\mathbf{A})=det{\left(s\left[\begin{matrix}1&0\\0&1\\\end{matrix}\right]-\left[\begin{matrix}-2&-2\\-1&-3\\\end{matrix}\right]\right)}=s^2+5s+4=0 det{(}s\mathbf{1}-\mathbf{A})=det{\left(s\left[\begin{matrix}1&0\\0&1\\\end{matrix}\right]-\left[\begin{matrix}-2&-2\\-1&-3\\\end{matrix}\right]\right)}=s^2+5s+4=0](https://esezam.okno.pw.edu.pl/filter/tex/pix.php/d07bdc6e52786b0888ed20961f680de6.gif)
![\left[\begin{matrix}-2&-2\\-1&-3\\\end{matrix}\right]\left[\begin{matrix}x_{11}\\x_{21}\\\end{matrix}\right]=-4\left[\begin{matrix}x_{11}\\x_{21}\\\end{matrix}\right] \left[\begin{matrix}-2&-2\\-1&-3\\\end{matrix}\right]\left[\begin{matrix}x_{11}\\x_{21}\\\end{matrix}\right]=-4\left[\begin{matrix}x_{11}\\x_{21}\\\end{matrix}\right]](https://esezam.okno.pw.edu.pl/filter/tex/pix.php/98da58619946ea857c1005bef1b98001.gif)
![\left[\begin{matrix}-2&-2\\-1&-3\\\end{matrix}\right]\left[\begin{matrix}x_{12}\\x_{22}\\\end{matrix}\right]=-1\left[\begin{matrix}x_{12}\\x_{22}\\\end{matrix}\right] \left[\begin{matrix}-2&-2\\-1&-3\\\end{matrix}\right]\left[\begin{matrix}x_{12}\\x_{22}\\\end{matrix}\right]=-1\left[\begin{matrix}x_{12}\\x_{22}\\\end{matrix}\right]](https://esezam.okno.pw.edu.pl/filter/tex/pix.php/c81dc8f12a0f80bda575a1e51e513543.gif)




![\mathbf{x}_1=\left[\begin{matrix}1\\1\\\end{matrix}\right],\mathrm{\ \ \ \ \ \ }\mathbf{x}_2=\left[\begin{matrix}2\\-1\\\end{matrix}\right] \mathbf{x}_1=\left[\begin{matrix}1\\1\\\end{matrix}\right],\mathrm{\ \ \ \ \ \ }\mathbf{x}_2=\left[\begin{matrix}2\\-1\\\end{matrix}\right]](https://esezam.okno.pw.edu.pl/filter/tex/pix.php/231f9723294c0a77728f203942e5e532.gif)








![\left[\begin{matrix}\frac{di_L}{dt}\\\frac{du_C}{dt}\\\end{matrix}\right]=\left[\begin{matrix}\frac{-R}{L}&\frac{-1}{L}\\\frac{1}{C}&0\\\end{matrix}\right]\left[\begin{matrix}i_L\\u_C\\\end{matrix}\right]+\left[\begin{matrix}\frac{1}{L}&\frac{R}{L}\\0&\frac{-1}{C}\\\end{matrix}\right]\left[\begin{matrix}e\\i\\\end{matrix}\right] \left[\begin{matrix}\frac{di_L}{dt}\\\frac{du_C}{dt}\\\end{matrix}\right]=\left[\begin{matrix}\frac{-R}{L}&\frac{-1}{L}\\\frac{1}{C}&0\\\end{matrix}\right]\left[\begin{matrix}i_L\\u_C\\\end{matrix}\right]+\left[\begin{matrix}\frac{1}{L}&\frac{R}{L}\\0&\frac{-1}{C}\\\end{matrix}\right]\left[\begin{matrix}e\\i\\\end{matrix}\right]](https://esezam.okno.pw.edu.pl/filter/tex/pix.php/fb3855516ee5b7e6948c420c55ca4acf.gif)
![{x}=\left[\begin{matrix}i_L\\u_C\\\end{matrix}\right] {x}=\left[\begin{matrix}i_L\\u_C\\\end{matrix}\right]](https://esezam.okno.pw.edu.pl/filter/tex/pix.php/956799a2d420b84bc8f1b0d6d561425a.gif)
![{u}=\left[\begin{matrix}{e}\\{i}\\\end{matrix}\right] {u}=\left[\begin{matrix}{e}\\{i}\\\end{matrix}\right]](https://esezam.okno.pw.edu.pl/filter/tex/pix.php/a2c54486997109a8e09b846d597839b0.gif)
![{A}=\left[\begin{matrix}-2&-1\\1&0\\\end{matrix}\right] {A}=\left[\begin{matrix}-2&-1\\1&0\\\end{matrix}\right]](https://esezam.okno.pw.edu.pl/filter/tex/pix.php/0cfbf4ba42c61c3c826679610dd3c15e.gif)
![det{(}s\mathbf{1}-{A})=det{\left(\left[\begin{matrix}s&0\\0&s\\\end{matrix}\right]-\left[\begin{matrix}-2&-1\\1&0\\\end{matrix}\right]\right)}=s^2+2s+1=0 det{(}s\mathbf{1}-{A})=det{\left(\left[\begin{matrix}s&0\\0&s\\\end{matrix}\right]-\left[\begin{matrix}-2&-1\\1&0\\\end{matrix}\right]\right)}=s^2+2s+1=0](https://esezam.okno.pw.edu.pl/filter/tex/pix.php/1980f239138794cc029e468ad7496c3c.gif)