Control parameters of PSO, CSA and GSA for optimization.
This chapter presents the implementation of stable, accurate, and wideband second-order microwave integrators (SOMIs). These SOMI designs are obtained by the use of various cascading combinations of transmission line sections and shunt stubs. In order to obtain the optimal values of the characteristic impedances of these line elements, the particle swarm optimization (PSO), cuckoo search algorithm (CSA) and gravitational search algorithm (GSA) are used to approximate the magnitude response of the ideal second-order integrator (SOI). Based on magnitude response, absolute magnitude error, phase response, convergence rate, pole-zero plot, and improvement graph, the performance measure criteria for the proposed SOMIs are performed. The results of the simulation and statistical analysis reveal that GSA exceeds the PSO and CSA in order to approximate the ideal SOI in all state-of-the-art eligible for wide-band microwave integrator. The designed SOMI is compact and suitable for applications covering ultra-wideband (UWB). The designed SOMI structure is also simulated on Advanced Design Software (ADS) in the form of a microstrip line on a dielectric constant 2.2 RT/Duroid substrate with a height of 0.762 mm. In the 3–15 GHz frequency range, the simulated magnitude result agrees well with the ideal one.
- Cuckoo search algorithm
- gravitational search algorithm
- particle swarm optimization
- line elements
- microwave integrator
Integration plays an important role in many scientific and engineering applications. An integrator is an electronic circuit which produces the output that is the integral of the input applied. Electronic analogue integrators are the basis of analog computers and charge amplifiers, which are performed in the continuous time domain. The integrator is widely used in analog computers, analog-to-digital converters and wave-shaping circuits. Initially, an RC integrator is a circuit that approximates the mathematical process of integration. A simple R-C integrator circuit is shown in Figure 1, in which a capacitor (C) in series with a resistor (R) and the source (). The output () of the circuit is taken across the capacitor (C).
Let i is the resulting current. Applying Kirchhoff’s voltage law to the circuit,
Multiplying throughout by C, we get
as , the term may be neglected
Integrating with respect to on both sides of Eq. (3)
From Eq. (2),
Eq. (7) shows that the output of an integrator circuit is the integral of the input signal. These analog integrators are limited for low frequency application. Thus, the researcher moved to design digital integrators. Digital integrator is a system that performs mathematical operations on a sampled discrete time signal to reduce or enhance certain aspects of that signal. It is commonly used for applications such as waveform shaping, coherent detection, edge detection, and accumulator analysis in biomedical engineering and signal processing. It is widely utilized in biomedical engineering and signal processing applications, for example, as waveform shaping, coherent detection, edge detection, and accumulator analysis. It is also used in radar applications such as the allocation of mobile satellites, enterprise networks, commercial television services and digital services . In order to design the wideband digital integrators, various methods were intended. Using the Newton-cotes integration rule and various digital integration techniques, the Recursive wideband digital integrators have been designed . For low-speed applications up to barely a few hundred MHz, the integrators are primarily designed and implemented. Therefore, to cover wideband applications such as radar and wireless communication, the design and implementation of integrators for high-frequency applications is necessary. The microwave integrator is essentially used to measure the time integral of the input signal at microwave frequencies (0.3–300 GHz). Using wideband integrators, the high-frequency active filters can be introduced, and these wideband integrators can also be used for industrial and real-time applications for ultra-wideband (frequency range 3.1–10.6 GHz) applications [3, 4, 5]. In the z-domain, Hsue et al. have introduced a first-order trapezoidal-rule microwave integrator using a chain-scattering transmission matrix with an operating frequency range from 1 to 10 GHz . By cascading equal-length transmission line sections in the z-domain, three microwave integrators and differentiators have been designed and implemented with different time constants [7, 8]. The second-order microwave integrator (SOMI) was designed by Tsai et al.  in the z-domain. A further first-order microwave integrator was designed by Gautam et al. using the ABCD transmission matrix with a bandwidth of 4 to 10 GHz . SOMI was designed by Gupta et al. in the z-domain over the 1.5–5.5 GHz frequency range . Another SOMI was designed by Gautam et al. in the z-domain over the frequency range 3–15 GHz . Nowadays, analyzers are gravitating towards the use of population-based meta-heuristic algorithms to optimize system coefficients in order to design complex or multi - modal systems [13, 14, 15, 16, 17]. This chapter introduces modern and compact SOMI designs that lead to wide bandwidth. These designed SOMIs are accomplished by cascading three transmission line sections and two single section stubs of equal length. By population-based meta heuristic algorithms, the optimum value of characteristic impedances of these line elements are obtained. A global cost-function solution is achieved by minimizing the error gap between the ideal second-order integrator (SOI) and the designed SOMIs. The design-1 SOMI approximates the ideal SOI over the 2.5 to 16 GHz frequency range, and the design-2 SOMI approximates the ideal SOI over the 3 to 15 GHz frequency range. These designed wideband SOMIs would operate with a wider frequency band to be used in mobile communication on a mobile network such as 4 G and 5 G (above 3 GHz) . All the simulated results are obtained by MATLAB and ADS. These simulated outcomes are formulated to be in close agreement with the ideal one. The novelty of these designed SOMIs exists in terms of wide bandwidth and miniaturization of hardware.
2. Problem formulation of SOMIs
To realize second-order microwave integrator in microwave range, consider a two-port network, which is illustrated in Figure 2.
Its scattering matrix is defined as .
where and are incident waves at port 1 and port 2, respectively, and and are reflected waves at port 1 and port 2, respectively . The chain scattering matrix of a two-port network can be established from the scattering matrix (S-matrix). The chain-scattering matrix of two-port network is defined as .
Then the chain-scattering matrix is
The formulation of SOMIs is employed with equal length line elements is cascading. The overall transfer function of a cascaded network can be established by multiplying the chain scattering matrices of the line elements. These line elements can be transmission line sections and stubs. Assume the length of all transmission line sections and stubs is , where represents the wavelength of the lines at the normalizing angular frequency or we can say that the electrical length of each section (stubs and transmission lines) is set to be at the normalizing frequency. The electrical length of line element is .
The frequency response of an ideal second order integrator (SOI) is given by
where ω represents the angular frequency in radians per second. Since the magnitude response of an open-circuited stub decreases with frequency, an open-circuited stub can be chosen to design an integrator. Assume the impedance of an open-circuited stub to be . Its chain scattering is .
where represents the reference characteristic impedance that is 50 Ω, is the characteristic impedance of an open circuited stub and represents the phase constant.
Assume is the angular frequency and is the propagation delay attributable to the length . Therefore, the term can be expressed by , which can be considered as a unit of delay, that is .
Subsequently, the chain scattering matrix of open-circuited stub is
and is given by
where . Likewise, the chain scattering matrix for the transmission line section in domain is given by .
where the reflection coefficient is given by
where represents the characteristic impedance of serial transmission line section. Similarly, the chain scattering matrix of a short-circuited stub is given by 
where the coefficient is given by
where represents the characteristic impedance of the short-circuited stub. Serial line sections and short-circuited stubs shunted with open circuited stub should be employed in the transmission line configuration to design an integrator by cascading. A cascaded connection of two-port network is equivalent to a single two-port network containing a product of matrices. Assume the SOI is composed of open-circuited stubs, short-circuited stubs and transmission line sections. The overall in generalized form is given by
where, denotes the number of open-circuited stubs, denotes the number of short-circuited stubs, denotes the number of transmission line section, is the total number of line elements and , , ..., Coefficients are functions of the reflection coefficient of each section of the transmission line and stub characteristic impedances (in terms of open-circuited stub coefficient k and short-circuited stub coefficient δ). In addition, to design the SOMI, serial transmission lines cascaded with shunt circuited stubs can be employed. Then the frequency-domain response of the transfer function, which has an integrator characteristic, is thus obtained. In order to design a wideband SOMI, the following two methods are used.
Design-I (SOMI using two open stubs and three transmission line sections): In this design, we may select two open-circuited stubs and three transmission line sections as shown in Figure 3, where . Overall of the design-1 SOMI is given by
Design-II (SOMI using one open stub, one short stub and three transmission line section): In this design, a short-circuited stub is cascaded with the transmission line sections and an open-circuited stub as shown in Figure 5, where . Overall of the design-2 SOMI is given by
The term is due to open-circuited stub, the term is due to short-circuited stub, and is the delay factor of the transmission line section. These line elements are transmission lines of equal length with a length of at an operating frequency of 12.5 GHz.
The next task be to achieve the optimum value of characteristic impedances of design-1 and design-2 SOMI line elements. The optimization algorithms are used to obtain these characteristic impedances of the line elements. In order to lower the cost function, the design of SOMI is considered as an approximation problem. The cost function differs between the response in magnitude of an ideal SOI and the SOMI designed. The cost function is formulated in the sense of least squares and can be expressed as
where is the cost function, in which is the frequency response of an ideal SOI.
3. Employed optimization methods
In infinite impulse response (IIR) systems, the error surface is generally non-quadratic and multimodal with respect to the system parameters. Minimization of such error fitness function using derivative-based search algorithm is difficult. This is due to the fact that the derivative-based search algorithm may not converge to the global minima and get stuck in local minima. Moreover, IIR systems are associated with the stability issues as the poles of the systems may lie outside the unit circle. Such techniques are found unfit to solve multi-objective, multi-modal complex problems and a fine tuning of algorithm parameters is required. To conquer these disadvantages, several practitioners rely on meta-heuristic algorithms, which are based on natural evolution. The meta-heuristic algorithms are nature inspired population-based search techniques which have the ability to serve a global optimal solution with high convergence by accumulating random search and selection principle. Therefore, intelligent search paradigms and optimization methods are adopted in this work for an optimal differentiator design (a multi-modal problem) in short computation time and with high accuracy. Three population-based heuristic search algorithms PSO, CSA, and GSA are employed in this section to diminish the cost function in order to find optimum values of characteristic impedances for designed wideband SOMIs of line elements. Table 1 displays the optimum set of algorithm control parameters for designed SOMIs. A brief description of all three algorithm is discussed below.
|Discovering rate of alien egg, Pa||—||0.25||—|
3.1 Particle swarm optimization
Particle swarm optimization was proposed by Kennedy and Eberhart in nineties based on swarm behavior. PSO has simplicity, high solution quality and superior convergence characteristics as compared to other algorithms. It is more efficient, easy to implement and flexible to control between global and local exploration of the search space [13, 14]. In PSO, every particle has a candidate solution and each candidate solution has its position and velocity. The position vector gives the required solution and the velocity vector gives the current position by which it reaches at the new position. After every iteration, the position and the velocity vectors are updated unless final coefficients are obtained. The velocity and position vector are updated according to the following equations:
where, is the velocity of the particle and is the weight factor. and are positive cognitive parameters while and are two random parameters has range from [0, 1]. is the global best position component at iteration and is the particular best value of particle. is the position vector of particle at iteration. Some parameters are selected to obtain the characteristic impedances of line elements, which are provided in Table 1.
3.2 Cuckoo search algorithm
The cuckoo search algorithm is developed by Yang and Deb in 2009 which is inspired by the concept of unique breeding behavior of cuckoo bird in combination with Lèvy (λ) fights. The theoretical concepts of CSA are well developed and tested in [15, 16, 17]. The single parameter setting in CSA is proven to be a crucial superiority factor as compared to other nature-based algorithms. As a result, the optimized results are executed in very less time. CSA is employed in this section for determining the optimum characteristic values of impedances of transmission line elements on account of minimizing the magnitude error of designed SOMI.
The species of cuckoo birds lay eggs in other bird nests, where the host birds either throw off the detected strange eggs or leave their nests and move into a new spot. The algorithm symbolizes each host nest to a capable solution for this design problem and assigns it a fitness value, as defined in Eq. (24). Furthermore, CSA starts to exchange the current fitness value with a better solution iteratively. The concept of Lèvy flights is then introduced in the process for exploration of new solutions, mathematically modelled using the Lèvy distribution, 1 < λ ≤ 3  with an infinite variance and infinite mean.
3.3 Gravitational search algorithm
This algorithm is based on Newton’s theory of gravitational force and was introduced by E Rashedi et al. in 2009 . Candidates in GSA are supposed to be objects with a given mass. In accordance with Newton’s Law, every candidate in a space has an attraction force with every other candidate. This force correlates inversely to the square of the distance between the candidates and is directly proportional to the product of their masses [20, 21]. Each candidate in the GSA has certain parameters: position of candidate’s mass, inertial mass of candidate, active gravitational mass and passive gravitational mass between candidates. The position of the mass of candidates resembles the solution of the problem devised . The gravitational mass and inertial mass are determined using a cost function of designed problem. The mathematical representation of GSA is considering candidates within a system and all candidates are randomly positioned in search space. The gravitational force on candidate from candidate at time in dimension is defined as .
where, represents gravitational constant at time, represents the passive gravitational mass related to candidate , represents the active gravitational mass related to candidate , represents the Euclidian distance between two candidates and and is a small constant. Then is calculated as .
where, represents the initial value of gravitation constant, is descending coefficient, the current iteration is represented by , and represents the maximum number of iterations.
Then the total force that acts on candidate in dimension is calculated as
where is a random number. As reported to the motion’s law, the acceleration of candidate is given as
where represents the mass of the candidate . The inertial mass and the gravitational mass are updated as follows .
where represents the cost function of the candidate .
The best and worst value are given by
Then the position and the velocity of candidates are given by .
where represents a random number. Finally, the position and velocity of the candidate are obtained.
|Pseudo code for GSA for the design of SOMI-GSA|
|Define and fitness function, Eq. (24)|
|Initialize population size of candidates, ng and other control parameters|
|Set upper and lower bounds, maximum iterations, generate population|
|While iteration, l increases, do|
|Compute fitness using Eq. (24)|
|Evaluate gravitational constant, and|
|Compute masses, gravitational forces, and acceleration using Eqs. (27)–(32). for each candidate|
|Updae the velocity and position of each candidate using Eqs. (35), (36)|
|New solutions are updated|
|Record the best solution|
|Characteristic impedances optimized|
4. Simulation results
In this section, simulation results are discussed and analyzed. All the simulation results are carried out in MATLAB environment. The same control parameters for PSO, CSA and GSA have been selected for appropriate comparison of optimization algorithms. The lower and upper limits of the optimized coefficients are set to be 10 and 150 for functional realizability. Absolute magnitude error (AME), phase response, pole-zero plot, convergence rate and improvement rate are taken into account in assessing the performance of the proposed SOMI magnitude response.
4.1 Design-1 [SOMI using two open stub and three serial line sections]
Figure 3 presents the configuration of the design-1 SOMI, and Table 2 displays the characteristic impedances of design-1 SOMI achieved with PSO, CSA and GSA. The algorithm steps have shown in section-3 by which these characteristic impedances are attained. The magnitude response, the AME, the phase response, the pole-zero plot, the fitness rate and the improvement graph are the parameters selected in Figure 4(a)–(f) for the frequency response analysis of the design-1 SOMI. Table 3 records the overall magnitude error and the phase error. Table 4 and Table 5 summarize the statistical analysis and the qualitatively analyzed data of magnitude error of the design-1 SOMI. It is observed that GSA-based results have least values in all aspects with respect to PSO and CSA. The pole-zero plot of the Design-1 SOMI for stability analysis in which all the poles and zeros lie inside the unit circle is shown in Figure 4(d). This plot asserts that the SOMI design-1 is stable. The convergence profile of all three algorithms in which GSA is clearly demonstrated to be faster than PSO and CSA in less execution time is shown in Figure 4(e), which is also summarized in Table 6. The percentage improvement comparison in magnitude error for the design-1 SOMI as CSA to PSO, GSA to PSO, and GSA to CSA and the percentage improvement comparison in phase errors summarized in Table 7 are shown in Figure 4(f). The design-1 SOMI GSA-based magnitude response closely matches that of the ideal one in the 2.5–16 GHz frequency range. On the basis of these observations, compared to the PSO and CSA, GSA has the lowest magnitude error and the highest convergence speed.
|Algorithm||Total magnitude error||Phase error|
|Algorithm||Iteration cycle||Execution time (s)|
|Execution time (s)|
|Algorithm||Magnitude PI (%)||Phase PI (%)|
4.2 Design-2 [SOMI using one open stub, one short stub and three serial lines]
The design-2 SOMI configuration is illustrated in Figure 5. The optimized characteristic impedances of design-2 SOMI obtained by PSO, CSA, and GSA. The algorithm steps have discussed in Section 3 by which these characteristic impedances are attained. Table 8 offers the optimized characteristic impedances of the line elements of Design-2 SOMI. For the frequency response analysis of the design-2 SOMI, the same parameters are selected as the design-1 SOMI as shown in Figure 6(a)–(f). In Table 9, the total magnitude and phase errors are listed. Table 10 and Table 11 summarize the statistical analysis and the qualitatively analyzed data of different characteristics of magnitude error of the design-2 SOMI. It is observed that GSA-based results have the least values in all aspects as compared to the PSO and CSA. The pole-zero plot of the Design-2 SOMI for stability analysis in which all the poles and zeros lie inside the unit circle is shown in Figure 6(d). This plot asserts that SOMI Design-2 is stable. The convergence profile of all three algorithms in which GSA is clearly demonstrated to be faster than PSO and CSA with less execution time is shown in Figure 6(e), which is also summarized in Table 12. The percentage improvement comparison in magnitude error for the design-2 SOMI as CSA to PSO, GSA to PSO, and GSA to CSA and the percentage improvement comparison in process error, summarized in Table 13, are shown in Figure 6(f). The Design-2 SOMI GSA-based magnitude response closely matches that of the ideal one in the 3.0–15 GHz frequency range. Based on the observations, compared to PSO and CSA, GSA has the least magnitude error and highest convergence speed.
|Algorithm||Total magnitude error||Phase error|
|Algorithm||Iteration cycle||Execution time (s)|
|Execution time (s)|
|Algorithm||Magnitude PI (%)||Phase PI (%)|
Figure 7 depicts the comparison of the designed SOMIs with the existing second order microwave integrator. Figure 7(a) and (b), respectively, display the comparative magnitude response and AME of the designed SOMI GSA-based optimization with an existing integrator. The corresponding qualitative and statistical magnitude error analysis of the existing SOMIs and the designed integrator is shown in Table 14. Figure 7 and Table 14 confirm that GSA-based designed SOMIs are more suitable in terms of magnitude, minimal magnitude error, and wide bandwidth, especially for high-frequency ranges.
The results of the proposed-2 SOMI GSA-based leads in performance to the proposed-1 SOMI counterpart. All simulation results of the design-2 SOMI approximates with the ideal one and have linear phase response in its frequency range 3.0–15 GHz. Therefore, the designed-2 SOMI is simulated on ADS by using microstrip lines. To simulate the designed-2 SOMI on ADS, RT/duroid substrate is selected with dielectric constant =2.2, height of substrate , loss tangent , and copper cladding is 35 um. The magnitude and phase response of the designed-2 SOMI GSA-based is illustrated in Figure 8(a) and (b) respectively.
The study focused on the design and analysis of compact, stable and wideband second order microwave integrators. The designs are obtained by cascading the line elements i.e. transmission line sections, open-circuited stub and short-circuited stub. The optimum values of line elements are obtained by applying the PSO, CSA and GSA by which the magnitude response of designed integrators approximate the ideal magnitude response. The results are simulated statistically on MATLAB, which affirms that GSA outperforms the PSO and CSA in all state-of-the-art in terms of magnitude response. Furthermore, the designed-2 SOMI GSA-based is also simulated on ADS using microstrip lines. With the exception of lower frequency range (under 3 GHz), the simulated magnitude response of integrator has consistency over the frequency range from 3 to 15 GHz with ideal one in both MATLAB and ADS environment, which makes it appropriate for ultra-wideband applications.
Oppenheim AV and Shafer RW: Discrete-Time Signal Processing. New Jersey: Prentice-Hall. 1989
Ngo NQ: A new approach for the design of wideband digital integrator and differentiator. IEEE Transactions on Circuits and Systems II: Express Briefs; 2006; 53; 936–940
Skolink MI: Introduction to Radar Systems. New York: McGrawHill. 1980
Pozar DM: Microwave Engineering, 3rd ed., Wiley, Singapore, 2005
Chiang-Hsue: Design and Implementation of filters using transfer function in the Z Domain. IEEE Trans on Microwave theory and Techniques. May 2001; 49; 5
Hsue CW, Tsai LC and Kan ST: Implementation of a trapezoidal rule microwave integrator. Microwave and Optical Technology Letters. 2006; 48; 822–825
Hsue C-W, Tsai L-C, Tsai Y-H.: Time-constant control of microwave integrators using transmission lines. IEEE Transactions on Microwave Theory and Techniques. 2006; 54; 1043–1047
Tsai L-C and Wu Y-T: Time-constant control analysis of microwave differentiators. IET Microwaves, Antennas & Propagation. 2008; 3; 1044–1050
Tsai LC and Fang HS: Design and implementation of second-order microwave integrators. Microwave and Optical Technology Letters. 2011; 53; 1983–1986
Gautam U, Upadhyay DK and Rawat TK: New designs of first order microwave integrator. IEEE International Conference on Signal Processing (SPIN), India; 2016 285–289
Gupta M and Upadhyay DK: New design of second order microwave integrator. IEEE International Conference on Innovative Mechanisms for Industry Applications (ICIMIA), India; 2017; 556–560
Gautam U, Rawat, TK: The optimal design and analysis of wideband second order microwave integrator. International Journal of Microwave and Wireless Technologies. 2019; 11; 3; 227–236
Kennedy J and Eberhart R: Particle Swarm optimization. Proceedings of IEEE International Conference Neural Network; 1995; 4;1942–1948
Shi JY and Eberhart RC: Empirical study of particle swarm optimization. Proceedings of the Congress on Evolutionary Computation (CEC99), Washington, DC, USA; 1999; 3; 1945–1950
Yang XS and Deb S: Cuckoo search via Lvy flights. Proceedings of world congress on nature and biologically inspired computing. USA IEEE Publications. 2009; 210–214
Yang XS and Deb S: Cuckoo search: recent advances and applications. Neural Computing & Applications. 2014; 24; 169–174
Kumar M and Rawat TK: Optimal fractional delay-IIR filter design using cuckoo search algorithm. ISA Transaction. 2015; 59; 39–54
Hong JS and Lancanter MJ: Microstrip Filters for RF/Microwave Applications. New York: Wiley. 2001
Rashedi E, Nezamabadi S and Saryazdi S: GSA: a gravitational search algorithm. Information Sciences. 2009; 179; 2232–2248
Newton Isaac. In Experimental Philosophy Particular Propositions are Inferred from The Phenomena and Afterwards Rendered General by Induction, 3rd ed.: Andrew Motte’s English Translation published. 1729
Rashedi E, Nezamabadi S and Saryazdi S: Filter modelling using gravitational search algorithm. Engineering Applications of Artificial Intelligence’s. 2011; 24; 117–122