Elements of Feedback Control Systems
Elements of Feedback Control Systems
Chapter 2
Elements of feedback control systems
Signals & Systems
Linearity
Laplace transform
21
22
greeting & talk: hello, how are your?, where are you from? have you eaten?
sleep in a room
D. S. Bernstein, Introducing signals, systems, and control in grades K through 12, IEEE Control Systems Magazine, April 2003, pp. 10-12.
23
Signals
a signal is a function of time, e.g.,
f is the force on some mass
24
Domain of a signal
domain of a signal domain of a signal: ts for which it is defined
some common domains:
all t, i.e., R
nonnegative t: t > 0
(here t = 0 just means some starting time of interest)
t in some interval: a t b
well usually study signals defined on all reals, or for nonnegative reals
25
26
0, t < 0
1, t 0
0, t < 0
t, t 0
r(t) =
a sinusoidal signal:
r(t) =
u(t) =
1, a t b
0, otherwise
u(t) = a cos(t + )
a, , are called signal parameters
27
u(t)2 dt
|u(t)| dt
max |u(t)|
t
lim
T T
average-absolute value:
T
0
1/2
u(t)2 dt
1 ZT
|u(t)| dt
lim
T T 0
for some signals these measures can be infinite, or undefined
28
29
Systems
definition:
a system transforms input signals into output signals
we concentrate on systems with one input and one output signal, i.e., single-input,
single-output (SISO) systems
2 10
Block diagrams
systems often denoted by block diagram:
u
G(s)
2 11
d
dt
2 12
integrator: y(t) =
t
a u( ) d
(a is often 0 or -)
y(t) = sgn(u(t)) =
u
1, u(t) 0
1, u(t) < 0
y
2 13
u1
u2
difference system: y(t) = u1(t) u2(t)
u1
u2
2 14
u1
u2
comparator system:
y(t) =
1, u1(t) u2(t)
1, u1(t) < u2(t)
u1
u2
2 15
Electronic realizations
systems described previously can be realized as electronic circuits with op-amps
scaling: y(t) = (1 + R2/R1 )u(t)
+
u
R1
y
R2
R
+
u1
R
+u
2
y
R
2 16
t
0 u( ) d
R
+
u
2 17
Linearity
a system F is linear if the following properties hold.
homogeneity: if u is any signal and a is any scalar,
F (au) = aF (u)
superposition: if u and u are any two signals,
F (u + u) = F u + F u
(watch out, just a few symbols here express a very complex meaning)
in words, linearity means scaling before or after the system is the same
y
y
u
u
a
a
F
F
summing before or after the system is the same
u1
u1
y
u2
F
u2
F
F
running average
difference
examples of nonlinear systems
sign detector
multiplier
comparator
2 18
2 19
Interconnections of systems
we can interconnect systems to form new systems, e.g.,
cascade or series: y = G(F u) = GF u
u
F
G
2 20
feedback: y = F (u Gy)
u
G
the minus sign is just a tradition, and often isnt there
in general, block diagrams are just a symbolic way to describe a connection of systems
2 21
when initial conditions are all zero, LCCODE systems are linear
an LCCODE gives an implicit description of a system; soon well be able to explicitly
express y in terms of u
example: first- and second-order system
dy
+ 2y = u
dt
dy
d2 y
+2 +y = u
dt2
dt
2 22
Mass-spring-damper system
it represents vehicle suspension system, building during earthquake, . . .
simple model of vehicle suspension system
y
m
k
b
u
my + by + ky = bu + ku
2 23
Laplace Transform
Idea: the Laplace transform converts integral and differential equations into algebraic
equations
2 24
F (s) = 0 f (t)est dt
for those s C for which the integral makes sense
s is called the complex variable or frequency variable with units sec1, i.e.,
frequency
st is unitless
2 25
Laplace transform of et
lets find Laplace transform of f (t) = et
F (s) =
=
t st
dt
0 ee
Z
(1s)t
dt
0 e
1 (1s)t
=
e
1s
0
1
=
s1
provided we can say e(1s)t 0 as t , which is true for Re s > 1 (why?)
that is to say,
the integral defining F makes sense for s C with Re s > 1
but the resulting formula for F makes sense for all s C except s = 1
in practice, these important details are not mentioned, and we simply say that the
Laplace transform of et is 1/(s 1), or
1
t
(e
)
=
L
s1
2 26
More examples
impulse function: f (t) = (t)
F (s) =
st
(t)e
dt
=1
1 st 1
F (s) =
= e =
s
s
0
integral makes sense for Re s > 0; formula makes sense for s 6= 0
st
dt
0 e
sinusoid function
first express f (t) = cos(t) as
1 jt 1 jt
f (t) = e + e
2
2
2 27
1 1 1 1
=
+
2 s j
2 s + j
s
= 2
s + 2
F (s) is valid for Re s > 0; but formula makes sense for s 6= j
F (s) =
F (s) = 2
s + 2
2 28
L(f + g) = F + G
3s
5
t
= 3(t) 2e
L1
s1
2 29
Re
note that the integration of F (s)est is calculated along the dash line in the splane
surprisingly, this formula is not widely used, except when we do numerical inverse of
Laplace transform
2 30
Time scaling
let f be a signal, define signal g by
g(t) = f (at)
where a > 0, then
G(s) = (1/a)F (s/a)
this expression makes sense because times are scaled by a, frequencies scaled by 1/a
1
1
=
(s/a) 1 s a
time-scaling formula requires a > 0, but this is true for all values of a
2 31
Derivative
if signal f (t) is differentiable for t 0, then
1
=
(f
)
=
s
L
s1
s1
1
s
1 = 2
L(sin t) = s 2
s +1
s +1
2 32
2 33
Integral
let g be the running integral of a signal f , i.e.,
g(t) =
t
0 f ( ) d
then
1
G(s) = F (s)
s
i.e., time-domain integral becomes division by frequency variable s
example: we know that the Laplace transform of the unit delta function
L((t)) = 1
and the integral of the unit delta function is the unit step function, u(t), hence,
L(u(t)) =
1
s
example: the integral of the unit step function is the unit ramp function, t,
1
(t)
=
L
s2
2 34
Time delay
let f be a signal and T > 0
define the signal g as
g(t) =
g is f , delayed by T seconds
0,
0t<T
f (t T ), t T
G(s) = 0 estg(t) dt
Z
= T estf (t T ) dt
Z
= 0 es( +T )f ( ) d
= esT F (s)
Z
where = t T
f
sT
2 35
f (t) =
1, a t b
0, otherwise
b
u(t b)
F (s) = e
bs
(1/s) e
eas ebs
(1/s) =
s
2 36
Multiplication by t
let f be a signal and define g = tf (t), then we have
G(s) = F (s)
i.e., multiplication by time-domain variable corresponds to differentiation in the
frequency domain
c.f., L(f ) = sF (s) f (0)
example: L(et) = 1/(s + 1)
L(tet) =
1
d 1
=
ds s + 1 (s + 1)2
k
k (k)
L(t f (t)) = (1) F (s)
k t
L(t e ) =
(k 1)!
(s + 1)k+1
2 37
Exponential scaling
let f be a signal and a a scalar, and define
g(t) = eatf (t),
then
st at
L(g(t)) = 0 e e f (t) dt
Z
= 0 e(sa)tf (t) dt
= F (s a)
Z
i.e., exponential scaling in the time domain corresponds to a shift (by a) in the
frequency domain
example: L(cos t) = s/(s2 + 1),
(s + 1)
s+1
=
L(e cos t) =
(s + 1)2 + 1 s2 + 2s + 2
t
2 38
Convolution
the convolution of signals f and g, denoted h = f g, is the signal
h(t) =
same as h(t) =
t
0 f (t
t
0 f ( )g(t
) d
)g( ) d
symmetric: f g = g f
2 39
L(cos(t)) =
s
,
2
2
s +
L(sin(t)) =
s2 + 2
2 40
Patterns
while the details of the Laplace transform pairs differ, we can see some interesting
symmetric patterns between
the time domain (i.e., signals) and
2 41
1/s
1/s2
1/(s + a)
/(s2 + 2)
s/(s2 + 2)
2 42
Convolution systems
consider convolution system with input u, output y satisfying
y = hu
Z
t
= 0 h( )u(t ) d
Z
t
= 0 h(t )u( ) d
in the frequency domain:
Y (s) = H(s)U (s)
H(s) is called the transfer function (TF) of the system
y
h
y
H
2 43
Properties
composition of convolution systems corresponds to
multiplication of transfer functions
convolution of impulse responses
u
A
u
BA
other properties:
we can manipulate block diagrams with transfer functions as if they were simple
gains
convolution systems commute with each other
2 44
Examples
first order LCCODE:
y + y = u, y(0) = 0
take Laplace transform to get
sY (s) y(0) + Y (s) = U (s)
after simple algebra, we obtain
1
U (s)
s+1
transfer function H(s) = 1/(s + 1); impulse response h(t) = et
delay: with T > 0,
0,
0t<T
y(t) =
u(t T ), t T
Y (s) =
transfer function H(s) = 1/s; impulse response h(t) = u(t) (unit step function)
2 45
u
esT
F
y
u
F
esT
2 46
Feedback connection
what is transfer function H from u to y?
u
G(s)
Y (s) =
U (s)
1 + G(s)
as if G were a simple scaling system
so we have
G(s)
H(s) =
1 + G(s)
in time domain, we have complicated integral equation
y(t) =
t
0 g(t
)(u( ) y( )) d
2 47
Mass-spring-damper system
it represents vehicle suspension system, building during earthquake, . . .
simple model of vehicle suspension system
y
m
k
b
u
my + by + ky = bu + ku
2 48
Impulse response
if u = , we have
y(t) =
t
0
h(t )( ) d = h(t)
so h is the output response when u = (hence h has the name impulse response)
u
y
H
practical problem: linear model often fails for very large input signals
2 49
Step response
the unit step response is the output when the input is a unit step, i.e., u(t) = 1
s(t) =
t
0 h( ) d
2 50
example: H(s) =
1
s+1
h(t) = et,
s(t) = 1 et
1
Step Response
Impulse Response
0.5
0
0
2
4
Time (sec)
0.5
0
0
2
4
Time (sec)
2 51
DC gain
if a convolution system is stable and
H(0) =
0 h( ) d
t
0 h( ) d
H(0)u
so H(0) gives the gain for static signals in terms of step response
Final value theorem:
given F (s) = L(f (t)), if limt f (t) exists, then
s0
example: H(s) =
1
s+1
s0
DC gain is equal to 1
2 52
Stability
we say that the convolution system is stable if impulse response h(t) 0 as t
if H(s) = L(h(t)) is rational, stability means poles are in the left-half s-plane (LHP)
example:
h(t) = e2t,
H(s) =
1
s+2
Impulse Response
1
0.9
0.8
0.7
Amplitude
0.6
0.5
0.4
0.3
0.2
0.1
0.5
1.5
Time (sec)
2.5
2 53
Summary
signals convey information. systems transform input signals into output signals
the Laplace transform converts integral and differential equations into algebraic
equations
a system is linear time-invariant (LTI) if it is linear (it satisfies homogeneity &
superposition) and time-invariant (it commutes with delays)
the impulse response is the output of the system when input is a unit delta
the step response is the output of the system when the input is a unit step