Question 6 of 6: Reading a Kaiser-window highpass design backwards
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
Paper format. National Exams,
May 2017 — 16-Elec-B1 Digital Signal Processing. Three hours, closed book;
one of two approved calculators plus one double-sided aid sheet of tables and
formulas. Six questions are printed and five constitute a complete paper,
each worth 12 points; the marking scheme published on page 1 gives the per-part
split. All six questions are solved here, because the set is a
study resource rather than a timed attempt.
Reference texts.
A. V. Oppenheim and R. W. Schafer, Discrete-Time Signal Processing,
3rd ed., Pearson, 2010 — the source of the z-transform tables, the DFT
property list and the Kaiser-window formulas reproduced on pages 7–9 of
this exam.
J. G. Proakis and D. G. Manolakis, Digital Signal Processing:
Principles, Algorithms and Applications, 4th ed., Pearson, 2007.
S. K. Mitra, Digital Signal Processing: A Computer-Based Approach,
4th ed., McGraw-Hill, 2011 — filter structures and classical IIR
approximations.
Question 6: Reading a Kaiser-window highpass design backwards (12 marks)
Given. A completed Kaiser-window highpass design whose outputs are quoted, together with the two Kaiser design formulas.
Quantity
Symbol
Value
Cut-off frequency
$\omega_c$
$0.525\pi$ rad/sample
Kaiser shape parameter
$\beta$
$5.653$
Window order
$M$
$96.58 \cong 97$
Band type
—
highpass (stopband below, passband above)
Find. (a) the tolerance $\delta$ and the two band-edge frequencies $\omega_s$ and $\omega_p$; (b) what must still be changed so the design meets its specification at every frequency including $\omega = \pi$.
Approach. Both Kaiser formulas are invertible: $\beta$ gives $A$, $A$ gives $\delta$, and $A$ together with $M$ gives the transition width $\Delta\omega$. The windowed prototype places its cut-off at the midpoint of the transition band, so $\omega_c$ and $\Delta\omega$ fix the two band edges. Part (b) is then a question about the parity of $M$.
Invert the $\beta$ formula to recover $A$. The quoted $\beta = 5.653$ is large, so try the top branch, which is valid for $A \gt 50$: $$\beta = 0.1102(A - 8.7) \;\Longrightarrow\; A = \frac{\beta}{0.1102} + 8.7 = \frac{5.653}{0.1102} + 8.7 = 60.0 \ \text{dB}.$$ The result satisfies $A \gt 50$, so the branch choice was consistent. Substituting back, $0.1102(60 - 8.7) = 5.653$, which confirms the round trip.
Convert the attenuation into a tolerance. By definition $A = -20\log_{10}\delta$, so $$\boxed{\,\delta = 10^{-A/20} = 10^{-60/20} = 10^{-3} = 0.001\,}$$ The Kaiser method produces essentially equal ripple in both bands, so this one number is both the passband tolerance and the stopband tolerance: the passband must stay within $1 \pm 0.001$ and the stopband must not exceed $0.001$, i.e. 60 dB of attenuation.
Invert the order formula to recover the transition width. Using the unrounded value $M = 96.58$ (the rounded 97 would introduce a small error), $$\Delta\omega = \frac{A - 8}{2.285\,M} = \frac{60 - 8}{2.285 \times 96.58} = \frac{52}{220.69} = 0.2356\ \text{rad/sample}.$$ Expressed in units of $\pi$ this is a round number, $$\boxed{\,\Delta\omega = 0.075\pi\ \text{rad/sample}\,}$$ and substituting it back reproduces $M = 96.58$ exactly.
Place the band edges. A windowed design puts the prototype cut-off in the middle of the transition band, not at its edge, so $\omega_c = \tfrac{1}{2}(\omega_s + \omega_p)$ and $\Delta\omega = |\omega_p - \omega_s|$. For a highpass the stopband is the lower band, so $\omega_s \lt \omega_c \lt \omega_p$: $$\omega_s = \omega_c - \frac{\Delta\omega}{2} = 0.525\pi - 0.0375\pi, \qquad \omega_p = \omega_c + \frac{\Delta\omega}{2} = 0.525\pi + 0.0375\pi.$$
Report the specification. Carrying out the arithmetic, $$\boxed{\,\omega_s = 0.4875\pi = 1.5315\ \text{rad/sample}, \qquad \omega_p = 0.5625\pi = 1.7671\ \text{rad/sample}\,}$$ so the complete specification reads: $|H(e^{j\omega})| \le 0.001$ for $0 \le \omega \le 0.4875\pi$, and $0.999 \le |H(e^{j\omega})| \le 1.001$ for $0.5625\pi \le \omega \le \pi$, with the transition band left unconstrained.
Part (b): diagnose the failure at $\omega = \pi$. The design returned $M = 97$, an odd value, so the impulse response has length $M + 1 = 98$, an even number of samples. A symmetric response of even length is a Type II FIR filter, and its axis of symmetry $M/2 = 48.5$ falls between two samples. That is not a cosmetic detail: for a Type II filter $$H(e^{j\pi}) = \sum_{n=0}^{M} h[n](-1)^{n} = 0 \quad\text{identically},$$ because $n$ and $M - n$ have opposite parity whenever $M$ is odd, so the symmetric pairs $h[n]$ and $h[M-n]$ enter the sum with opposite signs and cancel term by term.
State the remedy. A highpass filter must have gain close to unity at $\omega = \pi$, so a Type II response can never meet the specification there no matter how large $\beta$ or $M$ becomes. What is required is a change of parity: $$\boxed{\text{take } M \text{ EVEN} \;-\; \text{here } M = 98, \text{ length } 99, \text{ a Type I filter}}$$ Type I has an integer symmetry point $M/2 = 49$, imposes no structural zero at $z = -1$, and keeps the same $\beta = 5.653$.
Confirm the rest of the specification still holds. Increasing $M$ from 96.58 to 98 can only help, because the transition width the design achieves is inversely proportional to $M$: $\Delta\omega = 52/(2.285 \times 98) = 0.2322$ rad $= 0.0739\pi$, which is narrower than the required $0.075\pi$. The stopband attenuation is controlled by $\beta$ alone and is unchanged at 60 dB. So the single change $M = 97 \to 98$ fixes the neighbourhood of $\pi$ while leaving every other requirement satisfied with a small margin to spare.
Figure Q6.1 — the highpass tolerance scheme recovered in part (a): stopband 0 to 0.4875 pi with tolerance delta = 0.001, transition band of width 0.075 pi centred on the cut-off 0.525 pi, and passband 0.5625 pi to pi within 1 +/- 0.001.
Figure Q6.2 — why part (b) is a parity question. With M odd the length M+1 is even, the symmetry point M/2 is a half-integer, and the symmetric pairs cancel at w = pi, forcing a zero of H(z) at z = -1.
Quantity
Symbol
Value
Stopband attenuation
$A$
$60$ dB
Ripple tolerance (both bands)
$\delta$
$0.001$
Transition band width
$\Delta\omega$
$0.075\pi = 0.2356$ rad/sample
Stopband corner
$\omega_s$
$0.4875\pi = 1.5315$ rad/sample
Passband corner
$\omega_p$
$0.5625\pi = 1.7671$ rad/sample
Filter type as designed
$M = 97$ (odd)
Type II — forces $H(e^{j\pi}) = 0$, unusable for a highpass
Required change
$M = 98$ (even)
Type I, length 99, same $\beta = 5.653$; achieved $\Delta\omega = 0.0739\pi$