10.15587/1729-4061.2016.66015 THE ANALYSIS OF METHODS FOR DETERMINING DIRECTION OF ARRIVAL OF SIGNALS IN PROBLEMS OF SPACE-TIME

The analysis of methods for estimation of the signals direction of arrival based on the analysis of the directional characteristics of the antenna array according to the criterion of maximum and minimum level of the incoming signal is conducted. It is shown that the minimum methods are more accurate in comparison to the maximum methods. Considering the direct functional connection of the amplitude-amplitude-phase distribution of signals over the array elements with the value of signals direction of arrival, the technique allowing estimation of the field distribution vector together with an estimate of the array weight vector, which showed the computational complexity of the procedure, is presented. A number of superresolution techniques have been analyzed. The results of methods modeling with numerical estimates of resolution have shown that the resolution of these methods allows their implementation in the problems of space-time access under strictly limited conditions. The analysis of the signal arrival estimation algorithm using the sum-differential directional pattern with the definition of the direction finding by means of the equisignal zone is conducted. It is shown that this method seems to be one of the most constructive and applicable to the problems of space-time access.


Introduction
As international literary sources [1,2] claim, by the year 2020 and in future, 5G mobile communication systems will be able to provide mobile users with unlimited high-speed access to information at any place and any time. To achieve the set goal, a considerably large variety of applications and devices is needed and networks of mobile communication and broadband wireless access currently have them. Due to this fact, there emerged a necessity to implement longterm technological methods in 5G systems aimed at solving problems of mobile user access and issues of effective link resources utilization. The key long-term technological solutions implemented in the mobile communication systems of the 5th generation are [1,2]: -application of evolutional massive (multi-dimensional) multi-antenna MIMO technologies; -ability to effectively use the modes of dynamic 3D-beamforming. This will allow considerable increasing the signal power for remote users in high frequency bands and improving coverage in ultradense micro-and picocells; -application of micro-, pico-and femtocells in areas of ultradense user location, which decrease the load on macrocells, with the division of transmission of user traffic and control signals between macro-and microcells in different frequency ranges; -implementation of the full duplex in common bandwidth (transmission and reception are on the same frequencies); -decrease in the level of intra-system interference at the account of using new methods of frequency-space scheduling, coordination of base stations radiation in adjacent cells, improving reception methods; -new methods of multitone modulation, which provide the increase in spectral efficiency with respect to the OFDM technology.
All this explains the rationale and necessity to use and develop space-time access (STA) methods on the basis of implementation of multi-beam antenna arrays (MBAAs). Implementation of the STA methods will increase the effectiveness of radiofrequency resource utilization, link bandwidth and capacity of the channel resources access system.
The main range of problems among the STA tasks with the use of MBAAs includes determining signals direction of arrival from subscriber stations (SSs).
The problems of determining signals direction of arrival can be solved using different methods for estimation of antenna array directional characteristics, which have their own effectiveness. The analysis of the given methods concerning their effectiveness and practicability is a rather rational issue when solving the STA problems.

Literature review and problem statement
The problem of determining the direction of electromagnetic wave (EMW) arrival can both have an independent value (direction finding for radiation sources) and be a composite component of more general problems of radio communication, navigation, etc.
A great number of direction finding techniques have been developed: amplitude, phase, narrow-basis and wide-basis, with Rayleigh resolution, with superresolution (SR), etc.
Estimation methods for direction of arrival (DoA) are considered as a basis of many tasks of telecommunications, including problems of space-time access in mobile communication systems.
On the basis of DoA, in the STA problems it is possible to solve the tasks of setting the main lobe in the DP in the direction of a calling SS during a session. A lot of publications are devoted to the DoA. In [3] the authors consider the problem of estimation for nominal arrival directions of signals from a number of distributed sources using a nested array, which is formed on the basis of the whitening filter. In [4] a new algorithm of finding DoA estimations for several signals correlated in time is proposed. The proposed approach is based on the common diagonalized structure of a set of space-time correlated matrices. In the contrast to subspace methods based on DoA estimates, there is no need to estimate a noise or signal in subspace in their explicit form. In [5] there is a new method of DoA estimation for coherent signals on the basis of decorrelation using the algorithm, in which every line of a covariation matrix can be used to form a full rank. In [6] the authors are oriented at the problem of DoA estimation using the statistics of the second order. In [7] it is proposed to improve DoA resolution within the scope of the OMP (orthogonal multiplexing with search of match) iteration selection algorithm, where a process of searching for a new DoA estimation in the considered sample is performed at each iteration. In [8] the DoA estimation of several sources of wideband signals on the basis of the maximum likelihood method considering the influence of nonuniform noise has been studied; that has given the possibility to obtain the solution closed to the Cramer-Rao criterion under high signal-to-noise ratio. In [9] the detailed comparison of investigation of operational characteristics for the family of "superresolutional" methods is considered to estimate spatial spectrum of Gaussian sources of noise radiation arriving at the antenna array.
Within the studied problematics, the authors of this work have conducted a general analysis of methods for determining the direction of arrival of the SS signals. The goal of the analysis is selection of the effective method for determining signals arrival and its practicability in the STA problems in mobile communication systems.
Analyzing publications devoted to this subject [3][4][5][6][7][8][9], it is worth mentioning that the relevance of the given problem does not decrease because nowadays there is a deficiency in using frequency, time and code resources. There has emerged an urgent necessity for implementation of additional physical resources, which realize the growth of link throughput and the increase in productivity of a mobile communication system in general. In this aspect, implementation of spatial resource, which possesses a considerable effective potential, is considered to be perspective.

Goal and tasks of research
The goal of the work is to analyze methods for determining the signals direction of arrival on the basis of the implemented approach for AA field distribution estimation jointly with the estimation of weight vector of SS received signals.
In order to achieve the goal the following tasks have been set: 1. To analyze the effectiveness of antenna array directional characteristics on the basis of different estimation methods.
2. To analyze the technique for estimation of the vector of AA distribution field together with the estimation procedure of weight vector of signals received by SSs from different directions in the AA.
3. To analyze "superresolution" methods when determining the direction of SS signal arrival and estimation of their theoretical resolution.

Analysis of directional characteristics for antenna arrays in signal arrival estimation
All the methods for determining the signal direction of arrival into the observation point (position of the phase front of electromagnetic wave in the selected coordinate system) rely on certain use of the direction pattern (DP) of the receiving antenna, which is considered as a spatial filter. From the antenna theory it is known [10][11][12] that any antenna has the maximum and minimum of reception (transmission) from a certain azimuth. Using these features of the DP it is possible to find corresponding directions at the antenna output, which are a normal with respect to the phase wavefront, according to the voltage values of max U and min U . The specifics of wireless system implementation will make us interested in the directions of azimuthal angles b. The angles of EMW arrival in the vertical axis a are usually situated along the horizon and in this case their analysis is less informative.
Thus, if DC ( ) f b , where b is the azimuthal angle, it is needed to turn according to the change in this angle, then the voltage at the antenna output will accordingly change, what brings some useful information due to which the DoA can be determined. φ is the complex amplitude of the n-th element excitation; n z cosb is the length difference of paths that come to the origin point.
The reception of SS signals by an antenna of a base station is performed at the main lobe (ML) of DP. It is obvious, that the narrower the ML receiving antenna is, the lesser is the possibility of the impact from undesired signals and interference. The level of a useful signal sent in the direction of the reception point is also growing under narrowing of the ML in transmitting antenna. When the aperture size is L/λ, the minimum value of ML width is 0,5 L θ = λ . According to the specifics of DC f( ) b , two basic methods of DoA are practically used: -the method of maximum when the direction finding b is found by the value max U (where max U is, there is also the direction finding i b ), in addition, n n n kz cos φ = b is determined, n b is the direction finding of the n-th SS; -the method of minimum, when the value min U is also used.
Both mentioned methods have their advantages and disadvantages. We analyze them as follows.
The method of maximum is implemented rather easily using almost any antenna. Alignment of the reception antenna to a correspondent ( s b = b ) is nothing but the maximum principle implementation. However, as a rule, the DP position is more blurred in comparison to the position of DP minimum. The resolution, accuracy in the angle 0 b determining under this method is not high and it has an order of the one sixth of the ML width in the DP [13]: where θ is the ML width in the DP at the level of 0,7 or 0,5 under DP according to power. However, the accuracy can be fixed by changing the parameters (dimension) of the antenna. It is known [13] that the width of antenna ML is connected with its electric length a l and it is determined by the ratio where k A is the antenna utilization coefficient, Thus, for λ=0,15 m and aperture radius is l=2 m, we receive θ=8 o . Obviously, such accuracy may not always be enough for practical application.
The resolution of two signals adjacent toward the signal arrival is connected with the implementation of the method of maximum. Fig. 1 shows the DP graphs of two similar antennas in the orthogonal coordinates, which demonstrate the ability of resolution for two adjacent signals with angles of arrival b 1 and b 2 .
The resolution of two signals is considered as limited for the given antenna with the characteristics f(b) (Rayleigh), if the angle distance Δb=b 2 -b 1 , is so that DC f(b 2 ) and f(b 1 ) intersect at the level, which is less than 0,7 (for DC according to voltage) and 0,5 (for DC according to power).
Despite the enormous DC blur in the maximum position, this method is being widely used in practice due to its simplicity and consistency with the resolution tasks, as well as with other tasks (detection, recognition, optimal reception, etc.). However, this method turns out to be the most effective in the complex with the method of minimum (the equisignal-zone method). In practical work, multi-beam arrays (MBAs) are widely used. These antennas are virtually the analogue of the comb filter, while direction finding is easily determined according to the fact of signal entering into one of the beams (Fig. 2). Such an MBA is usually implemented using a beam-forming Batler matrix [14]. The resolution of SS signals received using the MBA is possible at the angle distances ∆b, that satisfy the Rayleigh criterion. Considering the specifics of mobile system traffic, there is a necessity to create a large number of beams in a narrow sector of directions according to the number of calling SSs. The more the number of beams is, the more complicated and dimensions arrays are needed, but together with this the resolution of the adjacent signals of SSs is improved. However, such an MBA for the problems of wireless communication will hardly be suitable because there is a great possibility that a signal with the arrival angle i b (Fig. 2, а) will enter into zero of DP. This can happen, for instance, under relocation of an SS. The MBA with the controlled DC is more constructive, when a user signal is accompanied by a beam from an ML of the DP. Such MBAs include antenna arrays with adaptation to changes in the signal-to-noise ratio (adaptive antenna array -AAA (Fig. 3)). The method of minimum assumes the presence of a DP antenna with zone of DC minimums. Moreover, any antenna of a specific polarization has at least one DP minimum [11]. In addition loop antennas are often used. These antennas are used for practical direction finding; however, they possess several disadvantages, which lead to a shift in DoA estimations due to antenna effect, enormous impact of surrounding objects, small electric length, etc. In practical work, more complex electronically-controlled constructions are more frequently used. Thus, the outputs of two antennas A 1 and A 2 , the load of which is added to the general load, can have outputs of the corresponding sum or differential directional pattern. The electric scheme of such a construction and the corresponding Σ and Δ of diagrams are shown in (Fig. 4). Further, we will pay our attention to the methods of determining the direction of SS signals arrival based on using antenna arrays with the electronically-controllable weight vector (WV) due to which it is possible to form the necessary DC f(b) (2).

Analysis of methods for determining signals direction of arrival from subscriber station with estimation of field distribution vector
When solving the tasks of the STA at the uplink of SS → BS the explicit knowledge about the direction of arrival b is required because this parameter is mandatory as it is included into the access algorithm as an argument. In the backward direction BS → SS, there is a need for knowledge of reversed direction finding because its direct value is impossible to obtain due to the fact that as a rule the subscriber station contains 1, seldom 2 antennas, what prevents creating the DP with high resolution. The estimation of determining direction finding from the SS is difficult due to the constant change of orientation of this station.
Thus, the STA problem of spatial selection and DoA (direction finding) should be solved using BS resources, where these problems are solved on the basis of N-element antenna array. We are interested in distribution of the vector of EMF electric component for the signal at the aperture of the N-element AA.
where H -(N×L) is a matrix; h a e ,a e ,...,a e φ φ φ   = ë û is the N -dimensional vector, which characterizes distribution and spatial structure of the field of the i-th called signal; a ki is a standardized coefficient of gain for the k-th antenna element (AE) in the direction of the i-th signal arrival b i ; ki φ is the phase shift conditioned by the delay of the i-th signal at the output of the k-th AE with regard to the point accepted as a phase center of the AA. According to (2) the desired estimation of the direction finding is T is the operation of transposition.
To begin with, let us consider the specifics of the space-time problem at the uplink, when the SS needs the delivery of a communication service and it transmits a ringing signal (RS). The BS receives the RS of the given SS, selects it among many RSs from other SSs and extraneous interference, which get into the ringing frequency band. BS spatial selection is performed by a filter coherent with the broadband signal structure as well as on the basis of the algorithm of adaptation, in particular, using the procedure for WV estimation of the adaptive antenna array [12].
At the downlink, the STA problem is different. Due to the fact that the SS resources are limited, the problem itself is solved at the transmission side, i.e. on the BS. The solution of such a problem results in determining and estimating the arrival direction of the RS and setting a narrow beam of the BS transmitting antenna in the backward direction. In this case, we perform DP synthesis according to the given direction of the RS arrival from the given SS.
Let us obtain the estimations of the N-dimensional vector (5) in the assumption of neglecting the mutual connection between antenna elements on the basis of reception of L-dimensional signals where ( ) x t S t ,S t ,...,S t   = ë û is the L-dimensional vector of received signals.
is the optimal i-dimensional weight vector, which determines amplitude-phase current distribution at the output of the AE to be estimated.
Using (6) we obtain the estimation of H matrix: where {} E ⋅ is the estimation on the criterion of the ratio maximum between the SS i-th signal and the total of interfering signals of other SSs and noise (maximum SNIR), a minimum mean square error (MMSE), etc.
According to the WV estimation w i , it is possible to select a desired signal and find DoA estimations for the ith signals accepted by the given AA [12].
For the quality criterion of the estimation minimum, which is determined as a mean-square deviation of the received RS from the reference signal for the i-th SS, the optimal value of the WV W opt is found from the solution of the Wiener-Hopf matrix equation: where yy R is the correlation matrix of SS received signals, including the i-th SS as well as noises (t) ν ; is the vector of mutual correlation of the reference x i (t) and input y(t) signals.
For the criterion of estimation maximum determined as the ratio of the i-th signal to the total of interference signals and noises the WV estimation is: where IN R is the correlation matrix of interference signals and noises; V is the control vector corresponding to the signal of the i-th SS; b is the normalizing coefficient.
In the perfect case when choosing the value of (i) opt W there is a suppression of all acting SS ringing signals except for the i-th one and towards the BS of the i-th SS the maximum of DP is set. This maximum can be obtained as a result of created amplitude-phase distribution (APD).
The given DP maximum can be presented as ( ) max F Σ b , which provides cumulative yield of the sum-differential block [15]. To obtain F Δ differential DP, the aperture is divided into two halves and the mutual subtraction of signals of both halves is performed.
The correlation matrix of the received signals R yy is directly connected with the vector weight of the AA W opt by the ratio: where W 0 is the reference vector representing the initial weight vector. During operation this vector changes its orientation in accordance with the change of the signal direction of arrival from the calling subscriber station and it is characterized by (5).
In operation, the AA provides the suppression of all input signals except for the j-th one: It is obvious that the signals suppressed at the j th output include signals from other called stations and interference concentrated at the given part of the frequency band. To solve the STA problem at all ringing SSs we should estimate all the vectors h j , j=1, 2,…, L at the BS. Thus, the AA of the BS will be used in the L-fold manner to provide the access and processing of all the SS signals.
Let us find the estimation ĵ h using the results of the technique given in [16].
The technique for construction of the estimation algorithm ĵ h , j=1, 2,…, L. 1. The correlation matrix R yy can be obtained on the basis of sample y i , j=1, 2,…, k-L.
where Λ i is the eigen vector, corresponding to the eigen value of the matrix R xx . The eigen values λ i must be arranged in the decreasing manner starting with λ max . Moreover, λ min corresponds to the level of thermal noise. 4. Calculate the estimations of the cumulants for the 4th order whitening of the input signals vector .
For such a diagonalization we can use the generalized Givens transformation.
The presented technique is quite difficult to calculate. However, some of these calculations are performed for the benefit of other tasks: ensuring detection, access, anti-jamming, etc. Due to the complexity, this technique does not have practical application. A more constructive approach is the method of equisignal zone, which is widely used in radiolocation.

Method for direction of arrival estimation for SS signal using equisignal reception zone
The problem of direction finding measurement b i using the method of search for an equisignal zone consists of two parallel or sequential procedures: 1. Finding the estimation of WVs w opt , which provides the maximum of reception over the ML of DC f Σ (b) max .
2. Finding the differential DC f ( ) ∆ b by dividing the aperture of the antenna array into two equal parts.
The method of direction finding via an equisignal reception zone is a common technique to determine the angular coordinates of radio sources in radiolocation. The antenna array of the amplitude monopulse direction finder consists of several (two) identical reception links, forming the directional pattern with a small deviation of the main beam from the antenna equisignal direction.
To determine the angular coordinates b i of the radiation source by the amplitude monopulse method we use a comparison of the relative (normalized) Р і powers of signals received simultaneously by two halves (links) of the antenna array, and further the estimation of the radiation source position is reduced to solving a system in the general case of nonlinear equations [13]: where f i (P i ) is the direction pattern over the power of the і-t h receiving link in the antenna system. Within the process of solution the i and j links of AAA must be identical.
It is possible to impose the monopulse direction finding algorithms more stringent requirements on the speed and increased accuracy of the reception of information on the angular coordinates of the SS signals. However, the main limitation of these methods is well-known: the presence of multipath propagation conditions or two or more sources with comparable capacities leads to large errors in the measured angular coordinates [14].
The method of superresolution (SR), where the problem of determining the direction of arrival is concomitant is considered as a separate issue.

Methods of "superresolution" in determining direction of SS signal arrival
It is worth noting that the successful SR solution requires a higher reception quality (signal-to-noise ratio), more accurate coordinates of the antenna elements (AEs) installation and the design of antennas, known data on the signal statistics and interference, etc. In addition, there are some unpleasant outcomes of decisions, when as a result of the SR method implementation there emerge false marks of the detection of non-existent targets, shifts in the estimates of the direction of arrival. However, superresolution problems are perspective in practice [16], and their application should be given a certain priority.
When implementing the superresolution methods in the spatial domain it is possible to determine the position of the spectral components (directions of arrival) of the analyzed (received) signals.
In works of J. Capon et al., the methods for superresolution of signals separated by an arbitrarily small distance have been developed on the basis of nonlinear spectral analysis techniques.
There is a large number of algorithms that implement the method of SR. They include [16]: 1) methods based on determining the positions of the pseudospectrum ( ) D b local maxima on the scanning in space: Capon, Borgiotti-Lagunas, "thermal noise".
These methods assume the maximum likelihood estimate of the spatial correlation matrix: Spatial pseudospectrum ( ) D b for n -element system array is determined as where ( ) is the directing vector for the m-th signal of the n -element AA.
The Capon's method is implemented in the form: The method of "thermal noise" is implemented as: The Borgiotti-Lagunas method is implemented as follows: 2) Methods based on the expansion of the spatial correlation matrix R over its eigenvectors and k V and keigenvalues.
The MUSIC method is implemented as: The method of EV (eigenvector) is implemented as follows: 3) Methods based on the separation of signal and noise spaces.
These methods include ROOT-MUSIC and ESPRIT methods focused on the use of linear equivalent AAs. The ROOT-MUSIC algorithm uses the orthogonality of its own signal and noise vectors, allowing finding the angles of arrival of all the sources in the field of observation within a single computational procedure.
Reverse value of the spatial pseudospectrum is represented as: The formulas for finding the polynomial roots according to its coefficient exist only for N<5. When N is bigger, it is necessary to apply numerical methods.
Let us study the methods for estimation of direction of signals arrival with superresolution on the basis of the simulation in the MATLAB package.
It follows from the analysis of graphs in Fig. 6, a (SNR=35 dB) that the non-adaptive beamforming method (Fourier) does not allow spatial separation of the signalsas seen, 3 out of 4 signals are in the same resolution element. Other methods have higher signal resolution characteristics and allow estimating the arrival angles of the signals with different accuracy. As it can be seen, the most high-precision method is the MME, but there are false peaks in the spectrum.
With a decrease in the SNR ratio up to 20 dB (Fig. 6, b) the resolution of methods significantly worsens, however, among the considered set of methods, the method of multisignal classification (MUSIC) is the smoothest and it provides a maximum of signals arrival powers from the considered directions. a b Fig. 6. Graphs of the radiation power on the angle of arrival of uncorrelated signals in superresolution methods: a -SNR=35 dB; b -SNR=20 dB To analyze the quality of DoA methods for "superresolution" and Rayleigh methods Fig. 7 shows the dependence of the threshold and objective function of the "superresolution" method in the case of 2 signals resolution with a signal-tonoise ratio of a signal h 2 =30 dB [17]. The objective function is a minimum standard deviation of DPA in the directions to the signal sources. As a threshold function we assume the degree of depth of signal source suppression from its level to the background of the acting noise. Obviously, the use of a superresolution method is effective in the case where the threshold function is above the decreasing objective function that involves receiving signals with an SNR of 30 dB. This is consistent with the results of accuracy estimations of methods presented in Fig. 6. The calculations show [18] that when SNR of one of the sources is more than 30 dB, the resolution angular distance of the two signals can be 0,1 ∆b ≅ θ. However, we should bear in mind that the implementation of the "superresolution" method requires rather high energy parameters of the signals that are calculated as 30 dB and above and that is not always possible to achieve in mobile communication lines.
Furthermore, accurate ratios of antenna design are required and there is no guarantee that false artifacts will not appear.
The analysis of the simulation results of the estimation for the spatial resolution of the signals arrival angles from 4 sources of radiation has shown: 1. The classic non-adaptive beamforming method is the easiest in implementation, however, it has a major drawback associated with the restriction of the Rayleigh resolution.
2. The estimation of sources spatial resolution according to Capon's method is superior to the classical method of beamforming, because it uses every available degree of freedom for the concentration of power produced in the direction of interest direction finding (16). Among the non-parametric methods of continuous analysis such as: "thermal noise", Borgiotti-Lagunas methods, the Capon's method has the lowest resolution (Fig. 6, a).
3. The method of multisignal classification (MUSIC) is the limiting case of the Capon's method when signalto-noise ratio tends to infinity, which is consistent with a higher resolution of the MUSIC in comparison to the Capon's algorithm. Analyzing the graphs shown in Fig. 6, a, b, we conclude that the MUSIC method is able to provide a resolution of the incoming signals under sufficient angular diversity of signal sources and at low signal-to-noise ratio values.
4. The simulation results have confirmed the statistical consistency of the proposed methods for determining the signals direction of arrival of the point radiation sources over the space-time sample at the output of the linear equidistant AA of the non-correlated sources.

Discussion of results of generalized analysis for selection of method to determine signals direction of arrival
Knowledge of direction finding, the directions of SS signals arrival is important practical and useful information that allows to improve the quality of space-time access: the SS signal support is provided under moving, the main lobe in the DP of AA is set in the direction of the calling station, the convergence of adaptive algorithms of the array is facilitated, where the direction finding is used as initial conditions.
In most of the studies, selection of the method to determine the signals direction of arrival especially from nearby sources is based on the method of direct inversion to the sample correlation matrix, which assumes the formation of the sample estimate for the field distribution vector in the AA (14). The study presents a corresponding technique to estimate the vector of AA field distribution, which implements an optimal solution and does not depend on the dispersion of the correlation matrix eigenvalues. However, it is shown that this approach could be hard to implement in practice due to the ill-conditioned inversion of the correlation matrix. The more constructive method is the equisignal zone, which is applied in radiolocation.
The estimation method of equisignal zone with the use of artificial neural networks (ANNs), which allows high-quality processing in real time, is highly constructive and accurate. In addition, together with this, the AANs method allows specifying the number of received signals and the direction of their arrival and that requires further research in this area.
In the above presented aspects, the conducted analysis is considered to be a continuation of previous studies. A definite advantage of this study is a synthesis of the methods for determining the signals direction of arrival on the basis of the known approaches and methods for estimating the directional characteristics of antenna arrays. The approach presented in the study will allow further definition and selection of the most effective methods to be implemented in the STA problems in mobile communication systems.

Conclusions
1. There is a large number of methods and algorithms for estimating the signals direction of arrival, but they all result in the methods for estimation of the maximum or minimum of antenna signal reception. In this case, the minimum methods are more accurate in comparison to the maximum methods, what can be justified due to the different steepness of directional patterns in the given areas.
2. The STA problem of spatial selection and estimation of signals direction of arrival is solved using BS resources and the N-element antenna array. Considering the direct functional connection of amplitude-phase distribution of signals over the AA elements with the value of arrival direction, the procedure and method showing the ability to jointly calculate the vector of distribution with the estimate of the weight vector in AA has been presented. However, the computational complexity of this method limits the scope of its application.
3. The methods for determining the SS signal arrival directions using "superresolution" are promising. Under ideal conditions, these methods make it possible to obtain arbitrarily accurate values of the direction finding; in addition, it is assumed that geometrical and electrodynamic parameters of the AA of the resolved signals themselves are accurately known. Obviously, in reality, these data have certain dispersion, and the multipath propagation nature in mobile lines leads to shifted estimates and the appearance of false marks. The computational complexity of the methods and the possibility of unexpected errors in the estimation of the direction find-ing does not allow the implementation of the method for the STA tasks.
4. The algorithm of signal arrival estimation using the sum-differential DP is quicker and easier to implement. Determination of the direction finding is thus carried out by means of the equisignal zone, which allows to obtain a twice higher accuracy in comparison to the Rayleigh resolution. This method seems to be one of the most meaningful in application to the STA problems.