Scientific Research

An Academic Publisher

**Synthesis of Antennas for Field and Polarization Control** ()

Keywords

Share and Cite:

*Journal of Electromagnetic Analysis and Applications*,

**9**, 97-112. doi: 10.4236/jemaa.2017.97009.

1. Introduction

The increasing usage of electromagnetic polarization in telecommunication systems is well known and documented, ranging from applications such as diversity and frequency reuse to electromagnetic compatibility and future 4 k/8 k resolution Digital Television Broadcasting [1] - [6] . The polarization of a wave is defined by the path described by the electric field vector at a fixed point in space as time varies, and the polarization of an antenna is defined as the polarization of the wave it radiates when operating in transmitting mode [6] . The possibility of generating 3D polarizations may lead to enhancements in the existing telecommunication systems or even the development of new configurations and applications.

This work presents a way of generating 3D polarizations and field distributions, where the electric field intensity vector can be oriented in any direction in space. Therefore, it is possible to generate 3D polarizations and/or fields with different shapes, as time varies, such as cylinders, spheres, and others. In addition, the orthogonality between 3D polarizations is also discussed. Antennas naturally generate a polarization that is orthogonal to the desired one, normally referred to as cross polarization, but the idea here is to generate the orthogonal polariza- tions in a controlled manner, similarly to the use of orthogonal polarizations in frequency reuse systems. Cross polarization is normally undesired, with the exception of the application introduced in [7] .

It was previously demonstrated analytically and numerically by the authors in [8] that 3D polarizations can be generated at an observation point inside a 3D antenna array, where the elements are facing inwards. Volumetric control of the intensity of electric fields and sound has been previously published [9] [10] [11] [12] , but as far as the authors know, generation and application of 3D polarizations and fields have not yet been investigated in the way herein presented.

It was demonstrated in [8] that ideal antennas can generate 3D polarizations, and here the concept is extended for pyramidal horn antennas, as well as spatially distributed electromagnetic fields. A new method introduced by the authors for optimum design of these antennas was also employed in order to ensure greater gain stability with respect to manufacturing variations [13] .

A brief summary of the concept discussed in [8] is presented for completeness, also illustrating different scenarios than the ones covered in [8] , and it is shown that antenna arrays can be used for spatial multiplexing in wireless communi- cation systems by controlling the fields in 3D space, in order to reduce interference among users. This might be critical in the actual trending of the evolution of cellular systems towards nano and phantom cells [14] . Another possible application for the field and polarization control herein discussed is the generation of torus-knotted polarizations and fields, where the latter can be encountered in nature around planets, stars and other celestial objects, as indicated historically in Radioastronomy as well as possibly to be further studied with state-of-the-art Radiotelescopes [15] , and may be of importance for plasma confinement and nuclear magnetic resonance devices [16] [17] .

2. Three-Dimensional Polarization with Ideal Dipoles (Near-Field) and Orthogonality

The antenna array shown in Figure 1 is used to generate 3D polarizations at the observation point O, and was derived from the antenna array configurations introduced by one the authors in [9] [10] . All the computer simulations employed the computer package MathCad [18] . Note that with this configuration, each

Figure 1. Array of 3 ideal dipole antennas in free space. Although not discussed in this work, further consideration can also be given to array configurations of moving elements, including ones in known paths.

antenna produces fields only in the radial direction. This was done to simplify the equations discussed in [8] , even though it is possible to change the configuration for any desired direction. If the amplitude and phase of the antennas are constant, the resulting fields obtained by superposition (it is assumed a linear, homogenous, time invariant medium) normally lead to an elliptical polarization. However, by controlling the variations of the amplitude and/or phase of the elements it is possible to rotate this ellipse with time, and hence yielding a 3D polarization.

The electric field generated by one antenna is frequency dependent, however the figure formed by the fields generated in the array of ideal dipole are frequency independent if the same frequency is used for the all antennas. This occurs because the proportionality between the electric fields generated by individual antennas is kept constant when the frequency is changed. Thus, for illustration purposes, without loss of generality, we use f = 100 MHz, I = 1 mA, h = λ/100 (the length of all 3 antennas placed at a distance equal to λ from point O). It is important to note that for the scenarios considered in this study there is no propagation, so the 3D polarization needs to be visualized as the polarization of the resulting electric field, not the wave. Alternate configurations that allow the 3D polarization to move within a desired path are currently being investigated (by determining different sets of feed currents for each desired location along the trajectory, for example), but the theory does not lose its importance nor generality as herein presented. By varying continuously the phase of the ideal dipole P_{2} in Figure 1, i.e. ψ_{2} in the interval 0 to 2π, with ψ_{1} = ψ_{3} = 0˚ respectively for the other two dipoles, the resulting polarization is shown in Figure 2. In general, 3 components (basis) of electric fields are necessary in order to generate a 3D polarization. In the case of the components being orthogonal, they form an orthogonal base; the reader is referred to [8] for further details.

3. Three-Dimensional Polarization with Pyramidal Horn Antennas (Far-Field)

In order to show that it is also possible to generate 3D polarizations with antennas other than the ideal dipoles, we employed the array of pyramidal horns shown in Figure 3. The geometry of each horn is depicted in Figure 4.

Figure 2. Cylindrical 3D polarization obtained with the configuration of Figure 1 with ψ_{1} = ψ_{3} = 0˚ and ψ_{2} varying continuously in the interval 0 to 2π.

Figure 3. Array of pyramidal horn antennas.

Figure 4. Geometry of a pyramidal horn antenna.

The electric field of each horn can be approximated by (far-field) [19]

(1)

(2)

where E_{0} is the absolute value of the intensity of the electric field at the aperture (assuming the mode TE_{01}). The electric field at the aperture is given by

(3)

where. The functions I_{1} and I_{2} are given by

(4)

(5)

where C(x) and S(x) are the Fresnel functions, and are given in [19] . Although the equations are not exact, the generation of 3D polarization from them can be demonstrated without loss of generality since a more exact approach would lead to different values in the radius of the sphere, which can always be compensated by adjusting properly the amplitude and/or phases of the element excitations (which in this example are the electric fields in the aperture of each antenna). Using electric fields distributions in the aperture of each horn as

(6)

(7)

(8)

where ζ is a variable controlled by the amplitude variation of the electric field in the aperture of each antenna, it can be shown that the total electrical field (far-field analysis) at the observation point O is given by

(9)

which in the time domain is

(10)

where C_{0} and φ_{c} are respectively the absolute value and phase of (valid for O at the origin)

(11)

The method introduced by the authors in [13] was employed to design each pyramidal horn (all equal to each other), and is herein omitted for simplicity. The optimum design ensures greater gain stability with respect to variations of the aperture dimensions [13] . The design gain was 18 dBi, frequency of 6 GHz and the waveguide WR137 (R_{1} = 19.94 cm, R_{2} = 18.2 cm, A = 18.46 cm, B = 14.19 cm). The antennas are away 10^{3}λ from the origin (far-field), and E_{0} = 1 V/m. The spherical polarization shown in Figure 5 is obtained by varying ζ continuously.

4. Spatial Multiplexing by Electric Field Control

In the previous sections the field control was done in only one point in the space; the reader is referred to [8] for a sensitivity analysis of degradations around that point. For spatial multiplexing, the fields should be controlled in many points in the space. This section shows the development of a method for controlling the electric field using three-dimensional array of antennas in multiple points in space, and how it can be applied in spatial multiplexing during the transmission of information in wireless digital communication systems.

The number of users in a given region can be increased with the new method, which is based on space division multiplexing, using the same frequency and time resources, with the least possible interference between users. Thus, the channel capacity can be increased. The interference can be further reduced between nearby users, increasing the signal to noise ratio (S/N) with the control of orthogonal polarizations.

4.1. Method for Electric Field Control

When m antennas are used for the control of n points in space, individual electric fields generated by each antenna are added to obtain the resulting field. The

Figure 5. Spherical polarization with pyramidal horns (far-field).

most direct method for this is to decompose the spherical components a_{r} and a_{θ} directions in cartesian coordinates a_{x}, a_{y} and a_{z}:

(12)

where a_{rx}, a_{ry} and a_{rz} respectively represent the decomposition of the components of the unit vector a_{r} in the directions a_{x}, a_{y} and a_{z}. Also, a_{θx}, a_{θy} and a_{θz} respectively represent the decomposition of the components of the unit vector a_{θ} in the directions a_{x}, a_{y} e a_{z}.

Considering an antenna located at point P_{1}, and an observation point in O_{1}, the equations below are obtained in cartesian coordinates:

(13)

(14)

(15)

where P_{1} is the position vector of antenna 1, S_{1} is the vector direction of antenna 1, and O_{1} is the position vector of the observation point 1. Superscript T in the equations was used for simplification, and indicates the transpose of the row vector. All these vectors are referenced to the origin of the cartesian system. P_{1} and O_{1} start from the origin and go to the specified points, and the distances represent the amplitudes. S_{1} starts from point P_{1} and is defined here as a unit vector, and its amplitude equals 1. The distance vector D_{11} between O_{1} and P_{1} is. The angle formed between the vector D_{11} and S_{1} is given by

(16)

The distance vector D_{11} has the same direction of the vector a_{r}_{11}. In general, establishing the variable u as being equal to x, y or z, a_{ru}_{11} can be determined as:

(17)

With a_{rx}_{11}, a_{ry}_{11}, and a_{rz}_{11} values obtained through (17), a_{θ}_{x}_{11}, a_{θ}_{y}_{11}, and a_{θ}_{z}_{11} can be determined as:

(18)

Using, the contribution of the excitation source can be separated of the contribution responsible for the position and direction of the antenna:

(19)

(20)

where φ_{Eu}_{11} is the argument of E_{u}_{11}. The coefficient C_{u}_{11} has the unit Ω/m. C_{u}_{11} and φ_{Eu}_{11} are calculated as

(21)

(22)

with

(23)

(24)

The phase and amplitude can be grouped in a single complex variable σ_{u}_{11} as shown below

(25)

(26)

Once the values of the length of the antenna and the frequency of the excitation source are determined in the vacuum, the coefficient σ_{u}_{11} will vary only with changes in the position and direction of the antenna, as well as the distance of this to the observation point. With m antennas and n observation points (m = 3n), Equations (25) and (26) can be generalized to

(27)

(28)

where i and l are integers number, with,. Also, u = x, y or z. Since n can take any integer value, the electric field control in space can exist in the infinite points. In full matrix notation, this can be written as:

(29)

Calling s the matrix of coefficients, I the current vector and E the vector fields, the values of the excitation sources of the antennas can be determined by:

(30)

where is the inverse matrix of s. In practice, numerical methods such as Gaussian elimination or LU decomposition are used because they are more efficient and accurate than the numerical computation of the inverse matrix [20] . Thus, the values of the amplitudes and phases of the current phasor are obtained with the desired values of E at each point, as well with the calculation of the complex coefficient matrix s.

4.2. Digital Signal Transmission and Spatial Multiplexing Using the Method for Electric Field Control

In this section the volumetric antennas array and fields control can be employed for the transmission of digital signals as well as to minimize the interference by spatial multiplexing. For the simulation, the antenna array presented in Figure 6 was used.

For the observation point O_{1}, the main antennas are 1, 2 and 3 because they contribute more significantly to the resulting electric field at this point, since the value of the coefficient C_{u} has a higher value. Similarly, for point O_{2}, the main antennas are 4, 5 and 6. In the configuration of Figure 6, in each observation point the antennas are oriented such that the θ angle calculated by (26) is 90˚.

Assuming, at an operating frequency of 100 MHz, the observation points are defined as and. The vectors positions and directions of the antennas are:,

, , ,

and;,

, , , and

.

Note that the mutual couplings between the transmitting antenna and the reception points are negligible, since the distances between the points of observations and the antennas are in far-fields. The amplitudes of the mutual impedances of the antennas in transmission and reception are considerably reduced if the distances between them are large [19] . Within this context, we obtain a matrix of complex coefficients s:

Figure 6. Array with 6 ideal antennas, working together for the electric field control of two observation points.

It was considered for the simulation digital transmissions ASK with binary signals with m_{a} < 1. Although other modulations can be employed, we adopted this modulation for being simpler than the others, given that the focus of the simulation is to illustrate the transmission of signals to two different receiving points, sustaining isolation between the points at a same frequency and polarization.

For the simulation, it was also considered as bit 0 an amplitude of electric field equals to 2 × 10^{−3} V/m, and bit 1 an amplitude equals to 4 × 10^{−3} V/m. In the absence of transmission, the amplitude must be zero. A circular polarization was

used in the plane, noting that all antennas need to be excited

to generate the desired electric field.

The simplest condition happens when there is no transmission at both observation points. In this configuration, the desired fields E_{0} and E_{1} are zero, and current feeds of the antennas are zero. For other conditions, Equation (30) was used to obtain the amplitudes and phases of the currents. Table 1 presents all results.

In the transition times of the desired signals, the currents must be changed in these moments. According to the table, all antennas must be fed, even when there is transmission at only one point. This is due to the fact that all antennas should work together to nullify the field in one of the reception points. It is also noted that the electric current amplitudes of main antennas for the point where

Table 1. Results of the antennas feed current for ASK transmission, achieving isolation between the two observations points.

there is no transmission is small compared with the current amplitudes of the main antennas for the point where there is transmission.

Thus, all antennas work together and are responsible for the transmission of all signals. As each antenna carries a small portion of the information, the fading can be mitigated, due to the path diversity.

For illustration, Figure 7 shows the electric fields obtained when the desired signal is bit 1 at point 1 and null at point 2.

Figure 7 shows that the electric fields were generated according to the desired

signals, that is, a circular polarization in the plane, with an

amplitude equals to 4 × 10^{−3} V/m for the point of observation O_{1} (a), and an amplitude equals zero for the point of observation O_{2} (b).

In this way, any electric field is generated at the observation point, regardless of the desired field at another point, resulting in a transmission isolated by spatial separation or spatial multiplexing.

To observe how the signal is degraded when users move away from the observation point used in the design, consider other example using the same conditions, except for, which is approximately 3 km for 100 MHz. In order to obtain the same desired signals at O_{1} e O_{2}, with the new current values obtained by (30), electric fields were simulated for different points

and.

The time-average amplitude of the resulting field at the observation point i is calculated as follows:

(31)

In possession of the value, the error of time-average amplitude of the resulting field regarding the value of the desired field is calculated as:

(a) (b)

Figure 7. Instantaneous electric fields resulting in points of observations O_{1} (a) and O_{2} (b), given in V/m, when the desired signals are bit 1 at O_{1} and null at O_{2}, and the antenna feed current shown in Table 1.

(32)

Instead of calculating the difference of electric field amplitude at each instant of time, (32) was used in order not to consider the propagation delays between observations points. As the polarizations are in different planes, and the observa- tions points are separated by multiple wavelengths, the Equation (32) becomes more appropriate to the intended results. Table 2 shows these results.

In Table 2, the positive or negative sign of the value of the error indicates that the time-average amplitude of the field is increased or reduced. It is evident that, when the desired signals are not zero at the points of observations, the errors are minor compared to the lack of transmission. In both cases, the errors are small, though the observations points used in the design are about 300 m away from the receptions points. Figure 8 shows the resulting electric fields for one case.

Despite that the O_{1n} and O_{2n} points are approximately 300 m from the points O_{1}

Table 2. Errors of time-average amplitude of the resulting electric field for the O_{1n} and O_{2n} points.

(a) (b)

Figure 8. Instantaneous electric fields resulting at points of observations O_{1n} (a) and O_{2n} (b), given in V/m, when the desired signals are null at O_{1} and bit 0 at O_{2}.

and O_{2}, the amplitudes of the electric fields generated and the polarization planes almost did not suffer degradation.

Considering the small decoupling between the receiving antennas due to the difference between the planes of polarization, the received power will be slightly smaller. Consequently, the largest error obtained in the example will have a small reduction.

In fact, spatial regions whose deviations are small are formed around the observation point O_{i}, and the isolation is preserved. Thus, the reception point could move within that region without suffer interference. The region depends on the threshold detection to identify the signals.

Although the distance between the new observations points and those used in the design are large, it represents only a percentage of the distance from the observation points O_{1} and O_{2} to the origin, which equals 3.000 m. Thus, as these points move away from their main antennas, the greater will be the region of reception, where users can move without any significant errors. This is because the amplitude of the electric fields suffers smooth transitions between the two observation points.

Ignoring the effects of fading, the theoretical perfect isolation can be maintained if the system knows the locations of all the reception points. So with the new position of the observation point, the matrix of complex coefficients s could be updated, if necessary and new values of currents are obtained using (30). The sensitivity of the signals with respect to various distances from the target points, as well as case examples with pyramidal horn antennas, is discussed in depth in [21] . Furthermore, this approach can be of substantial interest when considering the fact that cellular systems are moving towards micro and phantom cells [14] .

5. Torus-Knots Polarizations and Fields

The theory and possible applications related to torus-knotted electromagnetic fields is discussed in [16] [17] , which can be encountered in nature around celestial objects and are used for plasma confinement and with nuclear magnetic resonance devices. Here we only comment for completeness that antenna array configurations similar to the ones presented in the previous sections can also be employed to generate these distributions in a controlled and consistent fashion. Furthermore, we show that not only the torus-knotted fields but polarizations as well can be obtained with the procedure [21] . Figure 9 and Figure 10 show an example of torus-knots for an electromagnetic three-dimensional polarization distribution (as time varies). The synthesis procedure is herein omitted for simplicity, as well as the torus-knotted field distributions (at a given instant of time), and the reader is referred to [21] for more information.

6. Conclusions

This work extended substantially a procedure introduced by the authors in [8] to generate three-dimensional (3D) polarizations, where the electric field intensity

Figure 9. An illustrative example of torus-knotted polarization in the xy-plane, produced by an array of 3 ideal dipole antennas. Although not herein shown similar distributions for the electromagnetic fields at a given instant of time can also be obtained [21] .

Figure 10. An illustrative example of torus-knotted polarization (three-dimensional view of Figure 9). Although not herein shown similar distributions for the electromagnetic fields at a given instant of time can also be obtained [21] .

vector can be oriented in any direction in space. These polarizations were referred to as “3D” due to their distribution in space, although time is also inherent. It was shown analytically and numerically that it is possible to generate 3D polarizations with different shapes, as time varies, such as cylinders, spheres, and others. In contrast to [8] where only ideal dipoles were employed for generating 3D polarizations, here pyramidal horns were also used, not only for polarizations but also to generate controlled electromagnetic field distributions. In addition, the orthogonality between 3D polarizations was also mentioned. It is worth noting that for the scenarios considered in this study, there is no propagation, so the 3D polarization needs to be visualized as the polarization of the resulting electric field, not the wave. Alternate configurations that allow the 3D polarization to move within a desired path are currently under investigation (by determining different sets of feed currents for each desired location, for example), but the basic theory as presented here does not lose its generality. In fact, more accurate results can be obtained with the numerical modeling using commercial and/or in-house antenna codes, which will also allow determining the correct excitations to compensate for effects not accounted here, such as mutual coupling in smaller cells, near-field effects for non-ideal antennas, reflections, refractions, etc…

Although the main purpose was to introduce a mathematically consistent synthesis procedure for controlling the three-dimensionally distributed electromagnetic fields and polarizations, with applications yet to be adapted or discovered, an effort was made to illustrate that the technique can be used for digital transmission of information in wireless communication systems with spatial multiplexing. For each observation point, three antennas will be necessary. If more antennas are used, the fading will be smaller and the results will be closer to the desired. Volumetric regions are created with this array, and isolation could be preserved if the reception points are within these regions. The interference for spatial multiplexing can be further reduced if the method of electric field control is combined with the use of orthogonal 3D polarizations between close users, as in nano and phantom cellular systems [14] . Finally, another possible use of the theory was briefly mentioned for torus-knotted electromagnetic fields, as well as illustrated for three-dimensional polarization distributions. Torus-knotted fields can be found naturally around various celestial objects, as historically shown in Radioastronomy as well as possibly to be further investigated with advanced Radiotelescopes [15] , and can be of importance for plasma confinement and nuclear magnetic resonance devices [16] [17] . Measurement procedures for three-dimensional electromagnetic fields in space (at a given instant of time) and polarizations (as time varies) are subject of further discussion, as for mathematical multi-dimensional results.

Conflicts of Interest

The authors declare no conflicts of interest.

[1] | Stutzman, W.L. and Thiele, G.A. (2012) Antenna Theory and Design. 3rd Edition, John Wiley & Sons, Hoboken. |

[2] | Paul, C.A. (2006) Introduction to Electromagnetic Compatibility. 2nd Edition, John Wiley & Sons, Hoboken. |

[3] | Terada, M. (2005) Reflector Antennas. In: Chang, K., Ed., Wiley Encyclopedia of RF and Microwave Engineering, 5, 4450-4474. |

[4] |
Forenza, A. and Heath, R.W. (2006) Benefit of Pattern Diversity via Two-Element Array of Circular Patch Antennas in Indoor Clustered MIMO Channels. IEEE Transactions on Communications, 54, 943-954. https://doi.org/10.1109/TCOMM.2006.873978 |

[5] | Zooghby, A.E. (2005) Smart Antenna Engineering. Artech House, London. |

[6] | Stutzman, W.L. (1992) Polarization in Electromagnetic Systems. Artech House, London. |

[7] |
Terada, M.A.B. (2003) Novel Beam Scanning Technique for Circularly Polarized Reflector Antennas. Proceedings of the IEEE International Symposium on Antennas and Propagation, Columbus, 22-27 June 2003, 770-773. https://doi.org/10.1109/aps.2003.1219349 |

[8] | Pereira, L.P.S. and Terada, M.A.B. (2012) Three-Dimensional Electromagnetic Polarization. Brazilian Symposium of Microwaves and Optoelectronics (MOMAG 2012), 1, 1-5. |

[9] | Terada, M.A.B. (2004) Distributed Three-Dimensional Array Configuration for Volumetric-Controlled Field Coverage: A New Concept for High-Efficiency, Distance-Selective Wireless Local Area Network. URSI-B Symposium on Electromagnetic Theory (URSI_EMTS), Pisa, 23-27 May 2004, 975-977. |

[10] |
Terada, M.A.B. (2004) Distance-Selective Wireless Network with Enhanced Isolation Characteristics: A New Concept Based on Three-Dimensional Wave Synthesis for Volumetric-Controlled Field Coverage. Proceedings of the 2004 IEEE International Symposium on Antennas and Propagation, Monterey, 20-25 June 2004, 2, 1736-1739. https://doi.org/10.1109/APS.2004.1330532 |

[11] | Piestun, R. and Shamir, J. (2002) Synthesis of Three-Dimensional Light Fields and Applications. Proceedings of the IEEE, 90, 222-244. |

[12] | Berkhout, A.J. (1988) A Holographic Approach to Acoustic Control. Journal of the Audio Engineering Society, 36, 977-995. |

[13] |
Pereira L.P.S. and Terada, M.A.B. (2011) New Method for Optimum Design of Pyramidal Horn Antennas. Journal of Microwaves, Optoelectronics and Electromagnetic Applications, 10, 266-277. https://doi.org/10.1590/S2179-10742011000100025 |

[14] |
Wanderley, P. and Terada, M. (2012) Assessment of the Applicability of the Ikegami Propagation Model in Modern Wireless Communication Scenarios. Journal of Electromagnetic Waves and Applications, 26, 1483-1491. https://doi.org/10.1080/09205071.2012.703077 |

[15] |
Young, André, Terada, Marco A.B., de Villiers, Dirk I.L. and Davidson, David B. (2012) Assessment of the Sensitivity of the South African KAT-7 and MeerKAT/SKA Radio Telescope Reflector Antennas. IEEE 2012 International Conference on Electromagnetics in Advanced Applications, 1, 486-489. https://doi.org/10.1109/ICEAA.2012.6328676 |

[16] |
Arrayás, M. and Trueba, J.L. (2015) A Class of Non-Null Toroidal Electromagnetic Fields and Its Relation to the Model of Electromagnetic Knots. Journal of Physics A: Mathematical and Theoretical, 48, 1-14. https://doi.org/10.1088/1751-8113/48/2/025203 |

[17] |
Irvine, W.T.M. and Bouwmeester, D. (2008) Linked and Knotted Beams of Light. Nature Physics, 4, 716-720. https://doi.org/10.1038/nphys1056 |

[18] |
Mathcad (2017) www.ptc.com/engineering-math-software/mathcad |

[19] | Balanis, C.A. (2005) Antenna Theory: Analysis and Design. 3rd Edition, John Wiley & Sons, Hoboken. |

[20] | Chapra, S.C. and Canale, R.P. (2009) Numerical Methods for Engineers. 6th Edition, McGraw Hill Science/Engineering/Math, New York. |

[21] | Pereira, L.P.S. (2015) Synthesis of Antennas for Field and Polarization Control. Ph.D. Dissertation, University of Brasilia, Brasilia. |

Copyright © 2020 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.