1 Stability Analysis of closed loop system in z-plane
The stability of the following closed loop system
can be determined from the location of closed loop poles in z-plane which are the roots of the characteristic equation
1. For the system to be stable, the closed loop poles or the roots of the characteristic equation must lie within the unit circle in z-plane. Otherwise the system would be unstable.
2. If a simple pole lies at , the system becomes marginally stable. Similarly if a pair of complex conjugate poles lie on the marginally stable. Multiple poles at the same location on unit circle make the system unstable.
circle, the system is
Example 1: Determine the closed loop stability of the system shown in Figure 1 when K = 1.
Figure 1: Example 1
Solution:
Since
and
can be simplified as
We know that the characteristics equation is
Since
<1, thus the system is stable.
Three stability tests can be applied directly to the characteristic equation without solving for the roots. Schur-Cohn stability test Jury Stability test Routh stability coupled with bi-linear transformation. Other stability tests like Lyapunov stability analysis are applicable for state space system models which will be discussed later. Computation requirements in Jury test is simpler than SchurCohn when the co-efficients are real which is always true for physical systems.
1.1 Jury Stability Test
Assume that the characteristic equation is a polynomial in z as follows,
, where
Jury Table
where,
This system will be stable if:
Example 2: The characteristic equation is
Thus,
We will now check the stability conditions.
1.
First condition is satisfied.
2.
Second condition is satisfied.
3.
Third condition is satisfied.
Next we will construct the Jury Table.
Jury Table
Rest
of
the
elements
are
also
calculated
in
similar
fashion.
The
elements
are
One
can
see
All criteria are satisfied. Thus the system is stable.
Example
3: The
characteristic
equation
is
Thus
Stability conditions are:
1.
First condition is satisfied.
2.
Second condition is not satisfied.
Since one of the criteria is violated, we may stop the test here and conclude that the system is unstable. or presence of a root on the unit circle and in that case the system can at the most become marginally stable if rest of the conditions are satisfied.
indicates the
The stability range of a parameter can also be found from Jury's test which we will see in the next example.
Example 4: Consider the system shown in Figure 1. Find out the range of K for which the system is stable.
Solution:
The closed loop transfer function:
Characteristic equation: Since it is a second order system only 3 stability conditions will be there.
Combining all, the range of K is found to be 0 < K < 2.39 .
If K = 2.39 , system becomes critically stable. The characteristics equation becomes:
Sampling period T = 1 sec.
rad/sec
The above frequency is the frequency of sustained oscillation.
1.2 Singular Cases
The situation, when some or all of the elements of a row in the jury table are zero, indicates the presence of roots on the unit circle. This is referred to as a singular case. It can be avoided by expanding or contracting unit circle infinitesimally by an amount transformation is:
which
is equivalent to move the roots of
P(z) off
the unit circle. The
where is a very small number. When is positive the unit circle is expanded and when is negative the unit circle is contracted. The difference between the number of zeros found inside or outside the unit circle when the unit circle is expanded or contracted is the number of zeros on the unit circle.
Since
for both positive and negative , the transformation requires the coefficient of the
zn term
be
multiplied by
Example 5: The characteristic equation:
Thus,
We will now check the stability conditions.
1.
First condition is satisfied.
2.
Second condition is satisfied.
3.
Third condition is satisfied.
Jury Table
Since the element
b1 is zero, we know that some of the roots lie on the unit circle.
If we replace
z by (1 + ) z, the characteristic equation would become:
First Jury Table
three
stability
conditions
are
satisfied
when
0+.
and
Since,
when
0+,
, thus which implies that the roots which are not on the unit circle are actually inside it and the system is marginally stable. The roots of the characteristic equation are found out to be i and 0.25 which verifies our conclusion.
