Digital Filters

1. Digital Filters

2. x(t) A/D x[n] Computer y[n] D/A y(t) Example:

3. To get the current output y[n], which is sent to the D/A, we get the previous output y[n-1] (sent the last time to the D/A) and add to it the current input from the A/D. In a programming language like C++, the code would look something like this: yp = 0; // set “previous” y to zero while (true) { x = inp(ad); // get new sample y = yp+x; // perform filtering operation outp(y, da); // send filtered sample to D/A yp = y; // set “current” to “previous” // for next iteration }

4. Finding the z-transform transfer function:

5. Now what does this operation do to an input signal? Example: find the impulse and step responses to the previous filter. Solution: we could use z-transforms or convolution to find the output when the input is an impulse or a step, but instead, let us see how the output “evolves.” A table will be created with three columns: n, x[n] and y[n]. Each column will correspond to an instance in discrete-time. For the impulse input x[n]=d[n] we have

7. For the step input x[n]=u[n] we have

8. The plots of these outputs are as follows: y[n] for x[n]=d[n] n 1 2 3 4

9. y[n] for x[n]=u[n] 5 4 3 2 1 n 1 2 3 4

10. Example: Find the impulse response and the step response for the following transfer function: Solution: From the definition of the transfer function we have

11. Cross-multiplying and converting back to discrete-time domain we have In other words, the current output is equal to one-half times the previous output plus the current input.

13. For the step input x[n]=u[n] we have

14. The plots of these outputs are as follows: y[n] for x[n]=d[n] n 1 2 3 4

15. y[n] for x[n]=u[n] 2 1 n 1 2 3 4

16. The step and impulse functions for this transfer function look like those for an RC low-pass filter. Suppose that we wished to construct a digital filter corresponding to an RC low-pass filter? Can we convert an analog filter to a digital filter?

17. The Matched z-Transform In the matched z-transform digital filter design method we try to “match” the impulse response of the analog filter with that of the digital filter being designed. To match the impulse responses, we take the inverse Laplace transform of the analog filter H(s)h(t), then sample the impulse response h(t)h[n], then take the z-transform of the sampled impulse response to get the z-transform transfer function h[n]H(z).

18. Analog Prototype Digital Filter L-1 Z sample Once we have our z-transform transfer function H(z), we apply the definition of the transfer function to write our digital filter equations:

19. Example: Use the matched filter design method to design the digital equivalent of an integrator. Solution: The analog transfer function is The inverse Laplace transform is

20. We then sample the impulse response to get h[n]: Finally, we take the z-transform of the impulse response to get the digital filter transfer function.

21. Finally, we apply the definition of the z-transform transfer function to get the relationship between the input of the digital filter x[n] and the output of the digital filter y[n].

22. Applying this formula to an arbitrary input x[n], we have

23. We see that the output is the summationof the input. Thus the digital filter accurately represents the analog filter. The digital filter does not always accurately represent the analog filter as will be seen in the next example.

24. Example: Use the matched filter design method to design the following transfer function: where Wc = Ws/4, and Ws is the sampling frequency. Solution: First, we find the impulse response

25. Then we sample the impulse response: Then we take the z-transform of the (discrete-time) impulse response: