RU2613037C1 - Calculator of range rate of moving object - Google Patents

Calculator of range rate of moving object Download PDF

Info

Publication number
RU2613037C1
RU2613037C1 RU2015150365A RU2015150365A RU2613037C1 RU 2613037 C1 RU2613037 C1 RU 2613037C1 RU 2015150365 A RU2015150365 A RU 2015150365A RU 2015150365 A RU2015150365 A RU 2015150365A RU 2613037 C1 RU2613037 C1 RU 2613037C1
Authority
RU
Russia
Prior art keywords
unit
block
additional
inputs
outputs
Prior art date
Application number
RU2015150365A
Other languages
Russian (ru)
Inventor
Дмитрий Иванович Попов
Original Assignee
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Рязанский государственный радиотехнический университет"
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Рязанский государственный радиотехнический университет" filed Critical Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Рязанский государственный радиотехнический университет"
Priority to RU2015150365A priority Critical patent/RU2613037C1/en
Application granted granted Critical
Publication of RU2613037C1 publication Critical patent/RU2613037C1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems

Abstract

FIELD: measurement technology.
SUBSTANCE: calculator of radial velocity of the moving object comprises a delay unit, a complex conjugation block complex multiplication, averaging unit, a phase calculation unit correcting measurement limits, a multiplier, a key calculation unit module, the first memory unit, a control unit, a threshold unit, a second memory unit , the timing, the first and second two-channel keys, additional averaging unit, an additional delay unit, an additional block of complex conjugation and additional complex multiplication unit connected in a certain way and performing between-period correlation processing raw samples.
EFFECT: improving the accuracy of measuring the speed due to the smaller number of functional transformation and expansion of the range of uniquely measured radial velocity while maintaining a clear distance measurement.
11 dwg

Description

Изобретение относится к вычислительной технике и предназначено для вычисления на основе корреляционного принципа радиальной скорости движущегося объекта; может использоваться в автоматизированных системах управления воздушным движением для обнаружения и измерения скорости летательных аппаратов.The invention relates to computer technology and is intended for calculation based on the correlation principle of the radial velocity of a moving object; can be used in automated air traffic control systems to detect and measure the speed of aircraft.

Известен многоканальный неследящий фильтровой измеритель [1], каждый канал которого содержит последовательно соединенные согласованный фильтр и детектор, выходы каналов объединены решающим устройством. Однако данное устройство обладает невысокой точностью измерения.Known multi-channel non-tracking filter meter [1], each channel of which contains a matched filter and detector connected in series, the outputs of the channels are combined by a resolver. However, this device has a low measurement accuracy.

Известно также устройство обработки сигнала движущейся цели [2], содержащее последовательно включенные блоки задержки, умножитель комплексных чисел и вычитатель. Однако это устройство обладает низкой точностью и неоднозначностью измерения.A device for processing a signal of a moving target [2] is also known, which contains series-delayed blocks, a complex number multiplier, and a subtractor. However, this device has low accuracy and ambiguity of measurement.

Наиболее близким к изобретению является измеритель доплеровских сигналов [3], выбранный в качестве прототипа, содержащий блок задержки, блок комплексного сопряжения, блок комплексного умножения, блок усреднения, блок вычисления фазы, умножитель, ключ, блок вычисления модуля, первый блок памяти, блок управления, пороговый блок, второй блок памяти и синхрогенератор. Однако данное устройство обладает неоднозначностью и невысокой точностью измерения за счет наличия большого числа функциональных преобразований.Closest to the invention is a Doppler signal meter [3], selected as a prototype, comprising a delay unit, a complex conjugation unit, a complex multiplication unit, an averaging unit, a phase calculation unit, a multiplier, a key, a module calculation unit, a first memory unit, a control unit , a threshold block, a second memory block, and a clock generator. However, this device has ambiguity and low measurement accuracy due to the presence of a large number of functional transformations.

Задачей, решаемой в изобретении, является расширение диапазона однозначно измеряемой радиальной скорости и повышение точности измерения за счет меньшего числа функциональных преобразований при применении совместной обработки неэквидистантных импульсных радиосигналов.The problem solved in the invention is to expand the range of unambiguously measured radial velocity and increase the accuracy of the measurement due to the smaller number of functional transformations when using joint processing of non-equidistant pulsed radio signals.

Для решения поставленной задачи в вычислитель радиальной скорости движущегося объекта, содержащий блок задержки, блок комплексного сопряжения, блок комплексного умножения, блок усреднения, блок вычисления фазы, блок коррекции пределов измерения, умножитель, ключ, первый блок памяти, блок вычисления модуля, блок управления, пороговый блок, второй блок памяти и синхрогенератор, введены первый и второй двухканальные ключи, дополнительный блок усреднения, дополнительный блок задержки, дополнительный блок комплексного сопряжения и дополнительный блок комплексного умножения.To solve the problem in the calculator of the radial speed of a moving object, containing a delay unit, complex conjugation unit, complex multiplication unit, averaging unit, phase calculation unit, measurement limits correction unit, multiplier, key, first memory unit, module calculation unit, control unit, the threshold block, the second memory block and the sync generator, the first and second two-channel keys, an additional averaging block, an additional delay block, an additional complex pairing block, and additional Yelnia complex multiplication unit.

Дополнительные блоки, введенные в предлагаемый вычислитель, являются известными. Так, соединенные вместе блок задержки, блок комплексного сопряжения и блок комплексного умножения позволяют выделить доплеровский набег фазы за интервал между соседними импульсами. Однако неизвестно совместное применение блока задержки, блока комплексного сопряжения, блока комплексного умножения, первого и второго двухканальных ключей, блока управления, дополнительного блока задержки, дополнительного блока комплексного сопряжения и дополнительного блока комплексного умножения. Новыми являются связи первого и второго двухканальных ключей с блоком комплексного умножения и блоком управления, блока усреднения с первым двухканальным ключом и дополнительным блоком задержки, дополнительного блока усреднения со вторым двухканальным ключом и дополнительным блоком комплексного сопряжения, дополнительного блока комплексного умножения с дополнительным блоком задержки и дополнительным блоком комплексного сопряжения, дополнительного блока комплексного умножения с блоком вычисления фазы, блоком коррекции пределов измерения и блоком вычисления модуля, что обеспечивает расширение диапазона однозначно измеряемой радиальной скорости и повышение точности измерения за счет меньшего числа функциональных преобразований при применении совместной обработки радиоимпульсных сигналов. Связи между синхрогенератором и всеми блоками вычислителя радиальной скорости движущегося объекта обеспечивают согласованную обработку неэквидистантной последовательности радиосигналов.Additional blocks entered in the proposed calculator are known. Thus, the delay unit, the complex conjugation unit, and the complex multiplication unit connected together allow one to isolate the Doppler phase incursion for the interval between adjacent pulses. However, the joint use of the delay unit, the complex conjugation unit, the complex multiplication unit, the first and second two-channel keys, the control unit, the additional delay unit, the additional complex conjugation unit and the additional complex multiplication unit is not known. The connections of the first and second two-channel keys with the complex multiplication unit and the control unit, the averaging unit with the first two-channel key and the additional delay unit, the additional averaging unit with the second two-channel key and the additional complex conjugation unit, the additional complex multiplication unit with the additional delay unit and the additional are new complex conjugation unit, an additional complex multiplication unit with a phase calculation unit, a limit correction unit, and measurement and a module calculation unit, which provides an extension of the range of unambiguously measured radial velocity and increased measurement accuracy due to fewer functional transformations when using joint processing of radio-pulse signals. The connections between the sync generator and all the blocks of the calculator of the radial speed of a moving object provide a coordinated processing of non-equidistant sequence of radio signals.

Сравнение с техническими характеристиками, известными из опубликованных источников информации, показывает, что заявляемое решение обладает новизной и имеет изобретательский уровень.Comparison with the technical characteristics known from published sources of information shows that the claimed solution has novelty and has an inventive step.

Заявляемое решение носит технический характер, осуществимо, воспроизводимо и, следовательно, является промышленно применимым.The claimed solution is technical in nature, feasible, reproducible and, therefore, is industrially applicable.

На фиг. 1 представлена структурная электрическая схема вычислителя радиальной скорости движущегося объекта; на фиг. 2 - блока задержки; на фиг. 3 - блока комплексного сопряжения; на фиг. 4 - блока комплексного умножения; на фиг. 5 - блока усреднения; на фиг. 6 - блока вычисления фазы; на фиг. 7 - блока коррекции пределов измерения; на фиг. 8 - блока присвоения знака; на фиг. 9 - блока вычисления модуля, на фиг. 10 - двухканального ключа; на фиг. 11 - блока управления.In FIG. 1 is a structural electrical diagram of a calculator of the radial velocity of a moving object; in FIG. 2 - delay unit; in FIG. 3 - block complex conjugation; in FIG. 4 - block complex multiplication; in FIG. 5 - averaging unit; in FIG. 6 - phase calculation unit; in FIG. 7 - block correction limits of measurement; in FIG. 8 - character assignment unit; in FIG. 9 - module calculation unit, in FIG. 10 - two-channel key; in FIG. 11 - control unit.

Вычислитель радиальной скорости движущегося объекта (фиг. 1) содержит блок 1 задержки, блок 2 комплексного сопряжения, блок 3 комплексного умножения, блок 4 усреднения, блок 5 вычисления фазы, блок 6 коррекции пределов измерения, умножитель 7, ключ 8, первый блок 9 памяти, блок 10 управления, блок 11 вычисления модуля, пороговый блок 12, второй блок 13 памяти, синхрогенератор 14, первый двухканальный ключ 15, второй двухканальный ключ 16, дополнительный блок 17 усреднения, дополнительный блок 18 задержки, дополнительный блок 19 комплексного сопряжения и дополнительный блок 20 комплексного умножения, при этом выходы блока 1 задержки соединены с входами блока 2 комплексного сопряжения, выходы которого соединены с первыми входами блока 3 комплексного умножения, вторые входы которого объединены с входами блока 1 задержки, выход первого блока 9 памяти соединен с первым входом умножителя 7, выход которого соединен с входом ключа 8, выход порогового блока 12 соединен с управляющим входом ключа 8, первый вход порогового блока 12 соединен с выходом второго блока 13 памяти, выходы блока 3 комплексного умножения соединены с объединенными входами первого 15 и второго 16 двухканальных ключей, управляющие входы которых соединены соответственно с первым и вторым выходами блока 10 управления, выходы первого двухканального ключа 15 соединены с входами блока 4 усреднения, выходы которого соединены с входами дополнительного блока 18 задержки, выходы второго двухканального ключа 16 соединены с входами дополнительного блока 17 усреднения, выходы которого соединены с входами дополнительного блока 19 комплексного сопряжения, выходы дополнительного блока 18 задержки соединены с первыми входами дополнительного блока 20 комплексного умножения, вторые входы которого соединены с выходами дополнительного блока 19 комплексного сопряжения, выходы блока 20 комплексного умножения соединены с объединенными входами блока 5 вычисления фазы, вторым и третьим входами блока 6 коррекции пределов измерения и блока 11 вычисления модуля, выход блока 11 вычисления модуля соединен со вторым входом порогового блока 12, выход блока 6 коррекции пределов измерения соединен со вторым входом умножителя 7, выход синхрогенератора 14 соединен с синхровходами блока 1 задержки, блока 2 комплексного сопряжения, блока 3 комплексного умножения, блока 4 усреднения, блока 5 вычисления фазы, блока 6 коррекции пределов измерения, умножителя 7, ключа 8, первого блока 9 памяти, блока 11 вычисления модуля, порогового блока 12, второго блока 13 памяти, первого 15 и второго 16 двухканальных ключей, дополнительного блока 17 усреднения, дополнительного блока 18 задержки, дополнительного блока 19 комплексного сопряжения и дополнительного блока 20 комплексного умножения, причем входами вычислителя радиальной скорости движущегося объекта являются входы блока 1 задержки, а первым и вторым выходами - соответственно выход ключа 8 и выход порогового блока 12.The radial velocity calculator of a moving object (Fig. 1) contains a delay unit 1, complex conjugation unit 2, complex multiplication unit 3, averaging unit 4, phase calculation unit 5, measurement limits correction unit 6, multiplier 7, key 8, first memory unit 9 , control unit 10, module calculation unit 11, threshold unit 12, second memory unit 13, clock generator 14, first two-channel key 15, second two-channel key 16, additional averaging unit 17, additional delay unit 18, additional complex pairing and add-on unit 19 an additional complex multiplication unit 20, while the outputs of the delay unit 1 are connected to the inputs of the complex conjugation unit 2, the outputs of which are connected to the first inputs of the complex multiplication unit 3, the second inputs of which are combined with the inputs of the delay unit 1, the output of the first memory unit 9 is connected to the first input a multiplier 7, the output of which is connected to the input of the key 8, the output of the threshold block 12 is connected to the control input of the key 8, the first input of the threshold block 12 is connected to the output of the second memory block 13, the outputs of the complex multiplication block 3 I am connected to the combined inputs of the first 15 and second 16 two-channel keys, the control inputs of which are connected respectively to the first and second outputs of the control unit 10, the outputs of the first two-channel key 15 are connected to the inputs of the averaging unit 4, the outputs of which are connected to the inputs of the additional delay unit 18, the outputs the second two-channel key 16 is connected to the inputs of the additional block 17 averaging, the outputs of which are connected to the inputs of the additional block 19 complex pairing, the outputs of the additional block 18 back The arms are connected to the first inputs of the additional complex multiplication unit 20, the second inputs of which are connected to the outputs of the additional complex conjugation unit 19, the outputs of the complex multiplication unit 20 are connected to the combined inputs of the phase calculation unit 5, the second and third inputs of the measurement limits correction unit 6, and the calculation unit 11 module, the output of block 11 calculating the module is connected to the second input of the threshold block 12, the output of block 6 correction of the limits of measurement is connected to the second input of the multiplier 7, the output of the clock Ora 14 is connected to the clock inputs of the delay unit 1, complex conjugation unit 2, complex multiplication unit 3, averaging unit 4, phase calculation unit 5, measurement limits correction unit 6, multiplier 7, key 8, first memory unit 9, module calculation unit 11, threshold block 12, second memory block 13, first 15 and second 16 two-channel keys, additional averaging block 17, additional delay unit 18, additional complex conjugation block 19 and additional complex multiplication block 20, the inputs being calculated ator radial velocity of the moving object 1 are input delay unit and first and second outputs - key respectively output 8 and an output of the threshold unit 12.

Блок 1 задержки и дополнительный блок 18 задержки (фиг. 2) содержат две цифровые линии задержки 21, входами блоков задержки являются входы цифровых линий задержки 21, выходы которых являются выходами блоков задержки.The delay unit 1 and the additional delay unit 18 (Fig. 2) contain two digital delay lines 21, the inputs of the delay units are the inputs of the digital delay lines 21, the outputs of which are the outputs of the delay units.

Блок 2 комплексного сопряжения и дополнительный блок 19 комплексного сопряжения (фиг. 3) содержат инвертор 22, первый вход блока комплексного сопряжения является его первым выходом, вторым входом является вход инвертора, выход которого является вторым выходом блока комплексного сопряжения.The complex conjugation unit 2 and the additional complex conjugation unit 19 (Fig. 3) comprise an inverter 22, the first input of the complex conjugation block is its first output, the second input is the inverter input, the output of which is the second output of the complex conjugation block.

Блок 3 комплексного умножения и дополнительный блок 20 комплексного умножения (фиг. 4) содержат два канала (I, II), каждый из которых включает первый перемножитель 23, последовательно включенные второй перемножитель 24 и сумматор 25, выход первого перемножителя 23 одного канала соединен со вторым входом сумматора 25 другого канала, а первыми и вторыми входами блока комплексного умножения соответственно являются объединенные между собой первые входы первого и второго перемножителей 23, 24 каждого из каналов, объединенные вторые входы первых перемножителей 23 и объединенные вторые входы вторых перемножителей 24, а выходами блока комплексного умножения являются выходы сумматоров 25 каждого из каналов.The complex multiplication block 3 and the additional complex multiplication block 20 (Fig. 4) contain two channels (I, II), each of which includes the first multiplier 23, the second multiplier 24 and the adder 25 connected in series, the output of the first multiplier 23 of one channel is connected to the second the input of the adder 25 of the other channel, and the first and second inputs of the complex multiplication block, respectively, are the combined first inputs of the first and second multipliers 23, 24 of each channel, the combined second inputs of the first Ithel 23 and the combined second inputs of the second multipliers 24 and the output unit are complex multiplication outputs of adders 25 of each of the channels.

Блок 4 усреднения и дополнительный блок 17 усреднения (фиг. 5) содержат два канала (I, II), каждый из которых состоит из (N-3)/2 последовательно включенных цифровых линий задержки 26 и (N-3)/2 последовательно включенных сумматоров 27, входами блока усреднения являются объединенные входы первой линии задержки 26 и первого сумматора 27 каждого канала (I, II), а выход k-й [k=1 … (N-3)/2] линии задержки 26 соединен со вторым входом k-го [k=1 … (N-3)/2)] сумматора 27 каждого канала (I, II), выходами блока усреднения служат выходы [(N-3)/2]-x сумматоров 27.Block 4 averaging and additional block 17 averaging (Fig. 5) contain two channels (I, II), each of which consists of (N-3) / 2 series-connected digital delay lines 26 and (N-3) / 2 series-connected adders 27, the inputs of the averaging block are the combined inputs of the first delay line 26 and the first adder 27 of each channel (I, II), and the output of the k-th [k = 1 ... (N-3) / 2] delay line 26 is connected to the second input of the k-th [k = 1 ... (N-3) / 2)] adder 27 of each channel (I, II), the outputs of the averaging block are the outputs of the [(N-3) / 2] -x adders 27.

Блок 5 вычисления фазы (фиг. 6) содержит последовательно соединенные делитель 28 и функциональный преобразователь 29, входами блока вычисления фазы являются входы делителя 28, а выходом блока вычисления фазы является выход функционального преобразователя 29.The phase calculation unit 5 (Fig. 6) contains a series divider 28 and a functional converter 29, the inputs of the phase calculator are the inputs of the divider 28, and the output of the phase converter 29 is the output of the functional converter 29.

Блок 6 коррекции пределов измерения (фиг. 7) содержит последовательно включенные модульный блок 30, сумматор 31, блок 32 присвоения знака, первый ключ 33 и сумматор 34, при этом первый вход блока 6 коррекции пределов измерения через второй ключ 35 соединен со вторым входом сумматора 34, выход блока 36 памяти соединен со вторым входом сумматора 31, второй вход блока коррекции пределов измерения соединен с управляющими входами первого 33 и второго 35 ключей, второй вход блока 32 присвоения знака является третьим входом блока коррекции пределов измерения, выход сумматора 34 является его выходом.Block 6 correction of the limits of measurement (Fig. 7) contains sequentially connected modular block 30, the adder 31, the block 32 of the assignment of the sign, the first key 33 and the adder 34, while the first input of the block 6 correction of the limits of measurement through the second key 35 is connected to the second input of the adder 34, the output of the memory unit 36 is connected to the second input of the adder 31, the second input of the measurement range correction unit is connected to the control inputs of the first 33 and second 35 keys, the second input of the character assignment unit 32 is the third input of the measurement range correction unit, course of the adder 34 is its output.

Блок 32 присвоения знака (фиг. 8) содержит блоки 37, 40 умножения, блок 38 памяти и ограничитель 39, причем второй вход блока присвоения знака является первым входом блока 37 умножения, второй вход которого соединен с выходом блока 38 памяти, выход блока 37 умножения соединен с входом ограничителя 39, выход которого соединен с первым входом блока 40 умножения, второй вход которого является первым входом блока присвоения знака, выходом блока присвоения знака служит выход блока 40 умножения.The character assignment unit 32 (FIG. 8) contains multiplication units 37, 40, a memory unit 38 and a limiter 39, the second input of the character assignment unit being the first input of the multiplication unit 37, the second input of which is connected to the output of the memory unit 38, the output of the multiplication unit 37 connected to the input of the limiter 39, the output of which is connected to the first input of the multiplication unit 40, the second input of which is the first input of the character assignment unit, the output of the character assignment unit is the output of the multiplication unit 40.

Блок 11 вычисления модуля (фиг. 9) содержат два блока 41 умножения, сумматор 42 и блок 43 извлечения квадратного корня, входами блока вычисления модуля являются входы блоков 41 умножения, выходы которых соединены с первым и вторым входами сумматора 42, выход которого соединен с входом блока 43 извлечения квадратного корня, выход которого является выходом блока вычисления модуля.The module calculation unit 11 (Fig. 9) contains two multiplication units 41, an adder 42 and a square root extraction unit 43, the inputs of the module calculation unit are the inputs of the multiplication units 41, the outputs of which are connected to the first and second inputs of the adder 42, the output of which is connected to the input a square root extraction unit 43, the output of which is the output of a module calculation unit.

Первый 15 и второй 16 двухканальные ключи (фиг. 10) содержат два ключа 44, входами двухканальных ключей являются входы ключей 44, выходы которых являются выходами двухканальных ключей.The first 15 and second 16 two-channel keys (Fig. 10) contain two keys 44, the inputs of the two-channel keys are the inputs of the keys 44, the outputs of which are the outputs of the two-channel keys.

Блок 10 управления (фиг. 11) содержит триггер 45 и элемент НЕ 46, входом блока управления является вход триггера 45, выход которого соединен с входом элемента НЕ 46, первым выходом блока 10 управления является выход триггера 45, а вторым выходом - выход элемента НЕ 46.The control unit 10 (Fig. 11) contains a trigger 45 and an element HE 46, the input of the control unit is the input of the trigger 45, the output of which is connected to the input of the element HE 46, the first output of the control unit 10 is the output of the trigger 45, and the second output is the output of the element NOT 46.

Вычислитель радиальной скорости движущегося объекта работает следующим образом.The calculator of the radial velocity of a moving object works as follows.

В заявляемом вычислителе обрабатывается неэквидистантная когерентно-импульсная последовательность N радиосигналов с чередующимися периодами повторения T1 и Т2, причем Τ12=ΔT. При отражении радиосигналов от движущейся цели их несущие частоты в соответствующих периодах приобретают доплеровские сдвиги фазыIn the inventive calculator, a non-equidistant coherent-pulse sequence of N radio signals with alternating repetition periods T 1 and T 2 is processed, with Τ 12 = ΔT. When radio signals are reflected from a moving target, their carrier frequencies in the corresponding periods acquire Doppler phase shifts

Figure 00000001
,
Figure 00000002
,
Figure 00000003
,
Figure 00000001
,
Figure 00000002
,
Figure 00000003
,

где

Figure 00000004
- доплеровская частота, νr - радиальная скорость цели,
Figure 00000005
- несущая частота радиосигналов, с - скорость распространения радиоволн.Where
Figure 00000004
is the Doppler frequency, ν r is the radial velocity of the target,
Figure 00000005
- carrier frequency of radio signals, s - propagation velocity of radio waves.

Отраженные от цели радиосигналы поступают на вход приемника, в котором усиливаются, в квадратурных фазовых детекторах переносятся на видеочастоту, а затем подвергаются аналого-цифровому преобразованию (соответствующие блоки на фиг. 1 не показаны). На вход вычислителя в одном элементе разрешения по дальности поступают цифровые отсчеты комплексной огибающейThe radio signals reflected from the target are fed to the input of the receiver, in which they are amplified, are transferred to the video frequency in quadrature phase detectors, and then undergo analog-to-digital conversion (the corresponding blocks in Fig. 1 are not shown). At the input of the calculator in one element of range resolution, digital readings of the complex envelope are received

Figure 00000006
, k=1 … N,
Figure 00000006
, k = 1 ... N,

где u1k, u2k - цифровые коды действительной и мнимой частей отсчетов Uk.where u 1k , u 2k are the digital codes of the real and imaginary parts of the samples U k .

Входные отсчеты Uk вычислителя (фиг. 1) в блоке 1 задержки (фиг. 2) под управлением синхронизирующих импульсов, вырабатываемых синхрогенератором 14, поочередно задерживаются на интервалы T1 и T2, что обеспечивает синхронность последующего комплексного умножения отсчетов по дальности. Синхрогенератор 14 управляется импульсами синхронизатора радиолокатора (на фиг. 1 не показан), следующими поочередно с интервалами T1 и Т2. В блоке 2 комплексного сопряжения (фиг. 3) осуществляется комплексное сопряжение задержанного отсчета (

Figure 00000007
). Далее в блоке 3 комплексного умножения (фиг. 4) реализуется попарное умножение отсчетов в соответствии с алгоритмомThe input samples U k of the calculator (Fig. 1) in the delay unit 1 (Fig. 2) under the control of the synchronizing pulses generated by the sync generator 14 are alternately delayed by the intervals T 1 and T 2 , which ensures synchronization of the subsequent complex multiplication of the samples in range. The sync generator 14 is controlled by the pulses of the radar synchronizer (not shown in Fig. 1), following alternately at intervals T 1 and T 2 . In block 2 complex conjugation (Fig. 3) is a complex conjugation of the delayed reference (
Figure 00000007
) Further, in block 3 of complex multiplication (Fig. 4), pairwise multiplication of samples is implemented in accordance with the algorithm

Figure 00000008
Figure 00000008

Попарные произведения (корреляции)

Figure 00000009
раздельно для каждого интервала T1 и Т2 соответственно через первый 15 и второй 16 двухканальные ключи раздельно поступают в блок 4 усреднения и в дополнительный блок 17 усреднения (фиг. 5). Поочередная коммутация первого 15 и второго 16 двухканального ключей осуществляется импульсами соответственно с первого и второго выходов блока 10 управления, синхронизируемого также импульсами синхронизатора радиолокатора.Pairwise products (correlations)
Figure 00000009
separately for each interval T 1 and T 2, respectively, through the first 15 and second 16 two-channel keys are separately received in block averaging 4 and in an additional block 17 averaging (Fig. 5). Alternate switching of the first 15 and second 16 two-channel keys is carried out by pulses, respectively, from the first and second outputs of the control unit 10, synchronized also by the pulses of the radar synchronizer.

В блоке 4 усреднения (фиг. 5) с помощью линий задержки 26 на интервал T12 и сумматоров 27 в каждом элементе разрешения по дальности осуществляется скользящее вдоль азимута когерентное суммирование (накопление) соответствующих интервалу T1 попарных произведений (корреляций). В результате на выходе блока 4 усреднения при нечетном N образуется пропорциональная корреляционному моменту отсчетов, соответствующих интервалу T1, величинаIn block 4 averaging (Fig. 5), using delay lines 26 for the interval T 1 + T 2 and adders 27 in each range resolution element, coherent summation (accumulation) of pairwise products (correlations) corresponding to the interval T 1 is performed along the azimuth. As a result, at the output of block averaging 4 for odd N, a value proportional to the correlation moment of samples corresponding to the interval T 1 is formed , the quantity

Figure 00000010
Figure 00000010

В дополнительном блоке 17 усреднения (фиг. 5) осуществляется аналогичное суммирование соответствующих интервалу Т2 попарных корреляций, что приводит к образованию на его выходе пропорциональной корреляционному моменту отсчетов, соответствующих интервалу Т2, величиныIn the additional averaging block 17 (Fig. 5), a similar summation of pairwise correlations corresponding to the T 2 interval is performed, which leads to the formation at its output proportional to the correlation moment of the samples corresponding to the T 2 interval,

Figure 00000011
:
Figure 00000011
:

Величина Y1 на выходе блока 4 усреднения (фиг. 5) по времени предшествует величине Y2 на интервал Т2, что компенсируется соответствующей данному интервалу задержкой Y1 в дополнительном блоке 18 задержки (фиг. 2). В дополнительном блоке 19 комплексного сопряжения (фиг. 3) инвертируется знак мнимой части величины Y2.The value of Y 1 at the output of the averaging unit 4 (FIG. 5) in time precedes the value of Y 2 by the interval T 2 , which is compensated by the corresponding delay Y 1 in the additional delay unit 18 (FIG. 2). In the additional block 19 complex conjugation (Fig. 3) the sign of the imaginary part of the value Y 2 is inverted.

Величины Y1 и

Figure 00000012
одновременно поступают соответственно на первые и вторые входы дополнительного блока 20 комплексного умножения (фиг. 4), на выходе которого вычисляется величинаValues Y 1 and
Figure 00000012
simultaneously arrive, respectively, at the first and second inputs of the additional complex multiplication unit 20 (Fig. 4), at the output of which the value is calculated

Figure 00000013
Figure 00000013

Величины ν1 и ν2 поступают на соответствующие входы блока 5 вычисления фазы (фиг. 6), где на основе блока 28 деления и арктангенсного функционального преобразователя 29 вычисляется оценкаThe values ν 1 and ν 2 are supplied to the corresponding inputs of the phase calculation unit 5 (Fig. 6), where, based on the division unit 28 and the arc tangent functional converter 29, an estimate is calculated

Figure 00000014
.
Figure 00000014
.

Последующие преобразования оценки

Figure 00000015
происходят в блоке 6 коррекции пределов измерения (фиг. 7) и зависят от знака ν1. При ν1>0 открыт второй ключ 35, и оценка
Figure 00000015
через сумматор 34 непосредственно поступает на выход блока коррекции пределов измерения. При ν1<0 открыт первый ключ 33, а второй ключ 35 закрыт. При этом в модульном блоке 30 образуется
Figure 00000016
, вычитаемый в блоке 31 из величины π, поступающей от блока 36 памяти. Полученной разности в блоке 32 присваивается знак величины ν2.Subsequent Assessment Conversions
Figure 00000015
occur in block 6 correction of the limits of measurement (Fig. 7) and depend on the sign of ν 1 . For ν 1 > 0, the second key 35 is open, and the estimate
Figure 00000015
through the adder 34 directly enters the output of the correction unit of the measurement limits. When ν 1 <0, the first key 33 is open, and the second key 35 is closed. Thus in the modular block 30 is formed
Figure 00000016
subtracted in block 31 from the value of π coming from the block 36 of the memory. The resulting difference in block 32 is assigned the sign of ν 2 .

Блок 32 присвоения знака (фиг. 8) работает следующим образом. На второй вход блока присвоения знака поступает величина ν2 [соотношение (1)], где в блоке 37 умножения производится ее умножение на постоянный множитель из блока 38 памяти с целью масштабирования и дальнейшего ограничения в ограничителе 39 по уровню ±1. Таким образом, после ограничения величина на выходе ограничителя 39 имеет смысл знака величины ν2, который, поступая на первый вход блока 40 умножения, присваивается разности

Figure 00000017
, поступающей с выхода блока 31 на первый вход блока 32 присвоения знака, т.е. на второй вход блока 40 умножения.Block 32 character assignment (Fig. 8) works as follows. The value ν 2 [relation (1)] is supplied to the second input of the sign-assigning unit, where in the multiplication unit 37 it is multiplied by a constant factor from the memory unit 38 in order to scale and further limit the limiter 39 to a level of ± 1. Thus, after the restriction, the value at the output of the limiter 39 has the meaning of the sign of the quantity ν 2 , which, arriving at the first input of the multiplication unit 40, is assigned the difference
Figure 00000017
coming from the output of block 31 to the first input of the character assignment block 32, i.e. to the second input of the multiplication block 40.

Рассмотренные операции позволяют в блоке 5 вычисления фазы сначала найти оценку доплеровского сдвига фазы, находящуюся в интервале [-π/2, π/2], а затем при помощи блока 6 коррекции пределов измерения расширить пределы ее однозначного измерения до интервала [-π, π] в соответствии с алгоритмомThe operations considered allow, in block 5 of the phase calculation, first to find an estimate of the Doppler phase shift located in the interval [-π / 2, π / 2], and then using block 6 to correct the measurement limits, expand the limits of its unambiguous measurement to the interval [-π, π ] according to the algorithm

Figure 00000018
Figure 00000018

Умножитель 7 (фиг. 1) осуществляет умножение найденной оценки сдвига фазы

Figure 00000019
на весовой коэффициент а, хранящийся в первом блоке 9 памяти, что позволяет найти однозначную оценку радиальной скорости в соответствии с выражениемThe multiplier 7 (Fig. 1) multiplies the found phase shift estimate
Figure 00000019
by the weight coefficient a stored in the first memory block 9, which allows one to find an unambiguous estimate of the radial velocity in accordance with the expression

Figure 00000020
Figure 00000020

где

Figure 00000021
- весовой коэффициент.Where
Figure 00000021
- weight coefficient.

Для уменьшения вероятности работы вычислителя по шумам в нем исключается выдача полученной оценки на выход в отсутствие отраженного от цели сигнала. В блоке 11 вычисления модуля (фиг. 9) вычисляется величинаTo reduce the likelihood of the computer working by noise, it excludes the issuance of the obtained estimate for output in the absence of a signal reflected from the target. In block 11, the calculation of the module (Fig. 9) calculates the value

Figure 00000022
Figure 00000022

которая поступает на второй вход порогового блока 12, в котором сравнивается с пороговым уровнем z0, записанным во втором блоке 13 памяти. Если происходит превышение порогового уровня z0, то с выхода порогового блока 12 поступает сигнал разрешения на прохождение результата вычисления с выхода блока 7 умножения через ключ 8 на первый выход вычислителя радиальной скорости движущегося объекта. В противном случае ключ 8 разомкнут. Кроме того, сигнал с выхода порогового блока 12, являющегося вторым выходом вычислителя радиальной скорости движущегося объекта, может быть использован для отсчета других координат цели, например дальности.which is supplied to the second input of the threshold block 12, which compares with the threshold level z 0 recorded in the second memory block 13. If the threshold level z 0 is exceeded, then the output of the threshold block 12 receives a permission signal for passing the calculation result from the output of the multiplication block 7 through the key 8 to the first output of the calculator of the radial velocity of a moving object. Otherwise, key 8 is open. In addition, the signal from the output of the threshold block 12, which is the second output of the calculator of the radial speed of a moving object, can be used to read other coordinates of the target, for example range.

Синхронизация вычислителя радиальной скорости движущегося объекта осуществляется подачей на все блоки заявляемого вычислителя последовательности синхронизирующих импульсов, вырабатываемых синхронизатором 14 (фиг. 1) с периодом повторения tк, определяемым из условия требуемой разрешающей способности по дальности.The synchronization of the calculator of the radial speed of a moving object is carried out by applying to all blocks of the inventive calculator a sequence of synchronizing pulses generated by the synchronizer 14 (Fig. 1) with a repetition period t k determined from the condition of the required range resolution.

Достижение технического результата объясняется следующим образом. В известном устройстве (прототипе) исходные доплеровские сдвиги фазы ϕ1 и ϕ2, по которым вычисляется величина Δϕ=ϕ12, имеют интервал однозначного измерения [-π, π], что соответствует интервалу однозначного измерения доплеровской частоты [-1/2T1], 1/2T1]] (по величине большего периода Т1). В предложенном устройстве величина Δϕ измеряется непосредственно, что соответствует интервалу однозначности доплеровских частот [-1/2ΔT, 1/2ΔT]. При этом интервал однозначного измерения доплеровской частоты и, следовательно, радиальной скорости расширяется в T1/ΔT раз, что соответствует решению поставленной задачи изобретения. Если в соответствии с условием

Figure 00000023
и с учетом
Figure 00000024
для максимально возможной скорости движущегося объекта νr max выбрать интервал
Figure 00000025
, то во всем диапазоне реальных скоростей движущегося объекта может быть осуществлено их однозначное измерение. При этом сохраняется однозначность измерения дальности, которая обеспечивается соответствующим выбором меньшего периода повторения импульсов Т2.The achievement of the technical result is explained as follows. In the known device (prototype), the initial Doppler phase shifts ϕ 1 and ϕ 2 , by which the value Δϕ = ϕ 12 is calculated, have a single-valued measurement interval [-π, π], which corresponds to a single-valued measurement Doppler frequency interval [-1 / 2T 1 ], 1 / 2T 1 ]] (by the value of the larger period T 1 ). In the proposed device, the value of Δϕ is measured directly, which corresponds to the interval of uniqueness of Doppler frequencies [-1 / 2ΔT, 1 / 2ΔT]. Moreover, the interval of unambiguous measurement of the Doppler frequency and, consequently, the radial velocity expands by T 1 / ΔT times, which corresponds to the solution of the problem of the invention. If in accordance with the condition
Figure 00000023
and given
Figure 00000024
for the maximum possible speed of a moving object ν r max choose the interval
Figure 00000025
, then in the entire range of real speeds of a moving object, their unambiguous measurement can be carried out. This preserves the uniqueness of the range measurement, which is ensured by the appropriate choice of a shorter pulse repetition period T 2 .

Обусловленные функциональными преобразованиями погрешности раздельного вычисления величин ϕ1 и ϕ2 являются статистически независимыми. В результате погрешность (дисперсия) разности ϕ12=Δϕ удваивается. В предложенном вычислителе при непосредственном вычислении величины Δϕ такое удвоение отсутствует, что соответствует повышению точности измерения.The errors of the separate calculation of ϕ 1 and ϕ 2 due to functional transformations are statistically independent. As a result, the error (dispersion) of the difference ϕ 12 = Δϕ is doubled. In the proposed calculator, when directly calculating Δϕ, such a doubling is absent, which corresponds to an increase in the measurement accuracy.

Таким образом, вычислитель радиальной скорости движущегося объекта позволяет расширить диапазон однозначно измеряемой радиальной скорости и повысить точность измерения за счет меньшего числа функциональных преобразований при применении предлагаемой совместной обработки неэквидистантных когерентно-импульсных радиосигналов.Thus, the calculator of the radial velocity of a moving object allows you to expand the range of unambiguously measured radial velocity and improve the accuracy of the measurement due to the smaller number of functional transformations when using the proposed joint processing of non-equidistant coherent-pulse radio signals.

БиблиографияBibliography

1. Ширман Я.Д. и Манжос В.Н. Теория и техника обработки радиолокационной информации на фоне помех. - М.: Радио и связь. - 1981. - С. 204. - Рис. 14.2.1. Shirman Y.D. and Manzhos V.N. The theory and technique of processing radar information against the background of interference. - M .: Radio and communication. - 1981. - S. 204. - Fig. 14.2.

2. Патент №63-49193 (Япония), МПК G01S 13/52. Радиолокационное устройство для обнаружения движущейся цели / К.К. Тосиба. Опубл. 03.10.1988. - Изобретения стран мира. - 1989. - Выпуск 109. - №15. - С. 52.2. Patent No. 63-49193 (Japan), IPC G01S 13/52. Radar device for detecting a moving target / K.K. Toshiba. Publ. 10/03/1988. - Inventions of the countries of the world. - 1989. - Issue 109. - No. 15. - S. 52.

3. Патент №2017167 (Россия), МПК G01S 13/58. Обнаружитель-измеритель доплеровских сигналов / Д.И. Попов, C.B. Герасимов и Е.Н. Матаев. Опубл. 30.07.1994. - Изобретения. - 1994. - №14. - С. 121.3. Patent No. 2017167 (Russia), IPC G01S 13/58. Detector-meter of Doppler signals / D.I. Popov, C. B. Gerasimov and E.N. Mataev. Publ. 07/30/1994. - Inventions. - 1994. - No. 14. - S. 121.

Claims (1)

Вычислитель радиальной скорости движущегося объекта, содержащий блок задержки, блок комплексного сопряжения, блок комплексного умножения, блок усреднения, блок вычисления фазы, блок коррекции пределов измерения, умножитель, ключ, первый блок памяти, блок управления, блок вычисления модуля, пороговый блок, второй блок памяти и синхрогенератор, при этом выходы блока задержки соединены с входами блока комплексного сопряжения, выходы которого соединены с первыми входами блока комплексного умножения, вторые входы которого объединены с входами блока задержки, выход первого блока памяти соединен с первым входом умножителя, выход которого соединен с входом ключа, выход порогового блока соединен с управляющим входом ключа, первый вход порогового блока соединен с выходом второго блока памяти, выход синхрогенератора соединен с синхровходами блока задержки, блока комплексного сопряжения, блока комплексного умножения, блока усреднения, блока вычисления фазы, умножителя, ключа, блока вычисления модуля, первого и второго блоков памяти, блока коррекции пределов измерения и порогового блока, отличающийся тем, что введены первый и второй двухканальные ключи, дополнительный блок усреднения, дополнительный блок задержки, дополнительный блок комплексного сопряжения и дополнительный блок комплексного умножения, при этом выходы блока комплексного умножения соединены с объединенными входами первого и второго двухканальных ключей, управляющие входы которых соединены соответственно с первым и вторым выходами блока управления, выходы первого двухканального ключа соединены с входами блока усреднения, выходы которого соединены с входами дополнительного блока задержки, выходы второго двухканального ключа соединены с входами дополнительного блока усреднения, выходы которого соединены с входами дополнительного блока комплексного сопряжения, выходы дополнительного блока задержки соединены с первыми входами дополнительного блока комплексного умножения, вторые входы которого соединены с выходами дополнительного блока комплексного сопряжения, выходы дополнительного блока комплексного умножения соединены с входами блока вычисления фазы, блока вычисления модуля, вторым и третьим входами блока коррекции пределов измерения, выход которого соединен со вторым входом умножителя, выход блока вычисления модуля соединен со вторым входом порогового блока, выход синхрогенератора соединен с синхровходами первого и второго двухканальных ключей, дополнительного блока усреднения, дополнительного блока задержки, дополнительного блока комплексного сопряжения и дополнительного блока комплексного умножения, причем входами вычислителя радиальной скорости движущегося объекта являются входы блока задержки, а первым и вторым выходами - соответственно выходы ключа и порогового блока.A radial velocity calculator of a moving object, comprising a delay unit, a complex conjugation unit, a complex multiplication unit, an averaging unit, a phase calculation unit, a measurement range correction unit, a multiplier, a key, a first memory unit, a control unit, a module calculation unit, a threshold unit, a second unit memory and a clock generator, while the outputs of the delay unit are connected to the inputs of the complex conjugation unit, the outputs of which are connected to the first inputs of the complex multiplication unit, the second inputs of which are combined with the input by the delay block, the output of the first memory block is connected to the first input of the multiplier, the output of which is connected to the key input, the output of the threshold block is connected to the control input of the key, the first input of the threshold block is connected to the output of the second memory block, the output of the clock generator is connected to the sync inputs of the delay block, block complex conjugation, complex multiplication block, averaging block, phase calculation block, multiplier, key, module calculation block, first and second memory blocks, correction unit for measurement limits and threshold th block, characterized in that the first and second two-channel keys, an additional averaging block, an additional delay block, an additional complex conjugation block and an additional complex multiplication block are introduced, while the outputs of the complex multiplication block are connected to the combined inputs of the first and second two-channel keys, control inputs which are connected respectively to the first and second outputs of the control unit, the outputs of the first two-channel key are connected to the inputs of the averaging unit, the outputs of which are dined with the inputs of the additional delay unit, the outputs of the second two-channel key are connected to the inputs of the additional averaging unit, the outputs of which are connected to the inputs of the additional complex conjugation unit, the outputs of the additional delay unit are connected to the first inputs of the additional complex multiplication unit, the second inputs of which are connected to the outputs of the additional complex unit interfaces, outputs of the additional complex multiplication block are connected to the inputs of the phase calculation unit, the calculation unit module, the second and third inputs of the block for the correction of the limits of measurement, the output of which is connected to the second input of the multiplier, the output of the unit for calculating the module is connected to the second input of the threshold block, the output of the clock is connected to the clock inputs of the first and second two-channel keys, an additional averaging block, an additional delay block an additional block of complex conjugation and an additional block of complex multiplication, and the inputs of the block are the inputs of the calculator of the radial speed of a moving object delay, and first and second output -, respectively, and outputs a key block threshold.
RU2015150365A 2015-11-24 2015-11-24 Calculator of range rate of moving object RU2613037C1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2015150365A RU2613037C1 (en) 2015-11-24 2015-11-24 Calculator of range rate of moving object

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2015150365A RU2613037C1 (en) 2015-11-24 2015-11-24 Calculator of range rate of moving object

Publications (1)

Publication Number Publication Date
RU2613037C1 true RU2613037C1 (en) 2017-03-14

Family

ID=58458291

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2015150365A RU2613037C1 (en) 2015-11-24 2015-11-24 Calculator of range rate of moving object

Country Status (1)

Country Link
RU (1) RU2613037C1 (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4825213A (en) * 1981-05-15 1989-04-25 Grumman Aerospace Corporation Simultaneous triple aperture radar
RU2017167C1 (en) * 1991-07-01 1994-07-30 Рязанский Радиотехнический Институт Pulse doppler locator
US5646623A (en) * 1978-05-15 1997-07-08 Walters; Glenn A. Coherent, frequency multiplexed radar
US5784026A (en) * 1996-09-23 1998-07-21 Raytheon E-Systems, Inc. Radar detection of accelerating airborne targets
RU154313U1 (en) * 2014-12-08 2015-08-20 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Рязанский государственный радиотехнический университет" MOVING OBJECT SPEED CALCULATOR

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5646623A (en) * 1978-05-15 1997-07-08 Walters; Glenn A. Coherent, frequency multiplexed radar
US4825213A (en) * 1981-05-15 1989-04-25 Grumman Aerospace Corporation Simultaneous triple aperture radar
RU2017167C1 (en) * 1991-07-01 1994-07-30 Рязанский Радиотехнический Институт Pulse doppler locator
US5784026A (en) * 1996-09-23 1998-07-21 Raytheon E-Systems, Inc. Radar detection of accelerating airborne targets
RU154313U1 (en) * 2014-12-08 2015-08-20 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Рязанский государственный радиотехнический университет" MOVING OBJECT SPEED CALCULATOR

Similar Documents

Publication Publication Date Title
RU2507536C1 (en) Coherent pulsed signal measuring detector
RU2642418C1 (en) Interference reject filter
RU2560130C1 (en) Pulsed radio signal detection and measurement device
RU2634190C1 (en) Interference rejecting counter
RU2674468C1 (en) Interference rejection filter
RU154313U1 (en) MOVING OBJECT SPEED CALCULATOR
RU2582877C1 (en) Adaptive compensator of passive interference phase
RU157108U1 (en) PASSIVE INTERFERENCE PHASE COMPENSATION DEVICE
CN104062511A (en) Accurate phase measurement method of multipath multi-target echo signal
RU2629642C1 (en) Doppler speed calculator of object movement
RU149732U1 (en) PHASOMETER OF COHERENT RADIO PULSES
RU2613037C1 (en) Calculator of range rate of moving object
RU161949U1 (en) COMPUTER FOR AUTO COMPENSATION OF SHIFT PHASE SHIFTS
RU161877U1 (en) COAGER RADIO SIGNAL DETECTOR-METER
RU161582U1 (en) MOBILE OBJECT SPEED CALCULATOR
RU2600111C1 (en) Detecting and measuring device of coherent pulsed radio signals
RU2546988C1 (en) Pulsed radio signal detector-measuring device
RU166743U1 (en) OBJECT RADIAL SPEED COMPUTER
RU146461U1 (en) RADIO-PULSE SIGNAL PROCESSING DEVICE
RU2513656C2 (en) Phase meter of coherent-pulse signals
RU2547159C1 (en) Phase indicator of radio pulse signals
RU2660803C1 (en) Filter of noise notching
RU150201U1 (en) RADIAL SPEED MEASURER
RU2629710C1 (en) Phase meter of coherent non-equidistant pulses
RU2609438C1 (en) Coherent pulsed radio signals phase meter

Legal Events

Date Code Title Description
MM4A The patent is invalid due to non-payment of fees

Effective date: 20171125