Design of Infinite Impulse Response Filters Based on Multi-Objective Particle Swarm Optimization
Abstract
:1. Introduction
2. Methodologies and System Description
2.1. Low-Pass Filter
2.2. High-Pass Filter
2.3. Band-Pass Filter
2.4. Butterworth Filter
- Maximal Flatness
- 2.
- Transition Band Performance
- 3.
- Flexible Order Design
- is a complex frequency variable (i.e., Laplacian variable), is the angular frequency. is the cutoff frequency. is the order of the filter. G is the gain, usually set to 1, so that when , the value of is .
2.5. Chebyshev Filter
- Type I Chebyshev Filter
- 2.
- Type II Chebyshev Filter
- Equations (8) and (9) are only true when , is the ripple coefficient, is the cutoff frequency, and is the order of the filter. is a Chebyshev polynomial of type I, and is a Chebyshev polynomial of type II.
2.6. Particle Swarm Optimization (PSO) Algorithm
- is the velocity of particle at the next time step.
- is the inertia weight, which controls the degree of retention of particle speed.
- is the speed of particle i at the current time step.
- and are learning factors, which are called personal learning factors and social learning factors, respectively. They are used to adjust the tendency of particles to move to pbest and gbest.
- and are random numbers in the range [0, 1], used to introduce randomness.
- is the best position found so far for particle , and is the best position found by all particles in the swarm.
- is the position of particle i at the current time step, and is the position of particle i at the next time step.
- is the maximum value of the inertia weight, and is the minimum value of the inertia weight.
- is the current iteration number, and is the maximum number of iterations.
2.7. Multi-Objective Particle Swarm Optimization (MOPSO) Algorithm
- External Archive
- 2.
- Selection of Guide Particles
- 3.
- Fitness Assessment and Dominance Judgment
- 4.
- Archive Update Strategy
- 5.
- Distance Measure and Goodness of Fit Ranking
2.8. Pareto Efficiency
2.9. Particle Swarm Algorithm Combined with Pareto Efficiency Process
Define fitness_function(x):
Define objective_constraints(f1, f2, f3, …, fn):
Initialize_particles(num_particles, dim):
Update_particles(particles, velocities, p_best, g_best, w, c1, c2):
Pareto_sort(particles, fitness_values):
MOPSO_algorithm(num_particles, dim, num_iterations, w, c1, c2):
|
- Define fitness_function(x).
- Define objective_constraints(f1, f2, f3, …, fn).
- Initialize particles using Initialize_particles(num_particles, dim).
- For each iteration from 1 to num_iterations:
- Calculate fitness values for each particle using fitness_function.
- Update personal best positions (p_best) and fitness values (p_best_fitness), If current fitness values are better than personal best values.
- Use Pareto_sort to get Pareto front indices.
- Update global best position (g_best) to the best solution in Pareto front.
- Output the best particle position (g_best) and the best particle fitness values (g_best_fitness).
2.10. Experimental Procedure
- Parameter Setting: In the initial basic parameter settings, we first need to set the filter type to be optimized and the filter type, such as low-pass, high-pass, band-pass, and other filter types. The initial filter parameter settings also include filter sampling. Frequency, passband frequency range, stopband frequency range, and frequency response; MOPSO initial parameter settings include objective function, number of iterations, number of particles, weights, and other related parameters; this study uses Butterworth filter, and two types of Chebyshev Type I filters as prototypes, and the design parameters of these two types of low-pass filters, high-pass filters, and band-pass filters are optimized.
- MOPSO Optimization: The MOPSO optimization algorithm is used to form a Pareto boundary for each set of iterated solutions and find a set of relatively optimal trade-off solutions through Pareto efficiency. All relevant constraints are defined according to the type of filter to be designed. Each candidate solution should be a vector containing the values of all objective functions. An external archive is used to store the Pareto optimal solution set (Pareto front). The solution of the current particle is compared with the solution in the Pareto Archive, and the archive is updated to keep it containing only non-dominated solutions (Pareto optimal solutions) until the stopping condition is met (for example, the maximum number of iterations is reached, or the solution no longer changes significantly in the external archive). Return to the parameter-setting part if the ideal filtering effect is not achieved and readjust parameters.
- Design Filter: The obtained optimal design parameter solution is used to design a filter and draw the filter amplitude, phase, pole, and zero-point diagrams; then, a section of white noise is randomly generated and filtered through the designed filter. Moreover, the signal-to-noise ratio (SNR) before and after filtering is calculated to confirm whether the designed filter achieves the filtering effect.
3. Experimental Results
Results
4. Discussion
- Optimization of frequency response
- 2.
- Optimization of selectivity and bandwidth
Comparison with Related Literature
5. Conclusions
Author Contributions
Funding
Data Availability Statement
Conflicts of Interest
References
- Chauhan, R.S.; Mehra, R.; Shallu. ASIC design of IIR butterworth digital filter for electrocardiogram. In Proceedings of the 2017 8th International Conference on Computing, Communication and Networking Technologies (ICCCNT), Delhi, India, 3–5 July 2017; pp. 1–6. [Google Scholar]
- Podder, P.; Hasan, M.M.; Islam, M.R.; Sayeed, M. Design and implementation of Butterworth, Chebyshev-I and elliptic filter for speech signal analysis. arXiv 2020, arXiv:2002.03130. [Google Scholar] [CrossRef]
- Ardakani, A.; Leduc-Primeau, F.; Gross, W.J. Hardware implementation of FIR/IIR digital filters using integral stochastic computation. In Proceedings of the 2016 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Shanghai, China, 20–25 March 2016; pp. 6540–6544. [Google Scholar]
- Getu, B.N. Digital IIR Filter Design using Bilinear Transformation in MATLAB. In Proceedings of the 2020 International Conference on Communications, Computing, Cybersecurity, and Informatics (CCCI), Sharjah, United Arab Emirates, 3–5 November 2020; pp. 1–6. [Google Scholar]
- Agrawal, N.; Kumar, A.; Bajaj, V.; Singh, G. Design of digital IIR filter: A research survey. Appl. Acoust. 2021, 172, 107669. [Google Scholar] [CrossRef]
- Loubna, K.; Bachir, B.; Izeddine, Z. Optimal digital IIR filter design using ant colony optimization. In Proceedings of the 2018 4th International Conference on Optimization and Applications (ICOA), Mohammedia, Morocco, 26–27 April 2018; pp. 1–5. [Google Scholar]
- Agrawal, N.; Kumar, A.; Bajaj, V. A new method for designing of stable digital IIR filter using hybrid method. Circuits Syst. Signal Process. 2019, 38, 2187–2226. [Google Scholar] [CrossRef]
- Matei, R. Analytic design of directional and square-shaped 2D IIR filters based on digital prototypes. Multidimens. Syst. Signal Process. 2019, 30, 2021–2043. [Google Scholar] [CrossRef]
- Agrawal, N.; Kumar, A.; Bajaj, V. Optimized design of digital IIR filter using artificial bee colony algorithm. In Proceedings of the 2015 International Conference on Signal Processing, Computing and Control (ISPCC), Waknaghat, India, 24–26 September 2015; pp. 316–321. [Google Scholar]
- Sharaf, A.M.; El-Gammal, A.A.A. A novel discrete multi-objective Particle Swarm Optimization (MOPSO) of optimal shunt power filter. In Proceedings of the 2009 IEEE/PES Power Systems Conference and Exposition, Seattle, WA, USA, 15–18 March 2009; pp. 1–7. [Google Scholar]
- Aimi, H.; Suyama, K. Design of IIR filters by determining particle reallocation space using multi-swarm PSO. In Proceedings of the 2015 15th International Symposium on Communications and Information Technologies (ISCIT), Nara, Japan, 7–9 October 2015; pp. 29–32. [Google Scholar]
- Yamamoto, K.; Suyama, K. Active enumeration of local minima for IIR filter design using PSO. In Proceedings of the 2017 Asia-Pacific Signal and Information Processing Association Annual Summit and Conference (APSIPA ASC), Kuala Lumpur, Malaysia, 12–15 December 2017; pp. 910–917. [Google Scholar]
- Dhanarasi, G.; Kumar, P.S.; Krishna, B.T. Design of Infinite Impulse Response Filter Using Particle Swarm Optimization and its Invariants. In Proceedings of the 2022 International Conference on Smart and Sustainable Technologies in Energy and Power Sectors (SSTEPS), Mahendragarh, India, 7–11 November 2022; pp. 306–310. [Google Scholar]
- Takase, Y.; Suyama, K. A Diversification Strategy for IIR Filter Design Using PSO. In Proceedings of the 2018 Asia-Pacific Signal and Information Processing Association Annual Summit and Conference (APSIPA ASC), Honolulu, HI, USA, 12–15 November 2018; pp. 1365–1369. [Google Scholar]
- Antoniou, A. Digital Filters: Analysis, Design and Applications, 2nd ed.; McGraw-Hill Education (ISE Editions): New York, NY, USA, 1993. [Google Scholar]
- Wang, D.; Tan, D.; Liu, L. Particle swarm optimization algorithm: An overview. Soft Comput. 2018, 22, 387–408. [Google Scholar] [CrossRef]
- Federico, M.; Walczak, B. Particle swarm optimization (PSO). A Tutor. Chemom. Intell. Lab. Syst. 2015, 149, 153–165. [Google Scholar]
- Gad, A.G. Particle swarm optimization algorithm and its applications: A systematic review. Arch. Comput. Methods Eng. 2022, 29, 2531–2561. [Google Scholar] [CrossRef]
- Saptarshi, S.; Basak, S.; Peters, R.A. Particle Swarm Optimization: A survey of historical and recent developments with hybridization perspectives. Mach. Learn. Knowl. Extr. 2018, 1, 157–191. [Google Scholar] [CrossRef]
- Cheng, R.; Jin, Y. A social learning particle swarm optimization algorithm for scalable optimization. Inf. Sci. 2015, 291, 43–60. [Google Scholar] [CrossRef]
- Ganjehkaviri, A.; Jaafar, M.M.; Hosseini, S.E.; Barzegaravval, H. Genetic algorithm for optimization of energy systems: Solution uniqueness, accuracy, Pareto convergence and dimension reduction. Energy 2017, 119, 167–177. [Google Scholar] [CrossRef]
- Figueiredo, E.M.; Ludermir, T.B.; Bastos-Filho, C.J. Many objective particle swarm optimization. Inf. Sci. 2016, 374, 115–134. [Google Scholar] [CrossRef]
Parameter | Coefficient |
---|---|
Acceleration factor (c1, c2) | 1.5 |
Weights | 0.7 |
Number of particles | 30 |
Number of iterations | 100 |
Maximum/minimum speed | 5/−5 |
Search maximum/minimum boundary | 50/0 |
random number | [0, 1] |
Objective function | |
Type | Before Filtering (dB) | After Filtering (dB) | |
---|---|---|---|
Butterworth | Low-Pass | −12.77 | −9.92 |
High-Pass | −11.53 | −6.34 | |
Band-Pass | −12.39 | −8.56 | |
Chebyshev Type I | Low-Pass | −13.62 | −11.19 |
High-Pass | −12.63 | −6.63 | |
Band-Pass | −10.72 | −8.30 |
References | N | M | τ | fP | fs | R | P | Imax |
---|---|---|---|---|---|---|---|---|
Haruna Aimi et al. [11] | 8 | 6 | 5 | 0.175 | 0.25 | 0.92 | 90 | 5000 |
Yamamoto et al. [12] | 6 | 4 | 4 | 0.25 | 0.33 | 0.90 | 80 | 10,000 |
Haruna Aimi et al. [11] | [This Work] | |
---|---|---|
MMSE (×10−2) | 2.12 | 1.83 |
ME (×10−2) | 2.72 | 2.34 |
SD (×10−2) | 0.83 | 0.03 |
Yamamoto et al. [12] | [This Work] | |
---|---|---|
MMSE (×10−2) | 3.26 | 2.83 |
ME (×10−2) | 3.58 | 3.34 |
SD (×10−2) | 4.74 | 2.72 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Su, T.-J.; Zhuang, Q.-Y.; Lin, W.-H.; Hung, Y.-C.; Yang, W.-R.; Wang, S.-M. Design of Infinite Impulse Response Filters Based on Multi-Objective Particle Swarm Optimization. Signals 2024, 5, 526-541. https://doi.org/10.3390/signals5030029
Su T-J, Zhuang Q-Y, Lin W-H, Hung Y-C, Yang W-R, Wang S-M. Design of Infinite Impulse Response Filters Based on Multi-Objective Particle Swarm Optimization. Signals. 2024; 5(3):526-541. https://doi.org/10.3390/signals5030029
Chicago/Turabian StyleSu, Te-Jen, Qian-Yi Zhuang, Wei-Hong Lin, Ya-Chung Hung, Wen-Rong Yang, and Shih-Ming Wang. 2024. "Design of Infinite Impulse Response Filters Based on Multi-Objective Particle Swarm Optimization" Signals 5, no. 3: 526-541. https://doi.org/10.3390/signals5030029
APA StyleSu, T. -J., Zhuang, Q. -Y., Lin, W. -H., Hung, Y. -C., Yang, W. -R., & Wang, S. -M. (2024). Design of Infinite Impulse Response Filters Based on Multi-Objective Particle Swarm Optimization. Signals, 5(3), 526-541. https://doi.org/10.3390/signals5030029