Bistritz stability criterion


In signal processing and control theory, the Bistritz criterion is a simple method to determine whether a discrete linear time invariant system is stable proposed by Yuval Bistritz. Stability of a discrete LTI system requires that its characteristic polynomials
is a stable polynomial, where is said to be stable if all its zeros are inside the unit circle, viz.
where . The test determines whether is stable algebraically . The method also solves the full zero location problem. Namely, it can count the number of inside the unit-circle zeros , on the unit-circle zeros zeros and outside the unit-circle zeros for any real or complex polynomial.
The Bistritz test is the discrete equivalent of Routh criterion used to test stability of continuous LTI systems. This title was introduced soon after its presentation. It has been also recognized to be more efficient than previously available stability tests for discrete systems like the Schur–Cohn and the Jury test.
In the following, the focus is only on how to test stability of a real polynomial. However, as long as the basic recursion needed to test stability remains valid, ZL rules are also brought.

Algorithm

Consider as above and assume . Define its reciprocal polynomial
The algorithm assigns to a sequence of symmetric polynomials
created by a three-term polynomial recursion. Write out the polynomials by their coefficients,
symmetry means that
so that it is enough to calculate for each polynomial only about half of the coefficients. The recursion begins with two initial polynomials driven from the sum and difference of the tested polynomial and its reciprocal, then each subsequent polynomial of reduced degree is produced from the last two known polynomials.
Initiation:
Recursion: For do:

Stability condition

The successful completion of the sequence with the above recursion requires
. The expansion of these conditions into
are called normal conditions.
Normal conditions are necessary for stability. This means that, the tested polynomial can be declared as not stable as soon as a is observed. It also follows that the above recursion is broad enough for testing stability because the polynomial can be declared as not stable before a division by zero is encountered.
Theorem. If the sequence is not normal then is not stable.
If normal conditions hold then the complete sequence of symmetric polynomials is well defined. Let
denote the count of the number of sign variations in the indicated sequence. Then
is stable if and only if.
More generally, if normal condition hold then has no UC zeros,
OUC zeros and IUC zeros.
Violation of various necessary conditions for stability may be used advantageously as early indications that the polynomial is not stable. The polynomial can be declared not stable as soon as a , or a, or a change of sign in the sequence of 's is observed.

Example

Consider the polynomial, where is a real parameter.
Q1:For what values of the polynomial is stable?
Construct the sequence:
Use their values at z=1 to form
All the entries in the sequence are positive for -4 < K < 22 . Therefore D is stable for −4 < K < 22.
Q2: Find ZL for K = 33 Var =2 => 2 OUC, 1 IUC zeros.
Q3: Find ZL for K = -11 Var=1 => 1 OUC, 2 IUC zeros.

Comments

The test bears a remarkable similarity to the Routh test. This is best observed when the Routh test is arranged appropriately into a corresponding three-term polynomial recursion.
The Bistritz test uses three-term polynomial recursion that propagates polynomials with symmetry as opposed to previously available classical tests for discrete systems that propagate polynomials with no particular structure using a two-term recursion. It stimulated the discovery of more algorithms in the area of digital signal processing and discrete systems collectively called "immittance" or "split" algorithms that adopted this technique to more efficient counterparts to also other classical so called "scattering" algorithms. The Bistritz test forms the "immittance" counterpart of the "scattering" type classical tests of Schur–Cohn and Jury.