RU2530320C2 - Formation method of virtual signal receiving channels - Google Patents
Formation method of virtual signal receiving channels Download PDFInfo
- Publication number
- RU2530320C2 RU2530320C2 RU2013102799/08A RU2013102799A RU2530320C2 RU 2530320 C2 RU2530320 C2 RU 2530320C2 RU 2013102799/08 A RU2013102799/08 A RU 2013102799/08A RU 2013102799 A RU2013102799 A RU 2013102799A RU 2530320 C2 RU2530320 C2 RU 2530320C2
- Authority
- RU
- Russia
- Prior art keywords
- antenna array
- vector
- auxiliary
- sources
- auxiliary sources
- Prior art date
Links
Landscapes
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
Description
Изобретение относится к области радиотехники и может использоваться при проектировании и эксплуатации комплексов радиопеленгации или систем радиосвязи портативного, мобильного (бортового) и стационарного базирования.The invention relates to the field of radio engineering and can be used in the design and operation of radio direction finding systems or radio communication systems of portable, mobile (airborne) and stationary basing.
Известен способ оценки напряженности электромагнитного или акустического поля по сигналам элементов антенной решетки, расположенных вблизи искажающего поле рассеивателя [1], заключающийся в том, что на каждом элементе антенной решетки записывают сигнал на временном интервале [0,Т], формируют дискретный спектр напряженности поля, выполняя над записанным сигналом процедуру преобразования Фурье, для каждой спектральной компоненты находят вектор В коэффициентов интерполяционной модели поля, удовлетворяющий матрично-векторному уравнению U(N)=QB, левая часть которого представляет собой вектор U(N) напряженностей поля спектральной компоненты элементов антенной решетки, а правая часть представляет собой произведение матрицы антенной решетки Q, элементы которой зависят от частоты и расположения элементов антенной решетки, и искомого вектора B коэффициентов интерполяционной модели поля, определяют значение поля спектральной компоненты в произвольной точке плоскости антенной решетки как скалярное произведение найденного вектора В коэффициентов интерполяционной модели поля и вектора q, зависящего от частоты и положения этой оцениваемой точки.A known method for evaluating the strength of an electromagnetic or acoustic field from the signals of the elements of the antenna array located near the field-distorting diffuser [1], which consists in the fact that on each element of the antenna array a signal is recorded on the time interval [0, T], a discrete spectrum of field strength is formed, performing the Fourier transform procedure on the recorded signal, for each spectral component, find the vector В of the coefficients of the interpolation field model that satisfies the matrix-vector equation U (N) = QB, the left part of which is the vector U (N) of the field strengths of the spectral component of the elements of the antenna array, and the right part is the product of the matrix of the antenna array Q, the elements of which depend on the frequency and location of the elements of the antenna array, and the desired vector B of the coefficients of the field interpolation model, determine the value of the field of the spectral component at an arbitrary point in the plane of the antenna array as the scalar product of the found vector. In the coefficients of the interpolation model of the field and the vector q, which depends on the frequency and the position of the estimated point.
Однако описанный выше способ основан на необходимости решения систем линейных алгебраических уравнений (СЛАУ), в состав которых входят плохообусловленные матрицы, из-за чего он теряет устойчивость функционирования при наличии погрешностей в измерении комплексных амплитуд на элементах физически существующей антенной решетки, вызванных действиями шумов или несимметричностью каналовHowever, the method described above is based on the need to solve systems of linear algebraic equations (SLAEs), which include poorly conditioned matrices, because of which it loses stability when there are errors in the measurement of complex amplitudes on the elements of a physically existing antenna array caused by noise or asymmetry channels
Изобретение направлено на повышение устойчивости методов оценки напряженности электромагнитного или акустического поля (формирования «виртуальных» каналов приема сигналов) при их функционировании в условиях наличия шумов или несимметричности каналов антенной системы.The invention is aimed at improving the stability of methods for assessing the intensity of electromagnetic or acoustic fields (the formation of "virtual" channels for receiving signals) when they operate in the presence of noise or asymmetry of the channels of the antenna system.
Это достигается тем, что на каждом элементе антенной решетки записывают сигнал на временном интервале [0,Т], формируют дискретный спектр напряженности поля, выполняя над записанным временным сигналом процедуру преобразования Фурье, для каждой спектральной компоненты находят вектор В комплексных амплитуд вспомогательных источников поля с помощью квазирешения матрично-векторного уравнения , например, в виде . При этом - вектор напряженностей поля спектральной компоненты элементов антенной решетки, Q - матрица вспомогательных источников, элементы которой зависят от частоты и взаимного расположения элементов антенной решетки и вспомогательных источников, Н H - знак Эрмитового сопряжения. При этом число вспомогательных источников определяется как, например, как число наиболее значимых собственных чисел автокорреляционной матрицы принятых антенной решеткой сигналов, т.е. наибольших и отличающихся от остальных по величине не менее, чем на порядок. Элементы матрицы Q определяются как Qn,m=exp(i·k0·rn,m)/rn,m, где k0 - волновое число свободного пространства, - расстояние от n-го элемента антенной решетки m-го до вспомогательного источника.This is achieved by the fact that on each element of the antenna array a signal is recorded on the time interval [0, T], a discrete field strength spectrum is formed, performing the Fourier transform procedure on the recorded time signal, for each spectral component the vector B of the complex amplitudes of auxiliary field sources is found using quasisolutions of a matrix-vector equation , for example, in the form . Wherein is the vector of the field strengths of the spectral component of the elements of the antenna array, Q is the matrix of auxiliary sources, the elements of which depend on the frequency and relative position of the elements of the antenna array and auxiliary sources, and H H is the sign of the Hermite pairing. The number of auxiliary sources is determined, for example, as the number of the most significant eigenvalues of the autocorrelation matrix of signals received by the antenna array, i.e. the largest and different from the rest in size no less than an order of magnitude. Elements of the matrix Q are defined as Q n, m = exp (i · k 0 · r n, m ) / r n, m , where k 0 is the wave number of free space, is the distance from the nth element of the antenna array of the mth to the auxiliary source.
Формируют «виртуальный» канал приема сигналов (определяют значение поля спектральной компоненты в произвольной точке плоскости антенной решетки) как скалярное произведение вектора комплексных амплитуд вспомогательных источников и вектора , элементы которого зависят от частоты, положения формируемого «виртуального» канала приема сигналов в пространстве и числа вспомогательных источников. При этом произвольный m-й элемент вектора определяют как gm=exp(i·k0·rm)/rm, где k0 - волновое число свободного пространства, - расстояние от точки, в которой формируется «виртуальный» канал приема сигналов (определяется значение поля спектральной компоненты) до m-го вспомогательного источника.Form a "virtual" channel for receiving signals (determine the value of the field of the spectral component at an arbitrary point on the plane of the antenna array) as the scalar product of the vector complex amplitudes of auxiliary sources and vectors , the elements of which depend on the frequency, position of the generated "virtual" channel for receiving signals in space and the number of auxiliary sources. Moreover, an arbitrary mth element of the vector defined as g m = exp (i · k 0 · r m ) / r m , where k 0 is the wave number of free space, - the distance from the point at which the "virtual" channel for receiving signals is formed (the field value of the spectral component is determined) to the m-th auxiliary source.
Значение напряженности электромагнитного или акустического поля в точке плоскости антенной решетки, получаемое в конечном счете, трактуется как полезный сигнал с некоторого «виртуального» антенного элемента, размещенного в данной точке пространства. Таким образом, можно говорить о формировании «виртуальных» каналов приема сигналов, информация с которых может быть использована для повышения отношения сигнал / шум, уточнения координат источников радиоизлучения и т.д.The value of the intensity of the electromagnetic or acoustic field at a point on the plane of the antenna array, which is ultimately obtained, is interpreted as a useful signal from some “virtual” antenna element located at a given point in space. Thus, we can talk about the formation of “virtual” channels for receiving signals, information from which can be used to increase the signal-to-noise ratio, to clarify the coordinates of radio emission sources, etc.
Отличительной особенностью данного метода является отсутствие необходимости решения системы линейных алгебраических уравнений (СЛАУ) относительно неизвестных комплексных амплитуд вспомогательных источников, в состав которых входят плохообусловленные матрицы. В данном случае вместо решения СЛАУ используется процедура квазирешения. Кроме того, число используемых вспомогательных источников определяется как величина наиболее значимых собственных чисел автокорреляционной матрицы входных сигналов, т.е. наибольших и отличающихся от остальных по величине не менее, чем на порядок.A distinctive feature of this method is the absence of the need to solve a system of linear algebraic equations (SLAE) with respect to unknown complex amplitudes of auxiliary sources, which include poorly conditioned matrices. In this case, instead of solving the SLAE, the quasi-solution procedure is used. In addition, the number of auxiliary sources used is determined as the value of the most significant eigenvalues of the autocorrelation matrix of input signals, i.e. the largest and different from the rest in size no less than an order of magnitude.
На чертеже изображена блок-схема предлагаемого устройства для реализации способа. Устройство содержит N элементов антенной решетки, каждый из которых соединен с соответствующим приемником сигнала 1.1÷1.N. Выход каждого приемника сигнала 1.1÷1.N соединен со входом соответствующего блока формирования спектра 2.1÷2.N. Спектральные компоненты 1÷L с выходов блоков формирования спектров 2.1÷2.N поступают на соответствующие входы 1÷N блоков обработки сигнала спектральной компоненты 3.1÷3.L. Каждый блок обработки сигнала спектральной компоненты 3.1÷3.L содержит последовательно соединенные узел нахождения амплитуд вспомогательных источников 4 и узел нахождения оценки напряженности поля в точке плоскости антенной решетки 5.The drawing shows a block diagram of the proposed device for implementing the method. The device contains N elements of the antenna array, each of which is connected to the corresponding signal receiver 1.1 ÷ 1.N. The output of each signal receiver 1.1 ÷ 1.N is connected to the input of the corresponding spectrum forming unit 2.1 ÷ 2.N. The spectral components 1 ÷ L from the outputs of the spectral forming units 2.1 ÷ 2.N are supplied to the corresponding inputs 1 ÷ N of the signal processing units of the spectral component 3.1 ÷ 3.L. Each signal processing unit of the spectral component 3.1 ÷ 3.L contains a series-connected node for finding the amplitudes of auxiliary sources 4 and a node for estimating the field strength at a point on the plane of the antenna array 5.
В состав устройства также входят узел определения числа вспомогательных источников 6, блок формирования матриц вспомогательных источников 7, а также блок формирования векторов «виртуальных» каналов приема сигналов 8. Выход узла определения числа вспомогательных источников 6 соединен со входом блока формирования матриц вспомогательных источников 7. Выходы блока формирования матриц вспомогательных источников 7 соединены, соответственно, со входами узлов нахождения амплитуд вспомогательных источников 4. Выходы блока формирования векторов «виртуальных» каналов приема сигналов 8 соединены, соответственно, со входами узлов нахождения оценки напряженности поля в точках плоскости антенной решетки. Выходы узлов нахождения оценки напряженности поля в точках плоскости антенной решетки 5 являются выходами устройства.The device also includes a node for determining the number of auxiliary sources 6, a block for generating matrices of auxiliary sources 7, and also a block for generating vectors of “virtual” channels for receiving signals 8. The output of the node for determining the number of auxiliary sources 6 is connected to the input of the block for generating matrices for auxiliary sources 7. Outputs block forming the matrix of auxiliary sources 7 are connected, respectively, with the inputs of the nodes of finding the amplitudes of the auxiliary sources 4. The outputs of the vector forming unit in "virtual" channels receive signals 8 are connected, respectively, to the inputs of nodes finding a bound field intensity at points plane array antenna. The outputs of the nodes finding estimates of the field strength at the points of the plane of the antenna array 5 are the outputs of the device.
Способ осуществляется следующим образом.The method is as follows.
Сигнал с каждого элемента антенной решетки поступает на вход приемника сигнала 1.1÷1.N, где производятся процедуры фильтрации, переноса на видеочастоту, аналого-цифрового преобразования и т.д. Видеосигналы с выходов приемников сигнала 1.1÷1.N поступают на входы соответствующих блоков формирования спектра 2.1÷2.N, где записываются на временном интервале [0, Т]. Затем в блоках формирования спектра 2.1÷2.N формируют дискретный спектр напряженности поля, выполняя над записанным временным сигналом процедуру дискретного преобразования Фурье длиной L. Спектральные компоненты 1÷L с выходов блоков формирования спектра 2.1÷2.N поступают на соответствующие входы 1÷N блоков обработки сигнала спектральной компоненты 3.1÷3.L.The signal from each element of the antenna array is fed to the input of the signal receiver 1.1 ÷ 1.N, where filtering, transferring to the video frequency, analog-to-digital conversion, etc. are performed. The video signals from the outputs of the signal receivers 1.1 ÷ 1.N are fed to the inputs of the corresponding spectrum forming units 2.1 ÷ 2.N, where they are recorded on the time interval [0, T]. Then, a discrete spectrum of field strength is formed in the blocks of spectrum formation 2.1 ÷ 2.N, performing a discrete Fourier transform of length L on the recorded temporary signal. The spectral components 1 ÷ L from the outputs of the spectrum forming blocks 2.1 ÷ 2.N are sent to the corresponding inputs 1 ÷ N signal processing units of the spectral component 3.1 ÷ 3.L.
В блоке 6, на основе принятых антенной решеткой сигналов, производится определение числа необходимых вспомогательных источников. При этом значение величины используемых вспомогательных источников определяется как число наиболее значимых собственных чисел автокорреляционной матрицы принятых антенной решеткой сигналов, т.е. наибольших и отличающихся от остальных по величине не менее, чем на порядок. Далее, полученное значение используется для формирования матриц вспомогательных источников Q. Каждый элемент матрицы Q определяется как Qn,m=ехр(i·k0·rn,m)/rn,m, где - расстояние от n-го элемента антенной решетки m-го до вспомогательного источника.In block 6, based on the signals received by the antenna array, the number of necessary auxiliary sources is determined. In this case, the value of the value of the auxiliary sources used is determined as the number of the most significant eigenvalues of the autocorrelation matrix of the signals received by the antenna array, i.e. the largest and different from the rest in size no less than an order of magnitude. Further, the obtained value is used to form matrices of auxiliary sources Q. Each element of the matrix Q is defined as Q n, m = exp (i · k 0 · r n, m ) / r n, m , where is the distance from the nth element of the antenna array of the mth to the auxiliary source.
Для каждой спектральной компоненты в узле нахождения амплитуд вспомогательных источников 4 находят вектор данных амплитуд . Вектор находят с помощью процедуры квазирешения матрично-векторного уравнения , в виде . При этом - вектор напряженностей поля спектральной компоненты элементов антенной решетки, значения которого поступают с блоков формирования спектра 2.1÷2.L, Q - матрица вспомогательных источников, элементы которой зависят от частоты и взаимного расположения элементов антенной решетки и вспомогательных источников, H - знак Эрмитового сопряжения.For each spectral component in the node finding the amplitudes of the auxiliary sources 4 find the vector of these amplitudes . Vector using the quasi-solution of the matrix-vector equation , as . Wherein is the vector of the field strengths of the spectral component of the elements of the antenna array, the values of which come from the spectrum forming units 2.1 ÷ 2.L, Q is the matrix of auxiliary sources, the elements of which depend on the frequency and relative position of the elements of the antenna array and auxiliary sources, H is the sign of the Hermite pairing.
Значение вектора , определенное в узле нахождения амплитуд вспомогательных источников поля 4 поступает на первый вход узла нахождения оценки напряженности поля в точке плоскости антенной решетки 5. На второй вход узла оценки напряженности поля в точке, лежащей в плоскости антенной решетки 5, поступает значение вектора «виртуального» канала приема сигналов , которое формируется в блоке 8. Произвольный m-й элемент вектора определяют как gm=ехр(i·k0·rm)/rm, где k0 - волновое число свободного пространства, - расстояние от точки, в которой формируется «виртуальный» канал приема сигналов (определяется значение поля спектральной компоненты) до m-го вспомогательного источника. В узле оценки напряженности поля в точке плоскости антенной решетки значение поля спектральной компоненты определяется как скалярное произведение вектора комплексных амплитуд вспомогательных источников и вектора оцениваемой точки: .Vector value defined in the node for finding the amplitudes of auxiliary sources of field 4 is fed to the first input of the node for estimating the field strength at a point on the plane of the antenna array 5. At the second input of the node for estimating the field strength at a point lying in the plane of the antenna array 5, the value of the vector of the "virtual" channel receiving signals , which is formed in block 8. The arbitrary mth element of the vector defined as g m = exp (i · k 0 · r m ) / r m , where k 0 is the wave number of free space, - the distance from the point at which the "virtual" channel for receiving signals is formed (the field value of the spectral component is determined) to the m-th auxiliary source. In the node for estimating the field strength at a point on the plane of the antenna array, the value of the field of the spectral component is defined as the scalar product of the vector of complex amplitudes of auxiliary sources and vectors estimated point: .
Источник информацииThe source of information
1. RU, патент №2405165 C2, кл., G01S 3/00, 27.11.2010.1. RU, patent No. 2405165 C2, class., G01S 3/00, 11.27.2010.
Claims (5)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2013102799/08A RU2530320C2 (en) | 2013-01-22 | 2013-01-22 | Formation method of virtual signal receiving channels |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2013102799/08A RU2530320C2 (en) | 2013-01-22 | 2013-01-22 | Formation method of virtual signal receiving channels |
Publications (2)
Publication Number | Publication Date |
---|---|
RU2013102799A RU2013102799A (en) | 2014-07-27 |
RU2530320C2 true RU2530320C2 (en) | 2014-10-10 |
Family
ID=51264664
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2013102799/08A RU2530320C2 (en) | 2013-01-22 | 2013-01-22 | Formation method of virtual signal receiving channels |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2530320C2 (en) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2140238A (en) * | 1983-05-11 | 1984-11-21 | Racal Communications Equip | Direction finding |
DE4128191A1 (en) * | 1991-08-24 | 1993-02-25 | Telefunken Systemtechnik | Short-baseline interferometer for direction-finding by correlation |
US6239747B1 (en) * | 1999-03-11 | 2001-05-29 | Lucent Technologies Inc. | Antenna system and method for direction finding |
US6249262B1 (en) * | 1999-11-03 | 2001-06-19 | Motorola, Inc. | Switchable antenna for radio communication devices |
RU2184980C1 (en) * | 2001-11-12 | 2002-07-10 | Рембовский Анатолий Маркович | Procedure measuring intensity of electromagnetic field of radio signals and device for its implementation |
RU2201599C1 (en) * | 2002-06-27 | 2003-03-27 | Рембовский Анатолий Маркович | Method of direction finding of radio signals and direction finder for its realization |
RU2267134C2 (en) * | 2003-12-03 | 2005-12-27 | Федеральное государственное унитарное предприятие "Ростовский-на-Дону научно-исследовательский институт радиосвязи" | Mode of direction finding of radio signals and a direction finder of radio signals |
US7652626B1 (en) * | 2007-12-05 | 2010-01-26 | The United States Of America As Represented By The Secretary Of The Army | Co-linear antenna for direction finding |
RU2405165C2 (en) * | 2008-12-09 | 2010-11-27 | Закрытое акционерное общество "ИРКОС" | Evaluation method of intensity of electromagnetic or acoustic field as per signals of elements of antenna array, which are located near field distorting diffuser (versions) |
-
2013
- 2013-01-22 RU RU2013102799/08A patent/RU2530320C2/en not_active IP Right Cessation
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2140238A (en) * | 1983-05-11 | 1984-11-21 | Racal Communications Equip | Direction finding |
DE4128191A1 (en) * | 1991-08-24 | 1993-02-25 | Telefunken Systemtechnik | Short-baseline interferometer for direction-finding by correlation |
US6239747B1 (en) * | 1999-03-11 | 2001-05-29 | Lucent Technologies Inc. | Antenna system and method for direction finding |
US6249262B1 (en) * | 1999-11-03 | 2001-06-19 | Motorola, Inc. | Switchable antenna for radio communication devices |
RU2184980C1 (en) * | 2001-11-12 | 2002-07-10 | Рембовский Анатолий Маркович | Procedure measuring intensity of electromagnetic field of radio signals and device for its implementation |
RU2201599C1 (en) * | 2002-06-27 | 2003-03-27 | Рембовский Анатолий Маркович | Method of direction finding of radio signals and direction finder for its realization |
RU2267134C2 (en) * | 2003-12-03 | 2005-12-27 | Федеральное государственное унитарное предприятие "Ростовский-на-Дону научно-исследовательский институт радиосвязи" | Mode of direction finding of radio signals and a direction finder of radio signals |
US7652626B1 (en) * | 2007-12-05 | 2010-01-26 | The United States Of America As Represented By The Secretary Of The Army | Co-linear antenna for direction finding |
RU2405165C2 (en) * | 2008-12-09 | 2010-11-27 | Закрытое акционерное общество "ИРКОС" | Evaluation method of intensity of electromagnetic or acoustic field as per signals of elements of antenna array, which are located near field distorting diffuser (versions) |
Also Published As
Publication number | Publication date |
---|---|
RU2013102799A (en) | 2014-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10955542B2 (en) | Radar apparatus and direction-of-arrival estimation device | |
Northardt et al. | Spatial compressive sensing for direction-of-arrival estimation with bias mitigation via expected likelihood | |
JP2019132850A (en) | program | |
RU2593276C1 (en) | Method of selecting moving targets | |
JP2012163403A (en) | Electronic scanning type radar system, and method and program for estimating direction of received wave | |
US20160173241A1 (en) | Detecting and processing weak signals using an array of antennas | |
CN104007421B (en) | Rowland C passive radar TOA method of estimation based on total variation and compressed sensing | |
KR102183439B1 (en) | Method and apparatus for estimating direction of arrival using combined beamspace music and tma | |
JP2016225981A (en) | Estimation device, estimation method of multipath time delay and receiver | |
RU2016129647A (en) | METHOD FOR DETERMINING COORDINATES OF SPACE VEHICLE BY SIGNALS OF NAVIGATION SATELLITES AND DEVICE FOR DETERMINING COORDINATES OF SPACE EQUIPMENT BY SIGNALS OF NAVIGATION SATELLITES | |
RU2718737C1 (en) | Method of determining coordinates of radio-frequency radiation sources | |
RU2524401C1 (en) | Method for detection and spatial localisation of mobile objects | |
RU2549207C2 (en) | Device for detecting hydroacoustic noise signals based on quadrature receiver | |
CA2948539C (en) | Calibration device | |
RU2491569C2 (en) | Method of direction finding with increased resolution ability | |
RU2530320C2 (en) | Formation method of virtual signal receiving channels | |
JP2010223895A (en) | Radar system | |
JP2019144077A (en) | Azimuth estimation method and device | |
Wang et al. | An approach for sound sources localization and characterization using array of microphones | |
RU2631944C1 (en) | Device for determining directions to radio-frequency sources | |
KR101589269B1 (en) | Method and apparatus for estimating location | |
Harley et al. | Matched field processing localization with random sensor topologies | |
JP5599200B2 (en) | Target detection / tracking device | |
RU2572792C1 (en) | Method of integrating noisy sea object detection systems | |
RU2555194C1 (en) | Processing of hydroacoustic signal from noise-emitting object |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20160123 |