Вычислитель доплеровской скорости движения объекта
Владельцы патента RU 2629642:
Федеральное государственное бюджетное образовательное учреждение высшего образования "Рязанский государственный радиотехнический университет" (RU)
Изобретение относится к вычислительной технике. Технический результат заключается в повышении точности измерения скорости за счет меньшего числа функциональных преобразований и расширении диапазона однозначно измеряемой доплеровской скорости. Вычислитель доплеровской скорости движения объекта содержит блок задержки, блок комплексного сопряжения, блок комплексного умножения, блок усреднения, блок вычисления фазы, блок коррекции пределов измерения, умножитель, ключ, блок вычисления модуля, первый блок памяти, блок управления, пороговый блок, второй блок памяти, синхрогенератор, первый и второй двухканальные ключи, дополнительный блок усреднения, дополнительный блок задержки, дополнительный блок комплексного сопряжения и дополнительный блок комплексного умножения, соединенные определенным образом и осуществляющие межпериодную корреляционную обработку исходных отсчетов. 11 ил.
Изобретение относится к вычислительной технике и предназначено для вычисления на основе корреляционного принципа доплеровской (радиальной) скорости движения объекта; может использоваться в автоматизированных системах управления воздушным движением для обнаружения и измерения скорости летательных аппаратов.
Известен многоканальный неследящий фильтровой измеритель [1], каждый канал которого содержит последовательно соединенные согласованный фильтр и детектор, выходы каналов объединены решающим устройством. Однако данное устройство обладает невысокой точностью измерения.
Известно также устройство обработки сигнала движущейся цели [2], содержащее последовательно включенные блоки задержки, умножитель комплексных чисел и вычитатель. Однако это устройство обладает низкой точностью и неоднозначностью измерения.
Наиболее близким к изобретению является измеритель доплеровских сигналов [3], выбранный в качестве прототипа, содержащий блок задержки, блок комплексного сопряжения, блок комплексного умножения, блок усреднения, блок вычисления фазы, умножитель, ключ, блок вычисления модуля, первый блок памяти, блок управления, пороговый блок, второй блок памяти и синхрогенератор. Однако данное устройство обладает неоднозначностью и невысокой точностью измерения за счет наличия большого числа функциональных преобразований.
Задачей, решаемой в изобретении, является расширение диапазона однозначно измеряемой доплеровской скорости и повышение точности измерения за счет меньшего числа функциональных преобразований при применении совместной обработки неэквидистантных импульсных радиосигналов.
Для решения поставленной задачи в вычислитель доплеровской скорости движения объекта, содержащий блок задержки, блок комплексного сопряжения, блок комплексного умножения, блок усреднения, блок вычисления фазы, блок коррекции пределов измерения, умножитель, ключ, первый блок памяти, блок вычисления модуля, блок управления, пороговый блок, второй блок памяти и синхрогенератор, введены первый и второй двухканальные ключи, дополнительный блок усреднения, дополнительный блок задержки, дополнительный блок комплексного сопряжения и дополнительный блок комплексного умножения.
Дополнительные блоки, введенные в предлагаемый вычислитель, являются известными. Так, соединенные вместе блок задержки, блок комплексного сопряжения и блок комплексного умножения позволяют выделить доплеровский набег фазы за интервал между соседними импульсами. Однако неизвестно совместное применение блока задержки, блока комплексного сопряжения, блока комплексного умножения, первого и второго двухканальных ключей, блока управления, дополнительного блока задержки, дополнительного блока комплексного сопряжения и дополнительного блока комплексного умножения. Новыми являются связи первого и второго двухканальных ключей с блоком комплексного умножения и блоком управления, блока усреднения с первым двухканальным ключом и дополнительным блоком задержки, дополнительного блока усреднения со вторым двухканальным ключом и дополнительным блоком комплексного сопряжения, дополнительного блока комплексного умножения с дополнительным блоком задержки и дополнительным блоком комплексного сопряжения, дополнительного блока комплексного умножения с блоком вычисления фазы, блоком коррекции пределов измерения и блоком вычисления модуля. Связи между синхрогенератором и всеми блоками вычислителя доплеровской скорости движения объекта обеспечивают согласованную обработку неэквидистантных радиосигналов.
Сравнение с техническими характеристиками, известными из опубликованных источников информации, показывает, что заявляемое решение обладает новизной и имеет изобретательский уровень.
Заявляемое решение носит технический характер, осуществимо, воспроизводимо и, следовательно, является промышленно применимым.
На фиг. 1 представлена структурная электрическая схема вычислителя доплеровской скорости движения объекта; на фиг. 2 - блока задержки; на фиг. 3 - блока комплексного сопряжения; на фиг. 4 - блока комплексного умножения; на фиг. 5 - блока усреднения; на фиг. 6 - блока вычисления фазы; на фиг. 7 - блока коррекции пределов измерения; на фиг. 8 - блока присвоения знака; на фиг. 9 - блока вычисления модуля, на фиг. 10 - двухканального ключа; на фиг. 11 - блока управления.
Вычислитель доплеровской скорости движения объекта (фиг. 1) содержит блок 1 задержки, блок 2 комплексного сопряжения, блок 3 комплексного умножения, блок 4 усреднения, блок 5 вычисления фазы, блок 6 коррекции пределов измерения, умножитель 7, ключ 8, первый блок 9 памяти, блок 10 управления, блок 11 вычисления модуля, пороговый блок 12, второй блок 13 памяти, синхрогенератор 14, первый двухканальный ключ 15, второй двухканальный ключ 16, дополнительный блок 17 усреднения, дополнительный блок 18 задержки, дополнительный блок 19 комплексного сопряжения и дополнительный блок 20 комплексного умножения.
Блок 1 задержки и дополнительный блок 18 задержки (фиг. 2) содержат две цифровые линии задержки 21. Блок 2 комплексного сопряжения и дополнительный блок 19 комплексного сопряжения (фиг. 3) содержат инвертор 22. Блок 3 комплексного умножения и дополнительный блок 20 комплексного умножения (фиг. 4) содержат два канала (I, II), каждый из которых включает первый перемножитель 23, второй перемножитель 24 и сумматор 25. Блок 4 усреднения и дополнительный блок 17 усреднения (фиг. 5) содержат два канала (I, II), каждый из которых состоит из (N-2)/2 цифровых линий задержки 26 и (N-3)/2 сумматоров 27. Блок 5 вычисления фазы (фиг. 6) содержит делитель 28 и функциональный преобразователь 29. Блок 6 коррекции пределов измерения (фиг. 7) содержит модульный блок 30, сумматор 31, блок 32 присвоения знака, первый ключ 33, сумматор 34, второй ключ 35 и блок 36 памяти. Блок 32 присвоения знака (фиг. 8) содержит блоки 37, 40 умножения, блок 38 памяти и ограничитель 39. Блок 11 вычисления модуля (фиг. 9) содержат два блока 41 умножения, сумматор 42 и блок 43 извлечения квадратного корня. Первый 15 и второй 16 двухканальные ключи (фиг. 10) содержат два ключа 44. Блок 10 управления (фиг. 11) содержит триггер 45 и элемент НЕ 46.
Вычислитель доплеровской скорости движения объекта работает следующим образом.
В заявляемом вычислителе обрабатывается неэквидистантная когерентно-импульсная последовательность N радиосигналов с чередующимися периодами повторения T1 и Т2, причем T1-T2=ΔT. При отражении радиосигналов от движущейся цели их несущие частоты в соответствующих периодах приобретают доплеровские сдвиги фазы
ϕ1=2πƒдT1, ϕ2=2πƒдT2, Δϕ=ϕ1-ϕ2=2πƒдΔT,
где ƒд=2νrƒн/c - доплеровская частота, νr - радиальная скорость цели, ƒн - несущая частота радиосигналов, с - скорость распространения радиоволн.
Отраженные от цели радиосигналы поступают на вход приемника, в котором усиливаются, в квадратурных фазовых детекторах переносятся на видеочастоту, а затем подвергаются аналого-цифровому преобразованию (соответствующие блоки на фиг. 1 не показаны). На вход вычислителя в одном элементе разрешения по дальности поступают цифровые отсчеты комплексной огибающей
Uk=u1k+iu2k, k=1…N,
где u1k, u2k - цифровые коды действительной и мнимой частей отсчетов Uk.
Входные отсчеты Uk вычислителя (фиг. 1) в блоке 1 задержки (фиг. 2) под управлением синхронизирующих импульсов, вырабатываемых синхрогенератором 14, поочередно задерживаются на интервалы T1 и T2, что обеспечивает синхронность последующего комплексного умножения отсчетов по дальности. Синхрогенератор 14 управляется импульсами синхронизатора радиолокатора (на фиг. 1 не показан), следующими поочередно с интервалами Т1 и Т2. В блоке 2 комплексного сопряжения (фиг. 3) осуществляется комплексное сопряжение задержанного отсчета . Далее в блоке 3 комплексного умножения (фиг. 4) реализуется попарное умножение отсчетов в соответствии с алгоритмом
, k=2…N.
Попарные произведения (корреляции) раздельно для каждого интервала T1 и Т2 соответственно через первый 15 и второй 16 двухканальные ключи раздельно поступают в блок 4 усреднения и в дополнительный блок 17 усреднения (фиг. 5). Поочередная коммутация первого 15 и второго 16 двухканального ключей осуществляется импульсами соответственно с первого и второго выходов блока 10 управления, синхронизируемого также импульсами синхронизатора радиолокатора.
В блоке 4 усреднения (фиг. 5) с помощью линий задержки 26 на интервал Т1+Т2 и сумматоров 27 в каждом элементе разрешения по дальности осуществляется скользящее вдоль азимута когерентное суммирование (накопление) соответствующих интервалу T1 попарных произведений (корреляций). В результате на выходе блока 4 усреднения при нечетном N образуется пропорциональная корреляционному моменту отсчетов, соответствующих интервалу T1, величина
.
В дополнительном блоке 17 усреднения (фиг. 5) осуществляется аналогичное суммирование соответствующих интервалу Т2 попарных корреляций, что приводит к образованию на его выходе пропорциональной корреляционному моменту отсчетов, соответствующих интервалу Т2, величины
.
Величина Y1 на выходе блока 4 усреднения (фиг. 5) по времени предшествует величине Y2 на интервал Т2, что компенсируется соответствующей данному интервалу задержкой Y1 в дополнительном блоке 18 задержки (фиг. 2). В дополнительном блоке 19 комплексного сопряжения (фиг. 3) инвертируется знак мнимой части величины Y2.
Величины Y1 и одновременно поступают соответственно на первые и вторые входы дополнительного блока 20 комплексного умножения (фиг. 4), на выходе которого вычисляется величина
.
Величины ν1 и ν2 поступают на соответствующие входы блока 5 вычисления фазы (фиг. 6), где на основе блока 28 деления и арктангенсного функционального преобразователя 29 вычисляется оценка
.
Последующие преобразования оценки происходят в блоке 6 коррекции пределов измерения (фиг. 7) и зависят от знака величины
. При ν1>0 открыт второй ключ 35, и оценка
через сумматор 34 непосредственно поступает на выход блока коррекции пределов измерения. При ν1<0 открыт первый ключ 33, а второй ключ 35 закрыт. При этом в модульном блоке 30 образуется |argV|, вычитаемый в блоке 31 из величины π, поступающей от блока 36 памяти. Полученной разности в блоке 32 присваивается знак величины ν2.
Блок 32 присвоения знака (фиг. 8) работает следующим образом. На второй вход блока присвоения знака поступает величина , где в блоке 37 умножения производится ее умножение на постоянный множитель из блока 38 памяти с целью масштабирования и дальнейшего ограничения в ограничителе 39 по уровню ±1. Таким образом, после ограничения величина на выходе ограничителя 39 имеет смысл знака величины ν2, который, поступая на первый вход блока 40 умножения, присваивается разности π-|argV|, поступающей с выхода блока 31 на первый вход блока 32 присвоения знака, т.е. на второй вход блока 40 умножения.
Рассмотренные операции позволяют в блоке 5 вычисления фазы сначала найти оценку доплеровского сдвига фазы, находящуюся в интервале [-π/2, π/2], а затем при помощи блока 6 коррекции пределов измерения расширить пределы ее однозначного измерения до интервала [-π, π] в соответствии с алгоритмом
Умножитель 7 (фиг. 1) осуществляет умножение найденной оценки сдвига фазы на весовой коэффициент а, хранящийся в первом блоке 9 памяти, что позволяет найти однозначную оценку доплеровской (радиальной) скорости в соответствии с выражением
,
где - весовой коэффициент.
Для уменьшения вероятности работы вычислителя по шумам в нем исключается выдача полученной оценки на выход в отсутствие отраженного от цели сигнала. В блоке 11 вычисления модуля (фиг. 9) вычисляется величина
,
которая поступает на второй вход порогового блока 12, в котором сравнивается с пороговым уровнем z0, записанным во втором блоке 13 памяти. Если происходит превышение порогового уровня z0, то с выхода порогового блока 12 поступает сигнал разрешения на прохождение результата вычисления с выхода блока 7 умножения через ключ 8 на первый выход вычислителя доплеровской скорости движения объекта. В противном случае ключ 8 разомкнут. Кроме того, сигнал с выхода порогового блока 12, являющегося вторым выходом вычислителя доплеровской скорости движения объекта, может быть использован для отсчета других координат цели, например дальности.
Синхронизация вычислителя доплеровской скорости движения объекта осуществляется подачей на все блоки заявляемого вычислителя последовательности синхронизирующих импульсов, вырабатываемых синхронизатором 14 (фиг. 1) с периодом повторения tк, определяемым из условия требуемой разрешающей способности по дальности.
Достижение технического результата объясняется следующим образом. В известном устройстве (прототипе) исходные доплеровские сдвиги фазы ϕ1 и ϕ2, по которым вычисляется величина Δϕ=ϕ1-ϕ2, имеют интервал однозначного измерения [-π, π], что соответствует интервалу однозначного измерения доплеровской частоты [-1/2T1, 1/2T1] (по величине большего периода T1). В предложенном устройстве величина Δϕ измеряется непосредственно, что соответствует интервалу однозначности доплеровских частот [-1/2ΔT, 1/2ΔT]. При этом интервал однозначного измерения доплеровской частоты и, следовательно, доплеровской (радиальной) скорости расширяется в T1/ΔT раз, что соответствует решению поставленной задачи изобретения. Если в соответствии с условием ƒд≤1/2ΔТ и с учетом ƒд=2νrƒн/с для максимально возможной скорости движущегося объекта νr max выбрать интервал ΔT≤с/4νr maxƒн, то во всем диапазоне реальных скоростей движущегося объекта может быть осуществлено их однозначное измерение. При этом сохраняется однозначность измерения дальности, которая обеспечивается соответствующим выбором меньшего периода повторения импульсов Т2.
Обусловленные функциональными преобразованиями погрешности раздельного вычисления величин ϕ1 и ϕ2 являются статистически независимыми. В результате погрешность (дисперсия) разности ϕ1-ϕ2=Δϕ удваивается. В предложенном вычислителе при непосредственном вычислении величины Δϕ такое удвоение отсутствует, что соответствует повышению точности измерения.
Таким образом, вычислитель доплеровской скорости движения объекта позволяет расширить диапазон однозначно измеряемой доплеровской (радиальной) скорости и повысить точность измерения за счет меньшего числа функциональных преобразований при применении предлагаемой совместной обработки неэквидистантных когерентно-импульсных радиосигналов.
Источники информации
1. Ширман Я.Д. и Манжос В.Н. Теория и техника обработки радиолокационной информации на фоне помех. - М.: Радио и связь. - 1981. - С. 204. - Рис. 14.2.
2. Патент №63-49193 (Япония), МПК G01S 13/52. Радиолокационное устройство для обнаружения движущейся цели / К.К. Тосиба. Опубл. 03.10.1988. - Изобретения стран мира. - 1989. - Выпуск 109. - №15. - С. 52.
3. Патент №2017167 (Россия), МПК G01S 13/58. Обнаружитель-измеритель доплеровских сигналов / Д.И. Попов, С.В. Герасимов и Е.Н. Матаев. Опубл. 30.07.1994. - Изобретения. - 1994. - №14. - С. 121.
Вычислитель доплеровской скорости движения объекта, содержащий блок задержки, блок комплексного сопряжения, блок комплексного умножения, блок усреднения, блок вычисления фазы, блок коррекции пределов измерения, умножитель, ключ, первый блок памяти, блок управления, блок вычисления модуля, пороговый блок, второй блок памяти и синхрогенератор, при этом выходы блока задержки соединены с входами блока комплексного сопряжения, выходы которого соединены с первыми входами блока комплексного умножения, вторые входы которого объединены с входами блока задержки, выход первого блока памяти соединен с первым входом умножителя, выход которого соединен с входом ключа, выход порогового блока соединен с управляющим входом ключа, первый вход порогового блока соединен с выходом второго блока памяти, выход синхрогенератора соединен с синхровходами блока задержки, блока комплексного сопряжения, блока комплексного умножения, блока усреднения, блока вычисления фазы, умножителя, ключа, блока вычисления модуля, первого и второго блоков памяти, блока коррекции пределов измерения и порогового блока, отличающийся тем, что введены первый и второй двухканальные ключи, дополнительный блок усреднения, дополнительный блок задержки, дополнительный блок комплексного сопряжения и дополнительный блок комплексного умножения, при этом выходы блока комплексного умножения соединены с объединенными входами первого и второго двухканальных ключей, управляющие входы которых соединены соответственно с первым и вторым выходами блока управления, выходы первого двухканального ключа соединены с входами блока усреднения, выходы которого соединены с входами дополнительного блока задержки, выходы второго двухканального ключа соединены с входами дополнительного блока усреднения, выходы которого соединены с входами дополнительного блока комплексного сопряжения, выходы дополнительного блока задержки соединены с первыми входами дополнительного блока комплексного умножения, вторые входы которого соединены с выходами дополнительного блока комплексного сопряжения, выходы дополнительного блока комплексного умножения соединены с входами блока вычисления фазы, блока вычисления модуля, вторым и третьим входами блока коррекции пределов измерения, выход которого соединен со вторым входом умножителя, выход блока вычисления модуля соединен со вторым входом порогового блока, выход синхрогенератора соединен с синхровходами первого и второго двухканальных ключей, дополнительного блока усреднения, дополнительного блока задержки, дополнительного блока комплексного сопряжения и дополнительного блока комплексного умножения, причем входами вычислителя доплеровской скорости движения объекта являются входы блока задержки, а первым и вторым выходами - соответственно выходы ключа и порогового блока.