Browse papers
A

Section A: Long Answer Questions

Attempt all / any as specified.

4 questions
1long12 marks

(a) Define a discrete-time LTI system. State and prove the conditions on the impulse response h[n]h[n] for the system to be (i) causal and (ii) BIBO stable. [6]

(b) A discrete-time LTI system has impulse response h[n]=(12)nu[n]h[n] = \left(\tfrac{1}{2}\right)^n u[n]. Determine the output y[n]y[n] when the input is x[n]=u[n]u[n4]x[n] = u[n] - u[n-4] using linear convolution, and comment on whether the system is stable. [6]

(a) Discrete-time LTI system; causality and stability conditions

A discrete-time LTI system is a system that is both linear (obeys superposition: T{ax1[n]+bx2[n]}=aT{x1[n]}+bT{x2[n]}T\{a x_1[n]+b x_2[n]\}=a\,T\{x_1[n]\}+b\,T\{x_2[n]\}) and time-invariant (a shift in input produces an identical shift in output). Such a system is completely characterized by its impulse response h[n]=T{δ[n]}h[n]=T\{\delta[n]\}, and its output is the convolution

y[n]=x[n]h[n]=k=x[k]h[nk].y[n]=x[n]*h[n]=\sum_{k=-\infty}^{\infty}x[k]\,h[n-k].

(i) Causality. A system is causal if the output at time nn depends only on present and past inputs. Writing y[n]=kh[k]x[nk]y[n]=\sum_k h[k]\,x[n-k], the term h[k]x[nk]h[k]x[n-k] uses a future input when k<0k<0. To exclude all future inputs we require

h[n]=0for n<0.\boxed{h[n]=0 \quad \text{for } n<0.}

Proof: If h[k]=0h[k]=0 for k<0k<0, then y[n]=k0h[k]x[nk]y[n]=\sum_{k\ge 0} h[k]x[n-k] involves only x[n],x[n1],x[n],x[n-1],\dots (past/present). Conversely, if h[m]0h[m]\neq0 for some m<0m<0, then applying x[n]=δ[n]x[n]=\delta[n] gives y[m]=h[m]0y[m]=h[m]\neq0 at time m<0m<0, i.e. output appears before the input — non-causal.

(ii) BIBO stability. A system is BIBO stable if every bounded input produces a bounded output. The necessary and sufficient condition is absolute summability of the impulse response:

n=h[n]<.\boxed{\sum_{n=-\infty}^{\infty}|h[n]|<\infty.}

Proof (sufficiency): If x[n]M|x[n]|\le M, then y[n]=kh[k]x[nk]Mkh[k]<|y[n]|=\big|\sum_k h[k]x[n-k]\big|\le M\sum_k|h[k]|<\infty. (Necessity): Choose the bounded input x[n]=sgn(h[n])x[n]=\operatorname{sgn}(h[-n]); then y[0]=kh[k]y[0]=\sum_k|h[k]|, which is unbounded if the sum diverges.

(b) Output by linear convolution

Given h[n]=(12)nu[n]h[n]=\left(\tfrac12\right)^n u[n] and x[n]=u[n]u[n4]={1,1,1,1}x[n]=u[n]-u[n-4]=\{1,1,1,1\} for n=0,1,2,3n=0,1,2,3 (and 0 elsewhere).

y[n]=k=03x[k]h[nk]=k=max(0,n3)min(3,n)(12)nk,0n.y[n]=\sum_{k=0}^{3}x[k]\,h[n-k]=\sum_{k=\max(0,n-3)}^{\min(3,n)}\left(\tfrac12\right)^{n-k},\qquad 0\le n.

For 0n30\le n\le 3 (rising part):

y[n]=k=0n(12)nk=m=0n(12)m=2(12)n.y[n]=\sum_{k=0}^{n}\left(\tfrac12\right)^{n-k}=\sum_{m=0}^{n}\left(\tfrac12\right)^{m}=2-\left(\tfrac12\right)^{n}.

So y[0]=1,  y[1]=1.5,  y[2]=1.75,  y[3]=1.875.y[0]=1,\; y[1]=1.5,\; y[2]=1.75,\; y[3]=1.875.

For n4n\ge 4 (decaying tail): all four input samples (k=0..3k=0..3) contribute:

y[n]=k=03(12)nk=(12)nk=032k=(12)n(15)=15(12)n.y[n]=\sum_{k=0}^{3}\left(\tfrac12\right)^{n-k}=\left(\tfrac12\right)^{n}\sum_{k=0}^{3}2^{k}=\left(\tfrac12\right)^{n}(15)=15\left(\tfrac12\right)^{n}.

Check n=4n=4: 15/16=0.937515/16=0.9375.

Result:

y[n]={2(12)n,0n315(12)n,n40,n<0y[n]=\begin{cases}2-\left(\tfrac12\right)^{n}, & 0\le n\le 3\\[4pt] 15\left(\tfrac12\right)^{n}, & n\ge 4\\ 0,&n<0\end{cases}

Numerically: {1,1.5,1.75,1.875,0.9375,0.4688,}\{1,\,1.5,\,1.75,\,1.875,\,0.9375,\,0.4688,\dots\}.

Stability: n=0(12)n=1112=2<\sum_{n=0}^{\infty}\left|\left(\tfrac12\right)^{n}\right|=\dfrac{1}{1-\tfrac12}=2<\infty, so h[n]h[n] is absolutely summable and the system is BIBO stable.

discrete-time-signalsdiscrete-time-systemsconvolution
2long16 marks

Consider an LTI system described by the difference equation

y[n]34y[n1]+18y[n2]=x[n]+x[n1].y[n] - \tfrac{3}{4}y[n-1] + \tfrac{1}{8}y[n-2] = x[n] + x[n-1].

(a) Determine the system function H(z)H(z) and its region of convergence assuming the system is causal. [5]

(b) Plot the pole-zero diagram and comment on the stability of the system. [4]

(c) Find the unit-impulse response h[n]h[n] using the inverse Z-transform (partial fraction method). [5]

(d) Sketch the general shape of the magnitude response H(ejω)|H(e^{j\omega})| and identify whether the system behaves as a low-pass or high-pass filter. [2]

(a) System function H(z)H(z) and ROC

Taking the Z-transform of y[n]34y[n1]+18y[n2]=x[n]+x[n1]y[n]-\tfrac34 y[n-1]+\tfrac18 y[n-2]=x[n]+x[n-1]:

Y(z)(134z1+18z2)=X(z)(1+z1).Y(z)\Big(1-\tfrac34 z^{-1}+\tfrac18 z^{-2}\Big)=X(z)\big(1+z^{-1}\big). H(z)=1+z1134z1+18z2=1+z1(112z1)(114z1).\boxed{H(z)=\frac{1+z^{-1}}{1-\tfrac34 z^{-1}+\tfrac18 z^{-2}}=\frac{1+z^{-1}}{\left(1-\tfrac12 z^{-1}\right)\left(1-\tfrac14 z^{-1}\right)}.}

Poles at z=12z=\tfrac12 and z=14z=\tfrac14; zero at z=1z=-1 (and a trivial zero at z=0z=0). For a causal system the ROC is outside the outermost pole:

ROC: z>12.\text{ROC}:\ |z|>\tfrac12.

(b) Pole-zero plot and stability

  • Poles (×\times): z=0.5z=0.5 and z=0.25z=0.25 on the positive real axis, inside the unit circle.
  • Zero (\circ): z=1z=-1 on the unit circle (negative real axis).

Description of the diagram: draw the unit circle; mark two ×\times on the positive real axis at 0.250.25 and 0.50.5, and one \circ at 1-1. Since both poles lie inside the unit circle and the causal ROC z>0.5|z|>0.5 includes the unit circle, the system is BIBO stable.

(c) Impulse response by partial fractions

Expand H(z)H(z) in z1z^{-1}. Write

H(z)=1+z1(112z1)(114z1)=A112z1+B114z1.H(z)=\frac{1+z^{-1}}{\left(1-\tfrac12 z^{-1}\right)\left(1-\tfrac14 z^{-1}\right)}=\frac{A}{1-\tfrac12 z^{-1}}+\frac{B}{1-\tfrac14 z^{-1}}.

Multiply out and equate. Let u=z1u=z^{-1}:

A(114u)+B(112u)=1+u.A\left(1-\tfrac14 u\right)+B\left(1-\tfrac12 u\right)=1+u.

Constant: A+B=1A+B=1. Coefficient of uu: 14A12B=1-\tfrac14 A-\tfrac12 B=1. From these: 14A12(1A)=114A12=1A=6-\tfrac14 A-\tfrac12(1-A)=1\Rightarrow \tfrac14 A-\tfrac12=1\Rightarrow A=6, hence B=5B=-5.

Using 11az1anu[n]\dfrac{1}{1-a z^{-1}}\leftrightarrow a^{n}u[n] (causal):

h[n]=6(12)nu[n]5(14)nu[n].\boxed{h[n]=6\left(\tfrac12\right)^{n}u[n]-5\left(\tfrac14\right)^{n}u[n].}

Check: h[0]=65=1h[0]=6-5=1 ✓ (matches the constant numerator term); h[1]=6(0.5)5(0.25)=31.25=1.75h[1]=6(0.5)-5(0.25)=3-1.25=1.75.

(d) Magnitude response and filter type

At ω=0\omega=0 (z=1z=1): H(1)=1+110.75+0.125=20.375=5.33H(1)=\dfrac{1+1}{1-0.75+0.125}=\dfrac{2}{0.375}=5.33 (large gain at DC). At ω=π\omega=\pi (z=1z=-1): H(1)=11=0H(-1)=\dfrac{1-1}{\dots}=0 (zero at z=1z=-1 forces zero gain at ω=π\omega=\pi). Thus H(ejω)|H(e^{j\omega})| is maximum at ω=0\omega=0 and falls to zero at ω=π\omega=\pi — a monotonically decreasing curve. The system behaves as a low-pass filter.

z-transformfrequency-responsediscrete-time-systems
3long12 marks

(a) Explain the impulse invariance method and the bilinear transformation method of designing IIR digital filters. Clearly state the mapping equations and discuss the problem of frequency warping. [6]

(b) Design a first-order digital low-pass IIR filter using the bilinear transformation from the analog prototype Ha(s)=Ωcs+ΩcH_a(s) = \dfrac{\Omega_c}{s + \Omega_c} with a 3 dB cut-off frequency of ωc=0.2π\omega_c = 0.2\pi rad/sample. Assume T=1T = 1 s and give the resulting transfer function H(z)H(z). [6]

(a) Impulse invariance vs bilinear transformation

Impulse invariance method. The digital filter is designed so that its impulse response is a sampled version of the analog filter's impulse response: h[n]=Tsha(nTs)h[n]=T_s\,h_a(nT_s). If Ha(s)=kAkspkH_a(s)=\sum_k \dfrac{A_k}{s-p_k}, the mapping of each pole is

Akspk  TsAk1epkTsz1,s=pkz=epkTs.\frac{A_k}{s-p_k}\ \longrightarrow\ \frac{T_s\,A_k}{1-e^{p_k T_s}z^{-1}},\qquad s=p_k \mapsto z=e^{p_k T_s}.

The analog frequency axis maps linearly (Ωω=ΩTs\Omega \to \omega=\Omega T_s), but because sampling is periodic in frequency, aliasing occurs — it is suitable only for band-limited (low-pass / band-pass) analog prototypes, not high-pass.

Bilinear transformation method. It maps the entire jΩj\Omega axis onto the unit circle once, avoiding aliasing, using

s=2T1z11+z1z=1+(T/2)s1(T/2)s.\boxed{s=\frac{2}{T}\,\frac{1-z^{-1}}{1+z^{-1}}}\qquad\Leftrightarrow\qquad z=\frac{1+(T/2)s}{1-(T/2)s}.

The left-half ss-plane maps inside the unit circle, so stability is preserved. However the frequency axis is mapped non-linearly:

Ω=2Ttan ⁣(ω2)(frequency warping).\Omega=\frac{2}{T}\tan\!\left(\frac{\omega}{2}\right)\qquad\text{(frequency warping).}

This compresses high analog frequencies into ωπ\omega\le\pi. Frequency warping distorts the frequency axis, so critical band edges must be pre-warped: Ωc=2Ttan ⁣(ωc2)\Omega_c=\dfrac{2}{T}\tan\!\left(\dfrac{\omega_c}{2}\right) before applying the transform.

(b) First-order low-pass IIR design (bilinear)

Given Ha(s)=Ωcs+ΩcH_a(s)=\dfrac{\Omega_c}{s+\Omega_c}, digital cut-off ωc=0.2π\omega_c=0.2\pi, T=1T=1 s.

Step 1 — Pre-warp the cut-off:

Ωc=2Ttan ⁣(ωc2)=2tan(0.1π)=2tan(18)=2(0.3249)=0.6498 rad/s.\Omega_c=\frac{2}{T}\tan\!\left(\frac{\omega_c}{2}\right)=2\tan(0.1\pi)=2\tan(18^\circ)=2(0.3249)=0.6498\ \text{rad/s}.

Step 2 — Apply s=2T1z11+z1=21z11+z1s=\dfrac{2}{T}\dfrac{1-z^{-1}}{1+z^{-1}}=2\dfrac{1-z^{-1}}{1+z^{-1}}:

H(z)=Ωc2(1z1)1+z1+Ωc=Ωc(1+z1)2(1z1)+Ωc(1+z1).H(z)=\frac{\Omega_c}{\dfrac{2(1-z^{-1})}{1+z^{-1}}+\Omega_c}=\frac{\Omega_c(1+z^{-1})}{2(1-z^{-1})+\Omega_c(1+z^{-1})}.

Step 3 — Substitute Ωc=0.6498\Omega_c=0.6498: Numerator =0.6498(1+z1)=0.6498(1+z^{-1}). Denominator =(2+0.6498)+(2+0.6498)z1=2.64981.3502z1=(2+0.6498)+(-2+0.6498)z^{-1}=2.6498-1.3502z^{-1}.

H(z)=0.6498(1+z1)2.64981.3502z1.H(z)=\frac{0.6498(1+z^{-1})}{2.6498-1.3502z^{-1}}.

Normalize by dividing by 2.64982.6498:

H(z)=0.2452(1+z1)10.5096z1.\boxed{H(z)=\frac{0.2452\,(1+z^{-1})}{1-0.5096\,z^{-1}}.}

Check (DC gain): at z=1z=1, H(1)=0.2452(2)10.5096=0.49040.4904=1H(1)=\dfrac{0.2452(2)}{1-0.5096}=\dfrac{0.4904}{0.4904}=1 — unity gain at DC, confirming a correct low-pass design. The pole at z=0.5096z=0.5096 (inside the unit circle) confirms stability.

digital-filter-designiir-filtersbilinear-transform
4long12 marks

(a) Derive the decimation-in-time (DIT) radix-2 FFT algorithm for an NN-point DFT and draw the complete signal-flow graph (butterfly diagram) for N=8N = 8. [8]

(b) Compare the number of complex multiplications and additions required for direct DFT computation versus the radix-2 FFT for N=1024N = 1024, and comment on the computational savings. [4]

(a) Decimation-in-time (DIT) radix-2 FFT derivation

The NN-point DFT is X[k]=n=0N1x[n]WNnkX[k]=\sum_{n=0}^{N-1}x[n]W_N^{nk}, where WN=ej2π/NW_N=e^{-j2\pi/N}. Split the sum into even (n=2rn=2r) and odd (n=2r+1n=2r+1) indexed samples:

X[k]=r=0N/21x[2r]WN2rk+r=0N/21x[2r+1]WN(2r+1)k.X[k]=\sum_{r=0}^{N/2-1}x[2r]W_N^{2rk}+\sum_{r=0}^{N/2-1}x[2r+1]W_N^{(2r+1)k}.

Using WN2rk=WN/2rkW_N^{2rk}=W_{N/2}^{rk} and factoring WNkW_N^{k} from the odd sum:

X[k]=G[k]DFT of even part+WNkH[k]DFT of odd part,k=0,,N/21.\boxed{X[k]=\underbrace{G[k]}_{\text{DFT of even part}}+W_N^{k}\underbrace{H[k]}_{\text{DFT of odd part}},\quad k=0,\dots,N/2-1.}

Because G[k]G[k] and H[k]H[k] are periodic with period N/2N/2 and WNk+N/2=WNkW_N^{k+N/2}=-W_N^{k}:

X[k+N/2]=G[k]WNkH[k].X[k+N/2]=G[k]-W_N^{k}H[k].

This pair forms the butterfly: two outputs from two inputs using one complex multiply by the twiddle factor WNkW_N^{k}. Recursively splitting each N/2N/2-DFT gives log2N\log_2 N stages.

N=8N=8 signal-flow graph (described). Inputs are applied in bit-reversed order: x[0],x[4],x[2],x[6],x[1],x[5],x[3],x[7]x[0],x[4],x[2],x[6],x[1],x[5],x[3],x[7]. There are log28=3\log_2 8 = 3 stages, each with N/2=4N/2=4 butterflies (12 butterflies total).

  • Stage 1: four 2-point butterflies using twiddle W20=1W_2^0=1 on adjacent pairs.
  • Stage 2: combine into 4-point DFTs using W40,W41W_4^0,W_4^1.
  • Stage 3: combine into the 8-point DFT using W80,W81,W82,W83W_8^0,W_8^1,W_8^2,W_8^3. Outputs X[0]X[7]X[0]\dots X[7] emerge in natural order. Each butterfly: top output =a+Wb=a+W\,b, bottom output =aWb=a-W\,b.
Stage1      Stage2        Stage3
x[0]\ /--•--\        /----•----\  X[0]
x[4]/ \  W2  \      / W4   /     X[1]
x[2]\ /--•--/  W4  /      \ W8   X[2]
x[6]/ \        \  /        \     X[3]
x[1]...                         X[4]
   (bit-reversed input -> natural-order output)

(b) Computational comparison for N=1024N=1024

Direct DFT: N2N^2 complex multiplications and N(N1)N(N-1) complex additions.

  • Multiplications =10242=1,048,576=1024^2=1{,}048{,}576.
  • Additions =1024×1023=1,047,552=1024\times1023=1{,}047{,}552.

Radix-2 FFT: N2log2N\dfrac{N}{2}\log_2 N complex multiplications and Nlog2NN\log_2 N complex additions. With log21024=10\log_2 1024=10:

  • Multiplications =10242×10=5120=\dfrac{1024}{2}\times10=5120.
  • Additions =1024×10=10,240=1024\times10=10{,}240.

Savings (multiplications): 1,048,5765120205×\dfrac{1{,}048{,}576}{5120}\approx 205\times fewer multiplications. In general the speed-up factor is N2(N/2)log2N=2Nlog2N\dfrac{N^2}{(N/2)\log_2 N}=\dfrac{2N}{\log_2 N}, which grows with NN — the FFT makes large-NN spectral analysis computationally feasible.

fftdftradix-2
B

Section B: Short Answer Questions

Attempt all / any as specified.

8 questions
5short6 marks

Compute the 4-point DFT of the sequence x[n]={1, 2, 3, 4}x[n] = \{1,\ 2,\ 3,\ 4\} using the DFT definition. Show the magnitude and phase of each DFT coefficient.

4-point DFT of x[n]={1,2,3,4}x[n]=\{1,2,3,4\}

Using X[k]=n=03x[n]ej2πkn/4X[k]=\sum_{n=0}^{3}x[n]\,e^{-j2\pi kn/4} with W4=ejπ/2=jW_4=e^{-j\pi/2}=-j, so W40=1, W41=j, W42=1, W43=jW_4^{0}=1,\ W_4^{1}=-j,\ W_4^{2}=-1,\ W_4^{3}=j.

X[0]X[0] (sum): 1+2+3+4=101+2+3+4=10.

X[1]X[1]: 1+2(j)+3(1)+4(j)=(13)+j(2+4)=2+2j1+2(-j)+3(-1)+4(j)=(1-3)+j(-2+4)=-2+2j.

X[2]X[2]: 1+2(1)+3(1)+4(1)=12+34=21+2(-1)+3(1)+4(-1)=1-2+3-4=-2.

X[3]X[3]: 1+2(j)+3(1)+4(j)=(13)+j(24)=22j1+2(j)+3(-1)+4(-j)=(1-3)+j(2-4)=-2-2j.

X[k]={10, 2+2j, 2, 22j}\boxed{X[k]=\{\,10,\ -2+2j,\ -2,\ -2-2j\,\}}

Magnitude and phase:

kkX[k]X[k]X[k]\lvert X[k]\rvertX[k]\angle X[k]
01010101000^\circ
12+2j-2+2j222.8282\sqrt2\approx2.828135135^\circ
22-222180180^\circ
322j-2-2j222.8282\sqrt2\approx2.828135-135^\circ

Note the conjugate symmetry X[3]=X[1]X[3]=X[1]^{*}, as expected for a real input sequence.

dftcircular-convolution
6short6 marks

State the Nyquist sampling theorem. An analog signal xa(t)=3cos(2π1000t)+2cos(2π3000t)x_a(t) = 3\cos(2\pi \cdot 1000\,t) + 2\cos(2\pi \cdot 3000\,t) is sampled at fs=4000f_s = 4000 Hz. Determine which frequency components are aliased and write the reconstructed signal's frequency content.

Nyquist sampling theorem

A band-limited signal containing no frequency components above fmaxf_{max} can be perfectly reconstructed from its samples if it is sampled at a rate fs>2fmaxf_s>2f_{max}. The minimum rate 2fmax2f_{max} is the Nyquist rate, and fs/2f_s/2 is the Nyquist (folding) frequency. Frequencies above fs/2f_s/2 are aliased (folded back) into the band [0,fs/2][0,f_s/2].

Analysis of the given signal

xa(t)=3cos(2π1000t)+2cos(2π3000t)x_a(t)=3\cos(2\pi\cdot1000\,t)+2\cos(2\pi\cdot3000\,t), sampled at fs=4000f_s=4000 Hz \Rightarrow folding frequency =2000=2000 Hz.

  • 1000 Hz component: 1000<20001000<2000 Hz, so it is below the folding frequency — sampled correctly (not aliased). It appears at 1000 Hz.
  • 3000 Hz component: 3000>20003000>2000 Hz, so it is aliased. Its apparent (alias) frequency is
falias=ffs=30004000=1000 Hz.f_{alias}=|f-f_s|=|3000-4000|=1000\ \text{Hz}.

So the 3000 Hz tone masquerades as a 1000 Hz tone.

Reconstructed signal

Both components fall at 1000 Hz after sampling and add (same amplitude/phase form of cosine):

x^(t)=3cos(2π1000t)+2cos(2π1000t)=5cos(2π1000t).\hat{x}(t)=3\cos(2\pi\cdot1000\,t)+2\cos(2\pi\cdot1000\,t)=5\cos(2\pi\cdot1000\,t).

The reconstructed signal contains only a single 1000 Hz component of amplitude 5 — the 3000 Hz tone is irretrievably lost to aliasing. To avoid this, fsf_s must exceed 2×3000=60002\times3000=6000 Hz.

samplingaliasing
7short6 marks

Explain the window method of FIR filter design. Compare the rectangular, Hamming and Blackman windows in terms of main-lobe width and peak side-lobe level, and discuss their effect on the transition band and stop-band attenuation of the designed filter.

Window method of FIR filter design

Procedure. Start from the desired (ideal) frequency response Hd(ejω)H_d(e^{j\omega}) and obtain its ideal impulse response hd[n]h_d[n] by the inverse DTFT (e.g., for an ideal low-pass filter hd[n]=sin(ωcn)πnh_d[n]=\dfrac{\sin(\omega_c n)}{\pi n}). This is infinite and non-causal, so it is truncated and tapered by multiplying with a finite-length window w[n]w[n]:

h[n]=hd[n]w[n],0nN1.h[n]=h_d[n]\,w[n],\qquad 0\le n\le N-1.

Multiplication in time = convolution in frequency: H(ejω)=HdWH(e^{j\omega})=H_d * W. The window's spectrum has a main lobe (which smears the ideal transition, setting the transition-band width) and side lobes (which produce passband/stopband ripple, setting the stopband attenuation). Choosing the window trades transition width against attenuation.

Comparison of windows

WindowMain-lobe widthPeak side-lobe levelMin. stopband attenuation
Rectangular4π/N4\pi/N (narrowest)13-13 dB (worst)21\approx 21 dB
Hamming8π/N8\pi/N41-41 dB53\approx 53 dB
Blackman12π/N12\pi/N (widest)58-58 dB (best)74\approx 74 dB

Effect on the designed filter

  • Rectangular has the narrowest main lobe, giving the sharpest (narrowest) transition band, but its high 13-13 dB side lobes cause large ripple (Gibbs phenomenon) and poor stopband attenuation (~21 dB).
  • Hamming lowers side lobes to 41-41 dB, giving ~53 dB stopband attenuation at the cost of a wider transition band (about twice the rectangular).
  • Blackman has the lowest side lobes (58-58 dB) and the best stopband attenuation (~74 dB), but the widest transition band.

Trade-off: lower side-lobe level \Rightarrow greater stopband attenuation but wider main lobe \Rightarrow wider transition band. The window is chosen to meet the required attenuation, and NN is increased to sharpen the transition.

digital-filter-designfir-filterswindowing
8short6 marks

Distinguish between linear convolution and correlation of two discrete-time sequences. Compute the cross-correlation rxy[l]r_{xy}[l] of x[n]={1, 2, 1}x[n] = \{1,\ 2,\ 1\} and y[n]={1, 1, 2}y[n] = \{1,\ -1,\ 2\} for all lags ll.

Linear convolution vs correlation

ConvolutionCorrelation
Definitiony[n]=kx[k]h[nk]y[n]=\sum_k x[k]\,h[n-k]rxy[l]=nx[n]y[nl]r_{xy}[l]=\sum_n x[n]\,y[n-l]
Second sequenceFolded (time-reversed) then shiftedNot folded, only shifted
PurposeOutput of an LTI system to an inputMeasure of similarity / time-shift between two signals
SymmetryCommutative: xh=hxx*h=h*xNot commutative: rxy[l]=ryx[l]r_{xy}[l]=r_{yx}[-l]

In short, correlation is convolution with one sequence time-reversed: rxy[l]=x[l]y[l]r_{xy}[l]=x[l]*y[-l].

Cross-correlation of x[n]={1,2,1}x[n]=\{1,2,1\} and y[n]={1,1,2}y[n]=\{1,-1,2\}

Using rxy[l]=nx[n]y[nl]r_{xy}[l]=\sum_n x[n]\,y[n-l] (slide yy relative to xx). Both length-3 sequences (indices 0..2), so ll ranges from 2-2 to +2+2.

  • l=2l=-2: overlap x[0]y[2]=12=2x[0]y[2]=1\cdot2=2. rxy[2]=2\Rightarrow r_{xy}[-2]=2.
  • l=1l=-1: x[0]y[1]+x[1]y[2]=1(1)+2(2)=1+4=3x[0]y[1]+x[1]y[2]=1(-1)+2(2)=-1+4=3. 3\Rightarrow 3.
  • l=0l=0: x[0]y[0]+x[1]y[1]+x[2]y[2]=1(1)+2(1)+1(2)=12+2=1x[0]y[0]+x[1]y[1]+x[2]y[2]=1(1)+2(-1)+1(2)=1-2+2=1. 1\Rightarrow 1.
  • l=+1l=+1: x[1]y[0]+x[2]y[1]=2(1)+1(1)=21=1x[1]y[0]+x[2]y[1]=2(1)+1(-1)=2-1=1. 1\Rightarrow 1.
  • l=+2l=+2: x[2]y[0]=11=1x[2]y[0]=1\cdot1=1. 1\Rightarrow 1.
rxy[l]={2, 3, 1, 1, 1},l=2,1,0,1,2\boxed{r_{xy}[l]=\{\,2,\ 3,\ \underline{1},\ 1,\ 1\,\},\quad l=-2,-1,0,1,2}

where the underlined value is at lag l=0l=0.

convolutioncorrelation
9short6 marks

(a) Determine the Z-transform and the region of convergence of x[n]=anu[n]+bnu[n1]x[n] = a^n u[n] + b^n u[-n-1], stating the condition on aa and bb for the ROC to exist. [4]

(b) List two important properties of the region of convergence of the Z-transform. [2]

(a) Z-transform and ROC of x[n]=anu[n]+bnu[n1]x[n]=a^n u[n]+b^n u[-n-1]

Right-sided term anu[n]a^n u[n]:

n=0anzn=11az1,ROC: z>a.\sum_{n=0}^{\infty}a^n z^{-n}=\frac{1}{1-a z^{-1}},\qquad \text{ROC: } |z|>|a|.

Left-sided term bnu[n1]b^n u[-n-1] (nonzero for n1n\le -1):

n=1bnzn=11bz1,ROC: z<b.\sum_{n=-\infty}^{-1}b^n z^{-n}=-\frac{1}{1-b z^{-1}},\qquad \text{ROC: } |z|<|b|.

Therefore

X(z)=11az111bz1,ROC: a<z<b.\boxed{X(z)=\frac{1}{1-a z^{-1}}-\frac{1}{1-b z^{-1}},\qquad \text{ROC: } |a|<|z|<|b|.}

The ROC is the annulus between the two poles. Condition for existence: the two regions overlap only if a<b|a|<|b|. If ab|a|\ge|b| the ROC is empty and the Z-transform does not converge.

(b) Two important properties of the ROC

  1. The ROC is a connected annular region (r1<z<r2r_1<|z|<r_2) centered at the origin, and it contains no poles (poles lie on its boundary).
  2. Sided-ness determines the ROC: for a causal/right-sided sequence the ROC is the exterior of the outermost pole (z>rmax|z|>r_{max}); for an anti-causal/left-sided sequence it is the interior of the innermost pole (z<rmin|z|<r_{min}). (Also: an LTI system is stable iff its ROC includes the unit circle.)
z-transformregion-of-convergence
10short6 marks

For the FIR system y[n]=x[n]x[n1]y[n] = x[n] - x[n-1], determine the frequency response H(ejω)H(e^{j\omega}), and obtain expressions for its magnitude and phase response. Sketch H(ejω)|H(e^{j\omega})| over 0ωπ0 \le \omega \le \pi and state whether the system is linear phase.

Frequency response of y[n]=x[n]x[n1]y[n]=x[n]-x[n-1] (first difference)

Impulse response h[n]={1,1}h[n]=\{1,-1\}, so

H(ejω)=nh[n]ejωn=1ejω.H(e^{j\omega})=\sum_n h[n]e^{-j\omega n}=1-e^{-j\omega}.

Factor out ejω/2e^{-j\omega/2}:

H(ejω)=ejω/2(ejω/2ejω/2)=ejω/22jsin ⁣(ω2).H(e^{j\omega})=e^{-j\omega/2}\big(e^{j\omega/2}-e^{-j\omega/2}\big)=e^{-j\omega/2}\cdot 2j\sin\!\left(\tfrac{\omega}{2}\right).

Writing 2j=2ejπ/22j=2e^{j\pi/2}:

H(ejω)=2sin ⁣(ω2)ej(π2ω2).\boxed{H(e^{j\omega})=2\sin\!\left(\tfrac{\omega}{2}\right)\,e^{\,j\left(\frac{\pi}{2}-\frac{\omega}{2}\right)}.}

Magnitude response:

H(ejω)=2sin ⁣(ω2).|H(e^{j\omega})|=2\left|\sin\!\left(\tfrac{\omega}{2}\right)\right|.

At ω=0\omega=0: 00. At ω=π\omega=\pi: 22. It rises monotonically from 0 to 2 over 0ωπ0\le\omega\le\pi — a high-pass (differentiator-like) characteristic; the sketch is a half-sine rising from the origin to a peak of 2 at ω=π\omega=\pi.

Phase response:

H(ejω)=π2ω2.\angle H(e^{j\omega})=\frac{\pi}{2}-\frac{\omega}{2}.

Linear phase? The phase π2ω2\dfrac{\pi}{2}-\dfrac{\omega}{2} is an affine (linear) function of ω\omega with constant slope 12-\tfrac12. Hence the system has constant group delay τ=dHdω=12\tau=-\dfrac{d\angle H}{d\omega}=\tfrac12 sample and is a (generalized) linear-phase system — expected since h[n]={1,1}h[n]=\{1,-1\} is anti-symmetric (Type-IV linear-phase FIR).

frequency-responsediscrete-time-systems
11short6 marks

Classify the following signals/systems with justification: (a) Is x[n]=cos(0.3πn)x[n] = \cos(0.3\pi n) periodic? If so, find its fundamental period. (b) Is the system y[n]=nx[n]y[n] = n\,x[n] linear, time-invariant, and stable?

(a) Periodicity of x[n]=cos(0.3πn)x[n]=\cos(0.3\pi n)

A discrete-time sinusoid cos(ω0n)\cos(\omega_0 n) is periodic iff ω02π=kN\dfrac{\omega_0}{2\pi}=\dfrac{k}{N} is rational. Here ω0=0.3π\omega_0=0.3\pi, so

ω02π=0.3π2π=0.32=320.\frac{\omega_0}{2\pi}=\frac{0.3\pi}{2\pi}=\frac{0.3}{2}=\frac{3}{20}.

This is rational, so the signal is periodic. The fundamental period is the smallest integer NN with ω0N=2πk\omega_0 N=2\pi k:

N=2πk0.3π=20k3  k=3, N=20.N=\frac{2\pi k}{0.3\pi}=\frac{20k}{3}\ \Rightarrow\ k=3,\ \boxed{N=20}.

(b) Properties of the system y[n]=nx[n]y[n]=n\,x[n]

Linear? Yes. For inputs x1,x2x_1,x_2: T{ax1+bx2}=n(ax1[n]+bx2[n])=anx1[n]+bnx2[n]=ay1[n]+by2[n]T\{a x_1+b x_2\}=n(a x_1[n]+b x_2[n])=a\,n x_1[n]+b\,n x_2[n]=a\,y_1[n]+b\,y_2[n]. Superposition holds.

Time-invariant? No. Shift the input: response to x[nn0]x[n-n_0] is nx[nn0]n\,x[n-n_0], but a shifted output is y[nn0]=(nn0)x[nn0]y[n-n_0]=(n-n_0)x[n-n_0]. Since nx[nn0](nn0)x[nn0]n\,x[n-n_0]\neq(n-n_0)x[n-n_0], the system is time-varying (the gain factor nn depends explicitly on time).

Stable? No. The multiplier nn grows without bound. For the bounded input x[n]=1x[n]=1 (all n0n\ge0), y[n]=ny[n]=n\to\infty, so a bounded input produces an unbounded output — BIBO unstable.

Summary: x[n]=cos(0.3πn)x[n]=\cos(0.3\pi n) is periodic with N=20N=20; the system y[n]=nx[n]y[n]=n x[n] is linear but time-varying and unstable.

discrete-time-signals
12short4 marks

Explain how the Discrete Fourier Transform (DFT) is related to the Discrete-Time Fourier Transform (DTFT). State two important properties of the DFT (e.g., circular shift, Parseval's theorem) with their mathematical expressions.

DFT–DTFT relationship

For a length-NN sequence x[n]x[n], the DTFT is the continuous-frequency spectrum

X(ejω)=n=0N1x[n]ejωn.X(e^{j\omega})=\sum_{n=0}^{N-1}x[n]e^{-j\omega n}.

The NN-point DFT is obtained by sampling the DTFT uniformly at NN equally-spaced frequencies ωk=2πkN\omega_k=\dfrac{2\pi k}{N}:

X[k]=X(ejω)ω=2πk/N=n=0N1x[n]ej2πkn/N,k=0,,N1.\boxed{X[k]=X(e^{j\omega})\Big|_{\omega=2\pi k/N}=\sum_{n=0}^{N-1}x[n]e^{-j2\pi kn/N},\quad k=0,\dots,N-1.}

Thus the DFT is the frequency-sampled DTFT over one period [0,2π)[0,2\pi); equivalently, it corresponds to the DTFT of the periodically-extended sequence. (It also equals the Z-transform evaluated at the NN equally-spaced points z=ej2πk/Nz=e^{j2\pi k/N} on the unit circle.)

Two important DFT properties

1. Circular shift. A circular time shift by mm introduces a linear phase:

x[(nm)modN]  DFT  X[k]ej2πkm/N.x\big[(n-m)\bmod N\big]\ \xleftrightarrow{\ \text{DFT}\ }\ X[k]\,e^{-j2\pi km/N}.

2. Parseval's theorem. Energy is conserved between time and frequency domains:

n=0N1x[n]2=1Nk=0N1X[k]2.\sum_{n=0}^{N-1}|x[n]|^{2}=\frac{1}{N}\sum_{k=0}^{N-1}|X[k]|^{2}.

(Other valid properties: linearity, circular convolution x[n]h[n]X[k]H[k]x[n]\circledast h[n]\leftrightarrow X[k]H[k], and symmetry X[Nk]=X[k]X[N-k]=X^{*}[k] for real x[n]x[n].)

dftfft

Frequently asked questions

Where can I find the BE Computer Engineering (IOE, TU) Digital Signal Analysis and Processing (IOE, EX 701 / ENEX 416) question paper 2078?
The full BE Computer Engineering (IOE, TU) Digital Signal Analysis and Processing (IOE, EX 701 / ENEX 416) 2078 (regular) question paper is available free on Kekkei. You can read every question online and attempt the paper under timed exam conditions.
Does the Digital Signal Analysis and Processing (IOE, EX 701 / ENEX 416) 2078 paper come with solutions?
Yes. Every question on this Digital Signal Analysis and Processing (IOE, EX 701 / ENEX 416) past paper includes a step-by-step solution, plus instant AI feedback when you attempt it on Kekkei.
How many marks is the BE Computer Engineering (IOE, TU) Digital Signal Analysis and Processing (IOE, EX 701 / ENEX 416) 2078 paper?
The BE Computer Engineering (IOE, TU) Digital Signal Analysis and Processing (IOE, EX 701 / ENEX 416) 2078 paper carries 80 full marks and is meant to be completed in 180 minutes, across 12 questions.
Is practising this Digital Signal Analysis and Processing (IOE, EX 701 / ENEX 416) past paper free?
Yes — reading and attempting this Digital Signal Analysis and Processing (IOE, EX 701 / ENEX 416) past paper on Kekkei is completely free.