This utility program can be a handy tool for anyone working in the field of signal processing with interest in implementation of the project through C++ programming. Although this work can be easily implemented in Matlab using simple inbuilt functions, the real challenge of designing and understanding the concept of filter design will be missed out. In order to understand the depth of designing an FIR filter using window functions, coding in C++ is a must. Described here is a C++ implementation of finite impulse response (FIR) filters using Blackman window method.
The term filter is often used to denote systems that process the input in a specified way. In this context, filtering describes a signal-processing operation that allows signal enhancement, noise reduction or increased signal-to-noise ratio. Systems for the processing of discrete time signals are also called digital filters. Depending on the requirement and application, the analysis of a digital filter may be carried out in time domain, z-domain or frequency domain.
A common application of digital filters is to modify the frequency response is some specified way. An ideal low-pass filter is one of the common examples of digital filters which passes the frequencies up to a specified value and blocks other frequencies. One way to approximate the ideal low-pass response is to do truncation of the impulse response by multiplication with the window functions. Mathematically, it can be expressed as follows:
where hd(n) is an infinite-duration impulse response and w(n) is the window function. There are various windows available in signal processing literature. However, in this article we specifically discuss Blackman window. In the case of Blackman window, w(n) is given as:
Plot of eqn(2) can be seen graphically as shown in Fig. 1.
Here, we have developed C++ based software for designing FIR filters using Blackman window to get frequency response. The resultant software has the following salient features:
1. Welcome screen of Blackman window design
2. It prompts you to enter filter’s type
3. It prompts you to enter window length (N) from 1 to 200
4. It prompts you to enter the cutoff frequency. The software is designed to accept the values in radians between 0.2 and 2.6
5. The program outputs normalised filter response. This is a filter response which is designed to pass through a cutoff frequency of 1 radian
6. The program outputs actual filter response. This is a filter response which is designed to pass through a cutoff frequency specified by the user
The software is designed to follow the structure as shown in Fig. 2, where the details of each file are explained as given below.
1. MAIN.CPP. This file contains the menu-driven program for user inputs
2. FUNC.H. It is a header file containing code related to Blackman window coefficients calculation
3. BLACK.H. It is a header file containing code for supporting C++ graphics
4. GRAPH.H. It is a header file again to support graphic content in C++
Output responses of the filter
Now, the results obtained after running the software is presented for Blackman window for a low-pass FIR filter for different window length N to observe the effect of N on frequency response. All values of |H(jω)| are in decibels. A normalised frequency scale is used.
Low-pass filter. Simulation of low-pass filter using Blackman window is done through the software. Figs 3 to 5 are the outputs of the low-pass filter for different window length of N=49, 99 and 149, respectively.
From Figs 3 to 5, we have following observations when N increases from 49 to 149 at cutoff frequency (fc)=1.2:
1. The widths of the pass band and side lobes decrease
2. The slope of the transition band increases, that is, the band develops a sharper cutoff
3. As N increases, the number of ripples increase in the pass band
4. The amplitude and shape of the end-ripple remains the same but it gets shifted towards the edge and finally out of the pass band for larger values of N
5. As N increases, the ripples in the pass band get smoothed out. A key difference in the observed responses of the rectangular is the absence of the Gibbs phenomenon. This owes itself to the smoother shape of the Blackman window in the time-domain, with the notable absence of a discontinuity that would lead to the Gibbs phenomenon.
6. As N increases the widths of the lobes decrease in the stop band
In the next section we present the outputs of high-pass, band-reject and band-pass filters.
Frequency response of high-pass filter with N=45, fc=1.2 is shown in Fig. 6. Frequency response of band-reject filter with N=67, with lower cutoff frequency (fl)=1.2 and upper cutoff frequency (fu)=2.0, is shown in Fig. 7. Frequency response of band-pass filter with N=91, fl=1.2 and fu=2.0 is shown in Fig. 8.
The program was tested on Windows XP, 32-bit PC using Turbo C++ version 3. For running Turbo C++ on 64-bit PC with Windows 7 or Windows 8, a separate Turbo C ++ IDE is available for download. You need to copy the source code and relevant header files into bin folder before running the program.
Download source code: click here
Digital signal processing finds application in almost every field of engineering. Its impact on consumer electronics is evident by the advent of high-definition televisions and digital imaging. Its application in the field of biomedical and signal processing includes enhancement of the quality of signals under interest using various digital filters. These filters can also be used in applications like noise removal. The filter routines can be easily coded to hardware description language like Verilog or VHDL and can turn out to be a good project for students or designers working in real-time image and video processing.
The author is an M.Tech from ABV-Indian Institute of Information Technology and Management, Gwalior, and B.Tech from Jamia Millia Islamia in New Delhi