Способ оценки угловых параметров ионосферных сигналов



Способ оценки угловых параметров ионосферных сигналов
Способ оценки угловых параметров ионосферных сигналов
Способ оценки угловых параметров ионосферных сигналов
Способ оценки угловых параметров ионосферных сигналов
Способ оценки угловых параметров ионосферных сигналов
Способ оценки угловых параметров ионосферных сигналов
Способ оценки угловых параметров ионосферных сигналов
Способ оценки угловых параметров ионосферных сигналов
Способ оценки угловых параметров ионосферных сигналов
Способ оценки угловых параметров ионосферных сигналов
Способ оценки угловых параметров ионосферных сигналов
Способ оценки угловых параметров ионосферных сигналов

 


Владельцы патента RU 2518013:

Федеральное государственное автономное образовательное учреждение высшего профессионального образования "Балтийский федеральный университет имени Иммануила Канта" (БФУ им. И. Канта) (RU)

Изобретение относится к радиотехнике, а именно к области пеленгации. Достигаемый технический результат - расширение возможностей пеленгации, сокращение времени расчета угловых параметров многолучевого ионосферного сигнала. Технический результат достигается тем, что круговую антенную систему, расположенную на поверхности земли, дополняют линейной системой вибраторов, расположенных вдоль вертикали к поверхности земли. С помощью сформированной таким образом антенной системы (трехмерная антенная система), многоканального приемника, многоканального аналого-цифрового преобразователя (АЦП) и временного преобразования Фурье формируют пространственно-временной массив комплексных данных E n , m , отображающий значения напряженности поля в n точках трехмерного пространства (n - номер вибратора) и в m-е моменты времени, с интервалами 1-2 секунды (индекс m определяет номер временного среза данных на n вибраторах, m=1÷M+1). Количество временных срезов данных берется на единицу больше, чем количество лучей М. Затем осуществляют соответствующую математическую обработку, фильтруют однолучевые поля из совокупности полей ионосферного сигнала, формируют для каждого выделенного поля диаграммы направленности, сканируют диаграммой направленности в диапазоне оценочных максимумов и углов места и оценивают азимуты, углы места и амплитуды по максимуму диаграммы направленности для М лучей ионосферного сигнала. 6 ил.

 

Изобретение относится к радиотехнике, а именно к области пеленгации, и может быть использовано для измерения азимутов и углов места ионосферных сигналов в условиях приема нескольких интерферирующих лучей, в том числе при малых углах места. При приеме ионосферных сигналов в точку приема, как правило, приходят несколько лучей (1÷4), отраженных от разных слоев ионосферы, с азимутами, находящимися в области главного лепестка диаграммы направленности (Δα~5÷10 градусов). Вследствие интерференции суммарный главный лепесток диаграммы направленности антенной системы, состоящий из суммы диаграмм направленности отдельных лучей, существенно меняется, часто создавая ложные азимуты (пеленги), отличающиеся от истинных азимутов (пеленгов) на десятки градусов. В результате достоверность оценок азимутов, а также углов места ионосферных сигналов оказывается низкой. Кроме того, при пеленгации ионосферных сигналов существует проблема малых углов места. При приеме ионосферных сигналов на расстояниях более 2000 км угол места β оказывается в области значений от 0° до 15÷20°. В этих условиях пространственная база для определения угла места L=Rsin(β) становится малой и значительно возрастает дисперсия измеряемых углов места. Задачей пеленгации ионосферных сигналов является определение координат источника излучения. В связи с этим в способе пеленгации необходимо решение задачи выделения лучевой структуры сигнала и определение, по измеренной лучевой структуре, координат излучателя. При этом важна информация как об азимутах, так и об углах места многолучевого ионосферного сигнала.

Известны фазовые способы измерения азимута и угла места (пеленгации), осуществляемые путем измерений разности фаз между вибраторами антенной системы и оценки по этим измерениям азимутов и углов места (Патент RU №2263327, опубликован 27.10.2005 г.; Патент RU №2365931, опубликован 27.08.2009 г.; Патент RU №2429500, опубликован 20.09.2011; Патент RU №2450283, опубликован 10.05.2012 г.; Заявка на изобретение RU №2010143935, опубликована 10.05.2012 г.; Заявка RU №2003108306 от 25.03.2003 г., G01S 3/14, опубликована 10.10.2004 г.). Недостатком вышеуказанных способов является то, что используется только фазовая информации. Однозначное определение фазы волны на вибраторах возможно на интервале 0÷360°. Это требует малого, по сравнению с длиной волны, пространственного разнесения вибраторов | R ¯ n + 1 R ¯ n | < λ 2 . В условиях приема двух или нескольких близких по азимуту лучей фазовые способы оценки азимута и угла места являются неустойчивыми. При разности фаз между лучами ~180° они дают отклонения азимута на десятки градусов (ложные пеленги), что значительно снижает достоверность оценок азимутов и углов места ионосферных сигналов. Разработанные к настоящему времени фазовые методы оценки азимутов и углов места (пеленгации) предназначены для оценок этих параметров для одного луча на заданной частоте. Антенная система, расположенная в плоскости земной поверхности, не может решить проблему малых углов в оценке азимутов и углов места ионосферных сигналов.

Известны способы пеленгации по максимуму диаграммы направленности антенной системы (Патент RU №2144200, опубликован 10.01.2000 г.; Патент RU №2258241, опубликован 10.08.2005 г.; Патент RU №2419805, опубликован 27.05.2011 г.; Патент RU №2201599, опубликован 27.03.2003 г.; Патент RU №2004100714, опубликован 20.06.2005 г.). В этом случае используется временное преобразование Фурье для частотного выделения сигнала от отдельных вибраторов антенной системы и различные формы пространственной обработки данных (формирование диаграммы направленности). Наилучшим образом диаграмма направленности формируется при использовании пространственного преобразования Фурье. В результате создается двумерный комплексный угловой спектр (диаграмма направленности антенной системы). При использовании пространственного преобразования Фурье диаграмма направленности формируется по выражению (комплексный вид):

U ( α ' , β ' ) = 1 N n = 1 N E n e i K ¯ ( α ' , β ' ) R ¯ n ( 1 )

где:

E n - комплексная амплитуда в n-точке пространства с радиус-вектором R ¯ n ,

K ¯ ( α ' , β ' ) - оценочный волновой вектор ионосферного сигнала,

, β ' ' - оценочные значения азимута и угла места ионосферного сигнала.

Квадрат модуля нормированной диаграммы направленности (расчетная форма) определяется выражением

U 2 ( a ' , β ' ) = E n cos ( ψ n K ¯ ( α ' , β ' ) R ¯ n ) ¯ 2 + E n sin ( ψ n K ¯ ( α ' , β ' ) R ¯ n ) ¯ 2 E ¯ n 2 ,          (2)

где черта сверху означает суммирование по индексу «n».

Недостатком вышеуказанных способов, связанных с оценкой углов места по максимуму диаграммы направленности антенной системы, является неустойчивость решения при наличии многолучевой интерференции при разности фаз между лучами ~180°. При приеме двух или нескольких интерферирующих лучей ионосферного сигнала (различие в азимутах меньше ширины диаграммы направленности) суммарная диаграмма направленности за счет изменения разности фаз между лучами (в области ~180 градусов) в значительной степени подавляется. Боковые или задние лепестки суммарной диаграммы направленности в этих условиях становятся больше по величине, чем основной лепесток диаграммы направленности. В результате азимут (пеленг) в амплитудном способе пеленгации определяется по максимальному боковому или заднему лепестку диаграммы направленности антенной системы. Отклонения азимута (ложные пеленги) достигают десятков градусов. В результате достоверность оценок азимутов и углов места ионосферных сигналов оказывается низкой. Кроме того, невозможно обеспечить работоспособность пеленгатора в широком частотном диапазоне 2÷30 МГц. Требования точности измерения пеленга в низкочастотной части диапазона (3÷4 МГц) определяют базу антенной системы 2R≥150 м. Однако при количестве вибраторов ~16 в высокочастотной части диапазона (20÷30 МГц) диаграмма направленности не формируется в связи с большим пространственным разнесением вибраторов по сравнению с длиной волны, что ограничивает частотный диапазон сверху. Расположение вибраторов антенной системы на поверхности земли (на плоскости) не может решить проблему малых углов при пеленгации ионосферных сигналов.

Известны способы определения азимута и угла места (пеленгации), антенная система в которых состоит из двух взаимно перпендикулярных линейных эквидистантных антенных решеток (Патент RU №2192651, опубликован 05.10.2000 г.). Способ включает в себя прием сигнала с помощью антенной системы, многоканального приемника, преобразование аналоговых сигналов в каждом канале в цифровую форму, использование двумерного углового преобразования Фурье, которое создает диаграмму направленности антенной системы. Максимум диаграммы направленности позволяет оценить азимут (пеленг) и угол места. Недостатком этого способа пеленгации является зависимость точности измерения пеленга от взаимной ориентации антенной системы и пеленга. При пеленгации ионосферных сигналов в условиях приема двух или нескольких интерферирующих лучей могут возникать ложные пеленги, что значительно снижает достоверность оценок азимутов и углов места ионосферных сигналов. Расположение вибраторов антенной системы на поверхности земли (в плоскости) не может решить проблему малых углов при пеленгации ионосферных сигналов данным способом пеленгации.

Известны способы оценки азимутов и углов места (пеленгации), антенная система в которых состоит из ограниченного количества вибраторов (3÷5 вибратора) (Патент RU №2262119, опубликован 10.10.2005 г.; Патент RU №2253877, опубликован 10.06.2005 г.). Недостатками этого способа пеленгации является малая помехоустойчивость, вследствие отсутствия статистической обработки данных, наличие ложных пеленгов при приеме двух или нескольких близких по азимуту лучей ионосферного сигнала, невозможность решения проблемы малых углов места.

Наиболее близким (прототип) к предлагаемому способу оценки угловых параметров при приеме многолучевых ионосферных сигналов является «Способ пеленгации с учетом корреляционной взаимосвязи между лучами», патент RU №2305294, МПК G01S 3/16, опубликован 27.08.2007 г.

Согласно данному способу пеленгации (прототип) последовательность действий следующая. С помощью антенной системы, состоящей из N-вибраторов, расположенных равномерно по окружности радиуса R (пространственная база сигнала) 1) принимают ионосферные сигналы, 2) преобразуют их по частоте, 3) усиливают посредством многоканального приемника. Аналоговые сигналы на выходах многоканального приемника в каждом канале (от каждого вибратора) 4) преобразуют в цифровую форму посредством многоканального аналого-цифрового преобразователя (АЦП). 5) Определяют амплитуды En и фазы ψn принятого ионосферного сигнала в каждом канале (от каждого вибратора) с помощью временного преобразования Фурье. 6) Формируют суммарную четырехмерную (для двух лучей) диаграмму направленности с учетом коэффициента корреляции A 1 между лучами по выражениям

U ( α 1 ' , β 1 ' , α 2 ' , β 2 ' ) = U 1 E n * exp ( i k 1 ¯ ( α 1 ' , β 1 ' ) R ¯ n ) ¯ + U 2 E n * exp ( i k 2 ¯ ( α 2 ' , β 2 ' ) R ¯ n ) ¯ E n E n * ¯          (3)

где

U 1 = E n exp ( i k 1 ¯ ( α 1 ' , β 1 ' ) R ¯ n ) ¯ + E n exp ( i k 2 ¯ ( α 2 ' , β 2 ' ) R ¯ n ) ¯ A 1 1 A 1 A 1 * ¯  

U 2 = E n exp ( i k 2 ¯ ( α 2 ' , β 2 ' ) R ¯ n ) ¯ + E n exp ( i k 1 ¯ ( α 1 ' , β 1 ' ) R ¯ n ) ¯ A 1 * 1 A 1 A 1 * ¯  

A 1 = exp ( i ( k ¯ 2 ( α 2 ' , β 2 ' ) k ¯ 1 ( α 1 ' , β 1 ' ) ) R ¯ n ¯

E n - комплексная амплитуда сигнала, измеренная на n вибраторе,

R ¯ n - радиус-вектор, определяющий местоположение n вибратора,

k ¯ 1 ( α 1 ' , β 1 ' ) , k ¯ 2 ( α 2 ' , β 2 ' ) - волновые векторы первого и второго лучей ионосферного сигнала, зависящие от оценочных азимутов α 1 ' , α 2 ' и углов места β 1 ' , β 2 ' ,

A 1 - коэффициент корреляции между лучами.

Черта сверху означает суммирование по индексу n.

Знак * означает комплексное сопряжение,

Знак ∧ означает комплексную величину.

7) Производят сканирование диаграммой направленности в четырехмерном пространстве за счет изменения оценочных азимутов и углов места α 1 ' , α 2 ' и β 1 ' , β 2 ' в указанных угловых диапазонах с определенным шагом и запоминая значения параметров α 1 ' , α 2 ' , β 1 ' , β 2 ' , U 1 , U 2 , U ( α 1 ' , β 1 ' , α 2 ' , β 2 ' ) . Азимут и угол места каждого луча определяются по максимуму четырехмерной диаграммы направленности.

Недостатком способа пеленгации с учетом корреляционной взаимосвязи между лучами (прототипа) является его ограниченность, связанная с приемом только одного или двух лучей ионосферного сигнала, в то время как в ионосферном сигнале количество лучей может достигать четырех и более. Вторым недостатком этого способа является большая трудоемкость, связанная с необходимостью перебора (сканирования диаграммой направленности) двух азимутов и двух углов места ( α 1 ' , β 1 ' и α 2 ' , β 2 ' ) с шагом по азимуту Δα≤0.5° в диапазоне 0÷360 градусов и с шагом по углу места Δβ≤0.5° в диапазоне 0÷90 градусов в четырехмерном пространстве азимутов и углов места. При увеличении количества лучей в этом способе пеленгации время, требуемое для расчета параметров, увеличивается экспоненциально с увеличением количества лучей. Если Т - необходимое время расчета параметров одного луча в одной точке области определения и М - количество точек, то для N лучей требуется время ТТ=Т*M2N. Третьим недостатком является наличие проблемы малых углов места при данном способе пеленгации ионосферных сигналов, обусловленное плоской (двумерной) антенной системой.

Блок-схема технической реализации данного способа пеленгации (прототипа) представлена на фиг.1. Согласно способу пеленгации с учетом корреляционной взаимосвязи между лучами последовательность действий следующая.

1. Принимают сигналы с помощью круговой антенной системы,

состоящей из N вибраторов ( N N 1 + 1 ( 1 | A 1 | 2 ) , N 1 = 8 ) , расположенных на поверхности земли равномерно по окружности радиуса R ( R 2 π 3 | k ¯ ( α 01 , β 01 ) k ¯ 2 ( α 02 , β 02 ) | ) (блок 1). Количество вибраторов определяется количеством параметров ионосферного сигнала и заданным значением отношения сигнал/шум. Радиус антенной системы определяется размером пространственных интерференционных максимумов на поверхности земли.

2. С помощью многоканального приемника преобразуют сигналы от каждой антенны (в каждом канале) по частоте, усиливают и фильтруют (блок 2).

3. Преобразуют аналоговые сигналы на выходах многоканального приемника от каждого вибратора в цифровую форму с помощью многоканального аналого-цифрового преобразователя (АЦП) (блок 3).

4. Определяют амплитуды En и фазы ψn сигналов от каждого

вибратора (в каждом канале), например, с помощью временного преобразования Фурье (блок 4).

5.Задают точку в пространстве оценочных азимутов α 1 ' , α 2 ' и углов места β 1 ' , β 2 ' и определяют согласно (3) четырехмерную диаграмму направленности U ( α 1 ' , β 1 ' , α 2 ' , β 2 ' ) в заданной точке по выборке данных E n и запоминают ее значения и значения переменных U 1 , U 2 (Блок 5).

6. Повторяют действия (5) с другими значениями параметров α 1 ' , α 2 ' и β 1 ' , β 2 ' (из области их определения) с шагом по азимуту Δα≤0.5° в диапазоне 0÷360 градусов и с шагом по углу места Δβ≤0.5° в диапазоне 0÷90 градусов и создают поверхность в четырехмерном пространстве азимутов и углов места (четырехмерную диаграмму направленности) (блок 6).

7. Определяют максимальное значение четырехмерной диаграммы направленности и фиксируют параметры α 1 ' = α 01 , β 1 ' = β 01 , α 2 ' = α 02 , β 2 ' = β 02 , а также U 1 = U 01 , U 2 = U 02 (блок 7).

8. Определяют ранг поля (однолучевое или двулучевое) по условию U01/U02≥3 и оставляют решения, удовлетворяющие этому условию (блок 8).

Целью предлагаемого изобретения «Способ оценки угловых параметров ионосферных сигналов» (технический результат) является:

- расширение возможностей «Способа пеленгации с учетом корреляционной взаимосвязи» на случай приема М лучей (М=1÷4 и более),

- решение проблемы малых углов места при пеленгации ионосферных сигналов,

- сокращение времени расчета угловых параметров многолучевого ионосферного сигнала.

Технический результат достигается тем, что круговую антенную систему, расположенную на поверхности земли, дополняют линейной системой вибраторов, расположенных вдоль вертикали к поверхности земли (например, 5 вибраторов с интервалом между вибраторами 10 м). С помощью сформированной таким образом антенной системы (трехмерная антенная система), многоканального приемника, многоканального аналого-цифрового преобразователя (АЦП) и временного преобразования Фурье формируют пространственно-временной массив комплексных данных E n , m , отображающий значения напряженности поля в n точках трехмерного пространства (n - номер вибратора на поверхности земли и вдоль вертикальной линейки вибраторов) и в m моменты времени, с интервалами 1-2 секунды (индекс m определяет номер временного среза данных на n вибраторах, m=1÷М+1). Количество временных срезов данных берется на единицу больше, чем количество лучей М. Формируют корреляционную матрицу A , размерностью, равной количеству лучей ((М)*(М)), и правый вектор-столбец данных b ¯ = ( b 1 ÷ b м ) T , элементы которых усредняются по индексу n (по вибраторам, по пространству), решая векторное уравнение (4), определяют неопределенные коэффициенты α 1 ÷ α м ,

A α ¯ = b ¯ ,                                 (4)

где

А - корреляционная матрица с элементами E n , k E n , p * ¯ ,

m=1÷M, p=1÷M, n=1÷N, черта сверху означает суммирование по индексу n. b ¯ = ( b 1 ÷ b M ) T - вектор-столбец с координатами b m = E n , M + 1 E n , m * ¯ (черта сверху определяет суммирование по индексу n, по вибраторам, m=1÷M).

α ¯ = ( α 1 ÷ α M ) T - вектор-столбец неопределенных коэффициентов.

Формируют полином Y ( ω ) (5), частотная зависимость которого позволяет оценить значения доплеровских сдвигов частот ωm лучей ионосферного сигнала по М минимумам полинома.

Y ( ω ) = x M + α M x M 1 + α M 1 x M 2 + α 2 x + α 1 = 0 ( 5 )

Где x = exp ( i ω Δ t ) .

Определяют элементы импульсных последовательностей М фильтров для выделения отдельных однолучевых полей E n ( m ) .

C m ( m ) = α 1 ( m ) x * , C 2 ( m ) = α 2 ( m ) + C 1 ( m ) ) x * , , C m ( m ) = α m ( m ) + C m 1 ( m ) ) x * , , C M ( m ) = 1

Выделяют (фильтруют) однолучевые m поля E n ( m ) из совокупности М полей ионосферного сигнала по выражениям (6). Верхний индекс в скобках в выражении (6) определяет номер однолучевого поля.

E n ( 1 ) = E n ,1 C 1 ( 1 ) + E n ,2 C 2 ( 1 ) + + E n , M C M ( 1 ) = U 1 exp ( i K ¯ 1 R ¯ n ) E n ( M ) = E n ,1 C 1 ( M ) + E n ,2 C 2 ( M ) + + E n , M C M ( M ) = U M exp ( i K ¯ M R ¯ n ) ( 6 )

Формируют диаграммы направленности для выделенных М лучей ионосферного сигнала с помощью пространственного преобразования Фурье, сканируют ими по азимуту и углу места и по максимумам углового спектра оценивают азимуты, углы места и амплитуды лучей ионосферного сигнала α m ' , β m ' , U m ' .

Существенные отличия предлагаемого способа пеленгации от прототипа следующие.

- Вместо плоской (двумерной) антенной системы (прототип) создается трехмерная антенная система, с помощью которой решается проблема малых углов места при пеленгации ионосферных сигналов.

- Исключается операция создания четырехмерной диаграммы направленности и сканирование этой диаграммой в четырехмерном пространстве, которая значительно увеличивает время расчета параметров ионосферного сигнала и не позволяет реализовать способ (прототип) с большим количеством лучей.

- Пространственная информация, получаемая с помощью антенной системы, дополняется необходимой информацией во времени (срезы 10 данных) и проводится наиболее оптимальная совместная пространственно-временная обработка данных.

- Определяются элементы М импульсных последовательностей фильтров для выделения М однолучевых полей из суммы полей ионосферного сигнала и выделяются (фильтруются) однолучевые поля ионосферного сигнала.

- Вместо оценок параметров лучей по максимуму четырехмерной диаграммы направленности в прототипе, в предлагаемом способе пеленгации производятся оценки параметров лучей ионосферного сигнала по максимумам двумерных диаграмм направленностей выделенных М полей ионосферного сигнала.

Обоснованием предлагаемого способа оценки угловых параметров ионосферных сигналов является следующее.

Запишем сигнал на n вибраторе и в m точке по времени E n , m в комплексной форме в виде

E n , m = p = 1 M U p e i ( ω p t m K ¯ p ( α 0, p , β 0, p ) R ¯ n ) ( 6 )

где U p - комплексная амплитуда р луча ионосферного сигнала,

K ¯ p ( α 0, p , β 0, p ) - волновой вектор р луча ионосферного сигнала,

α 0, p , β 0, p - значения азимута и угла места р луча ионосферного сигнала,

R ¯ n - местоположение n вибратора,

tm - время отсчетов напряженности поля (m меняется в пределах от 1 до М+1, Р=М - количество лучей в ионосферном сигнале,

ωp - доплеровский сдвиг частоты р луча ионосферного сигнала.

Элементарные вибраторы антенной системы располагаются по окружности радиуса R и вдоль вертикали к поверхности земли Z.

Используя временную последовательность данных, исключим переменные, описывающие поля Р-1 лучей ( U p , K ¯ p ( α 0, p , β 0, p ) , ω p ) . В результате получим выражение (7) с неопределенными коэффициентами C 1 ÷ C M (М=Р), выражение в правой части которого описывает поле лишь одной волны (индекс р у переменных убран).

m = 1 M C m E n , m = U p e i ( ω t m K ¯ ( α 0 , β 0 ) R ¯ n ) ( 7 )

Неопределенные коэффициенты C 1 ÷ C M не зависят от номера вибратора n и номера среза данных m. Они зависят лишь от временного интервала Δt. Следовательно, они одинаковы для выражения (7) при разных значениях индекса m. На основании этого запишем

m = 1 M C m E n , m + 1 = U e i ( ω t m K ¯ ( α 0 , β 0 ) R ¯ n ) e i ω Δ t ( 8 )

Объединяя (7) и (8) и переобозначив коэффициенты, получим

m = 1 M a m E n , m = E n , m = E n , M + 1 ( 9 )

Взаимосвязь коэффициентов C m и a m определяется выражениями

C 1 = a 1 x * , C 2 = ( a 2 + C 1 ) x * , , C m = ( a m + C m 1 ) x * , ,

C M = a = M + 1 1, x * = exp ( i ω Δ t ) ( 10 )

Используя выражение (9) запишем функционал правдоподобия в виде

Δ ( a 1 a M ) = n N | m M a m E n , m E n , M + 1 | 2 ( 11 )

Дифференцируя выражение (11) по коэффициентам a m и приравнивая дифференциалы к нулю, можно получить матричное уравнение. В векторной форме оно имеет вид (4)

A a ¯ = b ¯ ,                                 (4)

где

A - корреляционная матрица с элементами E n , k E n , p * ¯ ,

m=1÷M, p=1÷M, n=1÷N, черта сверху означает суммирование по индексу n, b ¯ = ( b 1 ÷ b M ) T - вектор-столбец с координатами b m = E n , M + 1 E n , m * ¯ (черта сверху определяет суммирование по индексу n, по вибраторам, m=1÷М).

α ¯ = ( α 1 ÷ α M ) T - вектор-столбец неопределенных коэффициентов.

Решая векторное уравнение (4), можно получить значения коэффициентов a 1 ÷ a M вектора a ¯ . Используя коэффициенты a 1 ÷ a M , составим полином, зависящий от частоты ω (5).

Y ( ω ) = x M + a M x M 1 + a M 1 x M 2 + a 2 x + a 1 = 0 ( 5 )

Корни полинома М степени дают М решений x ( m ) = exp ( i ω m Δ t ) . Следовательно, по выражениям (10) можно определить коэффициенты C p ( m ) , которые представляют собой элементы импульсной последовательности фильтра для выделения однолучевого поля из суммы полей ионосферного сигнала. Используя коэффициенты C p ( m ) , можно выделить (отфильтровать) однолучевые поля E n m (выражения (6)).

E n ( 1 ) = E n ,1 C 1 ( 1 ) + E n ,2 C 2 ( 1 ) + + E n , M C M ( 1 ) = U 1 exp ( i K ¯ 1 R ¯ n ) E n ( M ) = E n ,1 C 1 ( M ) + E n ,2 C 2 ( M ) + + E n , M C M ( M ) = U M exp ( i K ¯ M R ¯ n ) ( 6 )

Для каждого однолучевого поля, выполняя пространственное преобразование Фурье, можно получить диаграмму направленности.

U m ( α m ' , β m ' ) = n N E n ( m ) exp ( i K ¯ m ( α m ' , β m ' ) R ¯ n ) , m = 1 ÷ M ( 12 )

По максимумам диаграмм направленности определяют значения азимутов αm, углов места βm и амплитуд U m , для М лучей ионосферного сигнала.

Дисперсия азимута Dα и дисперсия углов места Dβ для однолучевого поля (дисперсия Рао-Крамера) при трехмерной антенной системе определяются выражениями.

D α = 2 σ 2 | E | ( k 0 R cos ( β ) ) 2 N

D β = 6 σ 2 | E | 2 k 0 2 ( ( 3 R sin ( β ) ) 2 N + 2 ( Z cos ( β ' ) ) 2 M 1 ) ( 13 )

где σ2 - дисперсия шума,

| E | - модуль напряженности однолучевого поля,

K0=2π/λ - волновое число, λ - длина волны,

R -радиус антенной системы,

Z - координата, перпендикулярная поверхности земли,

N, M1 - количество вибраторов на окружности радиуса R и вдоль координаты Z.

Дисперсия азимута определяется только вибраторами, расположенными по окружности на земле (N). Дисперсия угла места зависит как от горизонтальной пространственной базы, так и от вертикальной базы антенной системы. При β=0 дисперсия угла места определяется выражением

D β | b = 0 = 3 σ 2 | E | 2 ( k 0 Z ) 2 M 1 ( 14 )

Таким образом, дисперсия угла места зависит от новой пространственной базы:

L = ( R sin ( β ) ) 2 N + 2 ( Z 3 cos ( β ' ) ) 2 M 1 ( 14 )

Эта база не обращается в нулевое значение при значении β, стремящемся к нулю. Следовательно, дисперсия трехмерной антенной системы не будет увеличиваться при приближении β к нулю. Задача малых углов места при пеленгации ионосферных сигналов оказывается решенной.

Сущность предлагаемого способа оценки угловых параметров ионосферных сигналов поясняется следующими рисунками и диаграммами.

На фиг.1 представлена структурная схема устройства, с помощью которого реализуется прототип «Способ пеленгации с учетом корреляционной взаимосвязи между лучами».

На фиг.2 представлена структурная схема устройства, с помощью которого реализуется предлагаемый способ «Способ оценки угловых параметров ионосферного сигнала». По сравнению с прототипом второй, третий, четвертый блоки идентичны блокам прототипа, вместо 1,5÷8 блоков прототипа предлагаются действия новых блоков 1,5÷12. С помощью этих блоков формируется пространственно-временной массив, создаются корреляционная матрица и правый вектор-столбец векторного уравнения. При решении векторного уравнения находятся значения коэффициентов, формируется полином и находятся его корни. Это дает возможность определить элементы импульсной последовательности фильтров и выделить (отфильтровать) однолучевые поля. Определяются диаграммы направленности для отдельных лучей, и по максимумам диаграмм направленностей оцениваются амплитуды и угловые характеристики.

На фиг.3 показано расположение антенн на поверхности земли, составляющих круговую антенную систему (16 вибраторов), дополненную вертикальной линейкой вибраторов (5 вибраторов, установленных через интервал 10 м) «Круг-Z». Радиус антенной системы равен 65 м.

На фиг.4 показаны угломестные диаграммы направленности антенной системы «Круг-Z» на частоте 10 МГц (линия 2), на частоте 15 МГц (линия 1) и угломестная диаграмма направленности плоской круговой антенной системы на частоте 10 МГц (линия 3). В области малых значений углов места плоская антенная система не имеет четко выраженного максимума, а следовательно, углы места определяются с большими погрешностями. Так, например, если плоская волна приходит под углом места 10°, тогда антенная система "Круг-Z" имеет максимум диаграммы направленности при β=10°. В то же время диаграмма направленности плоской антенной системы "Круг" имеет плоскую вершину от 0° до 20°. Это и приводит к проблеме малых углов при приеме ионосферных сигналов.

На фиг.5 показана дисперсия углов места антенной системы «Круг -Z» в зависимости от углов места плоской волны. Сплошная линия представляет расчет дисперсии углов места по выражению (13). Точками отмечены дисперсии углов места, полученных при ее оценке по 50-ти расчетам. Таким образом, добавление к плоской антенной системе типа "Круг" линейки вибраторов на высотах от 0 до 40 м вполне решает проблему малых углов при приеме ионосферных сигналов.

На фиг.6 показаны углы места многолучевого ионосферного сигнала, полученные с помощью восьмиканального пеленгатора, работающего по алгоритму предлагаемого способа пеленгации. С помощью пеленгатора ведется прием сигналов радиостанцией «Би-Би-Си» (Англия). f=6,02 МГц. Отмечается наличие как двухлучевой структуры ионосферного сигнала, так и трехлучевой структуры. Получить такие оценки углов места с помощью классических методов пеленгации не представляется возможным.

Структурная схема устройства, в котором реализуется способ оценки угловых параметров ионосферных сигналов, представлена на фиг.2.

Согласно предлагаемому способу последовательность действий следующая (фиг.2).

1. Принимают ионосферный сигнал с помощью круговой антенной системы, дополненной вертикальной по отношению к поверхности земли линейкой вибраторов. Создают трехмерную антенную систему (блок 1).

2. Преобразуют сигналы в каждом канале (от каждого вибратора) по частоте, усиливают и фильтруют с помощью многоканального приемника (блок 2).

3. Преобразуют аналоговые сигналы на выходах многоканального приемника в цифровую форму с помощью многоканального аналого-цифрового преобразователя (АЦП) (блок 3).

4. Получают частотный спектр сигнала в N каналах с помощью преобразования Фурье. Формируют пространственный массив комплексных данных E n , выделяя максимумы частотных спектров (блок 4).

5. Формируют пространственно-временной массив данных E n , m , повторяя действия 1÷4 через интервалы 1÷2 с. Количество временных срезов данных на единицу больше (m=1÷M+1), чем количество лучей М (блок 5).

6. Формируют корреляционную матрицу A , размерностью, равной количеству лучей ((М)*(М)), и правый вектор-столбец данных b ¯ = ( b 1 ÷ b M ) T , элементы которых усредняются по индексу n (по вибраторам, по пространству) (блок 6).

7. Определяют коэффициенты α ¯ = ( α 1 ÷ α M ) T при решении векторного уравнения (4) (блок 7). A a ¯ = b ¯ .                                 (4)

8. Формируют полином Y ( ω ) , степени М для оценки значений доплеровских сдвигов частот ионосферного сигнала по выражению (5) и оценки по минимумам полинома Y ( ω ) М доплеровских сдвигов частот ионосферного сигнала ωm (блок 8).

Y ( ω ) = x M + a M x M 1 + a M 1 x M 2 + a 2 x + a 1 = 0, x m = exp ( i ω m Δ t ) ( 5 )

9. Определяют элементы импульсных последовательностей фильтров C 1 ( k ) = a 1 x k * , C m k = ( a m + C m 1 k ) x k * , C m ( k ) = 1,

k=1÷M (блок 9).

10. Выделяют (фильтруют) однолучевые поля E n m из совокупности полей ионосферного сигнала по выражениям (6) (блок 10).

E n ( 1 ) = E n ,1 C 1 ( 1 ) + E n ,2 C 2 ( 1 ) + + E n , M C M ( 1 ) = U 1 exp ( i K ¯ 1 R ¯ n ) E n ( M ) = E n ,1 C 1 ( M ) + E n ,2 C 2 ( M ) + + E n , M C M ( M ) = U M exp ( i K ¯ M R ¯ n ) ( 6 )

11. Формируют для каждого выделенного поля E n m диаграммы направленности с помощью пространственного преобразования Фурье по выражению (12) и сканируют этой диаграммой направленности в диапазоне оценочных азимутов α'=0÷360° и углов места β'=0÷90°.

U m ( α m ' , β m ' ) = n N E n ( m ) exp ( i K ¯ m ( α m ' , β m ' ) R ¯ n ) , m = 1 ÷ M ( б л о к  11 ) ( 12 )

12. Оценивают азимуты α m ' = α 0, m , углы места β m ' = β 0, m и амплитуды U m по максимуму диаграммы направленности для М лучей ионосферного сигнала (блок 12).

Способ оценки угловых параметров ионосферных сигналов, включающий прием ионосферного сигнала с помощью круговой антенной системы, состоящей из N вибраторов, расположенных равномерно по окружности радиуса R, преобразование по частоте и усиление сигналов от каждого n вибратора с помощью многоканального приемника, преобразование аналоговых сигналов на выходах многоканального приемника в каждом n канале в цифровую форму с помощью многоканального аналого-цифрового преобразователя, определение амплитуды En и фазы ψn принятого ионосферного сигнала в каждом канале с помощью временного преобразования Фурье, отличающийся тем, что круговую антенную систему дополняют вертикальной линейкой вибраторов с номерами n=N+1÷N1, создавая трехмерную антенную систему, дополнительно с помощью антенной системы, многоканального приемника, многоканального АЦП, блока преобразования Фурье формируют пространственно-временной массив данных En,k, ψn,k, повторяя измерения через временной интервал Δt=1÷2 с, где индекс n=1÷N1 определяет номер вибратора, а индекс k=1÷K определяет номер временного среза данных, причем K=M+1, формируют с помощью пространственно-временного массива данных корреляционную матрицу, размерностью, равной количеству лучей М, с элементами , где индекс p=1÷K, и правый вектор-столбец , где T означает транспонирование вектора-столбца, элементы которых усредняют по индексу n (по вибраторам), решая векторное уравнение A ^ a ¯ = b ¯ , находят значения М коэффициентов вектора , равное количеству лучей, формируют полином степени, равной количеству лучей м, , где x ^ m = exp ( i ω m Δ t ) , и находят М решений (, m=1÷M) этого полинома, равное количеству лучей, определяют с помощью полученных решений элементы M импульсных последовательностей фильтров , равное количеству лучей, с помощью которых выделяют (фильтруют) М однолучевых полей ионосферного сигнала, равное количеству лучей, по которым формируют двумерные диаграммы направленности, оценивают по максимумам диаграмм направленностей азимуты αm и углы места βm М лучей ионосферного сигнала.



 

Похожие патенты:

Изобретение относится к радионавигации и может использоваться в радионавигационных системах для измерения угловых координат подвижных объектов как в азимутальной, так и в угломестной плоскостях относительно задаваемого наземным радиомаяком направления.

Использование: в устройстве для обнаружения сигналов и определения направления на их источник. Сущность: устройство для обнаружения сигналов и определения направления на их источник содержит выполненную определенным образом дискретную антенную решетку (ДАР), включающую N ненаправленных пассивных и М активно-пассивных электроакустических преобразователей, соответствующие им I каналы передачи информации, блок управления характеристикой направленности, блок вычисления относительных координат элементов ДАР, пороговое устройство, вычислитель порога принятия решения, индикатор, блок управления активно-пассивными элементами ДАР, а также формирователь характеристик направленности с временной задержкой сигналов.

Изобретение относится к области гидроакустики и может быть использовано для построения систем обнаружения зондирующих сигналов гидролокаторов, установленных на подвижном носителе.

Изобретение относится к области радиотехники и может быть использовано при решении задач радиопеленгации с помощью переносных (малогабаритных) средств в декаметровом и метровом диапазонах радиоволн.

Изобретение относится к измерительной технике, в частности к пеленгаторам. .

Изобретение относится к области радиолокации и может быть использовано в радиолокаторах поиска и слежения. .

Изобретение относится к радиотехнике и может быть использовано для обнаружения и пеленгации фазоманипулированных сигналов. .

Изобретение относится к вычислительной технике. .

Изобретение относится к тепловым выключателям, предназначенным для защиты электронных и электрических приборов от перегрева при неисправностях, позволяет повысить надежность термовыключателя.

Изобретение относится к конструктивному выполнению средств гидрофизических исследований и может быть использовано, например, при реализации систем акустической томографии или систем пассивного обнаружения шумящих объектов.

Изобретение относится к радиотехнике и может быть использовано при разработке систем для определения координат источника радиоизлучения (ИРИ), а также в пассивной радиолокации. Достигаемый технический результат - повышение точности оценки разности моментов приема сигналов источника радиоизлучения, в двух разнесенных приемных пунктах. Указанный результат достигается за счет того, что в заявленном способе осуществляют прием сигнала в двух разнесенных приемных пунктах, оцифровку напряжения с выхода антенны, обнаружение сигнала, в каждом приемном пункте, оценку разности моментов приема, включающем оценку задержки отраженного сигнала относительно прямого в каждом приемном пункте, оценку разности моментов приема отраженных сигналов, вычисление разности моментов приема прямых сигналов как суммы задержки между прямым и отраженным сигналом в первом приемном пункте и задержки между отраженными сигналами в первом и во втором приемном пункте, минус задержка между прямым и отраженным сигналом во втором приемном пункте. 4 ил.

Предлагаемое устройство относится к контрольно-поисковым средствам, а именно к устройствам обнаружения местоположения людей, оказавшихся под завалами, образовавшимися в результате стихийного (землетрясения, торнадо, цунами и др.) или иного бедствия, и поиска взрывчатых и наркотических веществ, и может быть использовано при техногенных авариях, природных катастрофах, террористических актах и при предотвращении опасных для населения акций. Технической задачей изобретения является повышение помехоустойчивости и достоверности приема и демодуляции сложных сигналов с фазовой манипуляцией путем подавления узкополосных помех. Устройство обнаружения людей под завалами и поиска взрывчатых и наркотических веществ содержит одетый на служебную собаку 1 ошейник 2, мобильный первичный преобразователь 3 и вторичный преобразователь 12. Первичный преобразователь 3 содержит тактильные сенсоры 4.1 и 4.2, коммутатор 5, усилитель 6, модулятор 7, радиопередатчик 8, источник 9 питания, световой 10 и звуковой 11 маячки, задающий генератор 18, фазовый манипулятор 19, триггер 17, однополярный вентиль 20, интегратор 21, пороговый блок 22, ключ 23, усилитель 24 мощности и передающую антенну 25. Вторичный преобразователь 12 содержит вибраторную антенну 26, рамочную антенну 27, усилители 28 и 29 высокой частоты, амплитудные детекторы 30 и 31, блок 32 деления, пороговый блок 33, ключ 15, демодуляторы 14 и 44, перемножители 34, 35, 38 и 39, узкополосные фильтры 36 и 40, фильтры 37 и 41 нижних частот, фазоинверторы 42 и 43, блок 45 вычитания и регистратор 16. 7 ил.

Изобретение относится к области радиотехники и может использоваться при проектировании и эксплуатации комплексов радиопеленгации или систем радиосвязи портативного, мобильного (бортового) и стационарного базирования. Технический результат - повышение устойчивости функционирования методов оценки напряженности электромагнитного или акустического поля Для этого на каждом элементе антенной решетки записывают интервал на временном интервале [0,Т], производят формирование дискретного спектра напряженности поля с использованием процедуры преобразования Фурье, при этом. для каждой из полученных спектральных компонент находят вектор комплексных амплитуд/вспомогательных источников как приближенное решение матрично-векторного уравнения с использованием процедуры квазирешения. Число вспомогательных источников определяется как число наиболее значимых собственных чисел автокорреляционной матрицы принятых антенной решеткой сигналов, т.е. наибольших и отличающихся от остальных по величине не менее, чем на порядок. Далее определяют значения поля спектральной компоненты в произвольной точке плоскости антенной решетки (формируют «виртуальный» канал приема сигналов) как скалярное произведение найденного вектора комплексных амплитуд вспомогательных источников и соответствующего вектора «виртуального» канала приема сигналов. 4 з.п. ф-лы, 1 ил.

Изобретение относится к системам управления безопасностью полетов. Достигаемый технический результат - повышение эффективности систем управления безопасностью полетов. Способ основан на введении в бортовое оборудование воздушных судов системного процессора, который объединен с бортовой системой объективного контроля, бортовой и наземной аппаратурой моноимпульсных вторичных радиолокаторов, работающих в дискретно-адресном режиме и модернизированных до режима «8», радиовещательного автоматического зависимого наблюдения, аппаратурой международной системы спасания терпящих бедствие КОСПАС-САРСАТ, устанавливаемой на навигационных искусственных спутниках Земли Глонасс-К1 и его модификациях. Системный процессор обеспечивает формирование признаков наличия или отсутствия целостности воздушных судов с последующей передачей их экипажу, центрам управления и спасания на основе указанного объединения, что позволяет осуществить в реальном масштабе времени упреждающие, коррективные действия, необходимые для поддержания требуемого уровня безопасности полетов от взлета до посадки, а также оперативное определение местоположения воздушного судна в случае аварийной посадки без использования аварийных бортовых устройств регистрации. 1 з.п. ф-лы, 3 ил.

Изобретение относится к измерительной технике, в частности к пеленгаторам. Достигаемый технический результат - возможность частотной и пространственной селекции источников сигналов. Технический результат достигается тем, что устройство для определения направления на источник сигнала, содержит первую магнитную антенну, ориентированную в направлении Север - Юг, вторую магнитную антенну, ориентированную в направлении Запад - Восток, электрическую антенну с круговой диаграммой направленности, шесть усилителей, десять аналого-цифровых преобразователей (АЦП), персональную электронно-вычислительную машину (ПЭВМ или микропроцессор), блок системы единого времени (GPS или Глонасс), блок связи с абонентами, четыре смесителя, десять управляемых фильтров, четыре коммутатора, пять цифроаналоговых преобразователей (ЦАП), четыре калибратора, формирователь, третью магнитную антенну с круговой диаграммой направленности, а также гониометр. Перечисленные средства выполнены и соединены между собой определенным образом. 1 ил.

Изобретение относится к измерительной технике, в частности к пеленгаторам. Достигаемый техническим результат - возможность частотной и пространственной селекции источников сигналов. Указанный результат достигается тем, что устройство для определения направления на источник сигнала содержит первую магнитную антенну, ориентированную в направлении Север-Юг, первый усилитель, вторую магнитную антенну, ориентированную в направлении Запад-Восток, второй усилитель, последовательно соединенные третью антенну с круговой диаграммой направленности, третий усилитель, а также первый, второй и третий аналого-цифровые преобразователи (АЦП), персональную электронно-вычислительную машину (ПЭВМ или микропроцессор), дополнительно содержит блок системы единого времени (GPS или Глонасс), блок связи с абонентами, первый коммутатор, второй коммутатор, первый управляемый фильтр, четвертый АЦП, третий коммутатор, четвертый коммутатор, второй управляемый фильтр, пятый АЦП, первый цифроаналоговый преобразователь (ЦАП), первый калибратор, второй ЦАП, второй калибратор, третий ЦАП, третий калибратор, четвертый ЦАП, формирователь, а также третий, четвертый и пятый управляемые фильтры, первый и второй смесители, а также гониометр. Перечисленные средства определенным образом выполнены и соединены между собой. 1 ил.

Изобретение относится к области радиолокации, в частности к юстировочным щитам. Юстировочный щит моделирует прямые и зеркально отраженные от земли радиосигналы, идущие от ракеты и цели на конечном участке наведения. Юстировочный щит находится в дальней зоне антенны радиопеленгатора и содержит лазерный и инфракрасный излучатели. Для имитации сигналов от приемоответчика ракеты и сигналов, отраженных от цели, щит снабжен генератором радиоимпульсов с синтезатором частот. Достигается повышение точности юстировки. 3 ил.

Изобретение относится к радиолокации и может использоваться в спутниковых радионавигационных системах и комплексах радиоэлектронного подавления. Достигаемый технический результат - обеспечение возможности определения наклонной дальности до прямолинейно движущейся цели в отсутствии априорных данных о скорости ее движения. Достижение указанного технического результата обеспечивается за счет того, что совместно обрабатывают два последовательных во времени измерения пеленгов (углов азимута на цель β1 и β2) и мощностей сигналов Pc1 и Рс2, принимаемых автономной угломерной системой для рассматриваемых моментов измерений t1, t2, учитывая, что базовая точка измерения соответствует геометрическому центру автономной угломерной системы, а линия Ц1ЦN - траектория движения цели - соответствует точкам Ц1 и Ц2 измерения пеленга на цель в моменты времени t1 и t2, предполагая, что цель движется прямолинейно, затем вычисляют величины: производят два последовательных измерения принятой частоты f1, f2 в моменты времени измерений t1, t2 и на основании производимых измерений определяют величины и Δβ21=β2-β1, затем определяют дальность до цели по формуле: R 2 = c T ( 1 − K f 12 ) ( K f 12 Q 21 + 1 ) cos ( Δ β 21 ) − ( Q 21 + K f 12 ) . 5 ил.

Изобретение относится к радиотехнике и может быть использовано для определения местоположения воздушных судов (объектов) по их радиоизлучениям. Достигаемый технический результат - повышение точности измерения угла места примерно на порядок с расширением области надежного определения дальности с 30 км до 70-80 км, что упрощает, делает более надежным и безопасным привод воздушного судна на необорудованные аэродромы, буровые площадки. Указанный результат достигается за счет того, что способ включает прием сигналов бортового передатчика с помощью антенн, образующих кольцевую решетку, располагаемую вблизи и параллельно земной поверхности, определение по принятым антеннами решетки сигналам азимута объекта, измерение комплексных амплитуд принятых сигналов и преобразование измерений в угловой спектр путем умножения на комплексно-сопряженные диаграммы направленности антенн, суммирования результатов умножения по совокупности антенн и определения квадрата модуля суммы. Одновременно сигналы принимают с помощью дополнительных антенн, располагаемых на нормали к плоскости решетки из ее центра, угловой спектр нормируют на сумму квадратов модулей диаграмм направленности антенн (ДНА) и определяют угол места объекта, как положение его максимума. При этом ДНА определяют в направлении полученного азимута и с учетом коэффициента отражения радиоволн от земной поверхности. 6 ил.

Изобретение относится к радиотехнике и может быть использовано в многопозиционных радиотехнических системах для определения координат источников радиоизлучения (ИРИ). Достигаемый технический результат - увеличение точности оценки временных задержек принимаемых сигналов ИРИ с одновременным уменьшением количества требуемых вычислений. Способ основан на том, что за счет предварительной обработки сигналов ИРИ после их ретрансляции реализуется однократное вычисление взаимокорреляционных функций для оценки временных задержек при распространении сигналов ИРИ. 1 ил.
Наверх