Motor
Motor
Motor
5. The unit-impulse response can be used to fully describe a linear, time-invariant system. Given the
impulse response, the system output is the sum of the products of the impulse response coefficients
and corresponding input samples, called the digital convolution sum.
6. BIBO is a type of stability in which a bounded input will produce a bounded output. A BIBO
system requires that the sum of the absolute impulse response coefficients be a finite number.
7. The digital convolution sum, which represents a DSP system, is evaluated in three ways: the
graphical method, evaluation of the formula, and the table method. The table method is found to
be most effective.
3.7 PROBLEMS
3.1. Sketch each of the following special digital sequences:
a. 5dðnÞ
b. 2dðn 5Þ
c. 5uðnÞ
d. 5uðn 2Þ
3.2. Calculate the first eight sample values and sketch each of the following sequences:
a. xðnÞ ¼ 0:5n uðnÞ
b. xðnÞ ¼ 5sinð0:2pnÞuðnÞ
c. xðnÞ ¼ 5cosð0:1pn þ 300 ÞuðnÞ
d. xðnÞ ¼ 5ð0:75Þn sinð0:1pnÞuðnÞ
3.3. Sketch each of the following special digital sequences:
a. 8dðnÞ
b. 3:5dðn 4Þ
c. 4:5uðnÞ
d. 6uðn 3Þ
3.4. Calculate the first eight sample values and sketch each of the following sequences:
a. xðnÞ ¼ 0:25n uðnÞ
b. xðnÞ ¼ 3sinð0:4pnÞuðnÞ
c. xðnÞ ¼ 6cosð0:2pn þ 300 ÞuðnÞ
d. xðnÞ ¼ 4ð0:5Þn sinð0:1pnÞuðnÞ
3.5. Sketch the following sequences:
a. xðnÞ ¼ 3dðn þ 2Þ 0:5dðnÞ þ 5dðn 1Þ 4dðn 5Þ
b. xðnÞ ¼ dðn þ 1Þ 2dðn 1Þ þ 5uðn 4Þ
3.6 Summary 81
3.6. Given the digital signals xðnÞ in Figures 3.24 and 3.25, write an expression for each digital
signal using the unit-impulse sequence and its shifted sequences.
x ( n)
4
3 2
2 1 1
1
1
n
1 0 1 2 3 4 5
FIGURE 3.24
The first digital signal in Problem 3.6.
x ( n)
n
1 0 1 2 3 4 5
1
FIGURE 3.25
The second digital signal in Problem 3.6.
x ( n)
4
3
2
1
n
1 0 1 2 3 4 5
FIGURE 3.26
The first digital signal in Problem 3.8.
82 CHAPTER 3 Digital Signals and Systems
x(n)
n
1 0 1 2 3 4 5
1
FIGURE 3.27
The second digital signal in Problem 3.8.
3.9. Assume that a DS processor with a sampling time interval of 0.01 second converts the
following analog signals xðtÞ to a digital signal xðnÞ; determine the digital sequence for each
of the analog signals.
a. xðtÞ ¼ e50t uðtÞ
b. xðtÞ ¼ 5sinð20ptÞuðtÞ
c. xðtÞ ¼ 10cosð40pt þ 300 ÞuðtÞ
d. xðtÞ ¼ 10e100t sinð15ptÞuðtÞ
3.10. Determine which of the following systems is a linear system.
a. yðnÞ ¼ 5xðnÞ þ 2x2 ðnÞ
b. yðnÞ ¼ xðn 1Þ þ 4xðnÞ
c. yðnÞ ¼ 4x3 ðn 1Þ 2xðnÞ
3.11. Assume that a DS processor with a sampling time interval of 0.005 second converts each of
the following analog signals xðtÞ to a digital signal xðnÞ; determine the digital sequence for
each of the analog signals.
a. xðtÞ ¼ e100t uðtÞ
b. xðtÞ ¼ 4sinð60ptÞuðtÞ
c. xðtÞ ¼ 7:5cosð20pt þ 600 ÞuðtÞ
d. xðtÞ ¼ 20e200t sinð60ptÞuðtÞ
3.12. Determine which of the following systems is a linear system.
a. yðnÞ ¼ 4xðnÞ þ 8x3 ðnÞ
b. yðnÞ ¼ xðn 3Þ þ 3xðnÞ
c. yðnÞ ¼ 5x2 ðn 1Þ 3xðnÞ
3.13. Determine which of the following linear systems is time invariant.
a. yðnÞ ¼ 5xðn 10Þ
b. yðnÞ ¼ 4xðn2 Þ
3.6 Summary 83
3.24. Determine the stability for each of the following linear systems.
PN k
a. yðnÞ ¼ k ¼ 0 ð1:5Þ xðn kÞ
PN k
b. yðnÞ ¼ k ¼ 0 ð0:5Þ xðn kÞ
X
N
yðnÞ ¼ xðkÞhðn kÞ
k ¼ N
N
ak y(n − k) = x(n) (2.26)
k=0
2.5 PROBLEMS
P2.1 Generate the following sequences using the basic MATLAB signal functions and the basic
MATLAB signal operations discussed in this chapter. Plot signal samples using the stem
function.
1. x1 (n) = 3δ(n + 2) + 2δ(n) − δ(n − 3) + 5δ(n − 7), −5 ≤ n ≤ 15.
5
2. x2 (n) = k=−5
e−|k| δ(n − 2k), −10 ≤ n ≤ 10.
3. x3 (n) = 10u(n) − 5u(n − 5) − 10u(n − 10) + 5u(n − 15).
4. x4 (n) = e0.1n [u(n + 20) − u(n − 10)].
5. x5 (n) = 5[cos(0.49πn) + cos(0.51πn)], −200 ≤ n ≤ 200. Comment on the waveform
shape.
6. x6 (n) = 2 sin(0.01πn) cos(0.5πn), −200 ≤ n ≤ 200. Comment on the waveform shape.
7. x7 (n) = e−0.05n sin(0.1πn + π/3), 0 ≤ n ≤ 100. Comment on the waveform shape.
8. x8 (n) = e0.01n sin(0.1πn), 0 ≤ n ≤ 100. Comment on the waveform shape.
P2.2 Generate the following random sequences and obtain their histogram using the hist
function with 100 bins. Use the bar function to plot each histogram.
1. x1 (n) is a random sequence whose samples are independent and uniformly distributed
over [0, 2] interval. Generate 100,000 samples.
2. x2 (n) is a Gaussian random sequence whose samples are independent with mean 10 and
variance 10. Generate 10,000 samples.
3. x3 (n) = x1 (n) + x1 (n − 1) where x1 (n) is the random sequence given in part 1 above.
Comment
on the shape of this histogram and explain the shape.
4
4. x4 (n) = k=1 yk (n) where each random sequence yk (n) is independent of others with
samples uniformly distributed over [−0.5, 0.5]. Comment on the shape of this histogram.
Copyright 2010 Cengage Learning. All Rights Reserved. May not be copied, scanned, or duplicated, in whole or in part. Due to electronic rights, some third party content may be suppressed from the eBook and/or eChapter(s).
Editorial review has deemed that any suppressed content does not materially affect the overall learning experience. Cengage Learning reserves the right to remove additional content at any time if subsequent rights restrictions require it.
54 Chapter 2 DISCRETE-TIME SIGNALS AND SYSTEMS
P2.3 Generate the following periodic sequences and plot their samples (using the stem function)
over the indicated number of periods.
1. x̃1 (n) = {. . . , −2, −1, 0, 1, 2, . . .}periodic . Plot 5 periods.
↑
2. x̃2 (n) = e0.1n [u(n) − u(n − 20]periodic . Plot 3 periods.
3. x̃3 (n) = sin(0.1πn)[u(n) − u(n − 10)]. Plot 4 periods.
4. x̃4 (n) = {. . . , 1, 2, 3, . . .}periodic + {. . . , 1, 2, 3, 4, . . .}periodic , 0 ≤ n ≤ 24. What is the
↑ ↑
period of x̃4 (n)?
P2.4 Let x(n) = {2, 4, −3, 1, −5, 4, 7}. Generate and plot the samples (use the stem function) of
↑
the following sequences.
1. x1 (n) = 2x(n − 3) + 3x(n + 4) − x(n)
2. x2 (n) = 4x(4 + n) + 5x(n + 5) + 2x(n)
3. x3 (n) = x(n + 3)x(n − 2) + x(1 − n)x(n + 1)
4. x4 (n) = 2e0.5n x(n) + cos (0.1πn) x (n + 2) , −10 ≤ n ≤ 10
P2.5 The complex exponential sequence ejω0 n or the sinusoidal sequence cos (ω0 n) are periodic if
ω0 K
the normalized frequency f0 = is a rational number; that is, f0 = , where K and N
2π N
are integers.
1. Prove the above result.
2. Generate exp(0.1πn), −100 ≤ n ≤ 100. Plot its real and imaginary parts using the stem
function. Is this sequence periodic? If it is, what is its fundamental period? From the
examination of the plot what interpretation can you give to the integers K and N above?
3. Generate and plot cos(0.1n), −20 ≤ n ≤ 20. Is this sequence periodic? What do you
conclude from the plot? If necessary examine the values of the sequence in MATLAB to
arrive at your answer.
P2.6 Using the evenodd function, decompose the following sequences into their even and odd
components. Plot these components using the stem function.
1. x1 (n) = {0, 1, 2, 3, 4, 5, 6, 7, 8, 9}.
↑
2. x2 (n) = e0.1n [u(n + 5) − u(n − 10)].
3. x3 (n) = cos(0.2πn + π/4), −20 ≤ n ≤ 20.
4. x4 (n) = e−0.05n sin(0.1πn + π/3), 0 ≤ n ≤ 100.
P2.7 A complex-valued sequence xe (n) is called conjugate-symmetric if xe (n) = x∗e (−n) and a
complex-valued sequence xo (n) is called conjugate-antisymmetric if xo (n) = −x∗o (−n).
Then, any arbitrary complex-valued sequence x(n) can be decomposed into
x(n) = xe (n) + xo (n) where xe (n) and xo (n) are given by
1 1
xe (n) = [x(n) + x∗ (−n)] and xo (n) = [x(n) − x∗ (−n)] (2.27)
2 2
respectively.
Copyright 2010 Cengage Learning. All Rights Reserved. May not be copied, scanned, or duplicated, in whole or in part. Due to electronic rights, some third party content may be suppressed from the eBook and/or eChapter(s).
Editorial review has deemed that any suppressed content does not materially affect the overall learning experience. Cengage Learning reserves the right to remove additional content at any time if subsequent rights restrictions require it.
Problems 55
1. Modify the evenodd function discussed in the text so that it accepts an arbitrary
sequence and decomposes it into its conjugate-symmetric and conjugate-antisymmetric
components by implementing (2.27).
2. Decompose the following sequence:
y(n) = x(nM )
y(n) = {. . . , −2, 3, 5, 8, . . .}
↑
to implement the above operation. Use the indexing mechanism of MATLAB with
careful attention to the origin of the time axis n = 0.
2. Generate x(n) = sin(0.125πn), − 50 ≤ n ≤ 50. Decimate x(n) by a factor of 4 to
generate y(n). Plot both x(n) and y(n) using subplot and comment on the results.
3. Repeat the above using x(n) = sin(0.5πn), − 50 ≤ n ≤ 50. Qualitatively discuss the
effect of down-sampling on signals.
P2.9 Using the conv_m function, determine the autocorrelation sequence rxx () and the
crosscorrelation sequence rxy () for the following sequences.
where k is the amount of delay in samples and α is its relative strength. We want to
estimate the delay using the correlation analysis.
Copyright 2010 Cengage Learning. All Rights Reserved. May not be copied, scanned, or duplicated, in whole or in part. Due to electronic rights, some third party content may be suppressed from the eBook and/or eChapter(s).
Editorial review has deemed that any suppressed content does not materially affect the overall learning experience. Cengage Learning reserves the right to remove additional content at any time if subsequent rights restrictions require it.
56 Chapter 2 DISCRETE-TIME SIGNALS AND SYSTEMS
1. Determine analytically the crosscorrelation ryx () in terms of the autocorrelation rxx ().
2. Let x(n) = cos(0.2πn) + 0.5 cos(0.6πn), α = 0.1, and k = 50. Generate 200 samples of
y(n) and determine its crosscorrelation. Can you obtain α and k by observing ryx ()?
Copyright 2010 Cengage Learning. All Rights Reserved. May not be copied, scanned, or duplicated, in whole or in part. Due to electronic rights, some third party content may be suppressed from the eBook and/or eChapter(s).
Editorial review has deemed that any suppressed content does not materially affect the overall learning experience. Cengage Learning reserves the right to remove additional content at any time if subsequent rights restrictions require it.
Problems 57
y = Hx
where linear shifts in h(n − k) for n = 0, . . . , Nh − 1 are arranged as rows in the matrix H.
This matrix has an interesting structure and is called a Toeplitz matrix. To investigate this
matrix, consider the sequences
function [y,H]=conv_tp(h,x)
% Linear Convolution using Toeplitz Matrix
% ----------------------------------------
% [y,H] = conv_tp(h,x)
% y = output sequence in column vector form
% H = Toeplitz matrix corresponding to sequence h so that y = Hx
% h = Impulse response sequence in column vector form
% x = input sequence in column vector form
Copyright 2010 Cengage Learning. All Rights Reserved. May not be copied, scanned, or duplicated, in whole or in part. Due to electronic rights, some third party content may be suppressed from the eBook and/or eChapter(s).
Editorial review has deemed that any suppressed content does not materially affect the overall learning experience. Cengage Learning reserves the right to remove additional content at any time if subsequent rights restrictions require it.
58 Chapter 2 DISCRETE-TIME SIGNALS AND SYSTEMS
1. Using the filter function, compute and plot the impulse response of the system over
0 ≤ n ≤ 100.
2. Determine the stability of the system from this impulse response.
3. If the input to this system is x(n) = [5 + 3 cos(0.2πn) + 4 sin(0.6πn)] u(n), determine the
response y(n) over 0 ≤ n ≤ 200 using the filter function.
P2.20 A “simple” digital differentiator is given by
which computes a backward first-order difference of the input sequence. Implement this
differentiator on the following sequences, and plot the results. Comment on the
appropriateness of this simple differentiator.
1. x(n) = 5 [u(n) − u(n − 20)]: a rectangular pulse
2. x(n) = n [u(n) − u(n − 10)] + (20 − n) [u(n − 10) − u(n − 20)]: a triangular pulse
πn
3. x(n) = sin [u(n) − u(n − 100)]: a sinusoidal pulse
25
Copyright 2010 Cengage Learning. All Rights Reserved. May not be copied, scanned, or duplicated, in whole or in part. Due to electronic rights, some third party content may be suppressed from the eBook and/or eChapter(s).
Editorial review has deemed that any suppressed content does not materially affect the overall learning experience. Cengage Learning reserves the right to remove additional content at any time if subsequent rights restrictions require it.