- Source: Jury stability criterion
In signal processing and control theory, the Jury stability criterion is a method of determining the stability of a discrete-time, linear system by analysis of the coefficients of its characteristic polynomial. It is the discrete time analogue of the Routh–Hurwitz stability criterion. The Jury stability criterion requires that the system poles are located inside the unit circle centered at the origin, while the Routh-Hurwitz stability criterion requires that the poles are in the left half of the complex plane. The Jury criterion is named after Eliahu Ibraham Jury.
Method
If the characteristic polynomial of the system is given by
f
(
z
)
=
a
n
+
a
n
−
1
z
1
+
a
n
−
2
z
2
+
⋯
+
a
1
z
n
−
1
+
a
0
z
n
{\displaystyle f(z)=a_{n}+a_{n-1}z^{1}+a_{n-2}z^{2}+\dots +a_{1}z^{n-1}+a_{0}z^{n}}
then the table is constructed as follows:
row
_
z
n
_
z
n
−
1
_
z
n
−
2
_
⋯
z
1
_
z
0
_
1
a
0
a
1
a
2
⋯
a
n
−
1
a
n
2
a
n
a
n
−
1
a
n
−
2
⋯
a
1
a
0
3
b
0
b
1
⋯
b
n
−
2
b
n
−
1
4
b
n
−
1
b
n
−
2
⋯
b
1
b
0
5
c
0
c
1
⋯
c
n
−
2
6
c
n
−
2
c
n
−
3
⋯
c
0
⋮
⋮
⋮
⋮
⋮
2
n
−
5
p
0
p
1
p
2
p
3
2
n
−
4
p
3
p
2
p
1
p
0
2
n
−
3
q
2
q
1
q
0
{\displaystyle {\begin{array}{lcccccc}{\underline {\text{row}}}&\ {\underline {z^{n}}}\ &\ {\underline {z^{n-1}}}\ &\ {\underline {z^{n-2}}}\ &\ \cdots \ &\ {\underline {z^{1}}}\ &\ {\underline {z^{0}}}\ \\[8pt]1&a_{0}&a_{1}&a_{2}&\cdots &a_{n-1}&a_{n}\\[4pt]2&a_{n}&a_{n-1}&a_{n-2}&\cdots &a_{1}&a_{0}\\[10pt]3&b_{0}&b_{1}&\cdots &b_{n-2}&b_{n-1}\\[4pt]4&b_{n-1}&b_{n-2}&\cdots &b_{1}&b_{0}\\[10pt]5&c_{0}&c_{1}&\cdots &c_{n-2}&\\[4pt]6&c_{n-2}&c_{n-3}&\cdots &c_{0}&&\\[10pt]\ \!\vdots &\vdots &\vdots &\vdots &\vdots &&\\[10pt]2n-5\quad &p_{0}&p_{1}&p_{2}&p_{3}&&\\[4pt]2n-4&p_{3}&p_{2}&p_{1}&p_{0}&&\\[10pt]2n-3&q_{2}&q_{1}&q_{0}&&&\end{array}}}
That is, the first row is constructed of the polynomial coefficients in order, and the second row is the first row in reverse order and conjugated.
The third row of the table is calculated by subtracting
a
n
a
0
{\displaystyle {\tfrac {a_{n}}{a_{0}}}}
times the second row from the first row, and the fourth row is the third row with the first n elements reversed (as the final element is zero).
a
0
a
1
…
a
n
−
1
a
n
a
n
a
n
−
1
…
a
1
a
0
a
0
−
a
n
a
n
a
0
a
1
−
a
n
−
1
a
n
a
0
…
a
n
−
1
−
a
1
a
n
a
0
0
a
n
−
1
−
a
1
a
n
a
0
…
a
1
−
a
n
−
1
a
n
a
0
a
0
−
a
n
a
n
a
0
0
{\displaystyle {\begin{array}{ccccc}a_{0}&a_{1}&\dots &a_{n-1}&\quad a_{n}\quad \\[4pt]a_{n}&a_{n-1}&\dots &a_{1}&a_{0}\\[4pt]a_{0}-a_{n}{\frac {a_{n}}{a_{0}}}&a_{1}-a_{n-1}{\frac {a_{n}}{a_{0}}}&\dots &a_{n-1}-a_{1}{\frac {a_{n}}{a_{0}}}&0\\[4pt]a_{n-1}-a_{1}{\frac {a_{n}}{a_{0}}}&\dots &a_{1}-a_{n-1}{\frac {a_{n}}{a_{0}}}&a_{0}-a_{n}{\frac {a_{n}}{a_{0}}}&0\end{array}}}
The expansion of the table is continued in this manner until a row containing only one non-zero element is reached.
Note the
a
n
a
0
{\displaystyle {\tfrac {a_{n}}{a_{0}}}}
is
a
n
{\displaystyle a_{n}}
for the 1st two rows. Then for 3rd and 4th row the coefficient changes (i.e.
b
n
−
1
b
0
{\displaystyle {\tfrac {b_{n-1}}{b_{0}}}}
). This can be viewed as the new polynomial which has one less degree and then continuing.
If
a
0
>
0
{\displaystyle {a_{0}}>0}
then for every value of
a
0
,
b
0
,
c
0
{\displaystyle a_{0},b_{0},c_{0}}
... that is negative, the polynomial has one root outside of the unit disc. This implies that the method can be stopped after the first negative value is found when checking for stability.
Sample implementation
This method is very easy to implement using dynamic arrays on a computer. It also tells whether all the modulus of the roots (complex and real) lie inside the unit disc. The vector v contains the real coefficients of the original polynomial in the order from highest degree to lowest degree.
See also
Liénard–Chipart criterion, another stability criterion derived from Routh-Hurwitz (for continuous-time systems)
References
For more details please check these references:
A note on the reduced Schur–Cohn criterion
Wikibooks on Control Systems - Jury's Test
For advanced resources:
Archived August 2, 2008, at the Wayback Machine
Benidir, M. (1996). "On the root distribution of general polynomials with respect to the unit circle". Signal Processing. 53 (1): 75–82. Bibcode:1996SigPr..53...75B. doi:10.1016/0165-1684(96)00077-1.
http://www.laas.fr/~henrion/Papers/lyap.ps.gz
For implementations:
http://www.ticalc.org/archives/files/fileinfo/426/42696.html (TI-83+/84+ graphing calculators)
Kata Kunci Pencarian:
7
6.6
7.7
6.9
6.535
5.5
6.1
7.136
7.433
No More Posts Available.
No more pages to load.
