NivaarExam PrepOfficial exam papers ↗

22-Elec-B1 Digital Signal Processing · May 2015

Question 4 of 6: Parallel and cascade realisations of an IIR system (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 4: Parallel and cascade realisations of an IIR system (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 causal system already written as the sum of a first-order and a second-order section.

Given data
First section$H_1(z)=\dfrac{1}{1-\tfrac12 z^{-1}}$
Second section$H_2(z)=\dfrac{1+2z^{-1}}{1+\tfrac34 z^{-1}-\tfrac14 z^{-2}}$
System classcausal, so the ROC is the exterior of the outermost pole

Find. (a) a stability verdict with reasons; (b) a parallel signal flow graph using direct form II for the second-order section; (c) a cascade signal flow graph of a first-order section and a transposed direct form II second-order section; (d) a comparison of multiplier and delay counts and a recommendation.

Approach. Factor the second-order denominator to expose all three poles, which settles part (a); draw the parallel structure directly from the given decomposition; combine the two fractions over a common denominator to obtain the single numerator needed for the cascade; then count multipliers and delays in each drawing.

  1. Locate all three poles (part a). The first section contributes a pole at $z=\tfrac12$. For the second, multiply the denominator by $z^{2}$: $$z^{2}+\tfrac34 z-\tfrac14=0 \quad\Longrightarrow\quad z=\frac{-\tfrac34\pm\sqrt{\tfrac{9}{16}+1}}{2}=\frac{-0.75\pm 1.25}{2},$$ giving $z=\tfrac14$ and $z=-1$. Equivalently, in $z^{-1}$, $$1+\tfrac34 z^{-1}-\tfrac14 z^{-2}=\left(1+z^{-1}\right)\left(1-\tfrac14 z^{-1}\right).$$
  2. Apply the stability test (part a). The complete pole set is $\left\{\tfrac12,\ \tfrac14,\ -1\right\}$. For a causal system the ROC is $|z| \gt \max|p| = 1$, which does not contain the unit circle because the pole at $z=-1$ lies exactly on it. Hence $$\boxed{\,\text{poles at } z=\tfrac12,\ \tfrac14,\ -1 \ \Longrightarrow\ \text{the system is NOT stable}\,}$$ The pole at $z=-1$ is a marginal, undamped mode: it contributes a term $(-1)^{n}$ to $h[n]$, which neither decays nor grows, so $h[n]$ is not absolutely summable and a bounded input at $\omega=\pi$ produces an unbounded output.

Part (b): parallel form. The system function is already a sum of two sections, so the parallel structure is drawn directly from the given expression. The first-order branch realises $$y_1[n]=x[n]+\tfrac12\,y_1[n-1],$$ and the second-order branch, in direct form II, splits into the recursive and non-recursive halves $$w[n]=x[n]-\tfrac34\,w[n-1]+\tfrac14\,w[n-2],\qquad y_2[n]=w[n]+2\,w[n-1],$$ with $y[n]=y_1[n]+y_2[n]$. Direct form II shares one delay chain between the two halves, which is why only two delays appear in that branch.

x[n]y1[n]z^-11/2first-order section H1(z) = 1/(1 - 0.5 z^-1)w[n]y2[n]z^-1z^-1-3/4+1/42second-order section, Direct Form II (b0 = 1 needs no multiplier)y[n]Parallel form: 4 multipliers, 3 delay elements
Figure 4.1 — parallel realisation. The first-order branch needs one multiplier and one delay; the direct-form-II second-order branch needs three multipliers and two delays, because the leading numerator coefficient b0 = 1 is a plain wire.

Part (c): cascade form. A cascade needs the two sections multiplied rather than added, so the two fractions must first be placed over a common denominator.

  1. Combine the two sections into one rational function. Using the factorisation from step 1, the common denominator is $$D(z)=\left(1-\tfrac12 z^{-1}\right)\left(1+z^{-1}\right)\left(1-\tfrac14 z^{-1}\right),$$ and the numerator is $$\begin{aligned} N(z)&=\left(1+\tfrac34 z^{-1}-\tfrac14 z^{-2}\right)+\left(1+2z^{-1}\right)\left(1-\tfrac12 z^{-1}\right)\\ &=\left(1+\tfrac34 z^{-1}-\tfrac14 z^{-2}\right)+\left(1+\tfrac32 z^{-1}-z^{-2}\right)\\ &=2+\tfrac94 z^{-1}-\tfrac54 z^{-2}. \end{aligned}$$
  2. Split the result into a first-order and a second-order factor. Keeping the numerator whole with the second-order denominator gives the cleanest pairing: $$\boxed{\,H(z)=\underbrace{\frac{1}{1-\tfrac12 z^{-1}}}_{\text{section A, 1st order}}\times\underbrace{\frac{2+\tfrac94 z^{-1}-\tfrac54 z^{-2}}{1+\tfrac34 z^{-1}-\tfrac14 z^{-2}}}_{\text{section B, 2nd order}}\,}$$ For reference the numerator's own roots are $z=\left(-9\pm\sqrt{241}\right)/16$, that is $z\approx 0.4078$ and $z\approx -1.5328$; they are irrational, so pairing the whole numerator with the second-order denominator avoids introducing ugly section coefficients for no benefit.
  3. Write the transposed direct form II recursion for section B. With $b_0=2$, $b_1=\tfrac94$, $b_2=-\tfrac54$ and $a_1=-\tfrac34$, $a_2=\tfrac14$ (using $H=\left(b_0+b_1z^{-1}+b_2z^{-2}\right)/\left(1-a_1z^{-1}-a_2z^{-2}\right)$), the transposed structure computes $$\begin{aligned} y[n] &= b_0\,v[n]+s_1[n-1],\\ s_1[n] &= b_1\,v[n]+a_1\,y[n]+s_2[n-1],\\ s_2[n] &= b_2\,v[n]+a_2\,y[n], \end{aligned}$$ where $v[n]$ is the output of section A. Both realisations were checked: running a unit impulse through each structure reproduces the impulse response of the original $H(z)$ sample for sample.
x[n]z^-11/2section A: 1/(1 - 0.5 z^-1)v[n]section B: transposed Direct Form II22.25-1.25z^-1z^-1y[n]-3/4+1/4Cascade form: 6 multipliers, 3 delay elementsevery coefficient of section B needs its own multiplier, including b0 = 2
Figure 4.2 — cascade realisation. Section A is the same first-order block; section B is a transposed direct form II, in which the delayed state variables climb from the lower adders into the upper ones and every one of the five coefficients needs its own multiplier.

Part (d): comparison. Counting a multiplier only for coefficients that are not $\pm 1$ (a unity coefficient is a plain wire), and one storage element per $z^{-1}$ block:

Multiplier and storage count
StructureSectionMultipliersDelays
Parallel (b)1st order: $\tfrac12$11
2nd order DF II: $-\tfrac34,\ \tfrac14,\ 2$ ($b_0=1$ is a wire)32
total43
Cascade (c)1st order: $\tfrac12$11
2nd order TDF II: $2,\ \tfrac94,\ -\tfrac54,\ -\tfrac34,\ \tfrac14$52
total63

Both structures are canonic in storage — three delays for a third-order system, which is the theoretical minimum — so storage does not separate them. The multiplier counts do: $$\boxed{\,\text{parallel: 4 multiplies, 3 delays}\quad\text{vs}\quad\text{cascade: 6 multiplies, 3 delays}\ \Rightarrow\ \text{the PARALLEL form is preferable}\,}$$ The saving arises because the parallel decomposition was handed to us with a unity leading numerator coefficient in each branch, whereas forming the common numerator for the cascade produced three non-trivial feed-forward coefficients. Two secondary arguments reinforce the same choice: in a parallel structure the round-off noise generated in one section is not filtered by the other sections, so noise does not accumulate along a chain, and the two branches can be evaluated concurrently in a pipelined processor. The one advantage the cascade retains is that its zeros are controlled explicitly section by section, which matters when the zero locations must be held accurately under coefficient quantisation.

Final results
PartQuantityResult
(a)Poles$z=\tfrac12$, $z=\tfrac14$, $z=-1$
(a)StabilityNOT stable — the pole at $z=-1$ lies on the unit circle, so the causal ROC $|z| \gt 1$ excludes $|z|=1$
(b)Parallel structure$H_1$ direct + $H_2$ in direct form II; 4 multipliers, 3 delays
(c)Cascade structure$\dfrac{1}{1-\tfrac12 z^{-1}}\times\dfrac{2+\tfrac94 z^{-1}-\tfrac54 z^{-2}}{1+\tfrac34 z^{-1}-\tfrac14 z^{-2}}$; 6 multipliers, 3 delays
(d)Preferred structureParallel — same storage, two fewer multiplications, and no round-off accumulation between sections