A lot of modern antenna systems are functioning by condition of a different series of various requirements to their radiation characteristics. One of such requirements is forming the deep minimums or zeros in the given points of near or far zones of antenna. This requirement corresponds to issues of the electromagnetic compatibility (EMC). Of course, all antenna systems are working at the condition that is characterized by the existence of various types of radio electronic systems and devices in its neighboring environment. Many engineering problems appear in this connection. The solution to these problems should provide the effective functioning neighbor antenna radiating systems that are placed closely. Therefore, the need for solving two interrelated problems of EMC appears. They are distribution of frequencies while working in the assigned frequency range, and a constructive improvement of the radio electronic system’s individual elements.
The EMC problems were stated and discussed since 80th of last century. In this regard, it should be noted pioneer monograph  , and many journal papers (see, for example     and references therein). The last papers appurtenant to this problem show that the development of the methods and means for providing the certain requirements to the field characteristics of antenna systems in the near and far zones is actual now. This leads to necessity of development of the modeling tools and respective software for solving the practical engineering problems related to EMC of antenna systems.
Some general recommendations how to improve EMC issues for radiating systems are set out in . Different technical solutions and measurements for improving the EMC requirements in the neighboring environment are discussed in papers   . Application of numerical techniques in order to provide the series of restrictions on the radiation characteristics of antennas on the whole   , to RP only   or to characteristics radiation in a near zone  gives possible to elaborate the general vision to EMC problem solving. However, a comprehensive solution to the EMC problem for radiating systems and antennas is far from complete. The generalized variational approach for solving the array’s synthesis problems, proposed in this paper, is one of the means that allows to solve this problem at a more fundamental level. The mathematical statement of synthesis problem including the EMC specifications for plane arrays and results of modeling are discussed in the presented paper.
The paper is organized as follows.
The generalization of usual synthesis problem for plane array allowing to prescribe the restrictions on the radiation characteristics in a near zone is proposed and explanations to terms of formulated functional are given in Section 2. The formula for RP of array (array factor) is discussed with introducing the generalized angular coordinates. In Section 3, the procedure of obtaining the nonlinear integral equation for optimal distribution of currents in array elements, corresponding to minimizing functional, is described. Reduction of explicit formula for components of electromagnetic (EM) field in a near zone of array is made in Section 4. Section 5 and Section 6 are devoted to methods of creating zeros (deep gaps) in RP of array and values of EM field in a near zone. The computational results of modeling for rectangular and hexagonal arrays are presented and discussed in Sections 7; Section 8 contains concluding remarks and proposal for the future investigations.
2. Variational Statement of Synthesis Problem
The functional aimed to take into account a series of requirements to both the RP and field distribution in a near zone is formulated as sum of several terms    , each of them ensures the best approximation to the array’s radiation characteristics 
function is prescribed amplitude RP (by the condition of problem, it is real and positive), characterizes the module of RP created (synthesized) by array, is a weight function that aimed to formulate specific demands to the module of , in particular, to create zeros.
The second term in (1) provides an opportunity to minimize the values of field’s amplitude in the fixed points of array’s near zone. Such an assignment is fulfilled by means of increasing the value in the above points in comparison with the rest of them. Value is amplitude of given field in a near zone, by this is field created by array. The value of additional optimizing parameter c is determined as .
The last term in (1) is applied to restrict the amplitude of currents in the array’s elements. This provides to solve the problem to make away with the effect of super-directivity    for array. The additional optimizing values are applied to restrict the amplitude of currents in the array’s radiators.
We use for RP (array factor) of rectangular array the next formula 
the RP of hexagonal array is determined similarly . are the currents in array’s radiators, are the RPs of separate radiators.
We will use below the generalized angular coordinates instead of usual . This provides ability to present the dependence of RP of array on the one generalized parameter that includes the geometrical size of array and length of wave or frequency simultaneously. Such coordinates are , where , , is wavenumber, is length of wave, values and correspond to distance between elements of array in Ox and Oy directions. It is assumed that the given RP F differs from zero within the angles .
3. Numerical Minimization of Functional
In order to find the minimum of (1), we apply the approach of the variation calculus that yields in receiving the formula for the gradient of functional. Using such formula, we pass to non-linear integral equation  with respect to unknown currents in the array’s elements.
The notations can be simplified in a great extent if to use the operator form as for synthesized RP. For this goal, the notion of adjoint operators and is applied too.
The nm-th component for gradient of (1) on the currents can be presented as
Here, operator B presents explicit relation for the field’s components in a near zone. It is determined similarly to A, the difference is that the additional multiplier, depending on the distance r to array is presented. Equating gradient to zero, we receive the non-linear vector-matrix equation for the optimal currents in array’s elements
In practice, we do not deal with solving (4) because in the iterative procedure of its solving there is necessary to calculate the actions of A, B, and operators. If one to act on (4) by operator A, we have the nonlinear integral equation of Hammerstein’s type with respect to RP f  , because the , are included into respective exponent.
The explicit relation for operator is 
In the last formula procedure of integration in limits and is replaced by twice integration over . This is result of use of new coordinates and instead of initial and . Determination of operator is realized similarly.
If one to use the subsequent application of A, , as well as B, operators and operator form (2) of RP, the numerical calculations are much simple. This allows to receive the nonlinear integral equation relatively RP and explicit formulas for the kernel of this equation .
4. Determination of Field Components in Near Zone
In order to carry out with explicit relation for the field components in a near zone using the currents in the radiators of array, formula for the electric vector potential 
is applied. Here
that is we will restrict ourselves to the value of order . In the case of plane array
so far as the array is placed in the plane.
For plane array, can be written as
Here represent the vector current distributions of array’s elements; by this radius is replaced by for all array’s radiators. Without loss of generality, one can assume the separation for currents onto two axes and , and to consider values of for respective components independently.
Following the above assumptions, we will accept the same relations for the E and H field components in a near zone. As result, we have
In such a way, we have obtained relations (9), (11) or (9), (12) for determination E and H field components. These subsequent formulas can be considered as the linear operator B. This is because of linear dependence of E and H on the distribution of currents . In such consideration, operator is determined by introducing the respective Hilbertian space, as it was done to obtain .
5. Creation of Zeros in the RP
We rewrite (2) as
Function is periodical, namely
that is important for consideration below.
Formulas (13)-(15) are used for presentation of RP using the Kotelnikov’s formula  and forming zeros in the given points . For this goal, we present (13) as
We use the greed
in the region of and variables.
Kotelnikov’s formula for (17) is
Substituting (18) into (13), we receive
Finally, the Kotelnikov’s formula for plane array is
, are the values of initial PR in the interpolation points. Of course,
Let us minimization point coincides with one of interpolation point, as example . We determine the function in the interpolation points as
The problem to synthesize zero is formulated as
where is RP with formed zero.
The condition of stationarity of (23) leads to formula for currents
In fact, values (24) are the coefficients of discrete Fourier transform of function . More simple formula for optimal currents is 
In the case, when point does not coincide with any interpolation point and , , optimal can be determined by the iterative process
where j is iteration number.
Iterative process (26) converges because 
6. Creation of Zeros in a Near Zone
Let us . The weight function is used to create zero in the certain point of near zone. As a rule, neighbor of this point is used, and integration within is replaced by integration within .
Minimizing (1), we minimize too, and we have optimal currents . The components of electric vector potential (9) are determined by known . After this, we determine the rest of components of EM field. The radiation power density (RPD) is determined as 
7. Computational Modeling
7.1. Creating Zeros in RP
Firstly, we show the computational results for the hexagonal array that has 127 elements. By this we assume known RPs of array radiators. The computations are carried out using free MathWorks software.
Results are shown for frequency equal to 11.99 GHz, distance between radiators in the linear subarrays is equal to 0.0078 m, and distance between the linear subarrays is equal to 0.0087 m. The prescribed amplitude RP in the given range (Figure 1). In Figure 2, Figure 3, the results that show forming zeros in the vicinity of the first sidelobe are shown for E-plane of array. Using the weight function in functional (1), we form zeros in points and .
Figure 1. Front view of the prescribed RP.
Figure 2. The initial synthesized amplitude RP.
Figure 3. The synthesized RP with zero in the area of first sidelobes, .
The level of the first sidelobe was reduced to −52.6 dB (Figure 3). The obtained results testify that such a way to create zeros has some limitation. Use of Kotelnikov’s procedure (25) or (26) is more complicate, but it can provide much small zeros.
7.2. Creating Zeros in a Near Zone
Secondly, the results for antenna with rectangular placement of radiator are considered. Array that consists of 121 radiators is examined, and distance d between them is equal to 0.01 m, that is square of array is equal to 0.01 m2; the work frequency is equal to 5.976 GHz, by this the condition is fulfilled that allows to form the main lobe of RP without several maximums (see  ). In Figure 4, values of the RPD are shown without forming zeros at the distance 0.5 m of array.
For this case, RPD changes from 0.16 W/m2 to 0.207 W/m2. The requirement to form zero was stated for the distance of array at the next points: , . To obtain minimal values of RPD in these points, the field in (1) was equated to zero; herewith values of the weight function were decreased on one order greater in the comparison of rest values of and angles.
Such modification of initial data allows to decrease the RPD from quantity of 0.1831 W/m2 to quantity of 0.0373 W/m2 that is 21.7% of its initial value. In the rest range of and values, the RPD increased not more than on 10.2% (see Figure 5). The obtained results testify that the field values and its structure outside the considered and for such r is close to constant value.
In Figure 6, the results related to minimization of the RPD in one point are shown. Similarly to the previous example, the RPR values diminish significantly.
If distance r to array’s plane increases the values of the RPD decrease quickly. The respective results are shown in Figure 7 for the distance that is equal to to array.
The presented results testify that the RPD diminishes in a great extent for such parameters of array; by this its values are more two order lower than for values corresponding to . For example, the maximum of RPD for reaches . It confirms that area of small values of the RPD is placed at distance equal about to 6.0m from the plane of array. This is important data for the creating zeros for the procedure of design of such
Figure 4. The initial distribution of the RPD at .
Figure 5. The RPD with forming zero in two points at .
Figure 6. The RPD with forming zero in one point at .
Figure 7. The RPD at the distance 10.0 m.
type of arrays. The form and values of the main and sidelobes of the array’s RP are varied few in comparison with the case when the synthesis problem is considered without the creating of zeros. The main lobe is wider at level of 0.5 power of the RP on 4.7%, and the level of sidelobes wherein increased on 3.7 dB.
One more requirement to increase the array’s performance consists of the ability to take into account the additional limitations on the amplitude and phase distribution of the currents in the array elements. This is can be carry out by the prescribing the weight multipliers instead of fixed t . For goal to decrease , the values should be increased. The numerical calculations demonstrate that such a way is suitable to avoid the superdirectivity of array’s elements to a large degree. The weight multipliers are prescribed so that the amplitude of currents is close to constant distribution.
The variational approach for solving the antenna synthesis problem by the given amplitude RP is modified for the goal to take into account the additional restrictions on the array radiation in the near and far zones. The optimizing functional is formularized in such a way that the deep zeros can be created in the points of near and far zones.
The results of the numerical calculations demonstrate that the level of formed zeros less than −80 dB in the area of far sidelobes of the RP can be reached. By this, the array’s RP in the main lobe varies insignificantly, and difference of the obtained and initial RP’s values does not exceed 3.7 dB. The minimization of the field at the distance about 10.0 m from array does not lead to the considerable variation of the array’s RP.
The future research in the topic under investigation is designated for consideration of vector character of array’s RP and non-transversal character of EM field in a near zone of array, taking into account the mutual coupling of array’s radiators, as well as at development of powerful numerical methods and respective software for modeling the array’s radiation properties in the real physical environment.
 Hansen, D. (2003) Review of EMC Main Aspects in Fast PLC Including Some History. 2003 IEEE International Symposium on Electromagnetic Compatibility (EMC ‘03), 1, 184-192.
 Lezynski, P. (2018) Random Modulation in Inverters with Respect to Electromagnetic Compatibility and Power Quality. IEEE Journal of Emerging and Selected Topics in Power Electronics, 6, 782-790.
 Lambrecht, N., Pues, H., De Zutter, D. and Ginste, D.V. (2017) Modeling of Contact Bounce in a Transient Electromagnetic Compatibility Test for the Analysis and Optimization of Nonlinear Devices. IEEE Transactions on Electromagnetic Compatibility, 59, 541-544.
 Liang, C.-H., Dang, X.-J., Wang, N. and Yuan, H.-B. (2011) Generalized Isolation between Antennas for EMC Problems in Complex EM Environments. IEEE Transactions on Electromagnetic Compatibility, 53, 645-652.
 Foged, L.J., Scialacqua, L., Mioc, F., Saccardi, F., Iversen, P.O., Shmidov, L., Braun, R., Quijano, J.L.A. and Vecchi, G. (2013) Echo Suppression by Spatial-Filtering Techniques in Advanced Planar and Spherical Near-Field Antenna Measurements [AMTA Corner]. IEEE Antennas and Propagation Magazine, 55, 235-242.
 Perez, J.R. and Basterrechea, J. (2007) Comparison of Different Heuristic Optimization Methods for Near-Field Antenna Measurements. IEEE Transactions on Antennas and Propagation, 55, 549-555.
 Bogdanov, F.G., Karkashadze, D.D., Jobava, R.G., Gheonjian, A.L., Yavolovskaya, E.A., Bondarenko, N.G. and Ullrich, C. (2010) Validation of Hybrid MoM Scheme with Included Equivalent Glass Antenna Model for Handling Automotive EMC Problems. IEEE Transactions on Electromagnetic Compatibility, 52, 164-172.
 Parrikar, R.P. and Gupta, K.C. (1998) Multiport Network Model for CAD of Electromagnetically Coupled Microstrip Patch Antennas. IEEE Transactions on Antennas and Propagation, 46, 475-483.
 Foged, L.J., Scialacqua, L., Saccardi, F., Quijano, J.L.A. and Vecchi, G. (2014) Application of the Dual-Equation Equivalent-Current Reconstruction to Electrically Large Structures by Fast Multipole Method Enhancement [AMTA Corner]. IEEE Antennas and Propagation Magazine, 56, 264-273.
 Ferrieres, X., Parmantier, J.-P., Bertuol, S. and Ruddle, A.R. (2004) Application of a Hybrid Finite Difference/Finite Volume Method to Solve an Automotive EMC Problem. IEEE Transactions on Electromagnetic Compatibility, 46, 624-634.
 Quijano, J.L.A. and Vecchi, G. (2010) Near-and Very Near-Field Accuracy in 3-D Source Reconstruction. IEEE Antennas and Wireless Propagation Letters, 9, 634-637.
 Schnattinger, G., Mauermayer, R.A.M. and Eibert, T.F. (2014) Monostatic Radar Cross Section Near-Field Far-Field Transformations by Multilevel Plane-Wave Decomposition. IEEE Transactions on Antennas and Propagation, 62, 4259-4268.
 Galindo-Israel, V. and Rahmat-Samii, Y. (1981) A New Look at Fresnel Field Computation Using the Jacobi-Bessel Series. IEEE Transactions on Antennas and Propagation, 29, 885-898.
 Yang, T.C. (2018) Performance Analysis of Superdirectivity of Circular Arrays and Implications for Sonar Systems. IEEE Journal of Oceanic Engineering, 43, 1-11.
 Diao, J. and Warnick, K.F. (2018) Practical Superdirectivity with Resonant Screened Apertures Motivated by a Poynting Streamlines Analysis. IEEE Transactions on Antennas and Propagation, 66, 432-437.
 Andriychuk, M.I. (2000) Investigation of Solution of the Nonlinear Synthesis Problem for the Waveguide Array. International Seminar/Workshop on Direct and Inverse Problems of Electromagnetic and Acoustic Wave Theory, Tbilisi, 47-51.
 Andriychuk, M.I. and Savenko, P.O. (2000) Synthesis of a Waveguide Array with Due Regard for the Mutual Coupling of Radiarors. International Conference on Mathematical Methods in Electromagnetic Theory, Kharkiv, 12-15 September 2000, Vol. 2, 604-606.
 Savenko, P.O. and Anokhin, V.J. (1997) Synthesis of Amplitude-Phase Distribution and Shape of Plane Antenna Aperture for a Given Power Pattern. IEEE Transactions on Antennas and Propagation, 45, 744-747.
 Andriychuk, M.I. and Voitovich, N.N. (2013) Antenna Synthesis According to Power Radiation Parttern with Condition of Norm Equality. International Seminar/Workshop on Direct and Inverse Problems of Electromagnetic and Acoustic Wave Theory, Lviv, September 23-26 2013, 137-140.