RU2010111583A

RU2010111583A - Direction Finding Method with Enhanced Resolution Ability

Info

Publication number: RU2010111583A
Application number: RU2010111583/07A
Authority: RU
Inventors: Николай Григорьевич Пархоменко; Николай Макарович Иванов; Валерий Николаевич Шевченко
Original assignee: Российская Федерация, от имени которой выступает Министерство обороны Российской Федерации
Priority date: 2010-03-25
Filing date: 2010-03-25
Publication date: 2013-01-20
Also published as: RU2491569C2

Abstract

1. Способ пеленгования с повышенной разрешающей способностью, заключающийся в том, что принимают многолучевой сигнал источника акустического или электромагнитного излучения антенной решеткой из N элементов, расположенных равномерно по окружности, формируют ансамбль сигналов, зависящих от времени и номера антенного элемента, синхронно преобразуют ансамбль принятых сигналов в цифровые сигналы, преобразуют цифровые сигналы в сигнал комплексной пространственной корреляционной матрицы, описывающий амплитуды и фазы взаимных сигналов, принятых элементами решетки, преобразуют сигнал пространственной корреляционной матрицы в сигналы собственных значений и собственных векторов, сравнивают сигналы собственных значений с порогом и при не превышении порога сигнал соответствующего собственного вектора идентифицируют как сигнал собственного вектора, принадлежащий шумовому подпространству, отличающийся тем, что преобразуют сигналы собственных векторов шумового подпространства в сигналы дискретных спектров собственных векторов, которые запоминают, формируют сигналы N-направленной комплексной фазирующей функции, зависящие от заданной частоты приема и описывающие возможные направления прихода сигнала от каждого потенциального источника, преобразуют сигналы N-направленной фазирующей функции в сигналы дискретных спектров фазирующей функции, которые запоминают, используя сигналы дискретных спектров собственных векторов и сигналы дискретных спектров фазирующей функции, формируют сигнал углового спектра, по максимумам сигнала углового спектра определяют азимут αи угол места βкаждого луча принятого мно1. The method of direction finding with increased resolution, which consists in the fact that they accept a multipath signal source of acoustic or electromagnetic radiation by an antenna array of N elements arranged uniformly around the circumference, form an ensemble of signals depending on the time and number of the antenna element, synchronously transform the ensemble of received signals into digital signals, convert digital signals into a signal of a complex spatial correlation matrix that describes the amplitudes and phases of the mutual signals, etc. taken by lattice elements, they transform the signal of the spatial correlation matrix into eigenvalues and eigenvectors, compare the eigenvalues with a threshold and, if the threshold is not exceeded, the signal of the corresponding eigenvector is identified as an eigenvector belonging to the noise subspace, characterized in that it converts eigenvectors the noise subspace into the signals of the discrete spectra of the eigenvectors that store, form the N-signals of a complex phasing function, depending on the given reception frequency and describing the possible directions of signal arrival from each potential source, transform the signals of the N-directional phasing function into the signals of the discrete spectra of the phasing function, which are stored using the signals of the discrete spectra of the eigenvectors and the signals of the discrete spectra of the phasing function , form the signal of the angular spectrum, from the maxima of the signal of the angular spectrum, determine the azimuth α and the elevation angle β of each ray received

Claims

1. The method of direction finding with increased resolution, which consists in the fact that they accept a multipath signal source of acoustic or electromagnetic radiation by an antenna array of N elements arranged uniformly around the circumference, form an ensemble of signals depending on the time and number of the antenna element, synchronously transform the ensemble of received signals into digital signals, convert digital signals into a signal of a complex spatial correlation matrix that describes the amplitudes and phases of the mutual signals, etc. taken by lattice elements, they transform the signal of the spatial correlation matrix into eigenvalues and eigenvectors, compare the eigenvalues with a threshold and, if the threshold is not exceeded, the signal of the corresponding eigenvector is identified as an eigenvector belonging to the noise subspace, characterized in that it converts eigenvectors of the noise subspace into the signals of the discrete spectra of eigenvectors that store, form N-signals of a complex phasing function, depending on the given reception frequency and describing the possible directions of signal arrival from each potential source, transform the signals of the N-directional phasing function into the signals of the discrete spectra of the phasing function, which are stored using the signals of the discrete spectra of the eigenvectors and the signals of the discrete spectra of the phasing function form the signal of the angular spectrum, the azimuth α ₀ and the elevation angle β _{0 of} each ray of the received multipath signal, the obtained two-dimensional bearings (α ₀ , β ₀ ) of the selected rays are displayed on a cartographic background.

2. The method according to claim 1, characterized in that the formation of signals of an N-directional complex phasing function is carried out according to

where

- the shift number of the phasing function in azimuth α; Q = M / N; M is the number of grid nodes in azimuth α; δ = Δ / Q is the azimuthal step α; Δ = 2π / N; β _ℓ - given grid nodes by elevation angle β;

- node number by angle β; L is the number of nodes; d (qδ + nΔ, β _ℓ ) is the complex radiation pattern of the nth element; r is the radius of the antenna array; f is the given reception frequency.

3. The method according to claim 1, characterized in that the formation of the signal of the angular spectrum is carried out according to the formula

, where F ^-1 {...} is the inverse DFT operator, α _m = mδ,

- grid node number in azimuth α, m = q + nQ, n = 0, ..., N-1,

is the current number of eigenvector signals, J is the number of eigenvector signals,

- signals of discrete spectra of eigenvectors,

- signals of discrete spectra of the phasing function.