02.10.2019 Views

UploadFile_6417

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Analysis of A/D Quantization Noise 549<br />

where R 0 is the constant term and R k ’s are the residues at the pole<br />

locations p k . This expansion can be computed using the residue function.<br />

Note that both poles and the corresponding residues are either real-valued<br />

or occur in complex-conjugate pairs. Then using (10.18a), we can show<br />

that (see [17] and also Problem P10.3)<br />

σq<br />

2<br />

σe<br />

2<br />

= R 2 0 +<br />

N−1<br />

∑<br />

k=1<br />

N−1<br />

∑<br />

l=1<br />

R k Rl<br />

∗<br />

1 − p k p ∗ l<br />

(10.22)<br />

The variance-gain expression in (10.22) is applicable for most practical<br />

filters since rarely do they have multiple poles. The approximate value of<br />

the variance-gain for IIR filters is given by<br />

σq<br />

2<br />

σe<br />

2<br />

≃<br />

K−1<br />

∑<br />

k=0<br />

|h(n)| 2 , K ≫ 1 (10.23)<br />

where K is chosen so that the impulse response values (magnitudewise)<br />

are almost zero beyond K samples. The following MATLAB function,<br />

VarGain, computes variance-gain using (10.19) or (10.22).<br />

function Gv = VarGain(b,a)<br />

% Computation of variance-gain for the output noise process<br />

% of digital filter described by b(z)/a(z)<br />

% Gv = VarGain(b,a)<br />

a0 = a(1); a = a/a0; b = b/a0; M = length(b); N = length(a);<br />

if N == 1<br />

% FIR Filter<br />

Gv = sum(b.*b);<br />

return<br />

else<br />

% IIR Filter<br />

[R,p,P] = residue(b,a);<br />

if length(P) > 1<br />

error(’*** Variance Gain Not computable ***’);<br />

elseif length(P) == 1<br />

Gv = P*P;<br />

else<br />

Gv = 0;<br />

end<br />

Rnum = R*R’; pden = 1-p*p’;<br />

H = Rnum./pden; Gv = Gv + real(sum(H(:)));<br />

end<br />

It should be noted that the actual output noise variance is obtained by<br />

multiplying the A/D quantization noise variance by the variance-gain.<br />

□ EXAMPLE 10.4 Consider an 8-order IIR filter with poles at p k = r e j2πk/8 , k =0,...,7. If r is<br />

close to 1, then the filter has 4 narrowband peaks. Determine the variance-gain<br />

for this filter when r =0.9 and r =0.99.<br />

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).<br />

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.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!