Open access peer-reviewed chapter

Analysis of Wideband Second-Order Microwave Integrators

By Usha Gautam and Tarun Kumar Rawat

Submitted: July 13th 2020Reviewed: October 29th 2020Published: January 11th 2021

DOI: 10.5772/intechopen.94843

Downloaded: 118


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

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 (Vin). The output (Vout) of the circuit is taken across the capacitor (C).

Figure 1.

R-C Integrator circuit.

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


Multiplying throughout by C, we get


as RCt, the term 0ti.dtmay be neglected


Integrating with respect to ton 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 [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.

Figure 2.

Two-port network.

Its scattering matrix is defined as [4].


where a1and a2are incident waves at port 1 and port 2, respectively, and b1and b2are reflected waves at port 1 and port 2, respectively [4]. 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 [5].


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 l=λο/4, 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 90οat the normalizing frequency. The electrical length of line element is θ=βl.

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 Zoc. Its chain scattering is [5].


where Z0represents the reference characteristic impedance that is 50 Ω, Zocis 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 l. Therefore, the term jtanβl=jtanωτcan be expressed by D1=ejωτ, which can be considered as a unit of delay, that is [5].


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


where k=Z02Zoc. Now we substitute z=D2in above Eq. (14) for designing purpose, which resembles a scaling by two on frequency axis. Then the chain scattering matrix of open-circuited stub is [6].


and S21is given by


where z=ejβl. Likewise, the chain scattering matrix for the transmission line section in Zdomain is given by [6].


where the reflection coefficient Γis given by


where ZTLrepresents the characteristic impedance of serial transmission line section. Similarly, the chain scattering matrix of a short-circuited stub is given by [8]


where the coefficient δis given by


where ZSCrepresents 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 Popen-circuited stubs, Qshort-circuited stubs and Rtransmission line sections. The overall S21in generalized form is given by


where, Pdenotes the number of open-circuited stubs, Qdenotes the number of short-circuited stubs, Rdenotes the number of transmission line section, Nis the total number of line elements N=P+Q+Rand αο, α1, α2,..., αNCoefficients 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 P=2,R=3. Overall S21of the design-1 SOMI is given by

Figure 3.

Configuration of the design-1 SOMI.


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 P=1,Q=1,R=3. Overall S21of the design-2 SOMI is given by

Figure 4.

Profile of design-1 SOMI (a) Magnitude response, (b) Magnitude error response, (c) Phase response, (d) Pole-zero plot, (e) Convergence profile, and (f) Improvement bae-graph.


The term 1+z1is due to open-circuited stub, the term1z1is due to short-circuited stub, and z3/2is the delay factor of the transmission line section. These line elements are transmission lines of equal length with a length of l=λ0/4at 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 Eω=HωS21is the cost function, in which Hω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.

Population sizenp=25nc=25ng=25
Maximum Iteration400400400
Discovering rate of alien egg, Pa0.25

Table 1.

Control parameters of PSO, CSA and GSA for optimization.

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, vikis the velocity of the ithparticle and ωis the weight factor. C1and C2are positive cognitive parameters whilerand1and rand2are two random parameters has range from [0, 1]. gbestkis the global best position component atkthiteration and pbestikis the particular best value of ithparticle. pbestikis the position vector of ithparticle at k thiteration. 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 [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 consideringNcandidates within a system and all candidates are randomly positioned in search space. The gravitational force on candidate ifrom candidate jat time tin dimension dis defined as [19].


where, Gtrepresents gravitational constant at timet, Mpitrepresents the passive gravitational mass related to candidate i, Majtrepresents the active gravitational mass related to candidate j, Rijtrepresents the Euclidian distance between two candidatesiand jand εis a small constant. Then Gtis calculated as [19].


where, G0represents the initial value of gravitation constant, αis descending coefficient, the current iteration is represented by itr, and maxitrrepresents the maximum number of iterations.

Then the total force that acts on candidate iin dimension dis calculated as


where randjis a random number. As reported to the motion’s law, the acceleration of candidate iis given as


where Maiirepresents the mass of the candidate i. The inertial mass and the gravitational mass are updated as follows [19].


where CFrepresents the cost function of the candidate i.


The best and worst value are given by


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


where randirepresents a random number. Finally, the position and velocity of the candidate are obtained.

Algorithm 1
Pseudo code for GSA for the design of SOMI-GSA
Define Hω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, N<400do
Compute fitness using Eq. (24)
Evaluate gravitational constant, G0and gbest
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 El+1<El
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.

Characteristic impedancePSOCSAGSA

Table 2.

Optimized characteristic impedances of the design-1 SOMI.

AlgorithmTotal magnitude errorPhase error

Table 3.

Comparison summary of total magnitude and phase error for design-1 SOMI.

AlgorithmMeanVarianceStandard deviation

Table 4.

Statistical data of magnitude error for design-1 SOMI.


Table 5.

Qualitatively data of magnitude error for design-1 SOMI.

AlgorithmIteration cycleExecution time (s)
per iteration
Execution time (s)
per cycle

Table 6.

Convergence profile for design-1 SOMI.

AlgorithmMagnitude PI (%)Phase PI (%)

Table 7.

Percentage improvement in magnitude and phase error for design-1 SOMI.

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.

Figure 5.

Configuration of the design-2 SOMI.

Characteristic impedancePSOCSAGSA

Table 8.

Optimized characteristic impedances of the design-2 SOMI.

Figure 6.

Profile of the design-2 SOMI (a) Magnitude response, (b) Magnitude error response, (c) Phase response, (d) Pole-zero plot, (e) Convergence profile, and (f) Improvement bar graph.

AlgorithmTotal magnitude errorPhase error

Table 9.

Comparison summary of total magnitude and phase error for design-2 SOMI.

AlgorithmMeanVarianceStandard deviation

Table 10.

Statistical data of magnitude error for design-2 SOMI.


Table 11.

Qualitatively data of magnitude error for design-2 SOMI.

AlgorithmIteration cycleExecution time (s)
per iteration
Execution time (s)
per cycle

Table 12.

Convergence profile for design-2 SOMI.

AlgorithmMagnitude PI (%)Phase PI (%)

Table 13.

Percentage improvement in magnitude and phase error for design-2 SOMI.

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.

Figure 7.

Comparison of the designed SOMI-GSA with existing integrator (a) Magnitude reponse and (b) Error response.

ReferenceTotal errorMax errorMin errorMean errorVarianceStandard deviation
2nd1–10 GHzGupta et al. [11]10.461.111.756×1040.1150.03080.1755
1–15 GHz[11]17.981.111.756×1040.1370.02780.1667
2nd2.5–16 GHzDesign-1 SOMI4.9950.2384.211×1040.0350.00240.0494
2nd3–15 GHzDesign-2 SOMI3.6570.1423.67×1040.0300.00120.0343

Table 14.

Comparison of magnitude error of the designed SOMIs with existing design method.

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 ϵr=2.2, height of substrate H=30mil0.762mm, loss tangent tanδ=0.001, 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.

Figure 8.

Using ADS software (a) Magnitude response and (b) Phase response.

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.

© 2021 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Usha Gautam and Tarun Kumar Rawat (January 11th 2021). Analysis of Wideband Second-Order Microwave Integrators, Innovations in Ultra-Wideband Technologies, Albert Sabban, IntechOpen, DOI: 10.5772/intechopen.94843. Available from:

chapter statistics

118total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Introductory Chapter: Ultra-Wideband Technologies

By Albert Sabban

Related Book

Advanced Radio Frequency Antennas for Modern Communication and Medical Systems

Edited by Albert Sabban

First chapter

Introductory Chapter: Novel Radio Frequency Antennas

By Albert Sabban

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us