NivaarExam PrepOfficial exam papers ↗

22-Elec-B1 Digital Signal Processing · May 2015

Question 2 of 6: Sampling a band-limited periodic signal (12 marks)

Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)

Notes on this paper

Paper format. National Exams, May 2015 — 07-Elec-B1 Digital Signal Processing. Three hours, closed book; candidates may use one approved Casio or Sharp calculator and one aid sheet written on both sides. Six questions are printed and five constitute a complete exam, each worth 12 marks (72 printed, 60 counted). Every question and sub-part is worked below, because this set is intended as a study resource rather than an exam script. The tables of z-transform pairs and DFT properties bound into the back of the paper are reproduced from Oppenheim & Schafer and are quoted where used.

Reference texts.

  • A. V. Oppenheim and R. W. Schafer, Discrete-Time Signal Processing, 3rd ed., Pearson — Ch. 2 (LTI systems and the DTFT), Ch. 3 (the z-transform, Tables 3.1 and 3.2), Ch. 4 (sampling), Ch. 6 (filter structures), Ch. 7 (IIR and FIR design), Ch. 8 (the DFT and circular convolution), Ch. 9 (the FFT). This is the text the printed appendix tables are taken from.
  • J. G. Proakis and D. G. Manolakis, Digital Signal Processing: Principles, Algorithms and Applications, 4th ed., Pearson — Ch. 5 (DFT), Ch. 8 (IIR design), Ch. 10 (FIR design).
  • A. V. Oppenheim and A. S. Willsky, Signals and Systems, 2nd ed., Prentice Hall — Ch. 3 (Fourier series), Ch. 7 (sampling).

Question 2: Sampling a band-limited periodic signal (12 marks)

Question text not reproduced: the examination questions are © Engineers and Geoscientists BC. Open the official past paper (linked at the top of this page) to read the question, then follow the worked solution below.

Given. A periodic continuous-time signal whose Fourier series is band-limited to nineteen harmonics, sampled uniformly at six samples per period.

Given data
Fundamental period of $x_c(t)$$T_p = 1\ \text{ms} = 10^{-3}\ \text{s}$
Fundamental frequency$f_0 = 1/T_p = 1\ \text{kHz}$
Harmonic content$a_k \ne 0$ only for $|k|\le 9$
Sample spacing$T = \tfrac16\times 10^{-3}\ \text{s}$
Sampling rate$f_s = 1/T = 6\ \text{kHz}$

Find. (a) whether the sampled sequence is periodic and with what fundamental period; (b) whether $f_s$ clears the Nyquist rate; (c) the discrete-time Fourier series coefficients of $x[n]$ expressed through the given $a_k$.

Approach. Substitute $t=nT$ into the given series; the ratio $T/T_p$ then decides the period, the highest harmonic decides the Nyquist rate, and grouping the exponentials that become indistinguishable after sampling produces the discrete-time coefficients.

  1. Sample the Fourier series. Putting $t=nT$ into the given expansion, $$x[n]=x_c(nT)=\sum_{k=-9}^{9}a_k\,e^{j2\pi k nT/10^{-3}} .$$ The whole problem hinges on the ratio of the two time scales, $$\frac{T}{T_p}=\frac{\tfrac16\times 10^{-3}}{10^{-3}}=\frac16 ,$$ so the exponent becomes $j2\pi k n/6$ and $$x[n]=\sum_{k=-9}^{9}a_k\,e^{j2\pi k n/6}.$$
  2. Read off the period (part a). Every exponential $e^{j2\pi kn/6}$ repeats when $n$ advances by $6$, because $e^{j2\pi k(n+6)/6}=e^{j2\pi kn/6}e^{j2\pi k}=e^{j2\pi kn/6}$ for integer $k$. No smaller advance works for all $k$ present, since $e^{j2\pi k n/6}$ with $k=1$ needs the full six samples. The sampled signal is therefore periodic, with $$\boxed{N=6\ \text{samples}} .$$ This is the general rule that sampling a periodic signal yields a periodic sequence exactly when $T/T_p$ is rational; here $T/T_p=1/6$ and six samples span precisely one continuous-time period.
  3. Compare the sampling rate with the Nyquist rate (part b). The highest frequency actually present is the ninth harmonic, $$f_{\max}=9f_0=9\times 1\ \text{kHz}=9\ \text{kHz},$$ so the Nyquist rate — twice the highest frequency — is $$f_{\text{Nyq}}=2f_{\max}=18\ \text{kHz}.$$ The actual sampling rate is $$f_s=\frac1T=\frac{1}{\tfrac16\times10^{-3}}=6\ \text{kHz}.$$ Since $6\ \text{kHz}$ is a factor of three below $18\ \text{kHz}$, $$\boxed{f_s = 6\ \text{kHz}\ \lt\ f_{\text{Nyq}} = 18\ \text{kHz}\ \Rightarrow\ T\ \text{is NOT small enough; the signal aliases}} $$ Only the harmonics with $|k|f_0 \lt f_s/2 = 3\ \text{kHz}$, that is $|k|\le 2$, survive unaliased; harmonics 3 through 9 all fold back into the base band.
f (kHz)Line spectrum of xc(t): 19 harmonics, spacing 1 kHz-9 kHz09 kHzhighest harmonic 9 kHz -> Nyquist rate 18 kHzf (kHz)After sampling at fs = 6 kHz: replicas every 6 kHz-fs/20fs/2 = 3 kHzfs = 6 kHz-fsred = shifted replicas; they overlap the baseband, so the sampling aliases
Figure 2.1 — the line spectrum of xc(t) (top) and the periodic replication produced by sampling at 6 kHz (bottom). Replicas spaced 6 kHz apart overlap the base band because the signal occupies +/-9 kHz, which is why aliasing occurs.

Part (c) asks for the discrete-time Fourier series. Because the discrete exponential $e^{j2\pi kn/6}$ depends on $k$ only through $k \bmod 6$, harmonics whose indices differ by a multiple of six are indistinguishable once sampled — that indistinguishability is the aliasing found in part (b), now expressed algebraically.

  1. Collect the aliased harmonics into six bins. Write $x[n]=\sum_{m=0}^{5}b_m e^{j2\pi mn/6}$ and gather every original coefficient whose index is congruent to $m$ modulo 6: $$\boxed{\,b_m=\sum_{\substack{|k|\le 9\\ k\equiv m\ (\mathrm{mod}\ 6)}} a_k ,\qquad m=0,1,\dots,5\,}$$
  2. Write the six sums out explicitly. Sweeping $k=-9,\dots,9$ and sorting by $k \bmod 6$ gives $$\begin{aligned} b_0 &= a_{-6}+a_{0}+a_{6}, &\qquad b_3 &= a_{-9}+a_{-3}+a_{3}+a_{9},\\ b_1 &= a_{-5}+a_{1}+a_{7}, &\qquad b_4 &= a_{-8}+a_{-2}+a_{4},\\ b_2 &= a_{-4}+a_{2}+a_{8}, &\qquad b_5 &= a_{-7}+a_{-1}+a_{5}. \end{aligned}$$ The nineteen given coefficients are accounted for exactly once each ($3+3+3+4+3+3=19$), which is a useful arithmetic check on the bookkeeping. Note that $b_3$ — the bin at $\omega=\pi$, the folding frequency — collects four terms rather than three, because both extreme harmonics $a_{\pm 9}$ land there.
Aliasing map: 19 harmonics fold onto N = 6 DFS binsbin mcontinuous-time harmonics k that fold onto itm = 0a(-6), a(0), a(6)m = 1a(-5), a(1), a(7)m = 2a(-4), a(2), a(8)m = 3a(-9), a(-3), a(3), a(9)m = 4a(-8), a(-2), a(4)m = 5a(-7), a(-1), a(5)each row is one DFS coefficient b(m) = sum of the listed a(k)
Figure 2.2 — the folding map. Each row is one discrete-time Fourier series coefficient; the red row at m = 3 is the folding frequency, where four continuous-time harmonics collapse together.

If instead the discrete Fourier series is written in the analysis form used in the appendix, $X[m]=\sum_{n=0}^{5}x[n]W_{6}^{mn}$, the two conventions differ only by the length factor, $X[m]=6\,b_m$.

Final results
PartQuantityResult
(a)Periodicity of $x[n]$periodic, fundamental period $N=6$ samples
(b)Highest frequency present$f_{\max}=9\ \text{kHz}$
(b)Nyquist rate$18\ \text{kHz}$
(b)Actual sampling rate$f_s = 6\ \text{kHz}$ — below Nyquist, so aliasing occurs
(b)Unaliased harmonics$|k|\le 2$ only
(c)DTFS coefficients$b_0=a_{-6}+a_0+a_6$; $b_1=a_{-5}+a_1+a_7$; $b_2=a_{-4}+a_2+a_8$; $b_3=a_{-9}+a_{-3}+a_3+a_9$; $b_4=a_{-8}+a_{-2}+a_4$; $b_5=a_{-7}+a_{-1}+a_5$