International Journal of Electrical and Computer Engineering (IJECE) Vol. No. October 2025, pp. ISSN: 2088-8708. DOI: 10. 11591/ijece. On design of a small-sized arrays for direction-of-arrivalestimation taking into account antenna gains Ilia Peshkov. Natalia Fortunova. Irina Zaitseva Department of Physics. Radio engineering and Electronics. Bunin Yelets State University. Yelets. Russia Article Info ABSTRACT Article history: In the paper a technique for designing antenna arrays composed of directional elements for direction-of-arrival (DOA) estimation is proposed. Especially this approach is applied for developing hybrid antenna arrays with increased accuracy which features digital spatial spectral estimation after preliminary analog beamforming. The earlier obtained explicit formula for calculating the CramyrAeRao lower bound (CRLB) which determines the relationship between the variance of the DOA-estimation and antenna elements' radiation patterns, array geometry, has been used. Main idea of the proposed technique is that it takes into account spatial pattern and gain of the antenna elements. The high gain unlike the number of the antenna elements or interelement distance is the most important factor which allows reducing the value of the DOA-estimation errors. A couple of the examples of calculating radiation patterns of antenna elements improving accuracy of DOA-estimation with super-resolution are provided in the paper. Proposed antenna arrays are modeled according to the method of moments (MoM). The values of the root mean square error after the DOA-estimation are It is shown that the resulting hybrid systems can reduce the error value in DOA-estimation with super-resolution. Received Aug 24, 2024 Revised Apr 14, 2025 Accepted Jul 3, 2025 Keywords: Cramer-Rao lower bound Direction-of-arrival estimation Hybrid antenna array Smart antenna Super-resolution This is an open access article under the CC BY-SA license. Corresponding Author: Ilia Peshkov Department of Physics. Radio engineering and Electronics. Institute of Mathematics. Nature Science and Technique. Bunin Yelets State University 39 Kommunarov street. Yelets 399770. Russia Email: ilvpeshkov@gmail. INTRODUCTION The state-of-art technologies actively use multiantenna wireless systems. The analysis of the spatial spectrum of signals underlies these devices . DOA-estimation is critical in applications requiring spatial awareness, minimal delay and exceptional reliability. For example, in 5G direction-of-arrival estimation contributes in beamforming by identifying user directions for optimal signal transmission . Detecting angles of incoming signals from obstacles/vehicles is utilized in automotive radar . elf-driving car. Low accuracy and resolution are the most serious problems of these systems using methods and algorithms for estimating spatial coordinates, which are limited by array aperture . Thus, the primary goal of the paper is how to improve the accuracy of DOA-estimation algorithms. However, the accuracy of the estimates cannot be increased in a easy way and has a lot of challenges. Today, there are several basic approaches regarding the accuracy of DOA-estimation. The most common way the researchers apply to improve accuracy is to enlarge the aperture with additional elements . The attitude is widespread in the millimeter range applications or military radars . Nevertheless, the complexity of the multi-antenna system architecture increases which results in a dramatic escalation of overall device costs . Journal homepage: http://ijece. Int J Elec & Comp Eng ISSN: 2088-8708 Another way improving the accuracy which has gained particular popularity is related to the design of an innovative DOA-estimation algorithm. Usually. It is tailored to particular use-case requirements . or instance, a multipath propagation channe. , a certain geometry of antenna array . uch as linear antenna array and estimation of signal parameters via rotational invariance technique (ESPRIT) . ), or for a signal waveform . o-called blind method. The main disadvantage of this approach is that new algorithms possess high computational complexity. For example, the well-known MUSIC algorithm has O(N. , at the same time ESPRIT and others O(N. depending on the number of antennas . In other words, the load increases a lot. This fact makes it difficult to implement for real-time applications. A similarly widespread technique increasing accuracy is based on the optimal arrangement of antenna elements . Several different criteria have been introduced such as random distribution of antenna coordinates . or complex coordinate expression of antennas . , or using arrays composed of omnidirectional radiators . , . The main disadvantage of the mentioned approaches is the consideration of ideal array patterns because omnidirectional antennas are absent in real applications. In other words, they are hardly implemented. The main idea of the proposed methodology for enhancing DOA estimation precision is the optimal placement of directional antennas, as well as their spatial patterns. The novelty of the proposed approach arises because that the directional spatial patterns of antenna elements are taken into account. So that the procedure for designing a dual-element array configuration for estimating spatial parameter estimation by means of manually forming the beam patterns and coordinates of the array elements to greatly reduce the In this way, the proposed approach highlights the gain of antenna elements enabling the sensor array to be developed for 360A range scanning by means of synthesizing the radiation patterns. Microstrip antennas can be designed using the method of moments by demonstrating the practical implementation of the presented method. Furthermore, the paper describes an example of creating an antenna array prototype for super-resolution DOA-estimation, which is computed using the closed-form analytical formulations of CramyrAeRao lower bound (CRLB) comparing to the optimization or neural networks which are black box solutions. The elaborated array consists of an analog beamforming scheme, which allows reducing the number of digital processing channels without decreasing the accuracy of DOA-estimation . Thus, the paper discusses the step-by-step process of antenna array design, starting from a concept based on analytical expressions and ending with the practical implementation on the basis of patch antennas. completed device is used to demonstrate the feasibility of the proposed approach. Especially, there is no unified methodology for designing that kind of devices which are built on sub-array preprocessing networks . or beam selection schemes . SYSTEM MODEL DOA-estimation formulation In this section the model of antenna systems widespread in direction finding is described. The arrays can be configured in different ways influencing on such properties as accuracy and resolution. Figure 1 shows such an array, which consists of a certain number of antenna elements. 2, y2, z. n, yn, z. (I1. E1, z. (I3. E3, z. (I4. E4, z. Figure 1. Antenna array view in the Cartesian system In the paper it is assumed that the signal is narrow-band and has the following spatial coordinates on azimuth and elevation I relative to the x, y, and z axes in the Cartesian system, respectively. Thus, the coordinates and I has to be estimated with maximum accuracy. The analytical model of the antenna array is expressed in the following manner . ycN yca. uE, yu. = . uE, yu. yce ycycoyc1 ycN U yciycA . uE, yu. yce ycycoycycA ] . On design of a small-sized arrays for direction-of-arrival-estimation taking A (Ilia Peshko. A ISSN: 2088-8708 2yuU where yco = . coycu , ycoyc , ycoyc ) = . inIcos, sinI ycycnycu yuE , cosI) is the spatial frequency specifying the oscillating yuI of the phase of the propagating wave in the x, y, z directions, ycycuycN = . cuycu , yycu , zycu )ycN is the position vector indicated to the nth antenna element and gn(. I) is the pattern of the nth antenna element. The antenna element output signals are represented by the complex vector . = ya. uEyco ) UI ycE. ycuEE. , where ycuE. = . , . , ycuycA . ]ycN is the vector of the dimension 1yN, which describes the antenna array output signals. = . , . , ycycA . ]ycN is the N-dimensional of the signals. ycuEE. = . , . , ycuycA . ]ycN is the noise vector. uEyco ) is the NCM dimensional matrix of the steering vectors . uE1 ). A , ycaE. uEycA )]. In practical applications, the spatial covariance matrix is estimated from a collection of K time samples . ycIC = Ocya ya ycuE. = yaCyc yuCyc yaCycya yaCycA yuCycA yaCycAya , yco=1 ycu ya where ycuE. is the N-dimensional vector at the k-th time sample at the digital antenna array output, the symbol AuEoAy describes the averaging over K samples, yaCycI is the signal eigenvector matrix, yaCycA is the noise eigenvector matrix, yuycI and yuycA is the eigenvalues matrices. Spatial spectrum by the method MUSIC is calculated as . ycEycAycOycIyaya . uc, yuE) = EEya yayca EE. uE) yca . uE)yaycA yaycA Reduction of DOA-estimation variance The variance of the bearings and I is assessed using the Cramer-Rao lower bound criterion, which is determined by the amount of noise, the positioning of the antenna elements in space and their radiation In this case, an arbitrary direction-of-arrival estimation algorithm cannot have values of variance below this limit, but only approaching it. Covariance error matrix for estimating the angular coordinates of signals with super-resolution in both azimuthal and elevation scanning planes can be written as . yaAycIycNycC = yu Eu . cNyc {[ 1 ]Oo[ yu4 Oe1 ] }] here yu1 = yayuEya ycEyaOu yayuE , yu2 = yayuEya ycEyaOu yayuo , yu3 = yayuoya ycEyaOu yayuE , yu4 = yayuoya ycEyaOu yayuo , yu = ycIyaya ycIOe1 yaycI. D DI are the matrices of steering vector differentiation a(. I) along the corresponding planes. K is the number of the Consequently, it is possible to reduce the values BSTO by influencing the coordinates of the antennas and the shape of their radiation patterns of DOA-estimator. Therefore, let us examine the CRLB expression given in . in greater detail. Previously, an explicit general formula of the CRLB was obtained for estimating the coordinates of the radiation source via Multisensor antenna arrays with directional sensors, also oriented in space in an arbitrary manner . ycycayc( yuE, yu. OO yua2 2ya 2 Oe1 Eu {(Ocycnyc yciycn2 yciyc2 . caycnA Oe ycaycA ) Ocycnyc . ciycnA yciyc Oe yciycA yciycn ) ) } . where ' denotes a derivative along or I depending on the scanning direction, i and j are the antenna elements indexes. The antenna array is supposed consisting of two antenna elements, then the equation . can be rewritten in the following form . uE, yu. = yua2 2ya {. ci12 yci22 . ca1A Oe yca2A )2 . ci2A yci1 Oe yci1A yci2 )2 )Oe1 } . As demonstrated in the equation's left-hand side . in curly brackets, the main physical factors determining the variance are the square of the gains of the sensor elements, just like the area occupied by the employed antenna array. Therefore, it is possible to increase the accuracy of DOA-estimates by changing one or both of these factors. On the contrary, the smallest possible number of antenna elements can be compensated by a larger gain, or by being located on a larger area. For example, if the gain of each antenna is increased twice, then the variance of DOA-estimates can be reduced radically. Int J Elec & Comp Eng. Vol. No. October 2025: 4642-4652 Int J Elec & Comp Eng ISSN: 2088-8708 Therefore, the proposed methodology for designing antenna arrays for DOA-estimation consists of synthesizing the directional patterns of individual elements g1 and g2 in order to minimize variance var(. I in the works familiar to the authors, which concern the optimization of the antenna array topology to increase accuracy, the gain and the radiation patterns of the array elements are not accounted for at all. Consider an example, two circular array configurations are studied: a 3-element and a 2-element array of directional antennas. The task is the variance of the estimates of the three-element array var(. I)3 must be equal to the variance of the two-element one var(. I)2. In this case the sum of the gains of the elements of the second array must be consistent to the sum of the gains of the three-element circular antenna array. The formula of the Cramer-Rao lower bound for the two-element antenna array must be remembered . and expression for the three-element array: I)3 OO 12 g 22 . '1 -a'2 )2 g 23 g12 . '1 -a'3 )2 g 22 g 23 . '2 -a'. '2 g1 -g1' g 2 )2 A ) } . '3 g1 -g1' g 3 )2 . '2 g 3 -g '3 g 2 )2 ) . According to the condition the following equation must be met: uE, yu. 2 = ycycayc. uE, yu. 3 In addition, it is assumed that two- and three-element circular arrays are located at the same radius from the center, as well as the noise power yua 2 and the number of samples K of the correlation matrices are the same for ycycayc( yuE, yu. 2 and ycycayc( yuE, yu. Thus, yua2 ci12 yci22 . ca1A Oeyca2A )2 . ci2A yci1 Oeyci1A yci2 )2 } = yua2 2ya yca yca yca = . ci12 yci22 . ca1A Oe yca2A )2 yci32 yci12 . ca1A Oe yca3A )2 yci22 yci32 . ca2A Oe yca3A )2 . yca = . ci2A yci1 Oe yci1A yci2 )2 . ci3A yci1 Oe yci1A yci3 )2 . ci2A yci3 Oe yci3A yci2 )2 ) . here and further yciycn2 and yciycn3 are the radiative characteristics of the elements of the considered two- and threlement arrays respectively. The equal terms, i. , the noise power and the number of samples of the correlation matrix can be eliminated, therefore it can be written that: yci12 yci22 . ca1A Oeyca2A )2 . ci2A yci1 Oeyci1A yci2 )2 = yci2 yci2. caA OeycaA )2 yci2 yci2. caA OeycaA )2 yci2yci2 . caA OeycaA )2A 1 2 1 3 1 1 2 3 2 ci2A yci1 Oeyci1A yci2 )2 . ci3A yci1 Oeyci1A yci3 )2 . ci2A yci3 Oeyci3A yci2 )2 . Consider only the denominators of the expression . yci12 yci22 . ci22 yci12 Oe yci12 yci22 )2 = yci13 yci23 yci33 yci13 yci23 yci33 . ci23 yci13 Oe yci13 yci23 )2 A A A A . ci33 yci13 Oe yci13 yci33 ) . ci23 yci13 Oe yci23 yci33 ) . Since it is clear that the product of the squares of the radiation patterns is much greater than the other terms, the simplification can be continued and: yci22 OO yci13 yci23 yci33 yci13 yci23 Now it is required to determine the gain of the pattern of the synthesized antenna array. It is hypothesized that the maximum value of the pattern of the rectangular antenna element is equal to 6, and the minimum value to Then it can be written that: 1222 OO 613 OO 13 After simple calculations it can be claimed that yci12 = 8. The value makes it possible to obtain the accuracy using a smaller number of antenna elements. Here, the missing elements are compensated by means of a higher gain, as it can be seen from the CRLB formula . , . On design of a small-sized arrays for direction-of-arrival-estimation taking A (Ilia Peshko. A ISSN: 2088-8708 It may be stated. the proposed approach allows designing antenna arrays having radiation patterns which minimize the variance in spatial coordinate estimates. On the contrary, it can be utilized to evaluate the accuracy of existing arrays by means of the closed form equations. The variance based on . for the spatial coordinate in azimuth equal to 90A can be calculated using . uE = 90Oo )2Oeyceycoycayc_ycycuycy OO uE = 90Oo )3 OO 2 yci2 2 yci2 yci2 yci2 yci2 yci2 OO 0. OO 0. As it can be seen from expressions . , according to the CRLB and the proposed methodology, the variance of the estimates using the dual-element array will promote smaller values in comparison with the trielement digital antenna array. With all this, there will only be mutual relationships between different geometries of antenna arrays, among which comparisons can be made and then the best one can be selected from the point of view of obtaining high accuracy of the DOA-estimation. Now, based on the obtained antenna gain value, a radiation pattern must be synthesized satisfying the requirements. In addition, the proposed hypothesis will be confirmed by fulfilling comparative modeling. The results are given below in the section 3. PERFORMANCE ANALYSIS Simulation study of the proposed methodology In the section the introduced strategy of designing an the dual-element array for DOA-estimation with super-resolution is will be approved via simulation study. A circular antenna array out of three directional elements is used for the comparison. In order to obtain the same accuracy as the three-element Each antenna of the dual-element array will have gain equal to yci12 = 8. 5 as it was proven in the previous part. The following antenna array structure is proposed for DOA-estimation, as represented in Figure 2. As can be seen from Figure 2, each antenna of the array must have radiation patterns which have to be synthesized and possess the gain yci12 = 8. 5 in order to minimize . The proposed patterns decreasing DOA-estimation errors are showed in Figure 3. Antenna Element 2 Antenna Element 1 DOA-estimation processor Figure 2. Proposed digital antenna array design for DOA-estimation with two antenna elements for 360A scanning on azimuth Figure 3. reveals the peak values are picked up in such a way that they coincide with the minimal of the pattern of the neighboring element. At the same time, the main part of the gain is focused in a wide azimuth range. A template radiation pattern is shown in Figure 3. that will be used for antenna synthesis. Further these radiation patterns are utilized to synthesize a linear antenna array, which is consists of ten directional elements with a uniform half-wave interelement distance. The linear antenna element shown in Figure 2 consists of directional antenna elements. Mathematically, the spatial pattern of each element is defined as . ya yciycu . uc, yuE) = 2yco . )yco . uE Oe 2yuUycu ycA ))yco , where D is the directivity, m controls D. Int J Elec & Comp Eng. Vol. No. October 2025: 4642-4652 Int J Elec & Comp Eng ISSN: 2088-8708 The radiation pattern of the linear array consisting of (N = . directional elements without taking into account mutual coupling can be obtained as . cuOe. ycc ycaycuyc. uE yuycu )) ya. uE) = OcycA ycu=1 yciycu . uE) UI ycaycu yce The linear array has inter-element spacing d, where each n-th element is weighted by amplitude aCo and phase Co. Their individual radiation patterns gCo() is characterized by . A genetic optimization algorithm is applied to synthesize the patterns of the entire linear array by means of MATLAB optimization toolbox . The relative error which is difference between the resulting field level and the intended pattern value across the M common sampling points is calculated using the formula . em = Eactual . )-Edesired . ), i = 1,2, . where E. denotes the electric field level of the linear array at a point with azimuth coordinate m. The value of least mean squares is used as the objective function according to the following expression . A , yuycA , yu1 . A , yuycA ) = ( OcycA yco=1 . ceyco | ) . ycA The values of phase n and amplitude excitation an of the linear array feeding network are optimized. Following optimization, the determined parameters . hase shifts Co and amplitude attenuations aC. appear in Table 1. Table 1. The values of phases and amplitudes for linear antenna array Antenna element index . Phase . , deg. Amplitude . As evident from Figure 3. the synthesized linear arrayAos radiation pattern closely matches the reference pattern. An obvious peak in the -5A to 5A azimuth is manifesting, as well as gain from -90A to 90A to achieve intended azimuthal coverage. In addition, it was possible to obtain the radiation pattern of the linear array that are higher than those of the single antenna, such as, for example, a rectangular microstrip. The achieved gain is the basis for reducing the error variance of DOA-estimation as it can be realized from . Thus, the steering vector yca. uE) for the two-element array is given by: uE) = [ yuU cuOe. ycc ycaycuyc. uE yuycu Oe yc yc ycaycuyc. uEOe 2 )yce yuI (OcycA ycu=1 yciycu . uE Oe ) UI ycaycu yce cuOe. ycc ycaycuyc. uE yuycu yc yc ycaycuyc. uE 2 )yce yuI (OcycA ycu=1 yciycu . uE ) UI ycaycu yce where r is the distance between antenna elements, n and an are given in Table 1. Figure 3. Radiation patterns of designed antenna elements . intended radiation patterns . ed and blue are indicated for corresponding antenna element. radiation patterns: blue Ae template, red Ae linear antenna array, yellow Ae common antenna On design of a small-sized arrays for direction-of-arrival-estimation taking A (Ilia Peshko. A ISSN: 2088-8708 Further a simulation is fulfilled in order to verify the hypothesis. The super-resolution algorithm MUSIC . is utilized to estimate the spatial coordinates of signal source. The correlation matrix . uses 100 samples with a 0 dB signal-to-noise ratio. The signal's azimuth coordinate varies from 0A to 180A, and we compute the RMSE at each angular position: ycIycAycIya = yaOe1 ocyayco=1. uE Oe yuECyco ) . Thus, the simulation process goes through several stages. At the first stage, the variance is set using . Then it is necessary to calculate the gain of the antenna pattern in the same way as . in order to obtain the acceptable level of variance from the previous step. The obtained values of g(. I) are substituted in . At the last step, the predefined signal coordinates are calculated according to AAE expressions . The process continues to obtain . Figure 4 represents circular arrays consisting of two (Figure 4. ) and three (Figure 4. ) directional elements respectively. Further they are compared with the array demonstrated in Figure 2. The singleelement radiation pattern of the circular arrays from Figure 4 is obtained from expression 8 for D = 6 db. The radiation patterns of circular arrays of two and three elements are shown in Figure 5() to 5. depicts the obtained flat-top pattern. Figure 5 reveals the circular antenna arrays from Figure 4. have radiation patterns in which local maxima coincide with local minima. Consequently, their product is At the same time, as evident from Figure 5. , the maximum is distributed along the azimuth in order to receive the signal uniformly over the entire range. In addition, the maximum here has a greater value in order to compensate for the minimum gain of the neighboring element. In Figure 6 several curves of RMSE of the MUSIC method are shown for the following types of digital antenna arrays: the two-element array from Figure 2, as well as the circular arrays depicted in Figure 4. Figure 6 demonstrates that the proposed digital array design method reduces RMSE across all azimuth scanning angles while maintaining accuracy comparable to . r better tha. the tri-element circular The RMSE of the estimates of the proposed dual-element antenna array having the synthesized patterns, which is shown in Figure 5(A), has a value of about 1A, which corresponds to the circular antenna array of the two elements (Figure 4. ) having radiation patterns as in Figure 5. Thus, the simulation confirmed the correctness of the proposed method of designing antenna arrays for DOA-estimation. At the same time, the error values obtained analytically from section 2 coincide with the graph in Figure 6. Figure 4. Circular antenna arrays consisting of . two and . three elements . Figure 5. The simulated radiation characteristics of the digital antenna array . depicted in Figure 4. , . illustrated in Figure 4. , and . presented in Figure 3 Int J Elec & Comp Eng. Vol. No. October 2025: 4642-4652 Int J Elec & Comp Eng ISSN: 2088-8708 Figure 6. Root mean square errors of the considered digital antenna arrays Simulation of DOA-estimator prototype Now in this section the approach explained in the section 2 is implemented by means of the Method of moments. The method allows obtaining antennas that are closest to real samples. The example of developing a microstrip antenna array with a synthesized radiation pattern which maximally reduces the variance of the DOA-estimation is described. Consider a two-element antenna array and its elements pattern shown in Figure 7. The circuit in Figure 7. has two antennas, i. AuAntenna 1Ay and AuAntenna 2Ay. An idealized perfect conductor (PEC) material is used as an emitter. Each antenna is represented in Figure 7. and must have a radiation pattern illustrated in Figure 5. This directional pattern can be achieved using several antenna elements, phase shifters and amplifiers that perform the phase shift n and amplitude excitation an. The synthesized radiation pattern of the antenna in Figure 7. is depicted in Figure 7. The optimization is fulfilled according to . The characteristics of the antenna arrays described in this section are obtained using the method of moments by AoMATLAB antenna toolboxAo . The calculated amplitude excitations and phase shifts for obtaining the radiation pattern are given in Table 2. Antenna 1 Antenna 2 Phaseshifters Analog Amplifiers Digital DOA-estimation processor . Figure 7. Two element DOA-estimator . sketch, . its antenna element, and . pattern of an element On design of a small-sized arrays for direction-of-arrival-estimation taking A (Ilia Peshko. A ISSN: 2088-8708 Table 2. The values of phases and amplitudes for antenna array depicted in Figure 7. Antenna element index . Phase . , deg. Amplitude . In addition, the dual- and tri-element circular antenna arrays are considered in Figure 5, which consist of rectangular patch antennas, which are shown in Figure 8. and its pattern in Figure 8. As it can be seen in Figure 8. , the patch element has the gain equal to 4. 7 dBi while the gain of the antenna from Figure 7. is approx. 7 dBi that can be viewed in Figure 7. At the same time, it is concentrated in a wide azimuth range. Thus, the hypothesis is that the designed antenna array will allow estimation direction-ofarrival coordinates with higher accuracy. The following figure shows the graphs of the pseudospectrum of the MUSIC method for the antenna arrays under consideration: the two-element array . lue curv. , the threlement circular . ed curv. and the two-element circular digital antenna array . ellow curv. The coordinate of the signal in the azimuthal plane is equal to 75A. SNR=5 dB. Thus, the steering vector . and the spatial correlation matrix . are calculated using the modified radiation patterns g() after the method-ofmoments. Figure 8. Electrodynamic model of the elements of the circular antenna array . rectangular patch antenna based on RO4730JXR dielectric and . its directional pattern As it can be seen from Figure 9, the two-element antenna array shown in Figure 7. and possessing the elements radiation patterns depicted in Figure 7. , allows obtaining spatial pseudospectra as sharp as the three-element circular antenna array. This is explained by the fact that the elements AuAntenna 1Ay and AuAntenna 2Ay have the gains which are greater than the elements of the circular antenna array, i. , 7 dBi (Figure 7. ) and 4. 7 dBi (Figure 8. ) respectively. Thus, as depicted in Figure 9, higher gain results in higher accuracy in DOA-estimation. Based on the procedure described in sections 2 and 3 the antenna array prototype is implemented and the results are consistent. The implemented scheme can be named as beamspace . , . Figure 9. The pseudospectrum of MUSIC Int J Elec & Comp Eng. Vol. No. October 2025: 4642-4652 Int J Elec & Comp Eng ISSN: 2088-8708 CONCLUSION The paper details a method for designing the digital antenna arrays for increasing DOA-estimation accuracy with super-resolution. Unlike the referenced methods, the proposed design approach accounts for the gain and directional patterns of individual antenna elements. Commonly, they focus on creating new DOA-estimation methods, adding more antennas, or optimizing array layouts, while neglecting antennaspecific properties by treating them as omnidirectional. It has been shown based on simulations and analytical expressions that the higher gain compensates for the lack of antennas. For example, the gain equaling to 8. 5 makes it possible to obtain the same DOA-estimation accuracy as an antenna array composed of more elements having the gain equaling to 6. Additionally, this technique has been utilized for designing a hybrid array architecture for directionof-arrival estimation with super-resolution. Modeling based on the Method-of-Moments incontestably demonstrated that the hybrid array with fewer antenna elements does not result in a decline in the DOAestimation accuracy. As a result, the reduced-element antenna array allows achieving the same level of accuracy in DOA estimation. In other words, the proposed approach can be implemented in practice as beam space hybrid construction. Thus, it can be argued that the considered technique can be used as a theoretical substantiation for the design of sensor arrays for spectral spatial processing. Thus, in the paper the methodology of the design flow of arrays for DOA-estimation with super-resolution has been represented, beginning from the analytical analysis on the closed-form equations up to the prototyping close to experimental solutions. Additionally, it handles the location of the antennas, the radiation patterns of the elements, the value of DOA-estimation errors simultaneously. All that allows moving from an abstract model to a final prototype vividly. REFERENCES