38
OK, I'll try to answer your questions:
Q1: the number of taps is not equal the to the filter order. In your example the filter length is 5, i.e. the filter extends over 5 input samples [$x(n), x(n-1), x(n-2), x(n-3), x(n-4)$]. The number of taps is the same as the filter length. In your case you have one tap equal to zero (the coefficient for $x(n-1)$), so ...
28
When choosing one of these 4 types of linear phase filters there are mainly 3 things to consider:
constraints on the zeros of $H(z)$ at $z=1$ and $z=-1$
integer/non-integer group delay
phase shift (apart from the linear phase)
For type I filters (odd number of taps, even symmetry) there are no constraints on the zeros at $z=1$ and $z=-1$, the phase shift ...
21
Digital filter design is a very large and mature topic and - as you've mentioned in your question - there is a lot of material available. What I want to try here is to get you started and to make the existing material more accessible. Instead of digital filters I should actually be talking about discrete-time filters because I will not consider coefficient ...
21
A linear phase filter will preserve the waveshape of the signal or component of the input signal (to the extent that's possible, given that some frequencies will be changed in amplitude by the action of the filter).
This could be important in several domains:
coherent signal processing and demodulation, where the waveshape is important because a ...
21
Let me add the following graphic to the great answers already given.
When a filter has linear phase, then all the frequencies within that signal will be delayed the same amount in time (as described mathematically in Fat32's answer).
Any signal can be decomposed (via Fourier Series) into separate frequency components. When the signal gets delayed through ...
20
FIR filters contain as many poles as they have zeros. but all of the poles are located at the origin, $z=0$.
because all of the poles are located inside the unit circle, the FIR filter is ostensibly stable.
this is probably not the FIR filter the OP is thinking about, but there is a class of FIR filters called Truncated IIR filters (TIIR) which may have a ...
answered Jan 7 '14 at 15:40
robert bristow-johnson
13.5k3 gold badges23 silver badges58 bronze badges
16
I agree that the windowing filter design method is not one of the most important design methods anymore, and it might indeed be the case that it is overrepresented in traditional textbooks, probably due to historical reasons.
However, I think that its use can be justified in certain situations. I do not agree that computational complexity is no issue ...
14
"Can you mathematically show that FIR filters have poles, because I'm not seeing it." – Jim Clay
can we assume this FIR is causal?
filter order is $N$. number of taps is $N+1$
the Finite Impulse Response: $ \quad h[n] = 0 \quad \forall \quad n>N, \ n<0$
transfer function of the FIR:
$$ \begin{align}
H(z) & = \sum_{n=-\infty}^{+\infty} h[n] ...
answered Jan 8 '14 at 5:12
robert bristow-johnson
13.5k3 gold badges23 silver badges58 bronze badges
14
Just to add to what's already been said, you can see this intuitively by looking at the following sinusoid with monotonically increasing frequency.
Shifting this signal to the right or left will change its phase. But note also that the phase change will be larger for higher frequencies, and smaller for lower frequencies. Or in other words, the phase ...
13
My favorite "Rule of thumb" for the order of a low-pass FIR filter is the "fred harris rule of thumb":
$N=[f_s/delta(f)]*[atten(dB)/22]$
where
delta(f) is the transition band, in same units of $f_s$
$f_s$ is the sample rate of the filter
atten(dB) is the target rejection in dB
For example if you have a transition band of 100 Hz in a system sampled at 1KHz,...
12
The essence and importance of linear phase property lies in the definition and the effect of group delay $$\tau(\omega) = - \frac {d\phi(\omega)}{d\omega}$$ on the applied signal $x[n]$, where $\phi(\omega)$ is the phase response of the filter; (phase of its frequency response).
Assume that a filter, with a fixed group delay of $n_0$ samples, is applied a ...
12
The answer to this question is already been explained clearly in the previous replies. Yet I wish to give it a try to present a mathematical interpretation of the same
Consider a linear time invariant System whose frequency response is governed by $H(w)$.
i.e if the input to this system is $e^{jw_{0}t}$ the output will be $H(w_{0})e^{jw_{0}t}$
Here $H(w_{...
11
Note that for stable IIR filters, the impulse response does approach zero as $n$ goes to infinity. It just never becomes exactly zero. However, the sum of the absolute values is finite. Just as an example, take the exponential impulse response
$$h[n]=a^nu[n],\qquad |a|<1\tag{1}$$
where $u[n]$ is the unit step function. The sum
$$\sum_{n=-\infty}^{\...
11
See How many taps does an FIR filter need?
In your case you'd need more than 1000 taps depending on the allowable ripple, as your cut-off frequency is less than fs/500.
Alternatives :
use an IIR, a simple order-1 DC removal filter could work great
Average your signal and subtract the average in order to remove the DC
Rick Lyons proposes a clever ...
10
You are correct. FFT based processing adds inherent latency to your system. However there are ways to tweak this.
Let's assume you have an FIR filter of length "N". This can be implement FFT-based using the standard overlap add or overlap save method, where the FFT length would be 2*N. Overall system latency will also be roughly 2*N: you need to accumulate ...
10
In reviewing fred harris Figures of Merit for various windows (Table 1 in this link) the Hamming is compared to the Hanning (Hann) at various values of $\alpha$ and from that it is clear that the Hanning would provide greater stopband rejection (The classic Hann is with $\alpha =2$ and from the table the side-lobe fall-off is -18 dB per octave). I provided ...
10
Note the difference between the zeros at $0.3 \pi$ and at $0.8 \pi$.
The first one is clearly a zero crossing, much like $abs(x)$ at $x=0$.
At $\theta = 0.8 \pi$, however, the curve is tangent to the horizontal axis, much like $x^2$ at $x=0$. So you have a doulbe zero here.
So your zeros are:
2 zeros at $z = e^{\pm j 0.3 \pi}$
2 double zeros at $z = e^{\...
8
FIR filters contain only zeros and no poles. If a filter contains poles, it is IIR. IIR filters are indeed afflicted with stability issues and must be handled with care.
EDIT:
After some further thought and some scribbling and google-ing, I think that I have an answer to this question of FIR poles that hopefully will be satisfactory to interested parties....
7
I would say that the answer to your question - if taken literally - is 'no', there is no general way to simply convert an FIR filter to an IIR filter.
I agree with RBJ that one way to approach the problem is to look at the FIR filter's impulse response and use a time domain method (such as Prony's method) to approximate that impulse response by an IIR ...
7
Windowed Sinc filters can be adaptively generated on the fly on processors barely powerful enough to run the associated FIR filter. Windowed Sinc filters can be generated in finite bounded time.
The generation of some simple windowed Sinc filters can be completely described (and inspected for malware, etc.) in a few lines of code, versus blind use of some ...
7
What you do in step 1 is simply truncate the infinite impulse response to approximate it by an FIR filter. If you use sufficiently many filter taps, the approximation becomes arbitrarily accurate. This means that the resulting FIR filter approximates the magnitude and the phase characteristic of the original IIR filter. So with this approach the phase will ...
6
If the Z-transform of the feedforward section is divisible by the Z-transform of the feedback section, the filter is FIR.
Consider your example: $y[n] = y[n-1] + x[n] - x[n-3]$.
The Z-transform is $\mathrm Y(z)- z^{-1}\mathrm Y(z) = \mathrm X(z) - z^{-3}\mathrm X(z)$, and the Z-transform of the response is $\mathrm H(z) = \mathrm Y(z)/\mathrm X(z) = (1 - z^...
6
To be precise the group delay of a linear phase FIR filter is $(N-1)/2$ samples, where $N$ is the filter length (i.e. the number of taps). The group delay is constant for all frequencies, because the filter has a linear phase, i.e. its impulse response is symmetrical (or asymmetric). A linear phase means that all frequency components of the input signal ...
6
In the general case you have
$$H(z)=\frac{P(z)}{Q(z)}$$
where $P(z)$ and $Q(z)$ are polynomials in $z$. If - as is the case in your example - $Q(z)$ just has one single term, $H(z)$ is definitely FIR, because you can simply divide each term of $P(z)$ by that respective power of $z$, and the number of terms of $H(z)$ equals the number of terms of $P(z)$.
...
6
Your questions still leave me wondering as to what you're actually designing. For software implementation on modern x86 CPUs, CICs make almost no sense, but they are extremely elegant in hardware.
These filter definitions are ridiculous if you're planning to use a FIR – a transition width of 1mHz means that a minimum phase equiripple filter [1, (5.75), p. ...
6
I'll show here one benefit of a windowed design and a trick to get the same benefit from Parks–McClellan.
For half-band, quarter-band etc. filters windowing retains the time-domain zeros of the scaled sinc function, which is the prototypical ideal low-pass filter. The zeros end up in the coefficients, reducing the computational cost of the filters. For a ...
6
Lets say we want to transmit a sequence of discrete data $\left\lbrace x[n] \right\rbrace$. But because we are living in analog world, the sequence must be modulated.
Call $T_s$ is symbol duration and use a set of orthonormal waveforms $\left\lbrace p_n(t) = p(t-nT_s), n \in \mathbb{Z} \right\rbrace$, (baseband) signal $x(t)$ can be written as
\begin{...
6
A moving average filter can be thought of as a type of low-pass filter that doesn't have any control over its bandwidth for a fixed number of taps.
For a finite impulse response (FIR) filter, the output signal $y[n]$ is given in terms of the input signal $x[n]$ and the filter taps $h[n]$:
$$
y[k] = \sum_{n=0}^{N-1}h[n]x[k-n].
$$
The filter length in this ...
6
There are several reasons why the two results don't match:
the coefficients of the FIR Hilbert transformer are wrong
the FIR Hilbert transformer is too short to even come close to the performance of the FFT-based implementation
the frequency of the input signal is too low for the FIR Hilbert transformer to perform properly. A FIR Hilbert transformer always ...
6
A causal first-order IIR filter is characterized by the following difference equation:
$$y[n]=b_0x[n]+b_1x[n-1]-a_1y[n-1]\tag{1}$$
with $x[n]$ the input signal, and $y[n]$ the output signal. The impulse response of that system can be computed via the $\mathcal{Z}$-transform or otherwise, and it turns out to be
$$h[n]=b_0\delta[n]+(-a_1)^{n-1}(b_1-b_0a_1)u[...
Only top voted, non community-wiki answers of a minimum length are eligible
