ELE 301: Signals and Systems: Prof. Paul Cuff
ELE 301: Signals and Systems: Prof. Paul Cuff
Fall 2011-12
Cuff (Lecture 3)
Fall 2011-12
1 / 55
Todays topics
Impulse response
Extended linearity
Response of a linear time-invariant (LTI) system
Convolution
Zero-input and zero-state responses of a system
Cuff (Lecture 3)
Fall 2011-12
2 / 55
Impulse Response
The impulse response of a linear system h (t) is the output of the system
at time t to an impulse at time . This can be written as
h = H( )
Care is required in interpreting this expression!
(t)
0
h(t, 0)
t
(t )
0
Cuff (Lecture 3)
H
h(t, )
Fall 2011-12
3 / 55
Cuff (Lecture 3)
Fall 2011-12
4 / 55
h(2) = H((2))
First, (2) is something like zero, so H(0) would be zero. Second, the
value of h(2) depends on the entire input waveform, not just the
value at t = 2.
H
(t)
0
(2)
2
Cuff (Lecture 3)
h(t, 0)
h(2, 0)
Fall 2011-12
5 / 55
Time-invariance
If H is time invariant, delaying the input and output both by a time
should produce the same response
h (t) = h(t ).
In this case, we dont need to worry about h because it is just h shifted in
time.
(t)
0
h(t)
t
(t )
0
Cuff (Lecture 3)
H
h(t )
0
ELE 301: Signals and Systems
t
Fall 2011-12
6 / 55
Cuff (Lecture 3)
Fall 2011-12
7 / 55
Extended Linearity
Summation: If yn = Sxn for all n, an integer from ( < n < ),
and an are constants
!
X
X
an yn = S
an x n
n
Cuff (Lecture 3)
Fall 2011-12
8 / 55
x
H
Cuff (Lecture 3)
Fall 2011-12
9 / 55
If the system obeys extended linearity we can interchange the order of the
system operator and the integration
Z
y (t) =
x( )H ( (t)) d.
Cuff (Lecture 3)
Fall 2011-12
10 / 55
Cuff (Lecture 3)
Fall 2011-12
11 / 55
Output
h(t)
(t)
0
(t )
0
x(t)
h(t )
t
(x()d)(t )
t
(x()d)h(t )
t
0
y(t)
x(t)
0
x(t) =
Cuff (Lecture 3)
t
x()(t )d
0
y(t) =
x()h(t )d
Fall 2011-12
12 / 55
System Equation
The System Equation relates the outputs of a system to its inputs.
Example from last time: the system described by the block diagram
+
-
Fall 2011-12
13 / 55
Solving the system equation tells us the output for a given input.
The output consists of two components:
The zero-input response, which is what the system does with no input
at all. This is due to initial conditions, such as energy stored in
capacitors and inductors.
x(t) = 0
0
Cuff (Lecture 3)
y(t)
0
Fall 2011-12
14 / 55
The zero-state response, which is the output of the system with all
initial conditions zero.
x(t)
0
y(t)
0
x(t) = 0
0
Cuff (Lecture 3)
y(t) = 0
Fall 2011-12
15 / 55
Example: Solve for the voltage across the capacitor y (t) for an arbitrary
input voltage x(t), given an initial value y (0) = Y0 .
i(t) R
x(t) +
Cuff (Lecture 3)
y(t)
Fall 2011-12
16 / 55
Cy 0 (t)
RCy 0 (t) + y (t) = x(t).
This is a first order LCCODE, which is linear with zero initial conditions.
First we solve for the homogeneous solution by setting the right side (the
input) to zero
RCy 0 (t) + y (t) = 0.
The solution to this is
y (t) = Ae t/RC
which can be verified by direct substitution.
Cuff (Lecture 3)
Fall 2011-12
17 / 55
1 t/RC
e
RC
Cuff (Lecture 3)
Fall 2011-12
18 / 55
Then
y (t) = A(t)e t/RC
Z t
1 /RC
= e t/RC
x( )
e
d + A(0)e t/RC
RC
0
Z t
1 (t )/RC
e
d + A(0)e t/RC
=
x( )
RC
0
At t = 0, y (0) = Y0 , so this gives A(0) = Y0
Z t
1 (t )/RC
y (t) =
x( )
e
d +
Y e t/RC
.
| 0 {z }
RC
|0
{z
} zeroinput response
zerostate response
Cuff (Lecture 3)
Fall 2011-12
19 / 55
Fall 2011-12
20 / 55
RC Circuit example
The impulse response of the RC circuit example is
h(t) =
1 t/RC
e
RC
Cuff (Lecture 3)
Example:
High energy photon detectors can be modeled as having a simple
exponential decay impulse response.
Doshi et al.: LSO PET detector
LE
1540
ts.
Coupler
ect LSOa
ghtguidea
V lens
era
er taper
Photomultiplier
Scintillating
Crystal
Average
Number of
Light
Energy
peak-tocrystals clearly
resolution
collection
resolved
!FWHN %# efficiency !%# valley ratio
13.0
19.9
27.2
35.0
19.5
Light 100.0
40.6
28.0
12.6
27.0
10.0
2.5
2.5
6.0
7.5
Light Fibers
9
7
7
6
9
Crystal
ergy resolution and light collection efficiency were measured with single
htguide elements.
Photon
METHODSDETECTOR CHARACTERIZATION
21 / 55
Energy spectra
Input: Photons
Timing resolution
E. Detector efficiency
A measure of the absolute detector efficiency was obtained. A 18F point source with known activity !68 " Ci# was
placed 15 cm away from the face of the detector module. The
actual gamma-ray flux impinging on the detector face was
calculated from the solid angle subtended by the detector
module at the source. The constant fraction discriminator
was set to eliminate electronic noise ($100 keV# and the full
energy spectrum was obtained for each crystal over a fixed
time. A background measurement without the 18F point
source was also obtained to subtract the LSO background
from the measurement. A lower energy window of 350 keV
was applied to all of the crystals and the number of counts
falling under the photopeak was calculated. The number of
counts detected was then divided by the total number of
gamma rays impinging on the detector module to obtain the
detector efficiency.
Output: Light
An image of the flood histogram from one detector module is shown in Fig. 6. All 81 crystals from the 9!9 LSO
array are clearly visible. An average peak-to-valley ratio of
3.5 was obtained over the central row of nine crystals. Not
all crystals are uniformly spaced in the flood histogram after
applying Anger logic. This may be a result of the nonuniform tapering of the optical fiber taper, the nonuniform packing of the reflectance powder between the crystals, or most
Fall 2011-12
22 / 55
Summary
Cuff (Lecture 3)
Fall 2011-12
23 / 55
=
=
Cuff (Lecture 3)
h(t)
y(t)
x(t)
h(t)
y(t)
Fall 2011-12
24 / 55
Since h (t) = 0 for t < , we can replace the upper limit of the integral
by t
Z t
y (t) =
x( )h (t) d.
Cuff (Lecture 3)
Fall 2011-12
25 / 55
h( ) [A cos(2f1 (t ) + )] d.
Z
h( ) sin(2f1 )d.
+A sin(2f1 t + )
Cuff (Lecture 3)
Fall 2011-12
26 / 55
h( ) cos(2f1 )d
h( ) sin(2f1 )d
Cuff (Lecture 3)
Fall 2011-12
27 / 55
Cuff (Lecture 3)
Fall 2011-12
28 / 55
Summary
Cuff (Lecture 3)
Fall 2011-12
29 / 55
Cuff (Lecture 3)
Fall 2011-12
30 / 55
Cuff (Lecture 3)
Fall 2011-12
31 / 55
Convolution Integral
The convolution of an input signal x(t) with and impulse response h(t) is
Z
y (t) =
x( )h(t ) d
= (x h)(t)
or
y = x h.
This is also often written as
y (t) = x(t) h(t)
which is potentially confusing, since the ts have different interpretations
on the left and right sides of the equation (your book does this).
Cuff (Lecture 3)
Fall 2011-12
32 / 55
Cuff (Lecture 3)
Fall 2011-12
>t
<t
Does not
contribute to y(t)
y(t)
x(t)
0
x(t) =
y(t) =
x()(t )d
33 / 55
x( )h(t )d
If x(t) is also causal, x(t) = 0 for t < 0, and the integral further simplifies
Z t
y (t) =
x( )h(t ) d.
0
Does not
contribute to y(t)
>t
<t
x(t)
0
x(t) =
y(t)
t
Cuff (Lecture 3)
x( )(t )d
y(t) =
Does not
contribute to y(t)
t
t
x( )h(t )d
Fall 2011-12
34 / 55
Graphical Interpretation
An increment in input x( ) (t)d produces an impulse response
x( )h (t)d . The output is the integral of all of these responses
Z
x( )h (t) d
y (t) =
h(t )
Cuff (Lecture 3)
h(t )
t
Fall 2011-12
35 / 55
x( )
h(t )
x( )h(t )
Cuff (Lecture 3)
Fall 2011-12
36 / 55
Simple Example
2
x()
1
0
-1
h()
1
0
-1
h()
2
1
0
-1
-1
Cuff (Lecture 3)
-1
-1
x()
t <0
1
0
h(t )
1
0
Cuff (Lecture 3)
x()
-1
1<t <2
h(t )
h(t )
Fall 2011-12
1
-1
0<t <1
1
x()
2<t <3
h(t )
1
y(t) = (x h)(t)
t >3
37 / 55
x()
1
x()
-1
h(t )
h(t )
1
-1
Fall 2011-12
38 / 55
y(t)
h(t)
Impulse response:
1.5
h(t)
0.5
0
0
tt
10
Fall 2011-12
39 / 55
x(t)
0.5
0
0
tt
10
y(t)
0.5
0
0
tt
10
Fall 2011-12
40 / 55
x(t)
1
u
0.5
0
0
y(t)
1
y
10
0.5
0
0
tt
10
Fall 2011-12
41 / 55
1
0
1
0
1
1
2
(t 1)
Cuff (Lecture 3)
1
0
0
1
1
0
(x h)(t)
h(t)
1
0
Fall 2011-12
42 / 55
Properties of Convolution
For any two functions f and g the convolution is
Z
(f g )(t) =
f ( )g (t ) d
= (g f )(t)
Fall 2011-12
43 / 55
x()
1
0
-1
xh
h(t )
-1
h()
hx
x()
1
0
-1
x(t )
-1
h()
1
y(t) = (x h)(t)
2
1
-1
Cuff (Lecture 3)
Fall 2011-12
44 / 55
Convolution is associative
Cuff (Lecture 3)
Fall 2011-12
45 / 55
Linearity
Convolution is also distributive,
f (g + h) = f g + f h
which is easily shown by writing out the convolution integral,
Z
f ( ) [g (t ) + h(t )] d
(f (g + h))(t) =
Z
Z
f ( )g (t ) d +
f ( )h(t ) d
=
= (f g )(t) + (f h)(t)
Fall 2011-12
46 / 55
Time-invariant
Cuff (Lecture 3)
Fall 2011-12
47 / 55
Fall 2011-12
48 / 55
Cuff (Lecture 3)
y
( f g)
Fall 2011-12
49 / 55
f
g
g
f
Cuff (Lecture 3)
Fall 2011-12
50 / 55
h(t)
t
Cuff (Lecture 3)
Fall 2011-12
51 / 55
h(t)
u(t)
d
dt
u(t)
0
d
dt
Cuff (Lecture 3)
h(t)
t
Fall 2011-12
52 / 55
Z
= e st
h( )e s d
y (t) = e st H(s)
Fall 2011-12
53 / 55
Fall 2011-12
54 / 55
y (t) = e jt H(j)
Cuff (Lecture 3)
Summary
Cuff (Lecture 3)
Fall 2011-12
55 / 55