<< Chapter < Page Chapter >> Page >

Fir systems

Recall the system representation of discrete-time LTI systems:
Image
Implementation of a discrete-time LTI system.
FIR systems are those that do not have poles/feedback, so the representation can be simplified:
Image
Implementation of an FIR discrete-time LTI system.
From the straightforward time domain expression of this kind of system ($y[n] = b_0 x[n]+ b_1 x[n-1]+ b_2 x[n-2]+ \cdots b_M x[n-M]$) we have the following transfer function:$\begin{align*} H(z)= \frac{Y(z)}{X(z)}&= b_0 + b_1 z^{-1} + b_2 z^{-2} + \cdots b_M z^{-M}\\&= z^{-M}(z-\zeta_1)(z-\zeta_2) \cdots (z-\zeta_M) \end{align*}$ In the z-domain, such systems have only zeros (not counting poles at $z=0$, or in the case of acausal systems, poles at $z=\infty$).

Fir filters

As FIR filters are FIR discrete-time LTI systems, their design is simply a matter of adjusting the $b$ coefficients of the time-domain representation (or equivalently, the location of the zeros in the z- domain). As with IIR systems, these systems are designed to perform to particular specifications in the frequency domain:
Image
CAPTION.
FIR filters must be of much higher order (i.e., have more non-zero coefficients in their time-domain representations) than IIR filters to achieve similar performance, but in exhange for this cost have several benefits over IIR filters. First, an FIR filter is GUARANTEED to be BIBO stable. Recalling that an LTI system is BIBO stable if and only if its impulse response is absolutely summable, the finite length of FIR impulse responses guarantees this condition. And second, unlike the non-linearity with IIR filters, FIR filters can be easily designed to have a generalized linear phase response (a phase response that is linear, with one caveat; see below).

In order to have a generalized linear phase response, FIR filters of length $N$ must be designed so that their impulse responses $h[n]$ are either even or odd symmetric about the center time value(s) ($(N-1)/2$ for odd $n$, $n=N/2$ and $n=N/2-1$ for even $N$). This design can be attained for $N$ that are either even or odd. We will consider when $N$ is odd, i.e. when $N=M+1$ where $M$ is an even number. For such a filter length, evenness implies that $h[n]=h[M-n]$. Suppose the impulse response is designed so that such is the case. Let's take a look at the system's frequency response, using that symmetry to simplify things:$\begin{align*} H(\omega)&= \sum_{n=0}^{M} h[n]e^{-j\omega n}\\&= \sum_{n=0}^{M/2-1} h[n]e^{-j\omega n} + h[M/2]e^{-j\omega M/2} + \sum_{n=M/2+1}^{M} h[n]e^{-j\omega n} \\&= \sum_{n=0}^{M/2-1} h[n]e^{-j\omega n} + h[M/2]e^{-j\omega M/2} + \sum_{n=M/2+1}^{M} h[M-n]e^{-j\omega n} \\&= \sum_{n=0}^{M/2-1} h[n]e^{-j\omega n} + h[M/2]e^{-j\omega M/2} + \sum_{r=0}^{M/2-1} h[r]e^{-j\omega (M-r)} \\&= h[M/2]e^{-j\omega M/2} + \sum_{n=0}^{M/2-1} h[n]\left(e^{-j\omega n} + e^{j\omega (n-M)} \right)\\&= \left( h[M/2] + \sum_{n=0}^{M/2-1} 2 h[n]\cos(\omega (n-M/2)) \right) e^{-j\omega M/2}\\&= A(\omega) e^{-j\omega M/2} \end{align*}$Take a close look at the $A(\omega)$ term. It is called the filter's amplitude, because its amplitude is the same as the amplitude of the frequency response: $|A(\omega)|=|H(\omega)|$. It is also purely real-valued. Thus the frequency response of the filter is a real number multiplied by an exponential function $e^{-j\omega M/2}$ that has linear phase. This means that the frequency response has a perfectly linear phase, except when the sign of $A(\omega)$ changes. However, since $A(\omega)$ will hover around the value of $1$ in the filter's pass-band, the filter will have linear phase in the pass-band.

Creating fir filters

For an FIR filter of odd-length ($N=M+1$, where $M$ is even) and even symmetry about $n=\frac{M}{2}$ ($h[n]=h[M-n]$) we have that the filter's frequency response is: $H(\omega)=A(\omega) e^{-j\omega M/2}$where $A(\omega)=h[M/2]~+~ \sum_{n=0}^{M/2-1} 2 \, h[n] \cos(\omega (n-M/2))$Because (generalized) linear phase is desirable, we can focus on meeting filter specifications in terms of modifying $A(\omega)$:
Image
FIR filter design amounts to creating an $A(\omega)$ to fit particular specifications.
The goal will be to meet the requirements with as small of a filter size as possible. To that end, there are two pieces of information that will be very useful for us. Both involve the "ripples" of $A(\omega)$ centered around $A(\omega)=1$ in the pass-band and $A(\omega)=0$ in the stop-band. It turns out that among all filters that meet some defined specifications on $A(\omega)$ (the error bounds in the pass- and stop-bands and the pass and stop frequencies), the filter of shortest length $M+1$ will: --be an equiripple filter , meaning that all of the ripple oscillations will be of the same deviation about $A(\omega)=1$ in the pass-band and $A(\omega)=0$ in the stop-band (in this sense, the "spread" equally the error out over the entire frequency range) --will touch theerror bounding box $\frac{M}{2}+2$ times Putting those properties together, we can see that there are two ways of going about FIR equiripple filter design. One way might be to specify the desired filter length $N=M+1$ (or, equivalently, filter order $M$), and then use the properties to create an equiripple filter based upon pass-band and stop-band frequencies; because of the first property, this filter will minimize the maximum deviation from $1$ in the pass-band and from $0$ in the stop-band. To design such a filter, one could use the command firpm in MATLAB. Another approach could be to specify these desired deviations, and then go about finding the filter of minimum length that will satisfy these deviations. To do that, one could use the firpmord command in MATLAB.

Equiripple filter of length $n=21$

Suppose we wanted to create a low-pass filter of length $N=21$ (order $M=20$), with a pass-band from $\omega=0$ to $\omega=.3\pi$ and a stop-band from $\omega=.35\pi$ to $\omega=\pi$. Among all possible FIR filters of length $N=21$, an FIR equiripple filter will have the smallest maximum deviation from the optimal frequency response magnitudes (of $|H(\omega)|=1$ in the pass-band and $|H(\omega)|=0$ in the stop-band). According to the equiripple filter properties, its frequency response will touch these error bounding boxes $\frac{M}{2}+2=12$ times.

Image
CAPTION.
Image
CAPTION.
CAPTION.
The equiripple filter is also designed, as described above, to have an impulse response that has even symmetry about its $n=10$ center time value. As a result, it will have a generalized linear phase:
Image
CAPTION.
Image
CAPTION.
CAPTION.
Finally, as this filter is an FIR filter, it will not have any poles (except for those at $z=0$ and, if it were acausal, at $z=\infty$). Since its order is $M=20$, it does have $20$ zeros:
Image
CAPTION.

Equiripple filter of length $n=101$

The filter in the first example, while optimal for its length, still has considerable deviation from the ideal frequency response of $1$ in the pass-band and $0$ in the stop-band. Having a more ideal response is straightforward, simply increase the length of the filter; as described above, the MATLAB command firpmord could be used to tell us the minimum length required to achieve a certain margin of tolerance from the ideal response. Suppose we had some requirement in mind, and that command told us a filter length of $101$ was required. The resulting equiripple filter of that length is below.

Image
CAPTION.
Image
CAPTION.
Image
CAPTION.
CAPTION.

Questions & Answers

I'm interested in biological psychology and cognitive psychology
Tanya Reply
what does preconceived mean
sammie Reply
physiological Psychology
Nwosu Reply
How can I develope my cognitive domain
Amanyire Reply
why is communication effective
Dakolo Reply
Communication is effective because it allows individuals to share ideas, thoughts, and information with others.
effective communication can lead to improved outcomes in various settings, including personal relationships, business environments, and educational settings. By communicating effectively, individuals can negotiate effectively, solve problems collaboratively, and work towards common goals.
it starts up serve and return practice/assessments.it helps find voice talking therapy also assessments through relaxed conversation.
miss
Every time someone flushes a toilet in the apartment building, the person begins to jumb back automatically after hearing the flush, before the water temperature changes. Identify the types of learning, if it is classical conditioning identify the NS, UCS, CS and CR. If it is operant conditioning, identify the type of consequence positive reinforcement, negative reinforcement or punishment
Wekolamo Reply
please i need answer
Wekolamo
because it helps many people around the world to understand how to interact with other people and understand them well, for example at work (job).
Manix Reply
Agreed 👍 There are many parts of our brains and behaviors, we really need to get to know. Blessings for everyone and happy Sunday!
ARC
A child is a member of community not society elucidate ?
JESSY Reply
Isn't practices worldwide, be it psychology, be it science. isn't much just a false belief of control over something the mind cannot truly comprehend?
Simon Reply
compare and contrast skinner's perspective on personality development on freud
namakula Reply
Skinner skipped the whole unconscious phenomenon and rather emphasized on classical conditioning
war
explain how nature and nurture affect the development and later the productivity of an individual.
Amesalu Reply
nature is an hereditary factor while nurture is an environmental factor which constitute an individual personality. so if an individual's parent has a deviant behavior and was also brought up in an deviant environment, observation of the behavior and the inborn trait we make the individual deviant.
Samuel
I am taking this course because I am hoping that I could somehow learn more about my chosen field of interest and due to the fact that being a PsyD really ignites my passion as an individual the more I hope to learn about developing and literally explore the complexity of my critical thinking skills
Zyryn Reply
good👍
Jonathan
and having a good philosophy of the world is like a sandwich and a peanut butter 👍
Jonathan
generally amnesi how long yrs memory loss
Kelu Reply
interpersonal relationships
Abdulfatai Reply
What would be the best educational aid(s) for gifted kids/savants?
Heidi Reply
treat them normal, if they want help then give them. that will make everyone happy
Saurabh
Got questions? Join the online conversation and get instant answers!
Jobilize.com Reply

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Discrete-time signals and systems. OpenStax CNX. Oct 07, 2015 Download for free at https://legacy.cnx.org/content/col11868/1.2
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Discrete-time signals and systems' conversation and receive update notifications?

Ask