Control parameters of PSO, CSA and GSA for optimization.

## Abstract

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.

### Keywords

- Cuckoo search algorithm
- gravitational search algorithm
- particle swarm optimization
- line elements
- microwave integrator

## 1. Introduction

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 (

Let i is the resulting current. Applying Kirchhoff’s voltage law to the circuit,

Multiplying throughout by C, we get

as

Integrating with respect to

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 [1]. 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 [2]. 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 [6]. 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. [9] 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 [10]. SOMI was designed by Gupta et al. in the z-domain over the 1.5–5.5 GHz frequency range [11]. Another SOMI was designed by Gautam et al. in the z-domain over the frequency range 3–15 GHz [12]. 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) [18]. 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 [4].

where

For network computation the chain-scattering matrix from S-matrix (on solving Eq. (8) and Eq. (9))

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

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

where

Assume

Subsequently, the chain scattering matrix of open-circuited stub is

where

and

where

where the reflection coefficient

where

where the coefficient

where

where,

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

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

The term

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

## 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.

Parameters | PSO | CSA | GSA |
---|---|---|---|

Population size | |||

Maximum Iteration | 400 | 400 | 400 |

2 | — | — | |

2 | — | — | |

0.01 | — | — | |

1.0 | — | — | |

0.1 | — | — | |

1.0 | — | — | |

Discovering rate of alien egg, Pa | — | 0.25 | — |

— | — | 100 | |

— | — | 20 | |

— | — | 0.001 | |

— | — | 2 | |

— | — | 1 |

### 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,

### 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 [17] 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 [19]. 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 [21]. The gravitational mass and inertial mass are determined using a cost function of designed problem. The mathematical representation of GSA is considering

where,

where,

Then the total force that acts on candidate

where

where

where

The best and worst value are given by

Then the position and the velocity of candidates are given by [19].

where

Algorithm 1 |

Pseudo code for GSA for the design of SOMI-GSA |

Define |

Initialize population size of candidates, ng and other control parameters |

Set upper and lower bounds, maximum iterations, generate population |

While iteration, l increases, |

Compute fitness using Eq. (24) |

Evaluate gravitational constant, |

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) |

For minimzation, |

If |

New solutions are updated |

EndIf |

EndWhile |

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 Figures 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.

Characteristic impedance | PSO | CSA | GSA |
---|---|---|---|

12 | 112.2 | 148.8 | |

49 | 17 | 10 | |

120 | 114 | 149 | |

14 | 11 | 43 | |

115 | 14 | 21 |

Algorithm | Total magnitude error | Phase error |
---|---|---|

PSO | 5.7746 | 1.0403 |

CSA | 5.5353 | 1.0031 |

GSA | 4.9954 | 1.0030 |

Algorithm | Mean | Variance | Standard deviation |
---|---|---|---|

PSO | 0.0410 | 0.0035 | 0.0588 |

CSA | 0.0393 | 0.0036 | 0.0603 |

GSA | 0.0354 | 0.0024 | 0.0494 |

Algorithm | Maximum | Minimum | Average |
---|---|---|---|

PSO | 0.2642 | 0.0410 | |

CSA | 0.3762 | 0.0393 | |

GSA | 0.2384 | 0.0354 |

Algorithm | Iteration cycle | Execution time (s) per iteration | Execution time (s) per cycle |
---|---|---|---|

PSO | 400 | 313.2654 | 0.6132 |

CSA | 400 | 291.3719 | 0.5547 |

GSA | 400 | 269.6370 | 0.4950 |

Algorithm | Magnitude PI (%) | Phase PI (%) |
---|---|---|

PSO | 4.1613 | 3.5758 |

CSA | 13.4935 | 3.5855 |

GSA | 9.7374 | 0.0099 |

### 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 Figures 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.

Characteristic impedance | PSO | CSA | GSA |
---|---|---|---|

116 | 26.87 | 141 | |

30.5 | 142.6 | 33 | |

145 | 84 | 138 | |

9.5 | 19 | 10 | |

110 | 138.09 | 116.89 |

Algorithm | Total magnitude error | Phase error |
---|---|---|

PSO | 3.9217 | 0.9861 |

CSA | 3.6920 | 1.0968 |

GSA | 3.6573 | 1.0595 |

Algorithm | Mean | Variance | Standard deviation |
---|---|---|---|

PSO | 0.0324 | 0.0237 | |

CSA | 0.0303 | 0.0013 | 0.0354 |

GSA | 0.0300 | 0.0012 | 0.0343 |

Algorithm | Maximum | Minimum | Average |
---|---|---|---|

PSO | 0.1048 | 0.0324 | |

CSA | 0.2097 | 0.0303 | |

GSA | 0.1425 | 0.0300 |

Algorithm | Iteration cycle | Execution time (s) per iteration | Execution time (s) per cycle |
---|---|---|---|

PSO | 400 | 301.5411 | 0.5918 |

CSA | 400 | 273.1862 | 0.5356 |

GSA | 400 | 244.6284 | 0.4802 |

Algorithm | Magnitude PI (%) | Phase PI (%) |
---|---|---|

PSO | 5.8571 | −11.22 |

CSA | 6.7419 | −7.443 |

GSA | 0.9398 | 3.400 |

Figure 7 depicts the comparison of the designed SOMIs with the existing second order microwave integrator. Figures 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.

Integrator Order | Frequency range | Reference | Total error | Max error | Min error | Mean error | Variance | Standard deviation |
---|---|---|---|---|---|---|---|---|

2nd | 1–10 GHz | Gupta et al. [11] | 10.46 | 1.11 | 0.115 | 0.0308 | 0.1755 | |

1–15 GHz | [11] | 17.98 | 1.11 | 0.137 | 0.0278 | 0.1667 | ||

2nd | 2.5–16 GHz | Design-1 SOMI | 4.995 | 0.238 | 0.035 | 0.0024 | 0.0494 | |

2nd | 3–15 GHz | Design-2 SOMI | 3.657 | 0.142 | 0.030 | 0.0012 | 0.0343 |

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

## 5. Conclusion

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.