Chapter 2-1 PDF
Chapter 2-1 PDF
Chapter 2-1 PDF
Settling time (ts):- The time at which the response of the system
reaches in to 98% or 95% of its final value.
• For 2% settling (98% reaching):
𝑦 𝑡𝑠 = 98% 𝑦 𝑡𝑓
= 0.98
• For 5% settling (95% reaching):
𝑦 𝑡𝑠 = 95% y 𝑡𝑓
= 0.95
Cont’d ...
Maximum Percent Overshoot (Mp):- It measures the deviation between the peak
value of the response and the final value.
𝑦 𝑡𝑝 − 𝑦(∞)
𝑀𝑝 = * 100%
𝑦(∞)
Cont’d ...
The s-plane expressions of the time domain specification requirements
are:
- Natural frequency (𝜔𝑛 ):
1.8
𝜔𝑛 ≥
𝑡𝑟
- Magnitude of the real part of the pole (𝜎):
4.6
𝜎 = 𝑅𝑒 𝑆𝑖 = 𝜉𝜔𝑛 ≥
𝑡𝑠
⇒ 2𝜋𝑘𝜔𝑛 = 1.8 𝜔𝑠
1.8
⇒ 𝜔𝑛 = 𝜔𝑠
2𝜋𝑘
0.29
⇒ 𝜔𝑛 = 𝜔𝑠
𝑘
𝜔𝑠
∴ 𝜔𝑛 =
3.5𝑘
In other words, the sample rate 𝜔𝑠 , should be 3.5𝑘 times faster than the natural
frequency 𝜔𝑛 . That means,
𝜔𝑠 = 3.5𝑘𝜔𝑛
2.2 Design by Emulation
A controller design is done as if the system is to be continuous and after a
sample period is selected, a discrete equivalent is computed and used in place of
the continuous design. This discrete controller may then be simulated and tested
in the discrete control loop and modifications made if necessary.
2.2.1 Discrete Equivalent Controllers
Design steps of continuous time lead compensator;
General transfer function of lead compensator:
𝑠+1 𝜏
𝐷 𝑠 = 𝐾𝑐 1 where 0 < 𝛼 < 1
𝑠+ 𝛼𝜏
Step 1:- From the given performance specifications of the system, determine the location
of the desired closed-loop poles.
𝑠𝑑 = −𝜉𝜔𝑛 ± 𝑗𝜔𝑛 1 − 𝜉 2
Step 2:- From the root locus of uncompensated system, check whether gain adjustment
alone can yield the requirement, if not calculate the angle deficiency (𝜙).
Cont’d…
Step 3:- Calculate the poles and zeros of compensator by assuming the zeros angle
contribution is equal to 90𝑜 .
Step 4:- Find the open-loop gain of the compensator from magnitude criteria.
Example 1: Consider the gain of the system:
4
𝐺𝐻(𝑠) = where H(s) is unity feedback
𝑠(𝑠+2)
iv) 𝑡𝑟 = 2𝑇
Cont’d…
Continuous time Lead Compensator Design
Step 1: Determining the location of the desired closed-loop poles:
• From the requirements, let 𝜉 = 0.44724, 𝜎 = 1.5
Then 𝜎 = 𝜉𝜔𝑛 = 1.5
𝜎 1.5
⟹ 𝜔𝑛 = = = 3.354
𝜉 0.44724
• The desired closed-loop poles will be
𝑠𝑑 = −𝜉𝜔𝑛 ± 𝑗𝜔𝑛 1 − 𝜉 2
= −𝜎 ± 𝑗𝜔𝑛 1 − 𝜉 2
= −1.5 ± 𝑗3.354 1 − 0.44724 2
= −1.5 ± 𝑗3
Step 2: Calculating the angle deficiency (𝜙):
• From angle criteria
𝜙 − 𝜃1 + 𝜃2 = −180°
⟹ 𝜃1 = 90° + tan−1 1.5 3 = 116.57°
⟹ 𝜃2 = tan−1 3 0.5 = 80.54°
Cont’d…
⟹ 𝜙 − 𝜃1 + 𝜃2 = −180°
⟹ 𝜙 − 116.57° + 80.54° = −180°
⟹ 𝜙 − 197.11° = −180°
⟹ 𝜙 = −180° + 197.11° = 17.11°
This angle deficiency should be corrected by the compensator. This means the angle
contributed by the compensator should be 17.11°.
⟹ 𝜙 = ∠𝐷 𝑠 |s=𝑠𝑑
Step 3: Calculating the poles and zeros of compensator:
𝑠+1 𝜏
𝐷 𝑠 = 𝐾𝑐 1 zero = −1 𝜏 , pole = −1 𝛼𝜏
𝑠+ 𝛼𝜏
• From our general design steps (refer to step 3), the angle contribution of the compensator zero
is equal to 90𝑜 .
⟹ 𝑅𝑒 𝑠𝑑 = 𝑧𝑒𝑟𝑜 ′ 𝑠 𝑜𝑓 𝐷 𝑠
⟹ −1 𝜏 = −1.5
⟹ 1 = 1.5
𝜏
• The pole of the compensator can be obtained from the above s-plane map by using a
trigonometric approach.
Cont’d… 1 1
−
𝛼𝜏 𝜏
⟹ tan 𝜙 =
3
1
−1.5
𝛼𝜏
⟹ tan 17.11° =
3
1
⟹ 3 tan 17.11° = − 1.5
𝛼𝜏
1
⟹ 1.5 + 3 tan 17.11° =
𝛼𝜏
1
⟹ = 2.4235
𝛼𝜏
Then
𝑠 + 1.5
𝐷 𝑠 = 𝐾𝑐
𝑠 + 2.4235
Step 4: Open-loop gain of the compensator
• It can be obtained from magnitude criteria as follows:
⟹ |𝐷 𝑠 𝐺𝐻 𝑠 |s=𝑠𝑑 = 1
Cont’d…
𝑠+1.5
⟹ |4𝐾𝑐 | =1
𝑠+2.4235 𝑠(𝑠+2) s=𝑠𝑑
(−1.5+𝑗3+1.5)
⟹ |4𝐾𝑐 | =1
−1.5+𝑗3+2.4235 (−1.5+𝑗3)(−1.5+𝑗3+2)
𝑗3
⟹ 4𝐾𝑐 =1
0.9235+𝑗3 −1.5+𝑗3 0.5+𝑗3
3
⟹ 4𝐾𝑐 =1
3.139∗3.354∗3.041
12𝐾𝑐
⟹ =1
3.139∗3.354∗3.041
⟹ 𝐾𝑐 = 2.668
⟹ 𝑇 = 0.268 𝑠𝑒𝑐
Cont’d…
• The discrete controller will have the following form:
𝑧 − 𝑧1
𝐷 𝑧 =𝐾
𝑧 − 𝑝1
Where
𝑧1 = 𝑒 𝑠𝑇 |𝑠=−1.5 𝑝1 = 𝑒 𝑠𝑇 |𝑠=−2.4235
= 𝑒 −1.5∗0.268 = 𝑒 −2.4235∗0.268
= 0.669 = 0.522
𝑧−0.669 𝑠+1.5
⟹ lim 𝐾 = lim 2.668
𝑧→1 𝑧−0.522 𝑠→0 𝑠+2.4235
1−0.669 2.668∗1.5
⟹𝐾 =
1−0.522 2.4235
0.331 2.668∗1.5
⟹𝐾 =
0.478 2.4235
2.668∗1.5∗0.478
⟹𝐾=
2.4235∗0.331
⟹ 𝐾 = 2.385
In general, the discrete equivalent of the designed continuous time lead compensator is
𝑧 − 0.669
𝐷 𝑧 = 2.385
𝑧 − 0.522
Implementing the Controller
Example 2: Give the difference equation that corresponds to the discrete transfer
function of the controller D(z) which can be used for implementation in a computer.
Solution: The difference equation of the discrete time compensator is
𝑈 𝑧 𝑧−0.669
⟹ = 2.385
𝐸 𝑧 𝑧−0.522
𝑈 𝑧 2.385𝑧−1.596
⟹ =
𝐸 𝑧 𝑧−0.522
a0 a n k b0 bn1k c0 cn2k
bk ck dk
an ak bn1 bk cn 2 ck
• The necessary and sufficient condition for stability is then:
b0 bn 1
a0 a n c0 cn 2
(1) (2)
P(1) 0 d 0 d n 3
(1) P(1) 0
n
m0 m2
Cont’d…
• Steps
i. From the given system description obtain the characteristic equation.
ii. Apply group 1 constraints, if any one fail, then the system is unstable, stop.
iii. If group 1 constraints are valid, start filling Jury’s table and check the
condition at each row. If condition fails, the system is unstable, stop.
iv. If all constraints are fulfilled, then the system is stable.
Example: for a system having the following characteristic equation, check
stability using Jury’s test.
P( z ) z 1.75 z 0.625 z 0.25
3 2
Solution:
Number of rows of Jury table: 2𝑛 − 3 = 2 ∗ 3 − 3 = 3
From the given characteristic equation: 𝑎3 = 1 , 𝑎2 = −1.75 , 𝑎1 = −0.625 , 𝑎0 = 0.25
Cont’d…
Jury table:
Row 𝑧0 𝑧1 𝑧2 𝑧3
1 𝑎0 𝑎1 𝑎2 𝑎3
2 𝑎3 𝑎2 𝑎1 𝑎0
3 𝑏0 𝑏1 𝑏2
Then
a0 a3 0.25 1
b0 (0.25 * 0.25) (1*1) 0.0625 1 0.9375
a3 a0 1 0.25
a0 a2 0.25 1.75
b1 {0.25 * (0.625)} (1.75 *1) 0.15625 1.75 1.59375
a3 a1 1 0.625
a0 a1 0.25 0.625
b2 {0.25 * (1.75)} (0.625 *1) 0.4375 0.625 0.1875
a3 a2 1 1.75
Cont’d…
Let us examine the condition for stability:
1) 𝑎0 < 𝑎3 ⟹ 0.25 < 1
i) 𝑃 𝑧 = 𝑧 4 − 𝑧 3 + 2𝑧 2 + 0.5
• The left half of the s-plane maps to a unit circle in the z-plane.
• The right half maps to a point out of the unit circle in z-plane.
Cont’d…
Generally, for a system given by transfer function:
n2
H (s) 2
s 2n s n2
ln r
ln 2 r 2
1
n ln 2 r 2
T
1 T
n ln r
Cont’d…
In summary, to get the specifications on acceptable pole locations in
the z-plane
• Estimate the desired 𝜔𝑛 , 𝜉, and 𝑀𝑝 from the continuous-time response
specifications. Compute σ = 𝜉𝜔𝑛 .
• Compute the radius r = 𝑒 −𝜉𝜔𝑛 𝑇 = 𝑒 −σ𝑇 .
• Obtain a plot of the z-plane showing lines of fixed damping and 𝜔𝑛 . The
MATLAB command zgrid will do this, plotting 𝜉 in steps of 0.1 from 0.1
to 0.9 and
𝑁𝜋
𝜔𝑛 = 𝑓𝑜𝑟 𝑖𝑛𝑡𝑒𝑔𝑒𝑟 𝑁 𝑓𝑟𝑜𝑚 1 𝑡𝑜 10.
10𝑇
• Mark the region of acceptable closed-loop pole locations on the plane.
The specification of steady state error:
consider the following discrete closed-loop system.
Proceeding as we for the continuous system, suppose the input r is a step, r 𝑘 = 1(𝑘), and the
disturbance w is zero. The transform of the error is computed using the same block-diagram
reduction tools that apply for continuous systems represented by their Laplace transforms,
except that now we use 𝐷 𝑧 and 𝐺(𝑧).
Doing this yields the transform of the error,
𝑅(𝑧)
𝐸 𝑧 =
1+𝐷 𝑧 𝐺 𝑧
If 𝐷 𝑧 𝐺(𝑧) has a pole at 𝑧 = 1, then the steady state error of Type 0 system for unit step
input is zero. Suppose there is a single pole at 𝑧 = 1. Then we have a Type 1 system and we
can compute the steady state error to a unit ramp input, r 𝑘𝑇 = 𝑘𝑇 1(𝑘𝑇) as follows.
Cont’d…
(ii) Type 1 system:
𝑇𝑧 1
𝐸 𝑧 = 2
𝑧−1 1+𝐷 𝑧 𝐺 𝑧
• Then, the steady-state error can be obtained by using final value theorem:
𝑒 ∞ = lim 𝑧 − 1 𝐸(𝑧)
𝑧→1
𝑇𝑧 1
= lim 𝑧 − 1
𝑧→1 𝑧−1 2 1+𝐷 𝑧 𝐺 𝑧
𝑇𝑧
= lim
𝑧→1 𝑧−1 (1+𝐷 𝑧 𝐺 𝑧 )
1
≜
𝐾𝑣
𝑧 − 𝑧1 𝑧 − 𝑧2 𝑧 − 𝑧3 … 𝑧 − 𝑧𝑛
𝐻 𝑧 =𝐾
𝑧 − 𝑝1 𝑧 − 𝑝2 𝑧 − 𝑝3 … 𝑧 − 𝑝𝑛
• Now suppose that 𝐻 𝑧 is the closed-loop transfer function that results from a Type 1 system, which implies
that the steady-state error of this system to a step input is zero.
𝐸 𝑧 =𝑅 𝑧 −𝑌 𝑧 but 𝑌 𝑧 = 𝐻 𝑧 𝑅(𝑧)
Cont’d…
Then
𝐸 𝑧 =𝑅 𝑧 −𝐻 𝑧 𝑅 𝑧
= [1 − 𝐻 𝑧 ]𝑅 𝑧
• The steady-state error of the system for unit step input can be obtained by using final value theorem as:
𝑒 ∞ = lim 𝑧 − 1 𝐸(𝑧)
𝑧→1
= lim 𝑧 − 1 [1 − 𝐻 𝑧 ]𝑅 𝑧
𝑧→1
𝑧
= lim 𝑧 − 1 1 − 𝐻 𝑧
𝑧→1 𝑧−1
= lim 1 − 𝐻 𝑧 𝑧
𝑧→1
= lim [1 − 𝐻 𝑧 ]
𝑧→1
• Since the system is a Type 1 system, its steady-state error for unit step input is zero. Then
𝑒 ∞ = lim [1 − 𝐻 𝑧 ]
𝑧→1
0 = lim [1 − 𝐻 𝑧 ]
𝑧→1
0=1−𝐻 1
∴ 𝐻 1 =1
Cont’d…
• In addition to this, the steady-state error of the system for a ramp input can be expressed as the
following:
𝑒 ∞ = lim 𝑧 − 1 𝐸(𝑧)
𝑧→1
= lim 𝑧 − 1 1 − 𝐻 𝑧 𝑅 𝑧
𝑧→1
𝑇𝑧
= lim 𝑧 − 1 1 − 𝐻 𝑧
𝑧→1 𝑧−1 2
𝑇𝑧
= lim (1 − 𝐻 𝑧 )
𝑧→1 𝑧−1
1
≜
𝐾𝑣
• Therefore (omitting a factor of z in the numerator, which makes no difference in the result)
1 1−𝐻 𝑧
= lim
𝑇𝐾𝑣 𝑧→1 𝑧 − 1
1−𝐻 1 1−1 0
= = =
1−1 1−1 0
Cont’d…
Since the limit is indeterminate, we can use L’Hopital’s rule
𝑑
1 1−𝐻 𝑧
= lim 𝑑𝑧
𝑇𝐾𝑣 𝑧→1 𝑑 (𝑧 − 1)
𝑑𝑧
𝑑
= lim − 𝐻(𝑧)
𝑧→1 𝑑𝑧
However, at 𝑧 = 1, we have
𝑑 1 𝑑 1 𝑑 𝑑
ln 𝐻(𝑧) = 𝐻 𝑧 = 𝐻 𝑧 = 𝐻(𝑧)
𝑑𝑧 𝐻(𝑧) 𝑑𝑧 𝐻(1) 𝑑𝑧 𝑑𝑧
so that
1 𝑑
= lim − ln 𝐻 𝑧
𝑇𝐾𝑣 𝑧→1 𝑑𝑧
Cont’d…
𝑑 (𝑧−𝑧𝑖 )
= lim − ln 𝐾
𝑧→1 𝑑𝑧 (𝑧−𝑝𝑖 )
𝑑
= lim − ln 𝑧 − 𝑧𝑖 − ln 𝑧 − 𝑝𝑖 + ln 𝐾
𝑧→1 𝑑𝑧
𝑑
= lim ln 𝑧 − 𝑝𝑖 − ln 𝑧 − 𝑧𝑖 − ln 𝐾
𝑧→1 𝑑𝑧
1 1
= lim −
𝑧→1 𝑧−𝑝𝑖 𝑧−𝑧𝑖
𝑛 𝑛
1 1
= −
1 − 𝑝𝑖 1 − 𝑧𝑖
𝑖=1 𝑖=1
Note:-
The farther the poles of the closed-loop system are from 𝑧 = 1, the larger the velocity
error constant and the smaller the errors.
Similarly, 𝐾𝑣 can be increased and the errors decreased by zeros close to 𝑧 = 1.
2.4.2 The Discrete Root Locus
The characteristic equation of the single-loop system is:
1+𝐺 𝑧 𝐻 𝑧 =0
The mechanics of drawing the root loci are exactly the same in the z-plane as in the s-plane.
The difference lies in the interpretation of the results because the pole locations in the z-plane
mean different things than pole locations in the s-plane when we come to interpret system
stability and dynamic response.
Angle and magnitude conditions
The characteristic equation can be written as:
𝐺 𝑧 𝐻 𝑧 = −1 ⇒ 𝐺 𝑧 𝐻 𝑧 ∠𝐺 𝑧 𝐻 𝑧 = 1∠180°
• Angle condition:
∠𝐺 𝑧 𝐻 𝑧 = ±180° 2𝑘 + 1 , 𝑘 = 0,1,2, …
• Magnitude condition:
𝐺 𝑧 𝐻 𝑧 =1
Cont’d…
• The value of z that fulfill both the angle and magnitude conditions are the roots of the
characteristic equation or the closed-loop poles.
• A plot of the points in the complex plane satisfying the angle condition alone is the
root locus. The roots of the characteristic equation corresponding to a given value of the
gain can be located on the root loci by the use of the magnitude condition.
General procedures for constructing root loci:
Step 1: Obtain the characteristic equation:
1+𝐺 𝑧 𝐻 𝑧 =0
and then rearrange this equation so that the parameter of interest such as given 𝐾 appears
as the multiplying factor in the form:
𝑧+𝑧1 𝑧+𝑧2 …(𝑧+𝑧𝑚 )
1+K =0 where 𝐾 > 0
𝑧+𝑝1 𝑧+2 …(𝑧+𝑝𝑛 )
Cont’d…
• From the factored form of the closed-loop transfer function, find the transfer function of
the open-loop pulse transfer function:
𝑚
𝑖=1 𝑧 + 𝑧𝑖
𝐺 𝑧 𝐻 𝑧 =𝐾 𝑛
𝑗=1 𝑧 + 𝑝𝑗
Then locate poles and zeros of the open-loop pulse transfer function in the z-plane.
Theorem 1: A root locus starts from an open-loop poles of the pulse transfer function at
K = 0 and terminates at a finite or infinite open-loop zero at 𝐾 = ∞.
Theorem 2: The number of branches of the root locus is max(m, n), where m and n are
the number of finite zeros and poles of 𝐺 𝑧 𝐻 𝑧 respectively.
Step 2: Determine the root loci on the real axis.
o Root loci on the real axis are determined by open-loop poles and zeros lying on it.
The conjugate complex poles and zeros of the open-loop pulse transfer function have no
effect on the location of the root loci on the real axis because the angle contribution of a
pair of conjugate complex poles or zeros is 360o on the real axis.
Cont’d…
o In constructing a root loci on the real axis, choose a test point on it.
If the total number of real poles and zeros to the right of this test point is odd, then
this point lies on a root locus.
Step 3: Determine the asymptotes of the root loci.
• If a test point z is located far from the origin, then the angles of all the complex poles
and zeros may be considered the same. One open-loop zero and one open-loop pole
then cancel each other. Therefore, the root loci for very large values of z must be
asymptotic to straight lines whose angles are given as follows:
±180° 2𝑘 + 1
A𝑛𝑔𝑙𝑒 𝑜𝑓 𝑎𝑠𝑦𝑚𝑝𝑡𝑜𝑡𝑒𝑠 = , 𝑘 = 0,1,2, …
𝑛−𝑚
j T 2 T
ze w j tan
T 2
Cont’d…
• Hence the frequency in s-plane and the frequency in w-plane are related by:
2 T
w tan
T 2
1+KGH(z) = 0
• The root locus of this system will be drawn using the standard steps of
continuous system.
Cont’d…
• Root locus method
• Example: Draw the root locus of the following system if 𝐷 𝑧 = 𝐾.
1 KG ( z ) 0
Where
z 1 1
G( z) 2
z s s 1
Cont’d…
• Using z-transform table and T = 1, the final characteristic equation becomes:
K (0.3679 z 0.2642)
1 KG ( z )
( z 1)( z 0.3679)
0.5
Imaginary Axis
-0.5
-1
-1.5
-5 -4 -3 -2 -1 0 1 2
Real Axis
Cont’d…
• Bode plot:- is another frequency response technique.
Graphical method which draws the magnitude and phase of the
closed-loop response as a function of frequency.
For discrete time systems two methods can be followed:
• Form the characteristic equation of the system in z.
• Change to w-plane and plot bode plot or
• Draw the bode plot in z-domain and convert it to s-domain using the
conversion formula.
2.6 Direct Design Method of Ragazzini
• Suppose we are given a discrete transfer function 𝐺 𝑧 of the plant (plus hold) and a desired
transfer function H(𝑧) between reference 𝑅(𝑧) and output 𝑌(𝑧). The structure is assumed to be
a unity feedback system and the design is to select the computer transfer function D(𝑧) to
realize 𝐻(𝑧).
• The overall transfer function is given by the formula:
𝐷(𝑧)𝐺(𝑧)
H z =
1 + 𝐷(𝑧)𝐺(𝑧)
This design calls for a 𝐷(𝑧) that will cancel the plant effects and add whatever is necessary to give the
desired result.
The problem is to discover and implement constraints on 𝐻(𝑧).
Consideration on constraints
i. The design must be causal.
If 𝐷 𝑧 is to be causal, then as 𝑧 → ∞ its transfer function is bounded and does not have a pole at
infinity. If 𝐺 𝑧 were to have a zero at infinity, then 𝐷 𝑧 would have a pole there. This pole should be
cancelled by 𝐻 𝑧 .
Thus the constraint for 𝐷 𝑧 to be causal is:
𝐻 𝑧 𝑚𝑢𝑠𝑡 ℎ𝑎𝑣𝑒 𝑎 𝑧𝑒𝑟𝑜 𝑎𝑡 𝑖𝑛𝑓𝑖𝑛𝑖𝑡𝑦 𝑜𝑓 𝑡ℎ𝑒 𝑠𝑎𝑚𝑒 𝑜𝑟𝑑𝑒𝑟 𝑎𝑠 𝑡ℎ𝑒 𝑧𝑒𝑟𝑜 𝑜𝑓 𝐺 𝑧 𝑎𝑡 𝑖𝑛𝑓𝑖𝑛𝑖𝑡𝑦.
ii. The roots of the characteristic equation of the closed-loop system are the roots
of equation:
1+𝐷 𝑧 𝐺 𝑧 =0
𝑐 𝑧 𝑏(𝑧)
• Let 𝐷 𝑧 = , 𝐺(𝑧) = 𝑤ℎ𝑒𝑟𝑒 𝑎, 𝑏, 𝑐, 𝑎𝑛𝑑 𝑑 𝑎𝑟𝑒 𝑝𝑜𝑙𝑦𝑛𝑜𝑚𝑖𝑎𝑙𝑠.
𝑑 𝑧 𝑎(𝑧)
• 𝐻 𝑧 𝑚𝑢𝑠𝑡 𝑐𝑜𝑛𝑡𝑎𝑖𝑛 𝑎𝑠 𝑧𝑒𝑟𝑜𝑠 𝑎𝑙𝑙 𝑡ℎ𝑒 𝑧𝑒𝑟𝑜𝑠 𝑜𝑓 𝐺 𝑧 𝑡ℎ𝑎𝑡 𝑎𝑟𝑒 𝑜𝑢𝑡𝑠𝑖𝑑𝑒 𝑡ℎ𝑒 𝑢𝑛𝑖𝑡 𝑐𝑖𝑟𝑐𝑙𝑒.
which implies
𝐻 1 =1
for (b): The velocity error constant requirement is that:
𝑇𝑧 1
𝑒 ∞ = lim (𝑧 − 1) 1 − 𝐻 𝑧 =
𝑧→1 (𝑧 − 1)2 𝐾𝑣
1 𝑇𝑧 1 − 𝐻 𝑧
= lim
𝐾𝑣 𝑧→1 𝑧−1
Cont’d…
We know that 1 − 𝐻 𝑧 is zero at 𝑧 = 1, so that to evaluate the limit, it is necessary to
use L’Hopital’s rule with the result:
𝑑
1 {𝑇𝑧 1 − 𝐻 𝑧 }
= 𝑑𝑧 |𝑧 = 1
𝐾𝑣 𝑑
(𝑧 − 1)
𝑑𝑧
𝑑
𝑇 1−𝐻 𝑧 −𝑇𝑧𝑑𝑧𝐻 𝑧
= |𝑧 = 1
1
𝑑
= 𝑇 1 − 𝐻(1) − 𝑇 𝐻 𝑧 |𝑧 = 1
𝑑𝑧
𝑑
= 𝑇 1−1 − 𝑇 𝐻 𝑧 |𝑧 = 1
𝑑𝑧
𝑑
= −𝑇 𝐻 𝑧 |𝑧 = 1
𝑑𝑧
Assignment #2
1. If the transfer function of the system is given as:
𝑃 𝑧 = 𝑧 3 + 0.5𝑧 2 − 2𝑘 + 1 𝑧 − 0.5
3. Draw the root locus of a negative feedback system with an open-loop transfer
function of
𝐾
𝐺𝐻 𝑧 =
𝑧(𝑧 + 3)2
4. Find the Bode plot of
𝑤−2
𝐺 𝑤 =
(𝑤 + 4)(𝑤 + 8)