2019handout1 Automatic Control Systems
2019handout1 Automatic Control Systems
2019handout1 Automatic Control Systems
In engineering the term control system is restricted to those systems whose major
function is to dynamically or actively command, direct, or regulate.
Example:
The figure shown below is an example of a control system which consists of a mirror
pivoted at one end and adjusted up and down with a screw at the other end. The
angle of reflected light is regulated by means of the screw.
The input is the stimulus or excitation applied to a control system from an external
energy source, usually in order to produce a specified response from the control
system.
The output is the actual response obtained from a control system. It may or may
not be equal to the specified response implied by the input.
The purpose of the control system usually identifies or defines the output and input.
Given the output and input, it is possible to identify or define the nature of the
system’s components.
A control system consists of subsystems and processes (or plants) assembled for
the purpose of obtaining a desired output with desired performance, given a
specified input.
disturbance 1 disturbance 2
input + + output
input controller process
or or
transducer or plant
reference + + controlled
summing summing variable
junction junction
disturbance 1 disturbance 2
input + + output
or input + controller + process + or
reference transducer or plant controlled
variable
output
transducer
or sensor
Computer-Controlled Systems
In many modern systems, the controller (or compensator) is a digital computer. The
advantage of using a computer is that many loops can be controlled or compensated
by the same computer through time sharing. Furthermore, any adjustments of the
compensator parameters required to yield a desired response can be made by
changes in software rather than hardware. The computer can also perform
supervisory functions, such as scheduling many required applications.
1.4 Feedback
Feedback is that property of a closed-loop control system which permits the output
(or some other controlled variable of the system) to be compared with the input to
the system (or an output to some other internally situated component or subsystem
of the system) so that the appropriate control action may be formed as some
function of the output and input.
For a system to be classified as a feedback control system, it is necessary that the
controlled variable be fed back and compared with the reference input. In addition,
the resulting error signal must actuate the control elements to change the output
so as to minimise the error.
Characteristics of Feedback
1. Increased accuracy. For example, the ability to faithfully reproduce the input.
2. Reduced sensitivity of the ratio of output to input to variations in system
characteristics.
3. Reduced effects of non-linearities and distortions.
4. Increased bandwidth. The bandwidth of a system is that range of frequencies
(of the input) over which the system will respond satisfactorily.
5. Tendency towards oscillation or instability.
Block Diagrams and Signal Flow Graphs are shorthand, graphical representations of
either the schematic diagram of a physical system, or the set of mathematical
equations characterising its parts.
Mathematical models, in the form of system equations, are employed when detailed
relationships are required. Every control system may theoretically be characterised
by mathematical equations. The solution to these equations represents the system’s
behaviour.
block
input output
The interior of the rectangle representing the block usually contains a description of
or the name of the element, or the symbol for the mathematical operation to be
performed on the input to yield the output. The arrows represent the direction of
unilateral information or signal flow.
Example:
+ - +
y y y
Any number of inputs may enter a summing point. The output is the algebraic sum
of the inputs.
In order to employ the same signal or variable as an input to more than one block or
summing point, a take-off point is used. This permits the signal to proceed
unaltered along several paths to several destinations.
x x take-off point
x x x
take-off x
point x
disturbance
actuating manipulated u
reference signal variable controlled
input control plant output
r + e=r±b elements m g2 c
± g1
b
forward path
primary
feedback
signal
feedback
elements
h
feedback path
The arrows of the closed-loop connecting one block with another represent the
direction of flow of control energy or information, and not the main source of
energy for the system.
Example: r = r(t)
1. The plant g2, also called the controlled system, is the body, process, or
machine, of which a particular quantity or condition is to be controlled.
2. The control elements g1, also called the controller, are the components
required to generate the appropriate control signal m applied to the plant.
7. The actuating signal e, also called the error or control action, is the
algebraic sum consisting of the reference input r plus or minus (usually minus)
the primary feedback b.
10. The forward path is the transmission path from the actuating signal e to the
controlled output c.
11. The feedback path is the transmission path from the controlled output c to
primary feedback signal b.
Negative feedback means that the summing point is a subtractor or e = r – b.
Positive feedback means that the summing point is an adder or e = r + b.
12. A Unity Feedback System is a feedback system in which the primary feedback
b is identically equal to the controlled output c.
Supplementary Terminology
13. A transducer is a device which converts one energy form into another.
Example: A potentiometer converts mechanical position into an electrical
voltage.
+
arm
v position position voltage
reference input input potentiometer output
voltage + v r
source r voltage
output
+
14. The command v is an input signal, usually equal to the reference input r. But
when the energy form of the command v is not the same as that of the primary
feedback b, a transducer is required between the command v and the reference
input r.
16. A stimulus or test input is any externally introduced input signal affecting the
controlled output c.
For example, the reference input r and a disturbance u are stimuli.
17. The time response of a system or element is the output as a function of time,
following the application of a prescribed input under specified operating
conditions.
Servomechanisms
A servomechanism is power-amplifying feedback control system in which the
controlled variable c is mechanical position, or a time derivative of position such as
velocity or acceleration.
Example:
A motor vehicle power-steering apparatus is a servomechanism.
The command input is the angular position of the steering wheel. A small rotational
torque applied to the steering wheel is amplified hydraulically, resulting in a force
adequate to modify the output, the angular position of the front wheels. Negative
feedback is necessary in order to return the control valve to the neutral position,
reducing the torque from the hydraulic amplifier to zero when the desired wheel
position has been achieved.
2. Laplace Transforms
Laplace Transforms make solving linear Ordinary Differential Equations (ODEs) and
related initial value problems, as well as systems of linear ODEs, much easier.
Applications abound: electrical networks, springs, mixing problems, signal
processing, and other areas of engineering and physics.
The key motivation for learning about Laplace Transforms is that the process of
solving an ODE is simplified to an algebraic problem (and transformations). This type
of mathematics that converts problems of calculus to algebraic problems is known as
operational calculus.
If f(t) is a function defined for all t ≥ 0, its Laplace Transform is the integral of
f(t) times e-st from t = 0 to ∞. It is a function of s, say, F(s), and is denoted by L(f);
thus
Notation
Original functions depend on t and their transforms on s. Original functions are
denoted by lowercase letters and their transforms by the same letters in
capital, so that F(s) denotes the transform of f(t), and Y(s) denotes the transform
of y(t), and so on.
The Transfer Function will allow separation of the input, system, and output into
three separate and distinct parts. The Transfer Function will also allow us to
algebraically combine mathematical representations of subsystems to yield a total
system representation.
where c(t) is the output, r(t) is the input, and the ai’s, bi’s, and the form of the
differential equation represent the system.
This equation is a purely algebraic expression. If we assume that all initial conditions
are zero, the equation reduces to
Now form the ratio of the output transform, C(s), divided by the input transform,
R(s):
The transfer function can be represented as a block diagram with the input on the
left, the output on the right, and the system transfer function inside the block.
Notice that the denominator of the transfer function is identical with the
characteristic polynomial of the differential equation. Also, we can find the output,
C(s) by using
C(s) = G(s)R(s)
Solution:
Taking the Laplace Transform of both sides, assuming zero initial conditions we have
Equivalent circuits for the electric networks may consist of passive linear
components: resistors, capacitors, and inductors.
In the analysis we combine electrical components into circuits, decide on the input
and output, and find the transfer function. Our guiding principles are Kirchhoff’s
laws. We sum voltages around loops or sum currents at nodes, depending on which
technique involves the least effort in algebraic manipulation, and then equate the
result to zero. From these relationships we can write the differential equations for
the circuit. Then we can take the Laplace Transforms of the differential equations
and finally solve for the Transfer Function.
Problem: Find the transfer function relating the capacitor voltage, VC(s), to the input
voltage, V(s) in the figure below.
RLC network
Solution
In any problem, the designer must first decide what the input and output should be.
From the question, we are to treat the capacitor voltage as the output and the
applied voltage as the input.
Summing the voltages around the loop, assuming zero initial conditions, yields the
integro-differential equation for this network as
Taking the Laplace Transform assuming zero initial conditions, rearranging terms,
and simplifying yields
(LCs2 + RCs + 1)VC(s) = V(s)
Z(s) = Ls + R + 1/Cs
3. Instead of writing the differential equation first and then taking the Laplace
Transform, we can draw the transformed circuit and obtain the Laplace transform
of the differential equation simply by applying Kirchhoff’s voltage law to the
transformed circuit.
Summary of Steps
1. Redraw the original network showing all time variables, such as v(t), i(t), and
vC(t), as Laplace transforms V(s), I(s), and VC(s), respectively.
2. Replace the component values with their impedance values. This replacement is
similar to the case of dc circuits, where we represent resistors with their
resistance values.
Transfer Function
E I
1/Z(s)
1
Z(s) =
1 / L1s + 1 / L2s + 1 / R1 + 1 / R2 + C1s + C2s
Example: Let it be desired to determine the equation relating the output voltage E2
to the input voltage E1 in the circuit below.
Solution:
Transfer Function
E1 R2(1 + R1C1s) E2
R1 + R2 + R1R2C1s
If the input impedance to the amplifier is high, then by Kirchhoff’s current law
Example: Find the transfer function, Vo(s) / Vi(s), for the circuit given below
Solution
Mechanical systems parallel electrical networks to such an extent that there are
analogies between electrical and mechanical components and variables. Mechanical
systems, like electrical networks, have three passive, linear components:
The spring and the mass are energy-storage elements. The viscous damper
dissipates energy. The two energy-storage elements are analogous to the two
electrical energy-storage elements, the inductor and capacitor. The energy dissipater
is analogous to electrical resistance.
The mechanical system requires just one differential equation, called the equation of
motion, to describe it.
1. We draw a free-body diagram, placing on the body all forces that act on the body
either in the direction of motion or opposite to it.
2. We use Newton’s law to form a differential equation of motion by summing the
forces and setting the sum equal to zero.
3. Assuming zero initial conditions, we take the Laplace transform of the differential
equation, separate the variables, and arrive at the transfer function.
Example: Find the transfer function, X(s)/F(s) for the system below
Place on the mass all forces felt by the mass. We assume the mass is travelling
toward the right. Thus, only the applied force points to the right; all other forces
impede the motion and act to oppose it. Hence, the spring, viscous damper, and the
force due to acceleration point to the left.
2. Now write the differential equation of motion using Newton’s law to sum to zero
all of the forces shown on the mass
block diagram
In order to work such a problem, we draw the free-body diagram for each point of
motion and then use superposition. For each free-body diagram we begin by holding
all other points of motion still and finding the forces acting on the body due only to
its own motion. Then we hold the body still and activate the other points of motion
one at a time, placing on the original body the forces created by the adjacent
motion.
Example: Find the transfer function, X2(s) =F(s), for the system below
Considering M1 only
Considering M2 only
where
block diagram
Example: Write, but do not solve, the equations of motion for the mechanical
network below
The system has three degrees of freedom, since each of the three masses can be
moved independently while the others are held still.
For M1
For M3
M1 has two springs, two viscous dampers, and mass associated with its motion.
There is one spring between M1 and M2 and one viscous damper between M1 and
M3.
We can solve these equations for any displacement, X1(s); X2(s); or X3(s), or transfer
function.
In determining inertia force the acceleration of a mass is always taken with respect
to the earth.
For series mechanical elements, the force f is equal to the summation of the forces
acting on each individual component, and each element undergoes the same
displacement.
Transfer Function
f x
1/Z(s)
or f= 1 * x = Zx
1 / K1 + 1 / K2 + 1 / C1s + 1 / C2s
Z(s) = 1
1 / K1 + 1 / K2 + 1 / C1s + 1 / C2s
Transfer Function
f x
1/Z(s)
A necessary condition for parallel elements is that the force be transmitted through
each element. Springs and dampers satisfy this condition because the force is the
same on both sides. However, this is not the case for a mass such as the one in (a)
in the figure below, because the difference in forces acting on both sides of a mass
is utilised in acceleration. Thus, a mass located between other elements cannot be in
parallel with them.
mechanical system
Writing the equations of motion for rotational systems is similar to writing them for
translational systems; the only difference is that the free-body diagram consists of
torques rather than forces.
We obtain the torques using superposition.
1. We rotate a body while holding all other points still and place on its free-body
diagram all torques due to the body’s own motion.
2. Then, holding the body still, we rotate adjacent points of motion one at a time
and add the torques due to the adjacent motion to the free-body diagram.
3. The process is repeated for each point of motion.
4. For each free-body diagram, these torques are summed and set equal to zero to
form the equations of motion.
Solution:
1. Obtain the schematic from the physical system. Even though torsion occurs
throughout the rod, we approximate the system by assuming that the torsion acts
like a spring concentrated at one particular point in the rod, with an inertia J1 to
the left and an inertia J2 to the right. We also assume that the damping inside the
flexible shaft is negligible.
There are two degrees of freedom, since each inertia can be rotated while the
other is held still. Hence, it will take two simultaneous equations to solve the
system.
(a) shows the torques on J1 if J2 is held still and J1 rotated. (b) shows the torques
on J1 if J1 is held still and J2 rotated. Finally, (c) is the final free-body diagram for
J1. The same process is repeated in Figure 2.24 for J2.
where
General Form:
Block Diagram
Example:
Electrical analogue
Parallel Analogue
A system can also be converted to an equivalent parallel analogue.
Example:
Block Diagrams
A Block Diagram of a system is a pictorial representation of the functions performed
by each component and of the flow of signals. Such a diagram depicts the
interrelationships that exist among the various components. A Block Diagram is
shorthand, graphical representation of a physical system, illustrating the functional
relationships among its components.
Block Diagrams permit evaluation of the contributions of the individual elements to
the overall performance of the system.
In a block diagram all system variables are linked to each other through functional
blocks. The functional block or simply block is a symbol for the mathematical
operation on the input signal to the block that produces the output. The transfer
functions of the components are usually entered in the corresponding blocks, which
are connected by arrows to indicate the direction of the flow of signals.
Definitions
1. G = direct transfer function = forward transfer function
2. H = feedback transfer function
3. GH = loop transfer function = open-loop transfer function
4. C/R = closed-loop transfer function = control ratio
5. E/R = actuating signal ration = error ratio
6. B/R = primary feedback ratio
In the following equations, the minus (-) sign refers to a positive feedback system,
and the plus (+) sign refers to a negative feedback system.
R E C
G
±
Any feedback system with only linear elements in the feedback loop can be put into
the form of a unity feedback system.
Example:
The characteristic equation for the unity feedback system, determined from 1±G = 0
is
DG ± NG = 0
Step 1: Put U = 0
Step 2: The system reduces to
Step 1:
Step 2:
Step 3:
Step 5:
Steps 1 and 2
We do not apply step 3 at this time, but go directly to step 4, moving takeoff point 1
beyond block G2 + G3
We may now rearrange summing points 1 and 2 and combine the cascade blocks in
the forward loop using Transformation 6, then, 1
Step 3:
Note that the same result could have been obtained after applying step 2 by moving
takeoff point 2 ahead of G2 + G3, instead of takeoff point 1 beyond G2 + G3.
Block G2 + G3 has the same effect on the control ratio C/R whether it directly follows
R or directly precedes C.
5. Signal-Flow Graphs
Example:
signal V(s) = R1(s)G1(s) - R2(s)G2(s) + R3(s)G3(s)
signal C2(s) = V(s)G5(s) = R1(s)G1(s)G5(s) - R2(s)G2(s)G5(s) + R3(s)G3(s)G5(s)
signal C3(s) = -V(s)G6(s) = - R1(s)G1(s)G6(s) + R2(s)G2(s)G6(s) - R3(s)G3(s)G6(s).
Not: in summing negative signals we associate the negative sign with the system
and not with a summing junction
X
-4 Z
3. Multiplication Rule
A cascade (series) connection of n-1 branches with transmission functions
A21, A32, A43 ... An(n01) can be replaced by a single branch with a new transmission
function equal to the product of the old ones.
10 -20 -200
reduces to
X Y Z X Z
4. The value of the variable represented by a node is equal to the sum of all the
signals entering the node.
5. The value of the variable represented by a node is transmitted through all the
branches leaving the node.
For the signal-flow graph above, y6 = a16y1 y7 = a17y1 y8 = a18y1
6. Parallel branches in the same direction connected between two nodes can be
replaced by a single branch with gain equal to the sum of the gains of the parallel
branches.
Example:
Example:
In simplifying a graph with self-loop, the loop can be replaced by a branch with gain
b*, and the node can be removed by adjusting the gain of all the paths through the
node by b*.
5.2 Definitions
4. A forward path is a path from the input node to the output node.
Example: X1 to X2 to X3 to X4 and X1 to X2 to X4 are forward paths.
7. The gain of a branch is the transmission function of that branch when the
transmission function is a multiplicative operator.
Example: A33 is the gain of the self-loop if A33 is a constant or transfer function.
9. The loop path is the product of the branch gains of the loop.
Example: The loop gain of the feedback loop from X2 to X3 and back to X2 is
A32A23.
In the canonical feedback system, if the signal-flow graph were to be drawn directly
from the equations, the “output node” would require an outgoing branch, contrary
to definition. This situation may be remedied by adding a branch with a transmission
function of unity entering a “dummy” node.
Example
The signal-flow graph of a linear feedback control system whose components are
specified by non-interacting transfer functions can be constructed by direct reference
to the block diagram of the system. Each variable of the block becomes a node
and each block becomes a branch.
Example:
The block diagram of the canonical feedback control system is given by
2. Arrange the m or n (whichever is larger) nodes from left to right. The nodes may
be arranged if the required loops later appear too cumbersome.
4. If the desired output node has outgoing branches, add a dummy node and a
unity-gain branch.
5. Rearrange the nodes and/or loops in the graph to achieve maximum pictorial
clarity.
Solution:
There are five variables: v1, v2, v3, i1 and i2. Using Kirchhoff’s Voltage and Current
Laws, we can write four independent equations. Proceeding from left to right in the
schematic, we have
Laying out the five nodes in the same order with v1 as an input node, and
connecting the nodes with the appropriate branches, we have
We can reduce complicated block diagrams to canonical form and obtain the the
following control ratio
Let us denote the ratio of the input variable to the output variable by T. For linear
feedback control systems, T = C/R.
∆ = 1 – (sum of loop gains) + (sum of all gain products of two nontouching loops)
- (sum of all gain products of three nontouching loops) + ...
∆1 = ∆ evaluated with all loops touching Pi eliminated
Two loops, paths, or a loop and a path are said to be nontouching if they have no
nodes in common.
∆ is called the signal-flow graph determinant or characteristic function, since
∆ = 0 is the system characteristic equation.
Example:
Determine C/R for the signal-flow graph of the following canonical feedback system
Solution:
There is only one forward path. Hence
P1 = G P2 = P3 = ... = 0
P11 = GH Pjk = 0, j ≠ 1, k ≠ 1
Then ∆ = 1 – P11 = 1 ± GH and ∆1 = 1 – 0 = 1
Finally,
Example:
Assume that two identical resistance networks are to be cascaded and used as the
control elements in the forward loop of a control system. The networks are voltage
dividers of the form
Signal-Flow Graph:
There are five variables: v1, v2, v3, i1 and i2. Using Kirchhoff’s Voltage and Current
Laws, we can write four independent equations. Proceeding from left to right in the
schematic, we have
We observe that the feedback branch –R3 in the above graph does not appear in the
signal-flow graph of the cascaded signal-flow graphs of the individual networks
connected from v2 to v’1:
This means that as a result of connecting the two networks, the second one loads
the first, changing the equation for v2 from
It is good general practice to calculate the gain of cascaded networks directly from
the combined signal-flow graph. Most practical control system components load each
other when connected in series.
Often, the easiest way to determine the control ratio of a complicated block diagram
is to translate the block diagram into a signal-flow graph. Takeoff points and
summing points must be separated by a unity-gain branch in the signal-flow graph.
If the elements of G and H of a canonical representation are desire, the direct
transfer function is
Example: Determine the control ratio C/R and the canonical diagram of the feedback
control system below
Example: Convert the cascaded, parallel, and feedback forms of the block diagrams
shown in (a), (b), and (c), respectively, into signal-flow graphs.
(a)
(b)
(c)
(a)
(b)
(c)
Solution:
Begin by drawing the signal nodes. Next, interconnect the nodes, showing the
direction of signal flow and identifying each transfer function.
signal nodes
signal-flow graph
Two approaches are available for the analysis and design of feedback control
systems:
6. Mason’s Rule
Mason’s Rule for reducing a signal-flow graph to a single transfer function requires
the application of one formula. The formula was derived by S. J. Mason when he
related the signal-flow graph to the simultaneous equations that can be written from
the graph (Mason, 1953).
6.1 Definitions
Loop Gain: The product of branch gains found by traversing a path that starts at a
node and ends at the same node, following the direction of the signal flow, without
passing through any other node more than once.
Example
Forward-path Gain: The product of gains found by traversing a path from the
input node to the output node of the signal-flow graph in the direction of signal flow.
There are two forward-path gains:
1. G1(s)G2(s)G3(s)G4(s)G5(s)(s)G7(s)
2. G1(s)G2(s)G3(s)G4(s)G6(s)G7(s)
Nontouching-loop Gain: The product of loop gains from nontouching loops taken
two, three, four, or more at a time.
The product of loop gain G2(s)H1(s) and loop gain G4(s)H2(s) is a nontouching-loop
gain taken two at a time.
In summary, all three of the nontouching-loop gains taken two at a time are
1. [G2(s)H1(s)] [G4(s)H2(s)]
2. [G2(s)H1(s)] [G4(s)G5(s)H3(s)]
3. [G2(s)H1(s)] [G4(s)G6(s)H3(s)]
eq1
where
k = number of forward paths
Tk = the kth forward-path gain
Example: Find the transfer function, C(s)/R(s), for the signal-flow graph below.
Solution
Thus, the combinations of nontouching loops taken two at a time are as follows:
Loop 1 and loop 2: G2(s)H1(s)G4(s)H2(s)
Loop 1 and loop 3: G2(s)H1(s)G7(s)H4(s)
Loop 2 and loop 3: G4(s)H2(s)G7(s)H4(s)
4. The nontouching loops taken three at a time are as follows:
Loops 1, 2, and 3: G2(s)H1(s)G4(s)H2(s)G7(s)H4(s)
We form ∆k by eliminating from ∆ the loop gains that touch the kth forward path:
∆1 = 1 - G7(s)H4(s) eq4
Since there is only one forward path, G(s) consists of only one term, rather than a
sum of terms, each coming from a forward path.
7.1 Introduction
The output response of a system is the sum of two responses: the forced
response and the natural response. Although many techniques, such as solving a
differential equation or taking the inverse Laplace transform, enable us to evaluate
this output response, these techniques are laborious and time-consuming.
Productivity is aided by analysis and design techniques that yield results in a
minimum of time. If the technique is so rapid that we feel we derive the desired
result by inspection, we sometimes use the attribute qualitative to describe the
method. The use of poles and zeros and their relationship to the time response of a
system is such a technique.
The concept of poles and zeros, fundamental to the analysis and design of control
systems, simplifies the evaluation of a system’s response.
Generally, poles are defined as the roots of the denominator of a transfer function.
Strictly speaking, the poles of a transfer function satisfy part (1) of the definition.
For example, the roots of the characteristic polynomial in the denominator are
values of s that make the transfer function infinite, so they are thus poles. However,
if a factor of the denominator can be cancelled by the same factor in the numerator,
the root of this factor no longer causes the transfer function to become infinite. In
control systems, we often refer to the root of the cancelled factor in the
denominator as a pole even though the transfer function will not be infinite at this
value.
Generally, zeros are defined as the roots of the polynomial of the numerator of a
Transfer Function.
Strictly speaking, the zeros of a transfer function satisfy part (1) of this definition.
For example, the roots of the numerator are values of s that make the transfer
function zero and are thus zeros. However, if a factor of the numerator can be
cancelled by the same factor in the denominator, the root of this factor no longer
causes the transfer function to become zero. In control systems, we often refer to
the root of the cancelled factor in the numerator as a zero even though the transfer
function will not be zero at this value.
To show the properties of the poles and zeros, let us find the unit step response of
the system. Multiplying the transfer function G(s) by a step function yields
where
Thus,
1. A pole of the input function generates the form of the forced response (that
is, the pole at the origin generated a step function at the output).
2. A pole of the transfer function generates the form of the natural response
(that is, the pole at -5 generated e-5t).
3. A pole on the real axis generates an exponential response of the form e-αt,
where -α is the pole location on the real axis. Thus, the farther to the left a pole is
on the negative real axis, the faster the exponential transient response will decay
to zero (again, the pole at -5 generated e-5t.
4. The zeros and poles generate the amplitudes for both the forced and natural
responses.
Example: Given the system below, write the output, c(t), in general terms. Specify
the forced and natural parts of the solution.
Thus,
Exercise:
A system has a transfer function
Write, by inspection, the output, c(t), in general terms if the input is a unit step.
Summary
A first-order system without zeros can be described by the transfer function shown
above. If the input is a unit step, where R(s) = 1/s, the Laplace transform of the
step response is C(s), where
C(s) = R(s)G(s) = a / [s(s + a)]
where the input pole at the origin generated the forced response cf(t) = 1, and the
system pole at -a, as shown above, generated the natural response cn(t) = -e-at .
When t = 1/a
or
Time Constant
We call 1/a the time constant of the response. The time constant can be
described as the time for e-at to decay to 37% of its initial value. Alternately, the
time constant is the time it takes for the step response to rise to 63% of its final
value.
The reciprocal of the time constant has the units (1/seconds), or frequency. Thus,
we can call the parameter a the exponential frequency.
Since the derivative of e-at is -a when t = 0, a is the initial rate of change of the
exponential at t = 0.
Thus, the time constant can be considered a transient response specification
for a first-order system, since it is related to the speed at which the system
responds to a step input.
The time constant can also be evaluated from the pole plot. Since the pole of the
transfer function is at -a, we can say the pole is located at the reciprocal of the time
constant, and the farther the pole from the imaginary axis, the faster the transient
response.
Rise Time Tr
Rise time is defined as the time for the waveform to go from 0.1 to 0.9 of its final
value.
Rise time is found by solving c(t) = cf(t) + cn(t) = 1 – e-at for the difference in
time at c(t) = 0.9 and c(t) = 0.1.
Hence,
Settling Time Ts
Settling time is defined as the time for the response to reach, and stay within, 2%
of its final value.
Letting c(t) = 0.98 and solving for time, t, in c(t) = cf(t) + cn(t) = 1 – e-at
we find the settling time to be
Ts = 4 / a
Consider a simple first-order system, G(s) = K /(s + a), whose step response is
If we can identify K and a from laboratory testing, we can obtain the transfer
function of the system.
Example: Assume the unit step response given by c(t) = cf(t) + cn(t) = 1 – e-at.
From the response, we measure the time constant, that is, the time for the
amplitude to reach 63% of its final value. Since the final value is about 0.72, the
time constant is evaluated where the curve reaches 0.63 x 0.72 = 0.45, or about
0.13 second.
Hence, a = 1/0.13 = 7.7.
The transfer function G(s) = 5/(s + 7) was used to generate the following response.
Overdamped Response
This function has a pole at the origin that comes from the unit step input and two
real poles that come from the system. The input pole at the origin generates the
constant forced response; each of the two system poles on the real axis generates
an exponential natural response whose exponential frequency is equal to the pole
location.
Hence, the output initially could have been written as
This response is called overdamped. We see that the poles tell us the form of the
response without the tedious calculation of the inverse Laplace transform.
Underdamped Response
This function has a pole at the origin that comes from the unit step input and two
complex poles that come from the system. We now compare the response of the
second-order system to the poles that generated it.
1. We will compare the pole location to the time function, and then we will compare
the pole location to the plot.
The poles that generate the natural response are at s = -1 ± j√8
Comparing these values to c(t) in the same figure, we see that the real part of the
pole matches the exponential decay frequency of the sinusoid’s amplitude, while
the imaginary part of the pole matches the frequency of the sinusoidal oscillation.
3. The steady-state response (unit step) was generated by the input pole located at
the origin. The type of response shown above is called an underdamped
response, one which approaches a steady-state value via a transient response
that is a damped oscillation.
Example: By inspection, write the form of the step response of the system below
Undamped Response
This function has a pole at the origin that comes from the unit step input and two
imaginary poles that come from the system. The input pole at the origin generates
the constant forced response, and the two system poles on the imaginary axis at ±j3
generate a sinusoidal natural response whose frequency is equal to the location of
the imaginary poles.
Hence, the output can be estimated as c(t) = K1 + K4 cos(3t – ϕ).
This type of response is called undamped. Note that the absence of a real part in
the pole pair corresponds to an exponential that does not decay. Mathematically, the
exponential is e-0t = 1.
This function has a pole at the origin that comes from the unit step input and two
multiple real poles that come from the system. The input pole at the origin
generates the constant forced response, and the two poles on the real axis at -3
generate a natural response consisting of an exponential and an exponential
multiplied by time, where the exponential frequency is equal to the location of the
real poles.
This type of response is called critically damped. Critically damped responses are
the fastest possible without the overshoot that is characteristic of the underdamped
response.
1. Overdamped Responses
Poles: Two real at -σ1, -σ2
Natural Response: Two exponentials with time constants equal to the reciprocal of
the pole locations, or c(t) = K1e-σ1t + K2e-σ2t
2. Underdamped Responses
Poles: Two complex at -σd ± jωd
Natural Response: Damped sinusoid with an exponential envelope whose time
constant is equal to the reciprocal of the pole’s real part. The radian frequency of
the sinusoid, the damped frequency of oscillation, is equal to the imaginary part
of the poles, or c(t) = Ae-σd*t cos(ωdt – ϕ)
3. Undamped Responses
Poles: Two imaginary at ±jω1
Natural Response: Undamped sinusoid with radian frequency equal to the
imaginary part of the poles, or c(t) = Acos(ω1t – ϕ)
author w chiyabu Page 74
4. Critically Damped Responses
Poles: Two real at -σ1
Natural Response: One term is an exponential whose time constant is equal to the
reciprocal of the pole location. Another term is the product of time, t, and an
exponential with time constant equal to the reciprocal of the pole location, or
c(t) = K1e-σ1t + K2te-σ1t
Exercise
For each of the following transfer functions, write, by inspection, the general form of
the step response:
questions answers
Natural Frequency, ωn
The natural frequency of a second-order system is the frequency of oscillation of the
system without damping. For example, the frequency of oscillation of a series RLC
circuit with the resistance shorted would be the natural frequency.
Damping Ratio, ζ
A second-order system’s underdamped step response is characterised by damped
oscillations. Our definition is derived from the need to quantitatively describe this
damped oscillation regardless of the time scale. Thus, a system whose transient
response goes through three cycles in a millisecond before reaching the steady state
would have the same measure as a system that went through three cycles in a
millennium before reaching the steady state. For example, the underdamped curve
has an associated measure that defines its shape. This measure remains the same
even if we change the time base from seconds to microseconds or to millennia.
A viable definition for this quantity is one that compares the exponential decay
frequency of the envelope to the natural frequency. This ratio is constant regardless
of the time scale of the response. Also, the reciprocal, which is proportional to the
ratio of the natural period to the exponential time constant, remains the same
regardless of the time base.
The general second-order system can be transformed to show the quantities ζ and
ωn.
Without damping, the poles would be on the jω-axis, and the response would be an
undamped sinusoid. For the poles to be purely imaginary, a = 0.
By definition, the natural frequency, ωn, is the frequency of oscillation of this system.
Since the poles of this system are on the jω-axis at ±j√b
ωn = √b or b = ω2n
Assuming an underdamped system, the complex poles have a real part, σ = -a/2.
The magnitude of this value is then the exponential decay frequency.
Hence,
Solution:
Comparing to
From
Solution:
From and
Exercise
Answers:
a. damping ratio, ζ = 0.3, natural frequency ωn = 20 system is underdamped
b. damping ratio, ζ = 1.5, natural frequency ωn = 30 system is overdamped
c. damping ratio, ζ = 1, natural frequency ωn = 15 system is critically damped
d. damping ratio, ζ = 0, natural frequency ωn = 25 system is undamped
Considering the step response for the above general second-order system, the
transform of the response, C(s), is the transform of the input times the transfer
function, or
A plot of this response appears below for various values of ζ, plotted along a time
axis normalised to the natural frequency.
Definitions
These definitions are valid for systems of order higher than 2, although analytical
expressions for these parameters cannot be found unless the response of the
higher-order system can be approximated as a second-order system.
Rise time, peak time, and settling time yield information about the speed of the
transient response. This information can help a designer determine if the speed and
the nature of the response do or do not degrade the performance of the system.
we get
or
Each value of n yields the time for local maxima or minima. Letting n = 0 yields t =
0, the first point on the curve that has zero slope. The first peak, which occurs at the
peak time, Tp, is found by letting n =1:
The term cmax is found by evaluating c(t) at the peak time, c(Tp).
and
we get
reaches and stays within ±2% of the steady-state value, cfinal. Settling Time is the
time it takes for the amplitude of the decaying sinusoid in c(t) to reach 0.02, or
Approximation for the Settling Time that will be used for all values of ζ
the rise time can be found. We first designate ωnt as the normalised time variable
and select a value for ζ. Using the computer, we solve for the values of ωnt that yield
c(t) = 0.9 and c(t) = 0.1. Subtracting the two values of ωnt yields the normalised
rise time, ωnTr, for that value of ζ.
Continuing in like fashion with other values of ζ, we obtain the results plotted in
below.
Solution:
ωn and ζ are calculated as 10 and 0.75, respectively. Now substitute ζ and ωn into
, , and
and find, respectively, that Tp = 0.475 second, %OS = 2.838, and Ts = 0.533
second. Using the table above, the normalised rise time is approximately 2.3
seconds. Dividing by ωn yields Tr = 0.23 second. This problem demonstrates that we
can find Tp, %OS, Ts, and Tr without the tedious task of taking an inverse Laplace
transform, plotting the output response, and taking measurements from the plot.
The pole plot for a general, underdamped second-order system is shown below.
From the Pythagorean Theorem, the radial distance from the origin to the pole is the
natural frequency, ωn, and the cosθ = ζ.
Comparing
and
with the pole location, we evaluate peak time and settling time in terms of the pole
location.
where ωd is the imaginary part of the pole and is called the damped frequency of
oscillation, and σd is the magnitude of the real part of the pole and is the
exponential damping frequency.
3. Since Damping Ratio, ζ = cosθ, radial lines are lines of constant damping ratio ζ.
Since percent overshoot is only a function of damping ratio ζ, radial lines are thus
lines of constant percent overshoot, %OS.
These concepts are depicted in the figure below, where lines of constant Peak
Time Tp, Settling Time Ts, and Percent Overshoot %OS are labelled on the s-
plane.
Lines of constant peak time, Tp, settling time, Ts, and percent overshoot, %OS.
Note: Ts2 < Ts1 ; Tp2 < Tp1; %OS1 < %OS2
The significance of the figure above can be better understood by examining the
actual step response of comparative systems.
Step responses of second-order underdamped systems as poles move with constant real part
As the poles move in a vertical direction, the frequency increases, but the envelope
remains the same since the real part of the pole is not changing. The figure shows a
constant exponential envelope, even though the sinusoidal response is changing
frequency. Since all curves fit under the same exponential decay curve, the settling
time is virtually the same for all waveforms. Note that as overshoot increases, the
rise time decreases.
Step responses of second-order underdamped systems as poles move with constant imaginary part
Let us move the poles to the right or left. Since the imaginary part is now constant,
movement of the poles yields constant frequency over the range of variation of the
real part. As the poles move to the left, the response damps out more rapidly, while
the frequency remains the same.
Notice that the peak time is the same for all waveforms because the imaginary part
remains the same.
Step responses of second-order underdamped systems as poles move with constant damping ratio
Moving the poles along a constant radial line yields does not change the percent
overshoot. Notice also that the responses look exactly alike, except for their speed.
The farther the poles are from the origin, the more rapid the response.
Example: Given the pole plot shown below, find ζ; ωn; Tp;%OS, and Ts
Solution:
The damping ratio is given by ζ = cosθ = cos[arctan(7/3)] = 0.394
The natural frequency, ωn, is the radial distance from the origin to the pole, or
ωn = √(72 + 32) = 7.616
The peak time is Tp = (π / ωd) = (π/7) = 0.449 s
The percent overshoot is
Example: Given the system shown below, find J and D to yield 20% overshoot and a
settling time of 2 seconds for a step input of torque T(t)
Solution:
First, the transfer function for the system is
Exercise: Find ζ; ωn; Ts; Tp; Tr, and %OS for a system whose transfer function is
Solution:
ζ = 0.421; ωn = 19; Ts = 0.5 s; Tp = 0.182 s; Tr = 0.079 s; and %OS = 23.3%: