Open access peer-reviewed chapter

Interference Pattern Representation on the Complex s-Plane

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: August 17th 2018Reviewed: September 3rd 2019Published: October 9th 2019

DOI: 10.5772/intechopen.89491

Downloaded: 521


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.


  • 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 iis 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.


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.


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


ARand ADare amplitudes, ωis the angular frequency, tis the time, and φRand φDare given by




xRand xDare the distances traveled by both beams and k=2πnλis the wavenumber: λis the wavelength and nis the refraction index.

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


Following, the irradiance ET2will be


Developing, we will obtain




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


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


If both beams have the same intensity IR=ID=Io, the total intensity will take the form


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


where Δx=xRxDis the length difference between the distances xRand xD. Basically, the irradiance is an interference pattern which is formed by two functions: enveloped and modulate. The enveloped function is fenv=2IoWm2and this function contains information from the optical source. The modulate function is given by fmod=1+cos4πnλΔxtand 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.


3. Complex function

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


Here, the angular frequency ωmwas proposed from the phase 4πnλΔxtand the interferometer system has not external perturbations.

To determinate the complex function Ins, we calculate the unilateral Laplace transform for our last expression


Substituting the trigonometric identity cosωmt=eiωmt+emt2into Eq. (14), the complex function can be estimated through




Solving the integrals, the complex function will be


Evaluating the limits,


Using the algebraic procedure, we obtain


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


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.


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


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


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


From our last results, the zeros are imaginary values


By the similar way,


To calculate the roots,


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 s1ands2and two poles s2ands3are over the imaginary axes and their locations depend on the angular frequency. The pole s1was 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 sis substituted by the term : iis the complex number and ωis the angular frequency. Such that, Eq. (20) takes the form


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


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


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


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


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


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


5. Intretrieval

As the Laplace transform is a linear transformation, during the transformation: IntInsand InsInt, 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 IntIns,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 InsInt.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


The integral complex is in the s-plane; their limits are γiTand γ+iT; the symbol L1·=12πilimTγiTγ+iT·estdsindicates the inverse Laplace transform. Substituting Eq. (21) into Eq. (33), the normalized interference pattern can be obtained by


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


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


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


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


Applying Table 1 , the solved inverse Laplace transform is

Laplace transformInverse Laplace transform
Time functionComplex functionComplex functionTime function

Table 1.

Fundamental Laplace transform [24].

Note: utis the Heaviside function.


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.


6. Numerical results and discussion

6.1 Results

To verify our proposal, we consider the next interference pattern


From Eq. (40), the enveloped fenvis a Gaussian function fenv=2Io=2et2and the modulate function is fmod=1+cos10t,where the angular frequency is ωm=10radiansse. If the interference pattern is normalized as Eq. (12), we obtain


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,


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


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


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 s1ands2and the poles s2ands3are 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 s1is 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


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


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


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ωm2. 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 s1is over the origin and it was generated by the direct component of the interference pattern.

  • s2and s3poles are over the imaginary axis and their position are ±m, where ωmis the angular frequency.

  • s1and s2zeros 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.


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.



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


Conflict of interest

The authors declare no conflict of interest.

© 2019 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

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 (October 9th 2019). Interference Pattern Representation on the Complex s-Plane, Advances in Complex Analysis and Applications, Francisco Bulnes and Olga Hachay, IntechOpen, DOI: 10.5772/intechopen.89491. Available from:

chapter statistics

521total 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

Inverse Scattering Source Problems

By Mozhgan “Nora” Entekhabi

Related Book

First chapter

Introductory Chapter: Frontier Research on Integral Equations and Recent Results

By Francisco Bulnes

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