способ определения скорости и направления движения носителя обзорной рлс
Классы МПК: | G01S13/58 для определения скорости или траектории движения; для определения знака направления движения |
Автор(ы): | Лихачев Владимир Павлович (RU), Нелюбов Дмитрий Михайлович (RU), Уткин Владимир Владимирович (RU) |
Патентообладатель(и): | Государственное образовательное учреждение высшего профессионального образования "Воронежское высшее военное авиационное инженерное училище (военный институт)" (RU) |
Приоритеты: |
подача заявки:
2007-12-20 публикация патента:
20.05.2010 |
Изобретение относится к области радиотехники и может быть использовано при создании многопозиционных пассивных радиолокационных станций (РЛС). Достигаемым техническим результатом изобретения является повышение точности и сокращение времени определения скорости и направления движения носителя обзорной РЛС в многопозиционной пассивной РЛС. Указанный результат достигается тем, что в j-м приемном пункте, где , j g, определяют пеленг на обзорную РЛС j и частоту сигнала обзорной РЛС f0j и передают их значения на ведущий приемный пункт, на ведущем приемном пункте определяют пеленг на обзорную РЛС g и разность частот Fgj=f0g-f0j, где f 0g - частота сигнала обзорной РЛС, измеренная ведущим приемным пунктом, рассчитывают угол между вектором скорости носителя обзорной РЛС и линией, соединяющей носитель обзорной РЛС и j-й приемный пункт, а также скорость носителя обзорной РЛС, по определенным вычислительным формулам. 7 ил., 1 табл.
Формула изобретения
Способ определения скорости и направления движения носителя обзорной РЛС, основанный на использовании многопозиционной пассивной РЛС, включающей J 3 приемных пунктов, g-й из которых, где , является ведущим, отличающийся тем, что в j-м приемном пункте, где , j g, определяют пеленг на обзорную РЛС j и частоту сигнала обзорной РЛС f0j и передают их значения на ведущий приемный пункт, на ведущем приемном пункте определяют пеленг на обзорную РЛС
k и разность частот Fgj=f0g-f0j, где f 0g - частота сигнала обзорной РЛС, измеренная ведущим приемным пунктом, рассчитывают угол между вектором скорости носителя обзорной РЛС и линией, соединяющей носитель обзорной РЛС и j-й приемный пункт по формуле
,
где - угол между вектором скорости носителя обзорной РЛС и линией, соединяющей носитель обзорной РЛС и g-й приемный пункт; gj= g- j - угол пересечения линий положения g-го и j-го приемных пунктов в точке истинного положения носителя обзорной РЛС;
g(j+1)= g- j+1 - угол пересечения линий положения g-го и (j+1)-го приемных пунктов в точке истинного положения носителя обзорной РЛС; знак «+» принимают, если j-й приемный пункт находится слева от g-го приемного пункта, знак «-» выбирают, если j-й приемный пункт находится справа от g-го приемного пункта; а скорость носителя обзорной РЛС определяют по формуле
где 0 - длина волны зондирующего сигнала обзорной РЛС.
Описание изобретения к патенту
Изобретение относится к области радиотехники, в частности к способам и технике пассивной радиолокации, и может быть использовано при создании многопозиционных пассивных РЛС.
Известны традиционные радиотехнические способы определения скорости носителя обзорной РЛС группировкой наземных приемных пунктов: триангуляционный (угломерный), разностно-дальномерный и их различные сочетания. Так, при триангуляционном способе [1] скорость носителя обзорной РЛС определяется как производная от расстояния, пройденного носителем обзорной РЛС за промежуток времени между несколькими оценками его местоположения. Данный способ отличается простотой реализации и минимально необходимым вычислительным ресурсом.
Разностно-дальномерный способ [2, стр.321-322], [3, стр.238] позволяет определить скорость носителя обзорной РЛС по результатам нескольких измерений разностей времени прихода сигналов обзорной РЛС на ведущий и ведомые приемные пункты.
Существенным недостатком способов-аналогов является необходимость завязки траектории движения носителя обзорной РЛС для последующей оценки его направления движения, что потребует значительного времени и многократных измерений его местоположения. Кроме того, при разностно-дальномерном способе оценивания скорости носителя при большой частоте следования импульсов обзорной РЛС появляется неоднозначность измерения его местоположения, которую устранить практически невозможно, а также имеются ограничения на максимальное расстояние между приемными пунктами, так как они должны одновременно находиться в створе диаграммы направленности антенны обзорной РЛС.
Наиболее близким к предлагаемому изобретению (прототипом) является способ определения скорости и направления движения носителя обзорной РЛС многопозиционной пассивной радиолокационной станцией (МПРЛС), описанный в [3, стр.382-387].
Сущность способа-прототипа поясняется фиг.1, где представлено взаимное расположение МПРЛС и носителя обзорной РЛС на плоскости. На фиг.1 обозначены: и Vrj - вектор скорости носителя обзорной РЛС и его радиальная скорость относительно j-го приемного пункта (ПП), где ; Rgj - расстояние между ведущим (g-м) и j-м ПП, которое известно и постоянно, , g j; 1gj, 2gj - углы пересечения линий положения g-го и j-го ПП в точке истинного положения носителя обзорной РЛС при первой и второй оценке его местоположения; D1g, D1 j и D2g, D2j - дальности от g-го и j-го ПП до носителя обзорной РЛС соответственно при первой и второй оценке его местоположения; - расстояние, которое преодолевает носитель обзорной РЛС за время
Тпер - время перемещения носителя обзорной РЛС из одного элемента разрешения по угловой координате g-го ПП в другой элемент разрешения; ППg - ширина диаграммы направленности (ШДН) антенны g-го ПП; V - тангенциальная скорость носителя обзорной РЛС, максимальная относительно одного из ПП, например g-го ПП; , и , - средние значения пеленгов на носитель обзорной РЛС, определенные g-м и j-м ПП соответственно при первой и второй оценке местоположения носителя обзорной РЛС.
Определение скорости и направления движения носителя обзорной РЛС в способе-прототипе основано на выполнении следующих операций (фиг.2):
1. Прием m пачек радиоимпульсов от обзорной РЛС j-м ПП в течение времени нахождения носителя обзорной РЛС в луче его диаграммы направленности
где Amj(t) - функция изменения огибающей m-й пачки радиоимпульсов, принятой j-м ПП; S(t) - функция модуляции единичного радиоимпульса длительностью u; Tu - период следования импульсов в пачке; nm(t) - закон изменения фазы единичного радиоимпульса в m-й пачке; f0 - частота зондирующего сигнала обзорной РЛС;
0 - длина волны зондирующего сигнала обзорной РЛС; - количество радиоимпульсов в пачке, обусловленное ШДН антенны, частотой следования импульсов и скоростью обзора пространства РЛС.
2. Измерение j-м ПП пеленга на работающую обзорную РЛС 1jm по каждой m-й принятой пачке импульсов обзорной РЛС и получение усредненного значения пеленга на обзорную РЛС
3. Передача усредненного значения пеленга на ведущий ПП.
4. Расчет значения угла пересечения линий положения g-го и j-го ПП в точке истинного положения носителя обзорной РЛС
5. Расчет значения дальности от носителя обзорной РЛС до ведущего ПП [4, стр.150]
6. Повторение операций п.п.1-5 при перемещении носителя обзорной РЛС в соседний элемент разрешения по угловой координате (луч диаграммы направленности) j-го ПП с целью определения второго среднего значения дальности от носителя обзорной РЛС до ведущего ПП
7. Расчет значения скорости носителя обзорной РЛС относительно ведущего ПП
Точность определения скорости носителя обзорной РЛС способом-прототипом определяется следующим выражением [5, стр.278]
где
- среднеквадратические ошибки (СКО) определения точек пересечения линий положения носителя обзорной РЛС g-го и j-го ПП при первой и второй оценке его местоположения соответственно; z - количество оценок местоположения носителя обзорной РЛС; D1g=D1g 1g, D1j=D1j 1j, D2g=D2g 2g, D2j=D2j 2j [7, стр.234] - СКО определения линий положения g-го и j-го ПП при первой и второй оценке местоположения носителя обзорной РЛС соответственно; , , , - СКО определения пеленгов на обзорную РЛС g-м и j-м ПП при первой и второй оценке его местоположения соответственно, определяемые по формуле [2, стр.186]
где qg, qj - отношения сигнал/шум по мощности для соответствующего ПП.
Значения qg, qj определяются по формуле [6, стр.118]
где [6, стр.236] энергия сигнала РЛС на входе соответствующих ПП, Ри - импульсная мощность обзорной РЛС; - длительность пачки радиоимпульсов обзорной РЛС;
РЛС - ШДН антенны обзорной РЛС; Тобз - время обзора пространства РЛС; GРЛС - коэффициент усиления антенны обзорной РЛС; GБЛg, GБЛj - коэффициенты, характеризующие уровень боковых лепестков диаграммы направленности антенны обзорной РЛС в направлении на соответствующий ПП; - скважность радиоимпульсов, излучаемых обзорной РЛС; N0=kбT°Kш [6, стр.246] - спектральная плотность шума на входе g-го и j-го ПП, kб=1,38·10-23 [Дж/град] - постоянная Больцмана, T°=300К -абсолютная температура в градусах Кельвина, Кш - коэффициент шума приемных устройств g-го и j-го ПП; - коэффициент, характеризующий потери энергии излученного сигнала РЛС в антенно-фидерном тракте ПП.
К недостаткам способа-прототипа относятся:
- низкая точность определения результирующей скорости носителя обзорной РЛС при малом количестве z;
- большое время, затрачиваемое на оценку направления движения носителя обзорной РЛС.
Техническим результатом изобретения является повышение точности и сокращение времени определения скорости и направления движения носителя обзорной РЛС в многопозиционной пассивной РЛС.
Указанный технический результат достигается тем, что в j-м приемном пункте определяют пеленг на обзорную РЛС j и частоту сигнала обзорной РЛС f0j и передают их значения на ведущий приемный пункт, на ведущем приемном пункте определяют пеленг на обзорную РЛС g и разность частот Fgj=f0g-f0j, где f 0g - частота сигнала обзорной РЛС, измеренная ведущим приемным пунктом, рассчитывают угол между вектором скорости носителя обзорной РЛС и линией, соединяющей носитель обзорной РЛС и j-й приемный пункт по формуле
где - угол между вектором скорости носителя обзорной РЛС и линией, соединяющей носитель обзорной РЛС и g-й приемный пункт; gj= g- j - угол пересечения линий положения g-го и j-го приемных пунктов в точке истинного положения носителя обзорной РЛС;
g(j+1)= g- j+1 - угол пересечения линий положения g-го и j+1-го приемных пунктов в точке истинного положения носителя обзорной РЛС; знак «+» принимается, если j-й приемный пункт находится слева от g-го приемного пункта, знак «-» выбирается, если j-й приемный пункт находится справа от g-го приемного пункта; а скорость носителя обзорной РЛС определяют по формуле
Сущность изобретения поясняется фиг.3, где представлено взаимное расположение приемных пунктов и носителя обзорной РЛС на плоскости для случая J=3.
На фиг.2 обозначены:
R31 и R32 - расстояния между ПП, которые известны и постоянны;
- вектор скорости носителя обзорной РЛС и радиальные скорости относительно соответствующего ПП;
1, 2 и 3 - углы между вектором скорости носителя обзорной РЛС и линиями, соединяющими носитель обзорной РЛС с соответствующим ПП;
1, 2, 3 - пеленги на обзорную РЛС, определенные соответствующими ПП;
31, 32 - углы пересечения линий положения ПП3, ПП1 и ПП3, ПП2 соответственно в точке истинного положения носителя обзорной РЛС.
В качестве ведущего ПП назначен ПП3. Принято, что носитель обзорной РЛС движется прямолинейно и равномерно относительно МПРЛС.
В предлагаемом способе выполняется следующая последовательность операций (фиг.4):
1. Прием первой пачки радиоимпульсов (1) от обзорной РЛС каждым ПП.
2. Определение каждым ПП пеленгов на работающую обзорную РЛС 1, 2, 3 и несущей частоты обзорной РЛС f01 , f02, f03.
3. Передача значений пеленгов 1, 2 и несущих частот сигнала обзорной РЛС f 01, f02 на ПП3.
4. Определение на ПП3 значений углов пересечения линий положения ПП1, ПП3 и ПП2, ПП3 в точке истинного положения носителя обзорной РЛС
5. Определение на ПП3 значений разностей частот сигнала (1) F31=f03-f01 и F32=f03-f02.
6. Расчет на ПП3 значений углов 1, 2 и 3
7. Определение на ПП3 скорости носителя обзорной РЛС
и усреднение полученных оценок
Таким образом, в предлагаемом способе определения скорости и направления движения носителя обзорной РЛС новыми существенными признаками изобретения являются:
1. Оценка несущих частот сигналов обзорной РЛС в каждом ПП.
2. Вычисление на ведущем ПП по первой пачке радиоимпульсов значений разностей частот сигнала обзорной РЛС.
3. Расчет углов 1, 2 и 3 (12)-(13) и оценка скорости носителя обзорной РЛС (14)-(15).
На фиг.5 представлен вариант построения МПРЛС, состоящей из J ПП, в которой реализован предлагаемый способ определения скорости и направления движения обзорной РЛС. Каждый ПП состоит из однотипного приемного устройства (ПрУ), в котором определяются несущая частота f01 f0J и пеленг на обзорную РЛС 1 J. Ведущий ПП (g-й ПП), кроме того, включает в себя: вычитающие устройства 1.1 1.j, осуществляющие вычисление значений Fg1 FgJ и g1 gJ, два делителя 2; J-1 блоков вычисления скорости носителя обзорной РЛС (БВСН РЛС 1 БВСН РЛС J-1); J-1 блоков вычисления направления движения носителя обзорной РЛС (углов между вектором скорости носителя обзорной РЛС и линией, соединяющей носитель обзорной РЛС и соответствующий ПП 1 J БВНД РЛС 1 БВНД РЛС J-1) и сумматор 5, реализующие вычисления по формулам (12)-(13) и (14)-(15).
Каждый БВСН РЛС включает в себя следующие устройства: три вычислителя синуса 3.g, 3.j, 3.j+1, осуществляющие определение синуса углов g, j и j+1; два вычитающих устройства 1.g и 1.j; два делителя 2.j и 2.j+1 и два умножителя 7.j и 7.j+1.
Каждый БВНД РЛС включает в себя: три вычитающих устройства 1.g, 1.j, 1.j+1; два делителя 2.j; два вычислителя синуса 3.j и 3.j+1, осуществляющие вычисление синусов углов gj и g(j+1), два вычислителя косинуса 4.j, 4.j+1, осуществляющие вычисление косинусов углов gj и g(j+1); три сумматора 5.g, 5.j. 5.j+1; вычислитель арктангенса 6.j; два умножителя 7.j, 7.j+1 и две схемы выбора знака 8.j.
Порядок определения скорости и направления движения носителя обзорной РЛС в предлагаемой МПРЛС заключается в следующем.
Сигнал обзорной РЛС поступает на вход всех приемных пунктов МПРЛС. В ПрУ каждого ПП вычисляются пеленги на обзорную РЛС 1 J и значения несущей частоты f01 f0J. Полученные значения пеленгов и несущих частот сигнала обзорной РЛС передаются на g-й ПП.
На ведущем ПП полученные значения пеленгов и несущих частот поступают на первые входы соответствующих устройств 1.1 1.J. На их вторые входы подаются значения f0g и g. В результате на выходах устройств 1.1 1.J формируются значения Fg1 FgJ и g1 gJ. Значение f0g также подается на вход устройства 2, на второй вход которого поступает постоянное значение скорости света (с). В результате на выходе устройства 2 формируется значение 0.
Значения Fg1 FgJ и g1 gJ поступают на соответствующие блоки вычисления направления движения носителя обзорной РЛС, где проводятся следующие операции.
В блоках 3.1 3.J и 4.1 4.J вычисляются значения синусов и косинусов углов g1 gJ соответственно. Далее в устройствах 1.j и 1.j+1 проводятся операции вычитания постоянного значения «1» из полученных значений косинусов углов gj и g(j+1). Полученные разности подаются на входы устройств 7.j, на вторые входы которых подается результат вычислений в блоке 2.j. С выходов устройств 7.j полученные произведения подаются на сумматоры 5.j, выходы которых подключены к устройству 2.j. Результат деления поступает на устройство 6.j, на выходе которого формируется значение угла g.
Значение угла g подается на вход устройств 5.j и 1.j. На вторые входы устройств подаются значения углов gj и g(j+1) с выходов устройств 8.j. Значение угла gj или g(j+1) поступает на выход соответствующего устройства 8.j только при наличии разрешающего сигнала определенного уровня «Л» или «П» (сигнал о местоположении соответствующего ПП слева или справа от ведущего ПП)
В итоге на выходе устройства 5.j формируется значение угла j или j+1, на выходе устройства 1.j - значение угла j+1 или j.
В результате вычислений на выходах БВНД РЛС формируются значения углов
1 J, которые характеризуют положение носителя обзорной РЛС относительно каждого ПП МПРЛС.
Значения углов g, j и j+1 поступают на устройства 3.j соответствующих БВСН РЛС. С выходов устройств 3.j значения синусов углов g, j и j+1 поступают на устройства 1.j. Полученные разности синусов поступают на устройства 2.j, на вторые входы которых поступают значения Fgj и Fg(j+1). Результаты деления подаются на устройства 7.j, на вторые входы которых подается значение 0 с выхода устройства 2. С выходов устройств 7.j значения скорости носителя обзорной РЛС и поступают на устройство 5. Сумма всех значений полученная во всех БВСН РЛС, с выхода устройства 5 поступает на вход устройства 2, на второй вход которого подается постоянное значение J-1. В результате на выходе устройства 2 формируется значение скорости носителя обзорной РЛС .
Все устройства, входящие в состав ведущего ПП и имеющие одинаковый номер, аналогичны по составу и принципу работы. Конструкция устройств 1 7 основана на использовании известных аналоговых и цифровых элементов и технических трудностей в реализации не представляет. Схема выбора знака 8 может состоять из мультиплексора с одним управляемым входом.
Оценка реализуемости и эффективности предлагаемого способа определения скорости и направления движения обзорной РЛС проводилась методом математического моделирования на ЭВМ. Операции 1-4 (фиг.4) известны, их реализация аналогична прототипу, а операции 5-7 были реализованы программно в среде MathCAD.
В результате моделирования была оценена точность определения скорости и направления движения носителя обзорной РЛС предлагаемым способом и степень влияния на нее величин, входящих в выражения (12)-(15).
В качестве исходных данных было принято, что носитель обзорной РЛС движется прямолинейно и равномерно относительно МПРЛС со скоростью ; значение Ри=103 Вт; длительность импульса и=10-6 с; период следования Т и=10-3 с, следовательно, Q=103; 0=0,03 м, значение РЛС в азимутальной и угломестной плоскости составляет 3°, тогда GРЛС=4000; Тобз=5 с, тогда пачки=42·10-3 с; дальность от ведущего ПП до носителя обзорной РЛС D3 принимала значения 50, 100 и 150 км; расстояния R31=R32=20 км; значение ПП в азимутальной и угломестной плоскости составляет 1,5°.
Были заданы 10000 реализаций углов пересечения линий положения 31 и 32 в точке истинного положения носителя обзорной РЛС и разностей частот сигнала обзорной РЛС F31, F32, распределенных по нормальному закону с математическим ожиданием 31= 3- 1, 32= 3- 2, , соответственно,
где Vr31=Vr3-Vr1, Vr32=Vr3-Vr2 и дисперсией, определяемой выражениями
где , , определяются по формуле (7); , , - СКО определения частоты сигнала обзорной РЛС соответствующим ПП, вычисляемая по формулам [2, стр.185]:
где q1, q2, q 3 - отношения сигнал/шум по мощности соответствующего ПП, определяемые по формуле (8).
Подставляя полученные в (7), (8) и (18) значения в (16) и (17) и проведя нелинейные преобразования сформированных случайных последовательностей были получены значения математического ожидания и дисперсии оценки скорости носителя обзорной РЛС . При заданных дальностях D3, равных 50, 100 и 150 км, значения составили 2,963, 3,471 и 3,929 м/с соответственно. Полученные значения были сопоставлены с потенциальной точностью оценки скорости движения носителя обзорной РЛС, определяемой по формуле [4, стр.140]
На фиг.6 представлены зависимости от времени наблюдения носителя обзорной РЛС t=zTпер для способа-прототипа. На фиг.7 представлены зависимости потенциальной от времени tп=m· пачки для заданных условий и значения , полученные в предлагаемом способе при значениях D3=50, 100 и 150 км. Из анализа графических зависимостей видно, что СКО определения скорости носителя обзорной РЛС в предлагаемом способе в отличие от прототипа практически не зависит от дальности до него.
Результаты сравнительного анализа точностей оценки скорости носителя обзорной РЛС способом-прототипом и предлагаемым способом для приведенных выше условий представлены в таблице.
Как видно из таблицы, в способе-прототипе для достижения точности определения скорости носителя обзорной РЛС, полученной в предлагаемом способе, необходимо выполнить не менее 6 оценок местоположения носителя обзорной РЛС, что составит в зависимости от дальности до носителя обзорной РЛС более 25 с. Следовательно, выигрыш по времени, возникающий в результате применения предлагаемого способа, составляет более 590 раз и возможен при приеме одной пачки импульсов от обзорной РЛС.
Таким образом, введение новых существенных признаков позволяет в МПРЛС достичь более высокой точности определения скорости и направления движения носителя обзорной РЛС при меньшем времени наблюдения по сравнению с прототипом.
Источники информации
1. Мамошин В.Р. Патент РФ № 2279105 от 02.08.2004 г. на изобретение «Комплексный способ определения координат и параметров траекторного движения авиационно-космических объектов, наблюдаемых группировкой станций слежения». БИ № 18 от 27.06.2006 г.
2. Смирнов Ю.А. Радиотехническая разведка. - М.: Воениздат, 2001. - 456 с.
3. Черняк B.C. Многопозиционная радиолокация. - М.: Радио и связь, 1993. - 416 с.
4. Караваев В.В., Сазонов В.В. Статистическая теория пассивной локации. - М.: Радио и связь, 1987. - 240 с.
5. Сосулин Ю.Г. Теоретические основы радиолокации и радионавигации: Учеб. пособие для вузов - М.: Радио и связь, 1992. - 304 с.
6. Ширман Я.Д. Теоретические основы радиолокации: Учеб. пособие для вузов. - М.: Сов. радио, 1970. - 560 с.
7. Казаринов Ю.М. Радиотехнические системы. Учеб. для вузов. - М.: Высшая школа. 1990. - 496 с.
Класс G01S13/58 для определения скорости или траектории движения; для определения знака направления движения