- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

International Journal of Antennas and Propagation

Volume 2013 (2013), Article ID 476878, 9 pages

http://dx.doi.org/10.1155/2013/476878

## Multiobjective-Optimized Design of a New UWB Antenna for UWB Applications

^{1}Department of Electromagnetism, University of Granada, Avenida Fuentenueva sn, 18071 Granada, Spain^{2}Department of Electric and Electronic Engineering, University of Stellenbosch, Matieland 7602, South Africa

Received 20 November 2012; Revised 6 March 2013; Accepted 7 March 2013

Academic Editor: Michael Yan Wah Chia

Copyright © 2013 C. Moreno de Jong van Coevorden et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

A multiobjective genetic algorithm has been applied to design a new printed, bow-tie antenna for ultrawideband applications, that is, ground penetrating radar, short range and high data rate communications, and so forth. The ultrawideband performance with respect to antenna impedance and gain is achieved by an optimized resistive loading profile and flare angle. A low-cost prototype is manufactured and numerical simulations are validated with measurements.

#### 1. Introduction

The antenna performance in ground penetrating radar (GPR) systems is measured by its ability to transmit and receive short pulses (on the order of a few nanoseconds), whose duration is a tradeoff between range resolution and penetration depth [1]. The transmission without distortion of these pulses needs ultrawideband (UWB), nondispersive antenna systems. Such ultrawideband antennas must exhibit a linear phase characteristic over the whole operating frequency band (so-called transient antenna [2]), apart from constant magnitude of the input impedance, polarization, and gain.

When an antenna is fed by a transient pulse (impulse antenna), the initial acceleration of the charges produces radiation from the feed point, which is the sole source of radiation until the traveling pulse reaches any antenna discontinuity. There, the current pulse that propagates along the antenna structure, is partially reflected, constituting secondary radiating sources [3]. Therefore, the time-domain antenna response can be divided into two parts: the main pulse and the ringing region. The main pulse results from the direct radiation of the excitation pulse at the feed point, while radiation originating at antenna discontinuities, gives rise to the ringing region. Depending on the relation between width of the pulse and size of the antenna, these regions could overlap.

A design goal for impulse GPR antennas is the removal or minimization of the ringing region. The main method to achieve this is to establish a proper distribution of resistive loads along the antenna to diminish the reflections of the current pulses at antenna discontinuities [4–6]. As a result, the input-impedance bandwidth of the antenna is increased [7]. However, the enhancement of the ultrawideband characteristics of an antenna, by loading its structure, reduces the antenna gain due to ohmic losses. The design of such antennas is a multiobjective engineering problem with opposite goals, where a tradeoff between bandwidth and gain needs to be found.

One kind of antenna widely used in GPR applications is the solid bow-tie antennas [4, 8] being simple to design and having ultrawideband impedance properties. However, in some cases, the use of the wire bow-tie antenna or strip bow-tie antenna may be advantageous, as they are more easily loaded with resistors than their solid counterparts and have adaptive properties [9, 10]. In [11], a microgenetic optimization algorithm (GA) was used to optimize the input impedance bandwidth of a thin-wire bow-tie antenna. The result was an antenna with a very high impedance bandwidth but a low gain. In an effort to increase the directivity of the antenna, a new thin-wire design was proposed in [12], where the front-to-back ratio and the broadband behavior of the input impedance were simultaneously optimized by means of a multiobjective GA. Improved directivity was achieved by bending the two arms of the antenna in a horn-like way, while the broadband impedance characteristics were accomplished by discrete resistive loading of the antenna wires. The result was greater broadside direction gain but at the expense of augmenting the antenna profile.

In this paper, a new optimized printed-strip bow-tie antenna is proposed. With the aim of using a multiobjective GA for the simultaneous optimization of three following antenna parameters: bandwidth, gain bandwidth, and gain in the broadside direction. This approach goes beyond simply seeking an antenna with high impedance bandwidth, as in the case described in [11]. Moreover, the possibility of using different load profiles on the different strips that compose the antenna geometry is allowed, so that it is possible to find better designs than in [11, 12], where the same load profile was used in all the antenna wires. The only condition imposed on the load profile is to have quadrangular symmetry. The optimization was carried out in two steps. First, a flat, thin-wire bow-tie antenna was designed by hybridizing an in-house multiobjective GA code with the method-of-moments code NEC [13] for modeling the thin-wire antenna response. Subsequently, in order to facilitate the construction of a prototype, the wire antenna was converted into a printed strip using the *same surface area* rule of thumb [14]. Finally, the equivalent printed antenna prototype was built and its measured parameter and gain were found to be in good agreement with the numerical results predicted by the commercial code CST Microwave Studio. (This two-step designing procedure is followed because the computation speed of the hybrid GA-NEC code is much faster than the GA-CST.)

The paper is organized as follows. Section 2 describes the parameters that characterize the antenna to be optimized, the specific strategy followed to produce a new UWB optimized design and the numerical results corresponding to that design. In Section 3, the parameter and antenna gain are measured and compared with numerical simulations. Conclusions are discussed in Section 4, where future directions of investigation are also outlined. Finally, in Appendix, a step-by-step description of the process followed to transform a thin-wire, bow-tie antenna to an equivalent printed antenna is given.

#### 2. Multiobjective Antenna Optimization

There are many strategies to handle multiobjective optimization problems: Vector Evaluated Genetic Algorithm [15], linear aggregation of objectives [16], or the most widely used, which is based on Pareto dominance [17]. In this paper, we will use this last strategy due to its high rate of success in electrical and electronic engineering applications [18, 19]. In particular, we chose the nondominated sorting genetic algorithm in the revised version of Deb et al. (NSGAII) because of its easy implementation, remarkable ability to reach the true Pareto front in mathematical test functions, and less computational complexity than other similar algorithms [17].

In this section, the system to be optimized is first described, setting what parameters can be tuned up and what parameters are kept fixed. Later, the basics of the optimization algorithm that is used in this work is outlined, ending with the presentation of the numerical results from the optimization process.

##### 2.1. Description of the Optimized System

The geometry of the thin-wire bow-tie antenna proposed for optimization is shown in Figure 1. It is made up by dipoles ( wires) forming a fanlike structure and fed at the center common point. Equally spaced gaps of width separated by a distance of are inserted in the wires. At these gaps, chip resistors are soldered in order to enhance the ultrawideband characteristics of the antenna to be optimized in the frequency range between 0.5 GHz and 3 GHz, which is a typical band in GPR applications. Four different resistors are located on each arm, but, since the maximum of the radiation pattern is required to be in the broadside direction, the distribution of resistors is forced to have quadrangular symmetry and therefore only different resistors will be considered in total for an antenna design. The number of dipoles were chosen to be because, although the UWB behavior of the antenna increases with the number of wires, no noticeable improvements are found over [20]; therefore, the load profile will consist in only 16 different resistors.

Among the parameters that define the antenna, the ones that will be allowed to vary during the optimization process are the exterior flare angle , keeping the wires equiangular spaced, and the values of the resistors at wire , and position on that wire, , where is the most external wire and the resistor placed in the vicinity of the feed point. The remaining parameters are kept fixed. The wires are chosen to have a radius of while their length is and the resistors, of body length (0805 SMD chip resistors), are equally spaced over the wires separated by a distance .

##### 2.2. Antenna Optimization Using a Multiobjective GA

GAs are optimization algorithms based on the theories of evolution and genetics [21]. GAs are iterative algorithms that consider a population of individuals, each individual representing a potential solution of the problem at hand, which is described by a set of genes. The process starts by evaluating the quality of individuals in the initial population by calculating a representative function, named *fitness* function, which is defined by the designer in terms of the expected performance of the optimized antenna. The best designs are selected and undergo the genetic operations of crossover and mutation, resulting in a new generation of individuals to be evaluated again. Crossover is usually accomplished by randomly selecting two individuals (called parents) and generating another two (called offsprings) by mixing the genes carried by the parents in some specified fashion. Mutation changes a gene with a certain probability within the allowed values in the design. This process is continued until a stop criterion is met.

The specific optimization algorithm that is used in the present work is the NSGAII. This algorithm is used to optimize, simultaneously, the parameter and gain () of the model proposed above, seeking to achieve the widest possible bandwidth in both magnitudes, with the highest value of the latter. As fitness functions to evaluate the performance of each individual, the following were chosen.:the width of the band where is below , calculated using as characteristic impedance, , being the mean value of the real part of in the frequency band of interest. With this, the algorithm searches antennas having input impedances with real parts near to and imaginary parts as close to 0 as possible. :the width of the band where the gain is between its maximum value and below it. :The maximum value of the Gain in the frequency band of operation.

The NSGAII starts with the creation of a population () (the exponent of means the iteration step of the optimization process) of individuals (antenna designs) with its parameters randomly chosen from a uniform distribution. In the present case, each antenna is represented by a real number, coding the flare angle , and an integer matrix of 16 rows of three integers which defines the value of each of the 16 different resistors of the structure as . After a preliminary study, the search space was reduced to the one defined by the discrete values and and the continuous interval , thereby speeding up the convergence of the optimization process. The main cycle of the iterative process (see Algorithm 1) starts with the evaluation of each individual of the population by computing and over the frequency band of interest by means of the method-of-moments-based frequency-domain code NEC. Then, the three fitness functions defined above are calculated for each individual and stored in a matrix . After the entire population is evaluated and the most promising solutions are selected, genetic operators (GO) are applied to the population to generate the new generation of potential solutions or offspring () from the previous generation ().

The specific genetic operators we employed were as follows. First, a binary tournament selection was used to identify the best individuals within the population. This operator randomly picks two individuals from the previous population and chooses the superior solution for a future crossover following the Pareto domination rules. As a mechanism to recombine the features of two individuals previously selected, a hybrid real-discrete crossover operator is chosen, where the angle is recombined via simulated binary crossover [22], and the resistors by means of a two point crossover [21]. The polynomial distribution introduced by Deb [23] is applied to mutate while to mutate the resistors the algorithm shown in Algorithm 2 is used. Among these operators, the specific multiobjective genetic operator is the tournament selection, in which the dominance in the pareto sense is implemented.

The values of parameter in Algorithm 2 define the variation range of each variable , where . On the other hand, after a short parametric study, it was found that the probability rate gives satisfactory results. This probability was chosen so that, after the individual was selected to mutate (according to the muting rate ), the value mutated resistor was in the neighborhood of the previous one.

Finally, to create the next generation (), an elitism operator was applied where the best individuals from were selected. This process was repeated until a proper set of solutions was found. For further details on the NSGAII algorithm, see [17].

##### 2.3. Optimization Results

To ensure the convergence of the Pareto front, five independent runs of the optimization algorithm were executed with a population of and crossover and mutation rates and , respectively. The achieved nondominated set of solutions appears in Figure 2, which displays designs with bandwidths from 1.6 : 1 to 5.7 : 1, gain bandwidths from 1.0 : 1 to 5.4 : 1, and maximum gains from to . From this set of solutions, the designer can choose the best antenna that matches a particular application. In this paper, a compromise design (marked in red in Figure 2) is chosen to construct a prototype and measure its electrical properties.

The parameters of the chosen antenna are while the load profile is shown in Table 1.

#### 3. Experimental Validation of the Antenna Performance

To illustrate the performance of the antenna design proposed in the previous section, two prototypes were fabricated and tested. One monopole over a ground plate fed with a coaxial connector and a dipole antenna fed by a balun. Both were made of metallic strips printed on FR-4 substrates (, ) of thickness ). A printed antenna model was chosen to facilitate its fabrication and testing in our laboratory facilities. Before the prototypes were built, numerical simulations were carried out by CST software, to determine how the antenna performance varied from the thin-wire model simulated by NEC and optimized by GA, to the printed strip model dealt with CST. The numerical results of all these experiments are presented in Appendix.

First, the monopole antenna shown in Figure 3 was built to measure its parameter. The antenna was printed on a FR-4 panel of dimension , situated perpendicularly to an aluminum ground plane and fed with a coaxial cable with its inner conductor attached to the antenna structure and its outer one to the ground plane. In order to know the minimum size of the ground plane, some numerical simulations were carried out in CST. The simulated antenna system was modeled as realistic as possible, including the FR-4 substrate, finite thickness ground plane and the coaxial connector feeding the antenna. These numerical simulations showed that, for calculations, a size of for the ground plane accurately approximates an infinite ground plane. The parameter was measured, in the frequency range from 500 MHz to 3 GHz, by placing the antenna under test inside an anechoic chamber, and using an HP 8510 C network analyzer. Figure 4 shows the measured return loss compared to the CST numerical results (the measured values were normalized to the mean value of the simulated antenna input resistance, in the frequency band considered, which was ). It can be seen that measurements and simulations match closely throughout the frequency band. Moreover, although it is not shown here, because the optimization was performed only up to 3 GHZ, it should be mentioned that the numerical return loss was below −10 db up to 8 GHz rendering a numerical bandwidth in greater than 10 : 1.

Next, we proceeded to measure the antenna gain pattern but, unlike what happened with the parameter, it was found, the CST simulations indicated that the size of the ground plane needed by the gain pattern of the monopole antenna to match that of the dipole antenna was on the order of several meters. Since this was an impractical size for measurements in the anechoic chamber available at the University of Stellenbosch, as a second prototype, the whole dipole antenna was designed and fabricated. It was printed on a FR-4 panel of dimension () and fed with a balun designed for 1.38 GHz (see Figure 5). Measurements were carried out just at the frequency of 1.38 GHz, while numerical simulations were performed over the entire frequency range. Figure 6 plots, with a symbol line, the numerically calculated gain in the broadside direction, versus frequency, showing a 3 dB gain bandwidth of 1 : 4.5. A comparison between the computed and measured gain patterns at 1.38 GHz is included in the inset of Figure 6, where a close matching between both results is noticed. The measured at 1.38 GHz using the balun was −8 dB.

Finally, since a good time-domain performance is a primary requirement when designing an UWB antenna for GPR applications, the time-domain behavior of the strip-dipole antenna has been investigated. The dipole antenna was excited at its center by a transient, voltage pulse with significant spectral components between 0.75 GHz and 3 GHz and its response was numerically calculated using CST. The shape of the voltage input signal, shown in a solid line in Figure 7, closely resembles that usually found in commercial GPR systems. It also has to be noted how the resistance loading has eliminated the multiple resonances associated to this type of antennas when unloaded [20]. These resonances are related to the presence of eigen modes which, in the case of loaded antennas, are shifted to lower frequencies, out of the frequency band of interest.

To quantify the level of late-time ringing in the antenna time-domain response, the fidelity between the input voltage and the output current at the feed point is calculated as the cross-correlation between both magnitudes using the expression [24]: where , , and are the cross-covariance and autocovariances of signals one and two, respectively. The fidelity compares the pulse shapes disregarding theirs amplitudes and time-delay factors. The value of the fidelity ranges from 0 (totally different signals) to 1 (the signals match perfectly). For the proposed antenna design, the fidelity factor computed for the input voltage and the current at the feed was found to be 0.937. The current pulse at the feed point is plot versus time in Figure 7, and it can be seen that its shape is quite similar to that of the voltage input pulse with a very low late-time ringing.

It is also of great interest to investigate the ability of an UWB antenna to preserve the waveform of the radiated fields in different directions. In the time-domain, and for transmission mode, the quality of the radiated signal waveform in relation to the input voltage applied to the antenna terminals is measured by the fidelity between the time integral of the transmitted field and the input voltage [24]. The time-domain electric fields radiated in different directions in the and planes were calculated using CST. In all cases the fields were computed at a distance of 0.5 m from the antenna feed point. They are plotted in Figures 8 and 9. The simulated fidelity factors between the input voltage and the time-integrated radiated fields are given in Table 2 for different observation angles. It can be seen that the transmitted electric field in the -plane presents high fidelity and its amplitude is almost constant for the different observation angles. Regarding the behavior of the electric field in the -plane, Figure 9 shows how, on the one hand, the amplitude of the field decays as the observation angle is increased and, on the other hand, the field has high fidelity for observation angles but it degrades dramatically for greater observation angles. These results are in accordance with the typical shape of dipole-antenna radiation patterns in the frequency domain, where the amplitude is almost constant in the -plane, while in the -plane varies from maximum values at broadside to minimum in the endfire direction [25].

#### 4. Conclusion

A multiobjective GA has been employed to the design of novel UWB, thin-wire bow-tie antenna for GPR applications. The maximum antenna gain and the impedance and gain bandwidths are simultaneously optimized. The advantage of using a multiobjective GA is that a vast set of design solutions is found in an optimization run instead of only one, allowing the engineer to pick the one most suitable to a particular application. An specific design has been chosen on the Pareto-optimal front and, following the *same surface rule*, an equivalent printed-strip model of the selected thin-wire antenna has been built and measured. The selected UWB antenna has been characterized in the frequency and time domains and good agreement was observed between numerical predictions and experimental results.

#### Appendix

#### Strip Model of the Thin-Wire, Bow-Tie Antenna

A printed antenna based on a thin-wire design requires the width of the printed strips to be chosen so that both models behave in the same manner. This is often done by using the *same surface area* rule [14]. In this appendix, we perform numerical simulations with CST to check whether this rule holds for this particular case.

In the thin-wire model, the area of each wire is , whereas in the printed antenna-model the area of each strip is , where is the strip width (the strip thickness is neglected). For the antenna design of radius of proposed in this paper (see Figure 10), the *same surface area* rule yields -wide strips for the printed model. Figure 10 shows the good agreement between the input impedance of the -wide, strip, dipole antenna and that of its -radius, thin-wire counterpart, in the considered frequency range. Both antennas were located in free space.

#### Acknowledgments

The work described in this paper and the research leading to these results has received funding from the Spain-South African joint research project HS2006-0012, the European Community's Seventh Framework Programme FP7/2007–2013 under Grant Agreement no. 205294 (HIRF SE project), from the Spanish National Projects TEC2010-20841-C04-04, CSD2008-00068, PHB2009-0067-PC, from the Junta de Andalucia Project P09-TIC-5327, and from the Granada Excellence Network of Innovation Laboratories (GENIL) project.

#### References

- D. J. Daniels,
*Surface Penetrating Radar*, IEE Radar, Navigation and Avionics Series, Institute of Electrical & Electronics Enginee, 1996. - P. R. Foster, “Performance of ultrawideband antennas,” in
*Ultrawideband Radar*, vol. 1631 of*Proceedings of SPIE*, pp. 134–145, January 1992. View at Scopus - R. Gómez-Martín, A. R. Bretones, and S. G. García, “Some thoughts about transient radiation by straight thin wires,”
*IEEE Antennas and Propagation Magazine*, vol. 41, no. 3, pp. 24–33, 1999. View at Publisher · View at Google Scholar · View at Scopus - K. L. Shlager, G. S. Smith, and J. G. Maloney, “Optimization of bow-tie antennas for pulse radiation,”
*IEEE Transactions on Antennas and Propagation*, vol. 42, no. 7, pp. 975–982, 1994. View at Publisher · View at Google Scholar · View at Scopus - T. Wu and R. King, “The cylindrical antenna with nonreflecting resistive loading,”
*IEEE Transactions on Antennas and Propagation*, vol. 13, no. 3, pp. 369–373, 1988. View at Google Scholar - M. Fernández-Pantoja, A. Monorchio, A. R. Bretones, and R. Gómez-Martín, “Direct GA-based optimisation of resistively loaded wire antennas in the time domain,”
*Electronics Letters*, vol. 36, no. 24, pp. 1988–1990, 2000. View at Publisher · View at Google Scholar · View at Scopus - S. Lee and K. Mei, “Analysis of zigzag antennas,”
*IEEE Transactions on Antennas and Propagation*, vol. 18, no. 6, pp. 760–764, 1970. View at Google Scholar - Y. Nishioka, O. Maeshima, T. Uno, and S. Adachi, “FDTD analysis of resistor-loaded bow-tie antennas covered with ferrite-coated conducting cavity for subsurface radar,”
*IEEE Transactions on Antennas and Propagation*, vol. 47, no. 6, pp. 970–977, 1999. View at Publisher · View at Google Scholar · View at Scopus - A. A. Lestari, A. G. Yarovoy, and L. P. Ligthart, “Numerical and experimental analysis of circular-end wire bow-tie antennas over a lossy ground,”
*IEEE Transactions on Antennas and Propagation*, vol. 52, no. 1, pp. 26–35, 2004. View at Publisher · View at Google Scholar · View at Scopus - A. A. Lestari, A. G. Yarovoy, and L. P. Ligthart, “Adaptive wire bow-tie antenna for GPR applications,”
*IEEE Transactions on Antennas and Propagation*, vol. 53, no. 5, pp. 1745–1754, 2005. View at Publisher · View at Google Scholar · View at Scopus - C. M. de Jong van Coevorden, A. R. Bretones, M. Fernández-Pantoja, F. J. García Ruiz, S. G. García, and R. Gómez-Martín, “GA design of a thin-wire bow-tie antenna for GPR applications,”
*IEEE Transactions on Geoscience and Remote Sensing*, vol. 44, no. 4, pp. 1004–1009, 2006. View at Publisher · View at Google Scholar · View at Scopus - C. Moreno de Jong van Coevorden, A. R. Bretones, M. Fernández-Pantoja et al., “Thin-wire antenna design for GPR applications using a multi-objective GA,”
*Near Surface Geophysics*, vol. 5, no. 1, pp. 23–28, 2007. View at Google Scholar · View at Scopus - G. Burke, A. Poggio, J. Logan, and J. Rockway, “NEC—numerical electromagnetics code for antennas and scattering,” in
*Antennas and Propagation Society International Symposium*, vol. 17, pp. 147–150, 1979. View at Publisher · View at Google Scholar - A. C. Ludwig, “Wire grid modeling of surfaces,”
*IEEE Transactions on Antennas and Propagation*, vol. 35, no. 9, pp. 1045–1048, 1987. View at Google Scholar · View at Scopus - J. D. Schaffer, “Multiple objective optimization with vector evaluated genetic algorithms,” in
*Proceedings of the 1st International Conference on Genetic Algorithms*, pp. 93–100, 1985. - P. Hajela and C. Y. Lin, “Genetic search strategies in multicriterion optimal design,”
*Structural Optimization*, vol. 4, no. 2, pp. 99–107, 1992. View at Publisher · View at Google Scholar · View at Scopus - K. Deb, A. Pratap, S. Agarwal, and T. Meyarivan, “A fast and elitist multiobjective genetic algorithm: NSGA-II,”
*IEEE Transactions on Evolutionary Computation*, vol. 6, no. 2, pp. 182–197, 2002. View at Publisher · View at Google Scholar · View at Scopus - S. Cui, A. Mohan, and D. S. Weile, “Pareto optimal design of absorbers using a parallel elitist nondominated sorting genetic algorithm and the finite element-boundary integral method,”
*IEEE Transactions on Antennas and Propagation*, vol. 53, no. 6, pp. 2099–2107, 2005. View at Publisher · View at Google Scholar · View at Scopus - Y. Kuwahara, “Multiobjective optimization design of Yagi-Uda antenna,”
*IEEE Transactions on Antennas and Propagation*, vol. 53, no. 6, pp. 1984–1992, 2005. View at Publisher · View at Google Scholar · View at Scopus - C. M. de Jong van Coevorden,
*Optimized antenna design using genetic algorithms [Ph.D. thesis]*, Electromagnetismo, Universidad de Granada, 2008. - K. A. de Jong,
*Analysis of the behavior a class of genetic adaptative systems [Ph.D. thesis]*, Dept. Comput. Commun. Sci., Univ. Michigan, Ann Arbor, Mitch, USA, 1975. - K. Deb and R. B. Agrawal, “Simulated binary crossover for continuous search space,”
*Complex Systems*, vol. 9, pp. 115–148, 1995. View at Google Scholar - K. Deb, “An efficient constraint handling method for genetic algorithms,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 186, no. 2–4, pp. 311–338, 2000. View at Google Scholar · View at Scopus - O. E. Allen, D. A. Hill, and A. R. Ondrejka, “Time-domain antenna characterizations,”
*IEEE Transactions on Electromagnetic Compatibility*, vol. 35, no. 3, pp. 339–345, 1993. View at Publisher · View at Google Scholar · View at Scopus - C. A. Balanis,
*Antenna Theory, Analisys and Design*, John Wiley & Sons, New York, NY, USA, 1982.