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

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

Иллюстрации

Показать все

Реферат

Изобретение относится к радиотехнике, а именно к области пеленгации, и может быть использовано для измерения азимутов и углов места ионосферных сигналов в условиях приема нескольких интерферирующих лучей, в том числе при малых углах места. При приеме ионосферных сигналов в точку приема, как правило, приходят несколько лучей (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) Определяют амплитуды E_n и фазы ψ_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 лучей требуется время Т_Т=Т*M^2N. Третьим недостатком является наличие проблемы малых углов места при данном способе пеленгации ионосферных сигналов, обусловленное плоской (двумерной) антенной системой.

Блок-схема технической реализации данного способа пеленгации (прототипа) представлена на фиг.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. Определяют амплитуды E_n и фазы ψ_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. Определяют ранг поля (однолучевое или двулучевое) по условию U₀₁/U₀₂≥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 ' .

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

- Вместо