Open access peer-reviewed chapter

Interference Pattern Representation on the Complex s-Plane

Written By

José Trinidad Guillen Bonilla, Alex Guillen Bonilla, Mario Alberto García Ramírez, Gustavo Adolfo Vega Gómez, Héctor Guillen Bonilla, María Susana Ruiz Palacio, Martín Javier Martínez Silva and Verónica María Bettancourt Rodriguez

Submitted: 17 August 2018 Reviewed: 03 September 2019 Published: 09 October 2019

DOI: 10.5772/intechopen.89491

From the Edited Volume

Advances in Complex Analysis and Applications

Edited by Francisco Bulnes and Olga Hachay

Chapter metrics overview

855 Chapter Downloads

View Full Metrics

Abstract

In this work, the normalized interference pattern produced by a coherence interferometer system was represented as a complex function. The Laplace transform was applied for the transformation. Poles and zeros were determined from this complex function, and then, its pole-zero map and its Bode diagram were proposed. Both graphical representations were implemented numerically. From our numerical results, pole location and zero location depend on the optical path difference (OPD), while the Bode diagram gives us information about the OPD parameter. Based on the results obtained from the graphical representations, the coherence interferometer systems, the low-coherence interferometer systems, the interferometric sensing systems, and the fiber optic sensors can be analyze on the complex s-plane.

Keywords

  • coherence interferometer system
  • Laplace transform
  • complex function
  • pole-zero map
  • Bode diagram
  • graphical representations

1. Introduction

Many coherence interferometers systems find practical applications for the physical parameter measurement, such as are temperature, strain, humidity, pressure, level, current, voltage, and vibration [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. Physical implementation and signal demodulation are very important for the good measurement. Many implementations are based on the Bragg gratings, fiber optics, vacuum, mirrors, crystals, polarizer, and their combinations [11, 12, 13, 14, 15]; whereas in the signal demodulation, has been applied commonly the Fourier transform [16, 17, 18, 19, 20]. This transform permits us to know all frequency components of any interference pattern, doing possible the signal demodulation for the interferometer systems.

The Laplace transform has many practical applications in topics such as control systems, electronic circuit analysis, mechanic systems, electric circuit system, pure mathematics, and communications. The linear transformation permits us to transform any time function into a complex function whose variable is s = + σ , where i is the complex operator, ω is the angular frequency, and σ is a real value. The complex function can represent in the complex s-plane, where their axes represent the real and imaginary parts of the complex variable s. This complex plane does feasible the study of dynamic systems, and some applications are the tuning closed-loop, stability, mathematical methods, fault detection, optimization, and filter design [21, 22, 23]. In addition, the s-plane permits graphical methods such as pole-zero map, Bode diagrams, root locus, polar plots, gain margin and phase margin, Nichols charts, and N circles [24].

In dynamic system analysis, pole-zero map and Bode diagrams are two graphical methods which have many practical applications. Both methods require a complex function, where the frequency response plays a very important role. In the pole-zero map, poles and zeros have been calculated from the complex function, and then, their locations are represented on the complex s-plane. It is usual to mark a zero location by a circle and a pole location a cross × [24]. In the Bode diagram, the magnitude and phase are calculated from the complex function, and then, both parameters are graphed. The graphic is logarithmic, and it shows the frequency response of our system under study.

Under our knowledge, the coherence interferometer system was not studied on the s-plane, and as a consequence, its interference pattern was not represented over the pole-zero map or Bode diagrams. In this work, the complex s-plane was used to represent the output signal of an interferometer system. Applying two graphical methods, such as pole-zero plot and Bode plot, the optical signal was represented. Numerically was verified that the pole location and the zero location depend directly on the optical path difference, while a Bode diagram shows the stability/instability of the interferometer.

Advertisement

2. Interference pattern

Figure 1 shows a schematic example of a Michelson interferometer [25]. The interferometer consists of a coherent source, an oscilloscope, a generator function, a beam splitter 50/50 and a PZT optical element. This interferometric system has two ways and its difference produces the optical path difference. The first one will be the reference. Its electrical field is

Figure 1.

Michelson interometer system.

E R t = A R e i ωt + φ R . E1

The second one is the signal measurement and its electrical field is

E D t = A D e i ωt + φ D . E2

A R and A D are amplitudes, ω is the angular frequency, t is the time, and φ R and φ D are given by

φ R = 2 k x R , E3

and

φ D = 2 k x D . E4

x R and x D are the distances traveled by both beams and k = 2 πn λ is the wavenumber: λ is the wavelength and n is the refraction index.

From Figure 1 , when the photodetector detects the total field, its signal is

E T t = A R e i ωt + φ R + A D e i ωt + φ D . E5

Following, the irradiance E T 2 will be

E T 2 = A R e i ωt + φ R + A D e i ωt + φ D A R e i ωt + φ R + A D e i ωt + φ D . E6

Developing, we will obtain

E T 2 = E R 2 + E D 2 + E R E D e i ωt + φ R e i ωt + φ D + e i ωt + φ D e i ωt + φ R E7

or

E T 2 = E R 2 + E D 2 + E R E D e i φ R φ D + e i φ R φ D . E8

Using the identity cos θ = e + e 2 , Eq. (8) takes the form

E T 2 = E R 2 + E D 2 + 2 E R E D cos φ R φ D . E9

In terms of intensity, the interferometer system produces the next interference pattern

I T = I R + I D + 2 I R I D cos φ R φ D . E10

If both beams have the same intensity I R = I D = I o , the total intensity will take the form

I T = 2 I o 1 + cos φ R φ D . E11

As seen in Eq. (11), the phase difference is due to the optical path difference between the two beams. Substituting Eqs. (3) and (4) into Eq. (11), the interference pattern in terms of intensity can be written as

I T t = 2 I o 1 + cos 4 πn λ Δ x t , E12

where Δ x = x R x D is the length difference between the distances x R and x D . Basically, the irradiance is an interference pattern which is formed by two functions: enveloped and modulate. The enveloped function is f env = 2 I o W m 2 and this function contains information from the optical source. The modulate function is given by f mod = 1 + cos 4 πn λ Δ x t and it contains information about the interference pattern. The modulate function consist of a constant (direct component) and a trigonometric function (cosine function) whose frequency depends on the optical path difference.

Advertisement

3. Complex function

Observing Figure 1 and Eq. (11), the phase difference φ R φ D is a time-varying function, and as a consequence, the phase 4 πn λ Δ x t is also a time-varying function. In this case, the instantaneous output voltage (or current) of our photodector is proportional to the normalized optical intensity I T t I o , where I o is the LASER intensity W m 2 [25]. Mathematically, the normalized interference pattern can be written as

I n t = I T t I o = 2 1 + cos ω m t . E13

Here, the angular frequency ω m was proposed from the phase 4 πn λ Δ x t and the interferometer system has not external perturbations.

To determinate the complex function I n s , we calculate the unilateral Laplace transform for our last expression

I n s = 0 I n t e st dt = 2 0 e st dt + 2 0 cos ω m t e st dt . E14

Substituting the trigonometric identity cos ω m t = e i ω m t + e m t 2 into Eq. (14), the complex function can be estimated through

I n s = 2 0 e st dt + 2 0 e i ω m t + e m t 2 e st dt , E15

or

I n s = 2 0 e st dt + 0 e i ω m t e st dt + 0 e m t e st dt . E16

Solving the integrals, the complex function will be

I n s = 2 s e st 0 + e i ω m t + s t i ω m t + s 0 + e m + s t m + s 0 . E17

Evaluating the limits,

I n s = 2 s + 1 s i ω m + 1 s + m . E18

Using the algebraic procedure, we obtain

I n s = 2 s + s + m + s i ω m s 2 i 2 ω m 2 = 2 s + 2 s s 2 + ω m 2 . E19

As seen in Eq. (19), the first term was produced by the direct component and the second term was produced by cosine function. Now, let us represent the complex function as

I n s = 4 s 2 + 2 ω m 2 s s 2 + ω m 2 . E20

Because the Laplace transform was used for the transformation, the normalized interference pattern can be studied in the time domain and on a complex s-plane. It is possible since both Eqs. (13) and (20) contain the same information.

Advertisement

4. Graphical representation

In mathematics and engineering, the s-plane is the complex plane which Laplace transform is graphed. It is a mathematical domain where, instead of view processes in the time domain modeled with time-based functions, they are viewed as equations in the frequency domain. Then, the function I n s can be graphed using the pole-zero map and the Bode diagrams. These graphical representations provide a basis for determining important system response characteristics.

4.1 Pole-zero plot

In general, the poles and zeros of a complex function may be complex, and the system dynamics may be represented graphically by plotting their locations on the complex s-plane, whose axes represent the real and imaginary parts of the complex variable s. Such graphics are known as pole-zero plots. It is usual to mark a zero location by a circle and a pole location a cross × . In this study, it is convenient to factor the polynomials in the numerator and denominator and to write the complex function in terms of those factors

I n s = P N s P D s = 4 s 2 + 2 ω m 2 s s 2 + ω m 2 , E21

where the numerator and denominator polynomials, P N s and P D s , have real coefficient defined by the system’s characteristic. To calculate the zeros, we require

P N s = 0 = 4 s 2 + 2 ω m 2 . E22

Solving last polynomial function, the roots (zeros) are localized at

s 2 = ω m 2 2 s = ω m 2 2 . E23

From our last results, the zeros are imaginary values

s 1 = i ω m 2 s 2 = i ω m 2 . E24

By the similar way,

P D s = 0 = s s 2 + ω m 2 . E25

To calculate the roots,

s 1 = 0 s 2 = i ω m s 3 = i ω m . E26

Using our previous results presented at Eq. (24) and Eq. (26), we represent a pole-zero plot for the interference pattern, see Figure 2 .

Figure 2.

Polo-zero map obtained from the interference pattern.

From Figure 2 , the interference pattern produces two zeros and three poles. Both zeros s 1 and s 2 and two poles s 2 and s 3 are over the imaginary axes and their locations depend on the angular frequency. The pole s 1 was obtained by the direct component; our normalized interference pattern and the location are over the origin.

4.2 Bode diagram

Based on the system theory and system graphic representation, the complex interference pattern can be represented through the Bode diagram. The graphical representation permits us to graph the frequency response of our interferometer system. It combines a Bode magnitude plot, expressing the magnitude (decibels) of the frequency response, and a Bode phase plot, expressing the phase shift.

As was mentioned, the complex interference pattern can be represented through the Bode diagram. To represent it, the term s is substituted by the term : i is the complex number and ω is the angular frequency. Such that, Eq. (20) takes the form

I n = I n ω = 2 2 2 + ω m 2 2 + ω m 2 . E27

The magnitude (in decibels) of the transference function above is given by decibels gain expression:

A vdB = 20 log I n . E28

Substituting Eq. (27) into Eq. (28), the magnitude will be

A vdB = 20 log 2 2 2 + ω m 2 2 + ω m 2 . E29

Applying the logarithm rules, Eq. (29) can express as

A vdB = 20 log 2 + 20 log 2 2 + ω m 2 20 log
20 log 2 + ω m 2 . E30

To determine the phase, Eq. (27) will express as

I n = 4 ω 2 + 2 ω m 2 i ω 3 + ω ω m 2 . E31

Here, i 2 = 1 was used. Last expression can be written as

I n = i 4 ω 2 + 2 ω m 2 ω 3 + ω ω m 2 . E32

From Eq. (32), the phase can also be determined.

Advertisement

5. I n t retrieval

As the Laplace transform is a linear transformation, during the transformation: I n t I n s and I n s I n t , the information is not lost and then the interference pattern can be studied in the time domain and on the complex s-plane. In Section 2, it was explained the transformation I n t I n s , and their poles and zeros were graphed over the pole-zero map. In addition, we developed interference pattern on the frequency plane, being possible to implement the Bode plot. Following, we recover the time function from the complex function I n s I n t . This section is didactic since the objective is to verify that the complex function’s information can also be represented in the time domain.

To recover the time function, we calculate the inverse Laplace transform through

I n t = 1 2 πi lim T γ iT γ + iT I n s e st ds = L 1 I n s . E33

The integral complex is in the s-plane; their limits are γ iT and γ + iT ; the symbol L 1 · = 1 2 πi lim T γ iT γ + iT · e st ds indicates the inverse Laplace transform. Substituting Eq. (21) into Eq. (33), the normalized interference pattern can be obtained by

I n t = 1 2 πi lim T γ iT γ + iT 4 s 2 + 2 ω m 2 s s 2 + ω m 2 e st ds = L 1 4 s 2 + 2 ω m 2 s s 2 + ω m 2 . E34

Applying the partial fraction, Eq. (34) can be expressed as

I n t = 1 2 πi lim T γ iT γ + iT A s + Bs + C s 2 + ω m 2 e st ds = L 1 A s + Bs + C s 2 + ω m 2 . E35

Here, A, B, and C are constants. To calculate the constant, we use next equality

4 s 2 + 2 ω m 2 s s 2 + ω m 2 = A s + Bs + C s 2 + ω m 2 4 s 2 + 2 ω m 2 = A s 2 + A ω m 2 + B s 2 + Cs E36

Using Eq. (36), we obtain the next equation system and their solutions as

A + B s 2 = 4 s 2 Cs = 0 A ω m 2 = 2 ω m 2 A = 2 B = 2 C = 0 . E37

Substituting all constants into Eq. (37), the time function will be

I n t = L 1 2 s + L 1 2 s s 2 + ω m 2 . E38

Applying Table 1 , the solved inverse Laplace transform is

Laplace transform Inverse Laplace transform
Time function Complex function Complex function Time function
f t = ku t F s = k s F s = k s f t = ku t
f t = tu t F s = 1 s 2 F s = k s f t = tu t
f t = t n u t F s = n ! s n + 1 F s = n ! s n + 1 f t = t n u t
f t = cos ωt u t F s = s s 2 + ω 2 F s = s s 2 + ω 2 f t = cos ω u t
f t = sen ωt u t F s = ω s 2 + ω 2 F s = ω s 2 + ω 2 f t = sen ωt u t
f t = cosh ωt u t F s = s s 2 ω 2 F s = s s 2 ω 2 f t = cosh ω u t
f t = senh ωt u t F s = ω s 2 ω 2 F s = ω s 2 ω 2 f t = senh ωt u t

Table 1.

Fundamental Laplace transform [24].

Note: u t is the Heaviside function.

I n t = 2 1 + cos ω m t . E39

Observing both Eq. (13) and Eq. (39), we recover the time function from the complex function. Thus, we confirm that the complex s-plane permits us to study the interferometer system through the complex s-plane, using the pole-zero map and Bode diagrams.

Advertisement

6. Numerical results and discussion

6.1 Results

To verify our proposal, we consider the next interference pattern

I T t = 2 e t 2 1 + cos 10 t . E40

From Eq. (40), the enveloped f env is a Gaussian function f env = 2 I o = 2 e t 2 and the modulate function is f mod = 1 + cos 10 t , where the angular frequency is ω m = 10 radians se . If the interference pattern is normalized as Eq. (12), we obtain

I n t = I T t I o = 2 1 + cos 10 t . E41

Figure 3 shows the interference pattern and the modulate function.

Figure 3.

(a) The simulated interference pattern. (b) The modulate function.

Calculating the Laplace transform,

I n s = P N s P D s = 4 s 2 + 200 s s 2 + 100 . E42

Using Expressions (24) and (42), the zeros are localized at the points

s 1 = i 10 2 s 2 = i 10 2 . E43

Now, using Eqs. (26) and (42), the poles are

s 1 = 0 s 2 = i 10 s 3 = i 10 . E44

Finally, its pole-zero map can observe in Figure 4 .

Figure 4.

Pole-zero plot determined from the complex modulate function (42).

As seen in Figure 4 , the zeros s 1 and s 2 and the poles s 2 and s 3 are over the imaginary axis. Their positions depend on the angular frequency, and therefore, their positions change due to the variations of the optical path difference. The pole s 1 is over the origin (of the complex s-plane), and it was generated by the direct component of our interference pattern.

To generate the Bode diagram, we consider the next complex function

I n s = 4 2 + 200 2 + 100 . E45

Combining Eqs. (28) and (45), the magnitude (in decibels) is

A vdB = 20 log 4 2 + 200 2 + 100 , E46

where s = was used. Applying the logarithm rules, the magnitude can calculate as

A vdB = 20 log 4 2 + 200 20 log 20 log 2 + 100 . E47

Using the Scientific MatLab software, we represent its Bode plot, see Figure 5 .

Figure 5.

Bode diagram obtained from the interference pattern.

Observing Figure 5 , the magnitude has a small variation between the intervals of 100 to 100.7 and from 101 to 102 while the phase is −900. These results confirm the integrative action indicated by Eq. (47). The interferometric system produces two asymptotics for the magnitude. First asymptotic is negative, its location is at the point 100.7 and the phase has transition from −90o to 900. Second asymptotic is positive, its location is 101 and the phase has transition from 900 to −90o. Last interval is between 100.7 and 101. In this case, the magnitude has small variation again but the phase is constant to 90o. These results confirm that the interferometer system will have integrative action and derivative action.

6.2 Discussion

Here, a complex function was obtained from the interference pattern produced by a coherence interferometer. Considering the pole-zero map, poles and zeros depend directly on the optical path difference of an interferometer. The interference pattern generates three poles and two zeros. A pole is over the origin and two poles are over the points ± i ω m . The zeros are over the points ± i ω m 2 . Now, considering the Bode plot, the interferometer can act as an integrator and as a derivator since the phase can take the value of −90o or 90o, see Figure 5 . Both graphical representations permit to know the optical path difference through the angular frequency and its dynamic response.

From our analysis and results, it is possible to infer a few key point of our novel method.

  • The complex s-plane permits us to study the interferometric systems.

  • The interference pattern can represent as a complex function whose poles are three and zeros are two.

  • Pole-zero map gives information about the optical path difference.

  • The pole s 1 is over the origin and it was generated by the direct component of the interference pattern.

  • s 2 and s 3 poles are over the imaginary axis and their position are ± m , where ω m is the angular frequency.

  • s 1 and s 2 zeros are over the imaginary axes and their locations depend on the angular frequency, see Figure 4 .

  • Bode diagram gives us information about the dynamic response of any interference pattern.

  • Based on the Bode diagram (Phase information), the interferometer can act as an integrative action and as a derivative action.

Based on our results, the interference pattern can be studied by both graphical methods. Those graphical representations can be applied to low-coherence interferometric systems, optical fiber sensors, communication systems, and optical source characterization.

Advertisement

7. Conclusion

In this work, applying the Laplace transform and inverse Laplace transform, we confirm that the interference pattern produced by an Interferometer, can study in the time domain and on the complex s-plane. The pole-zero plot and the Bode diagram were obtained from the complex interference pattern. Both graphical representations give us information about the interferometer. The optical Path Difference (OPD) information can measure through the pole-zero map and the behavior of interferometer can understand through the Bode diagram. Therefore, the interferometers can be studied on the complex s-plane, being possible measures physical parameters when those interferometers were disturbed. Also, the signal demodulation can be implemented for the quasi-distributed fiber sensor when the local sensors are interferometers. Some measurable parameters are temperature, string, displacement, voltage and pressure.

If a low-coherence interferometer is studied on the s-plane, then the fringe visibility and the magnitude of coherence grade can be measured.

Advertisement

Acknowledgments

Authors thank the Mexico’s National Council for Science and Technology (CONACyT) and university of Guadalajara for the support.

Advertisement

Conflict of interest

The authors declare no conflict of interest.

References

  1. 1. Zhao N, Lin Q, Jiang Z, Yao K, Tian B, Fang X, et al. High temperature high sensitivity multipoint sensing system based on three cascade Mach-Zehnder interferometers. Sensors. 2018;18(8):2688. DOI: 10.3390/s18082688
  2. 2. Jia X, Liu Z, Deng Z, Deng W, Wzng Z, Zhen Z. Dynamic absolute distance measurement by frequency sweeping interferometry based Doppler beat frequency tracking model. Optics Communication. 2019;430:163-169
  3. 3. Peng J, Lyu D, Huang Q, Qu Y, Wang W, Sun T, et al. Dielectric film based optical fiber sensor using Fabry-Perot resonator structure. Optics Communication. 2019;430:63-69
  4. 4. Vigneswaran D, Ayyanar VN, Sharma M, Sumahí M, Mani Rajan MS, Porsezian K. Salinity sensor using photonic crystral fiber. Sensors and Actuors A: Physical. 2018;269(1):22-28. DOI: 10.1016/j.sna2017.10.052
  5. 5. Kamenev O, Kulchin YN, Petrov YS, Khiznyak RV, Romashko RV. Fiber-optic seismometer on the basis of Mach-Zehnder interfometer. Sensors and Actuors A: Physical. 2016;244:133-137. DOI: 10.1016/j.sna.2016.04.006
  6. 6. Li L, Xia L, Xie Z, Liu D. All-fiber Mach-Zehnder interferometers for sensing applications. Optic Express. 2012;20(10):11109-11120. DOI: 10.1364/OE.20.011109
  7. 7. Yu Q, Zhou X. Pressure sensor based on the fiber-optic extrinsic Fabry-Perot interferometer. Photonic Sensors. 2011;1(1):72-83. DOI: 10.1007/s13320-010-0017-9
  8. 8. Yeo TL, Sun T, Grattan KTV. Fibre-optic sensor technologies for humidity and moisture measurement. Sensors and Actuators: Physical. 2008;144(2):280-295. DOI: 10.1016/j.sna.2008.01.017
  9. 9. Dong C, Li K, Jiang Y, Arola D, Zhang D. Evaluation of thermal expansion coefficient carbon fiber reinforced composites using electronic speckle interferometry. Optics Express. 2018;26(1):531. DOI: 10.1364/OE.26.000531
  10. 10. Wang S, Gao Z, Li G, Feng Z, Feng Q. Continual mechanical vibration trajectory tracking based on electro-optical heterodyne interferometric. Optics Express. 2014;22(7):7799. DOI: 10.1364/OE.22.007799
  11. 11. Wan X, Ge J, Chen Z. Development of stable monolithic wide-field Michelson interferometers. Applied Optics. 2011;50(21):4105-4114. DOI: 10-1364/AO.50.004105
  12. 12. Hassan MA, Martin H, Jiang X. Development of a spatially dispersed short-coherence interferometry sensor using diffraction grating orders: Publisher’s note. Applied Optics. 2018;57(1):5. DOI: 10.1364/AO.57.000005
  13. 13. Guillen Bonilla JT, Guillen Bonilla A, Rodríguez Betancourtt VM, Guillen Bonilla H, Casillas Zamora A. A theoretical study and numerical simulation of a quasi-distributed sensor based on the low-finesse Fabry-Perot interferometer: Frequency-division multiplexing. Sensors. 2017;17(4):859. DOI: 10.3390/s17040859
  14. 14. Liang Y, Zhao M, Wu Z, Morthier G. Investigation of grating-assisted trimodal interferometer biosensors based on a polymer platform. Sensors. 2018;18(5):1502. DOI: 10.3390/s.18052502
  15. 15. Guillen Bonilla JT, Guillen Bonilla H, Casillas Zamora A, Vega Gómez GA, Franco Rodríguez NE, Guillen Bonilla A, et al. Twin-grating fiber optic sensors applied on wavelength-division multiplexing and its numerical simulation. In: Rao SP, editor. Book of Numerical Simulation in Engineering and Science. 1st ed. IntechOpen. DOI: 10.5772/intechopen.75586
  16. 16. de SAuza JC, Oliveira ME, Dos Santos PAM. Brach-cut algorithm for optical phase unwrapping. Optics Letters. 2015;40(15):3456-3459. DOI: 10.1364OL.40.003456
  17. 17. Mizuno T, Kitoh T, Oguma M, Inoue Y, Shibata T, Hiroshi T. Mach-Zehnder interferometer with uniform wavelength period. Optics Letters. 2004;29(5):454-456. DOI: 10.1364/OL.29.000454
  18. 18. Miridonov SV, Shlyagin M, Tentori D. Twin-grating fiber optic sensor demodulation. Optics Communication. 2001;191(3–6):253-362. DOI: 10.1016/S0030-4018(01)01160-9
  19. 19. Kin JA, Kim JW, Kang CS, Jin J, Eom TB. Interferometric profile scanning system for measuring large planar mirror surface based on single-interferogram analysis using Fourier transform. Measurements. 2018;118:113-119. DOI: 10.1016/j.measurement.2018.01.023
  20. 20. Perea J, Libbey B, Nehmetallah G. Multiaxis heterodyne vibrometer for simultaneous observation of 5 degrees of dynamic freedom from a single beam. Optics Letters. 2018;43(13):3120-2123. DOI: 10.1364/OL.43.003120
  21. 21. Davoodi M, Mezkin N, Khorasani K. A single dynamic observer-based module for design of simultaneous fault detection, isolation and tracking control scheme. International Journal of Control. 2018;91(3):508-523. DOI: 10.1080/00207179.2017.1286041
  22. 22. Ordóñez Hurtado RH, Crisostomi E, Shorten RN. An assessment on the use of stationary vehicles to support cooperative positioning systems. International Journal of Control. 2018;91(3):608-621. DOI: 10.1080/00207179.2017.1286537
  23. 23. Xu S, Sun G, Li Z. Finite frequency vibration suppression for space flexible structures in tip position control. International Journal of Control, Automation and Systems. 2018;16(3):1021-1029. DOI: 10.1007/s12555-0160343-9
  24. 24. Wolovich WA. Automatic Control Systems: Basic Analysis and Design. 1st ed. Sounders College Pub.; 1994. 450 p. ISBN: 0030237734, 9780030237737
  25. 25. Andrade Barbosa FA, Nader G, Tokio Higuti R, Kitano C, Nelli Silva E. A simple interferometric method to measure the calibration factor and displacement amplification in piezoelectric flextensional actuators. Revista Controle & Automação. 2010;21(6):577-587

Written By

José Trinidad Guillen Bonilla, Alex Guillen Bonilla, Mario Alberto García Ramírez, Gustavo Adolfo Vega Gómez, Héctor Guillen Bonilla, María Susana Ruiz Palacio, Martín Javier Martínez Silva and Verónica María Bettancourt Rodriguez

Submitted: 17 August 2018 Reviewed: 03 September 2019 Published: 09 October 2019