Файл: Microcontroller based applied digital control (D. Ibrahim, 2006).pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 14.06.2025

Просмотров: 3277

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

200 System Stability

or

(z2 − 1.5z + 0.5) − (z − 0.2)(2z − 1.5) = 0,

which gives

z2 − 0.4z − 0.2 = 0

and the roots are at

z = −0.290 and z = 0.689.

4. The value of k at the breakaway points can be calculated from

k

= −

1

F (z) z=−0.290,0.689

which gives k = 0.12 and k = 2.08. The root locus of the system is shown in Figure 8.7. It is clear from this plot that the system is always stable since all poles are inside the unit circle for all values of k.

Lines of constant damping factor and constant angular frequency are plotted on the same axis in Figure 8.8.

Assuming that T = 1 s, ωn > 0.6 if the roots are on the left-hand side of the constant angular frequency line ωn = 0.2π/ T . The damping factor will be greater than 0.6 if the roots are below the constant damping ratio line ζ = 0.6. A point satisfying these properties has been chosen

Figure 8.7 Root locus for Example 8.10


NYQUIST CRITERION

201

Figure 8.8 Root locus with lines of constant damping factor and natural frequency

and shown in Figure 8.9. The roots at this point are given as s1,2 = 0.55 ± j 0.32. The value of k can now be calculated as

k

1

= − F (z)

0.55

j 0.32

z

=

±

which gives k = 0.377.

8.5 NYQUIST CRITERION

The Nyquist criterion is one of the widely used stability analysis techniques in the s-plane, based on the frequency response of the system. To determine the frequency response of a continuous system transfer function G(s), we replace s by j ω and use the transfer function G( j ω). In the s-plane, the Nyquist criterion is based on the plot of the magnitude |GH ( j ω)| against the angle GH ( j ω) as ω is varied.

In a similar manner, the frequency response of a transfer function G(z) in the z-plane can be obtained by making the substitution z = e j ωT . The Nyquist plot in the z-plane can then be obtained by plotting the magnitude of |GH (z)|z=e j ωT against the angle GH (z)|z=e j ωT as ω is varied. The criterion is then

Z = N + P,


202 System Stability

Figure 8.9 Point for ζ > 0.6 and ωn > 0.6

where N is the number of clockwise circles around the point −1, P the number of poles of GH(z) that are outside the unit circle, and Z the number of zeros of GH(z) that are outside the unit circle.

For a stable system, Z must be equal to zero, and hence the number of anticlockwise circles around the point –1 must be equal to the number of poles of GH(z).

If GH(z)has no poles outside the unit circle then the criterion becomes simple and for stability the Nyquist plot must not encircle the point −1.

An example is given below.

Example 8.11

The transfer function of a closed-loop sampled data system is given by

G(z)

,

1 + GH(z)

where

GH(z) =

0.4

.

(z

0.5)(z

0.2)

Determine the stability of this system using the Nyquist criterion. Assume that T = 1 s.


NYQUIST CRITERION

203

Solution

Setting z = e j ωT = cos ωT + j sin ωT = cos ω + j sin ω,

G(z)|z=e j ωT =

0.4

(cos ω

+

j sin ω

0.5)(cos ω

+

j sin ω

0.2)

or

0.4

G(z)|z=e j ωT = (cos2 ω − sin2 ω − 0.7 cos ω + 0.1) + j (2 sin ω cos ω − 0.7 sin ω) .

This has magnitude

0.4

|G(z)| =

(cos2 ω − sin2 ω − 0.7 cos ω + 0.1)2 + (2 sin ω cos ω − 0.7 sin ω)2

and phase

G(z)

=

tan−1

2 sin ω cos ω − 0.7 sin ω

.

cos2 ω − sin2 ω − 0.7 cos ω + 0.1

Table 8.2 lists the variation of the magnitude of G(z) with the phase angle. The Nyquist plot for this example is shown in Figure 8.10. Since N = 0 and P = 0, the closed-loop system has no poles outside the unit circle in the x -plane and the system is stable.

The Nyquist diagram can also be plotted by transforming the system into the w-plane and then using the standard s-plane Nyquist criterion. An example is given below.

Example 8.12

The open-loop transfer function of a unity feedback sampled data system is given by

G(z) =

z

.

(z

1)(z

0.4)

Derive expressions for the magnitude and the phase of |G(z)| by transforming the system into the w-plane.

Solution

The w transformation is defined as

z = 1 + w

1 − w

which gives

G(w)

=

+

(1 + w)/(1 − w)

=

1 + w

,

((1

w/1

w)

1)((1

+

w/1

w)

0.4)

2w(0.6

+

1.4w)

or

1 + w

G(w) = 1.2w + 2.8w2 .


Table 8.2 Magnitude and phase of G(z)

w

|G(z)|

G(z)

0

1.0000E+000

0

1.0472E−001

9.8753E−001

−1.9430E+001

2.0944E−001

9.5248E−001

−3.8460E+001

3.1416E−001

9.0081E−001

−5.6779E+001

4.1888E−001

8.3961E−001

−7.4209E+001

5.2360E−001

7.7510E−001

−9.0690E+001

6.2832E−001

7.1166E−001

−1.0625E+002

7.3304E−001

6.5193E−001

−1.2096E+002

8.3776E−001

5.9719E−001

−1.3492E+002

9.4248E−001

5.4789E−001

−1.4820E+002

1.0472E+000

5.0395E−001

−1.6089E+002

1.1519E+000

4.6505E−001

−1.7308E+002

1.2566E+000

4.3075E−001

1.7518E+002

1.3614E+000

4.0058E−001

1.6384E+002

1.4661E+000

3.7408E−001

1.5283E+002

1.5708E+000

3.5082E−001

1.4213E+002

1.6755E+000

3.3044E−001

1.3168E+002

1.7802E+000

3.1259E−001

1.2147E+002

1.8850E+000

2.9698E−001

1.1146E+002

1.9897E+000

2.8337E−001

1.0162E+002

2.0944E+000

2.7154E−001

9.1945E+001

2.1991E+000

2.6130E−001

8.2401E+001

2.3038E+000

2.5250E−001

7.2974E+001

2.4086E+000

2.4501E−001

6.3646E+001

2.5133E+000

2.3872E−001

5.4404E+001

2.6180E+000

2.3354E−001

4.5232E+001

2.7227E+000

2.2939E−001

3.6118E+001

2.8274E+000

2.2622E−001

2.7050E+001

2.9322E+000

2.2399E−001

1.8015E+001

3.0369E+000

2.2266E−001

9.0018E+000

Figure 8.10 Nyquist plot for Example 8.11