Take the 2-minute tour ×
Signal Processing Stack Exchange is a question and answer site for practitioners of the art and science of signal, image and video processing. It's 100% free, no registration required.

I would like to know whether this can be done and if then how would this feat be acomplished in Matlab given a filter calculated with Remez?

h = remez(...);

If this is not possible then I would appreciate an explanation whether this can be done approximately (approximations of the zeros and poles) or whether this is simply mathematically pointless because a broken rational complex polynomial representation does not exist.

share|improve this question

2 Answers 2

up vote 2 down vote accepted

If your impulse response is $h(n)$ then the $\mathcal{Z}$-transform of this impulse response (i.e. the transfer function of the filter) is

$$H(z)=\sum_{n=0}^{N-1}h(n)z^{-n}$$

where $N$ is the filter length, i.e. the number of coefficients.

Note that since it is an FIR filter, all its poles are at the origin of the $z$-plane. The zeros are of course the zeros of the polynomial with coefficients $h(n)$.

share|improve this answer
    
ŵell explained, it's ringing a bell. thanks –  Raffael May 13 at 18:52
    
are all zeros located within the unit-circle? –  Raffael May 13 at 18:55
1  
@Raffael: No, the Remez algorithm designs linear-phase FIR filters, and they have all their zeros either on the unit circle (in the stopband of the filter) or mirror-imaged at the unit circle (in the passband). I.e. for each zero inside the circle there must be one with inverse magnitude outside the unit circle. –  Matt L. May 13 at 18:58

Well, remez function in MATLAB returns your FIR filter coefficients b (of a numerator). Having them you can obtain zeros, poles and gain by calling:

[z, p, k] = tf2zpk(b, 1);

What's more, you can use the following tool to visualize your filter:

 fvtool(b, 1);

In R you can use the function polezero:

 polezero(b, 1); 

to get the plot.

share|improve this answer
    
+1: thanks for the matlab formula –  Raffael May 13 at 18:53

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.