Способ дискриминаторного измерения задержки отраженного импульса

 

Изобретение относится к радиолокационной технике и может быть использовано для обнаружения и измерения расстояний до разного рода подвижных и неподвижных объектов при реализации в РЛС зондирующего составного ЛЧМ импульса с псевдослучайными нвчальными частотами и фазами. Способ может быть использован в радиолокации и гидролокации, в тех областях и сферах деятельности, где необходимо измерять расстояния до недоступных объектов. Временной дискриминатор - это устройство, вычисляющее разность между измеряемой и ожидаемой величиной задержки отраженного сигнала.

Достигаемый технический результатсостоит в том, что предварительно осуществляется измерение реальной нелинейной временной дискриминаторной характеристики (ДХ) составного ЛЧМ фильтра сжатия, ее квантование с заданным шагом в пределах одного такта дискретизации входного сигнала, вычисление тангенсов наклона реальной ДХ в точках квантования, формирование линейной эталонной ДХ и запоминание этих характеристик и значений тангенсов; затем в процессе штатной работы осуществляется измерение по каждой цели реальной дискриминаторной оценки задержки (ДОЗ), вычисление квантованного значения этой ДОЗ, проводится ассоциативный выбор по реальной ДОЗ соответствующего эталонного значения ДОЗ, его уточнение с использованием значения тангенса наклона ДХ и расчет по уточненному значению ДОЗ искомого абсолютного значения задержки; при этом обеспечивается вычисление ДОЗ, не содержащей нелинейных ошибок, ДХ становится линейной и обеспечивается вычисление полного значения измеряемой задержки, а также пикового значения амплитуды сигнала с повышенной методической точностью.

Известны способы измерения дальности до цели, использующие дискриминаторные методы измерения времени прихода отраженного от цели сигнала (см. Л1. Радиоавтоматика: Учеб. пособие для студ. вузов спец. «Радиотехника» / В.А. Бесекерский, А.А. Елисеев, А.В. Небылов и др.; Под ред. В.А. Бесекерского. - М.: Высш. шк., 1985. - 271 с.). Основным элементом такого дискриминатора является преобразователь временного интервала между отраженным от цели импульсом и зондирующим импульсом в двоичный код, который используется для вычисления дальности до цели и замыкания петли обратной связи в системе автоматического сопровождения по дальности (АСД). Такой дискриминатор может успешно работать лишь в следящей системе АСД.

Известны способы и устройстваслежения за временной задержкой сигнала (см. Л2._СПОСОБ СЛЕЖЕНИЯ ЗА ВРЕМЕННОЙ ЗАДЕРЖКОЙ СИГНАЛА И УСТРОЙСТВО ДЛЯ ЕГО РЕАЛИЗАЦИИ, Патент РФ №2157593, H04L 7/02, Н04В 1/10, опубл. 10.10.2000, авторы Гармонов А.В.; Манелис В.Б.; Савинков А.Ю.). Данное изобретение используется в цифровых системах связи. Временное слежение в системах цифровой радиосвязи строится на основе формирования оценки временной ошибки синхронизации и последующей ее коррекции. Для формирования оценки временной ошибки обычно используется временной дискриминатор, оценивающий производную корреляционной функции полезного сигнала в заданные моменты времени. Такой дискриминаторможет использоваться только в системах слежения за задержкой с обратной связью.

Известны дискриминаторы, вычисляющие временное смещение измеряемого сигнала относительно ожидаемого времени по двум или трем его амплитудным отсчетам в дискретные моменты времени (см. Л3. "Применение цифровой обработки сигналов" под ред. Э. Оппенгейма, М, 1980 г., см. Л4. "Цифровой дискриминатор задержки", патент РФ №2059254, G01R 23/02, опубл. 27.04.1996, автор Наумов С.В.).

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

При реализации сложного составного ЛЧМ сигнала с псевдослучайными начальными частотами и начальными фазами дискриминаторная характеристика (ДХ), имеет нелинейный характер, вид которого изображен на фиг. 1 и 2. Этот вид ДХ является следствием неопределенной зависимости изменения амплитуд сжатого импульса во временных отсчетах (тактах дискретизации) при смещении входного сигнала в пределах одного такта дискретизации. При этом следует отметить два существенных фактора:

- максимальное пиковое значение амплитуды сжатого импульса всегда располагается в такте дискретизации, номер которого соответствует целочисленному значению кода задержки Nз, выраженному в числе тактов дискретизации сигнала на входе фильтра;

- для составного ЛЧМ сигнала с неизменными основными параметрами: длительностью ЗИ, девиацией частоты, частотой дискретизации и числом элементарных импульсов в ЗИ, единожды измеренная ДХ в силу цифрового характера аппаратуры обработки сохраняет свою форму неизменной на протяжении всего времени эксплуатации, причем эта форма остается неизменной даже при ограничении длительности ОИ, которое вызывается уменьшением дальности до цели (см. фиг. 1 и 2).

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

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

Поставленная цель достигается способом ассоциативного выбора по реально измеренному значению дискриминаторной оценки задержки (ДОЗ) соответствующего значения ДОЗ из эталонной ДХи уточнение этого значения путем дополнительных вычислений.

Для этого в способе-прототипе, предназначенном для оценки положения пика отраженного сигнала по дальности в режимах обнаружения и сопровождения и заключающемся в вычислении разности модулей двух амплитуд и делении ее на сумму этих модулей, для достижения поставленной цели согласно настоящему изобретениюреализуется предварительный этап, в процессе которого вычисляется реальная временная дискриминаторная характеристика (ДХ) составного ЛЧМ фильтра сжатия путем первоначального задания целочисленной задержки N_з отраженного импульса и последующего изменения задержки с малым шагомр=1/P квантования задержки от нуля до единицы интерваладискретизации сигнала Т_д, с вычислением безразмерных дискриминаторных оценок задержки (ДОЗ) dT₀ пo формуле:

где p - номер точки (шага) квантования на интервале дискретизации Т_д,m - номер отсчета локального максимума (ЛМ) в амплитудном рельефе сжатого сигнала на выходе фильтра сжатияот начала приема, причем (см. фиг. 4) знак плюс перед dT_0p ставится тогда, когда соседняя с максимумом амплитуда слева (A_m-1) больше соседней амплитуды справа (А_m+1), иначе - ставится знак минус; кроме того вычисляются тангенсы углов наклона ДХ tgα_p в точках квантования интервала дискретизации по формуле:

измеренные заранее значения реальной временной ДХ преобразуются в порядковые номера (адреса), по которым в постоянное или оперативное запоминающее устройство (ЗУ) заносятся и там хранятся реальные значения dT_0p, тангенсы углов наклона ДXtgα_p и точные (эталонные) значения dt_ЭТ ДХ; причем эталонная ДХ представляет собой прямую, идущую под углом 45° от значения -0,5Т_Д к значению +0,5Т_Д.

В процессе штатной работы в режиме сопровождения (см. фиг. 5) по формуле (1) и по вновь измеренным значениям амплитуд А_m, А_m±1, A_m±1 вычисляется ДОЗdT_n, от которой выбирается ближайшее слева квантованное значение dT_nкв, где n - номер шага квантования. Квантованное значение этой ДОЗdT_nкв является адресом, по которому из ЗУ выбирается соответствующее истинное значение ДОЗdt; при этом в режиме сопровождения проводится вычисление уточненного значения ДОЗ dt_уточн с использованием заданной дробной части Δ_C заданной задержки N_з по формулам:

где dN_з - остаток заданной задержки; dN_з.кв - квантованное значение остатка заданной задержки, выраженное в долях интервала дискретизации T_Д; значение Δ_C может поступать в составе целеуказания.

В режиме обнаружения (см. фиг. 6)проводится уточнение значения ДОЗ с использованием остатка измеряемой задержки Δс учетомтангенса угла наклона ДХ по формулам:

при этом измерение абсолютного значения задержки N_Зизм в режимах обнаружения и сопровождения осуществляется путем экстраполяции полного измеренного значения задержки m уточненным значением ДОЗ dt_уточн по формуле:

где m - номер отсчета, в котором находится амплитуда ЛМ сжатого ОИ.

Исходя из линейной зависимости амплитуды сжатого импульса в пределах одного такта дискретизации (см фиг. 3) пиковое значение амплитуды вычисляется по формуле:

где А_лм - амплитуда ЛМ, А_пик - искомое значение амплитуды пика.

При этом обеспечивается вычисление дискриминаторной оценки задержки, не содержащей нелинейных ошибок, ДХ становится линейной и обеспечивается вычисление полного значения измеряемой задержки и пиковой амплитуды сигнала с повышенной методической точностью;

Для наглядного представления вида временных дискриминаторных характеристик на фиг. 1 и 2 показаны реальные и идеальные (эталонные) ДХ составного ЛЧМ фильтра сжатия с указанными параметрами отраженного импульса с полной длительностью Т_и=40 мс и при таком же зондирующем импульсе отраженного импульса с укороченной длительностью Ти=4 мс из-за задержки подальности. Видно, что при обеих длительностях сигналов их дискриминаторные характеристики идентичны.

На фиг. 3 показан характер изменения амплитуды главного лепестка (ГЛ) сжатого импульса в пределах одного такта дискретизации. Видно, что ГЛ имеет вид равнобедренного треугольника. Пик сигнала совпадает с целочисленным номером такта дискретизации Т_д (номер смещения n=0).

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

На фиг. 5 изображен график, поясняющий процедуру вычисления дискриминаторной характеристики и тангенсов углов наклона на предварительном этапе и вычисления оценки задержки dt_уточн, в режиме сопровождения, на котором толстыми линиями со стрелками условно отображается порядок действий при измерении и вычислении времени задержки ОИ. По горизонтальной оси отложены значения заданных задержек Тз за вычетом целой части (Nз - ц.ч. Nз)Т_д, т.е. отложены дробные части задержки с шагом квантования, например, 0,1Т_Д, а именно - (0; 0,1; 0,2; … 0,9; 1)Т_д. По вертикальной оси отложены дискриминаторные оценки измеряемой задержки в пределах от -0,5Т_Д до +0,5Т_Д.

На фиг. 6 изображен аналогичный график, поясняющий процедуру вычисления дискриминаторной оценки задержки dt_уточн, в режиме обнаружения.

На фиг. 7 в качестве примера одного из возможных вариантов реализации этого способа приведенаупрощенная блок-схема дискриминатора.

На фиг. 7 обозначено:

1 - коммутатор амплитуд;

2 - блок сравнения амплитуд;

3 - вычитатель амплитуд;

4 - сумматор амплитуд;

5 - делитель;

6 - квантователь;

7 - шифратор;

8 - запоминающее устройство (ЗУ);

9 - вычитатель;

10 - делитель;

11 - сумматор;

12 - делитель;

13 - сумматор-вычитатель;

14 - вычитатель;

15 - умножитель.

Работает дискриминатор следующим образом.

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

Итак, заранеепутем моделирования на отдельном компьютере или на данном устройстве реализуется этап предварительного вычисления реальной дискриминаторной характеристики. Этот этап заключается в том, что на вход фильтра сжатия подается отраженный импульс с задержкой Nз, которая сначала задается целочисленной величиной, а затем изменяется с малым шагом квантования 1/Р такта дискретизации Т_д, например, с шагом 1/10, как показано на горизонтальной оси фиг. 5. После каждого цикла задания задержки с выхода фильтра сжатия снимаются отсчеты амплитуды локального максимума (ЛМ) Аm и соседних амплитуд Аm-1, Аm+1. В блоке сравнения 2, осуществляется сравнение амплитуд Аm-1 и Аm+1, и из них определяется наибольшая амплитуда, которая через коммутатор 1 по команде от блока 2 вместе с амплитудой ЛМ Am передается в вычитатель3 и в сумматор 4, с выходов которых коды поступают в блок 5 деления. Результату деления, вычисленному по формуле (1), приписывается знак плюс или минус, который вырабатывается в блоке сравнения, как поясняется на фиг. 4, по принципу: если Аm-1>Аm+1, то знак плюс, если Аm-1≤Аm+1, то знак минус. Результат деления - безразмерная дискриминаторная оценка задержки dТ₀реальной ДХ, заноситсяв ЗУ8в соответствии с квантованным в блоке 6 значением p(N_з - ц.ч. N_з)/P (р=0÷1), и сформированным шифратором 7 кодом адреса.

Примечание. Для облегчения визуального восприятия алгоритма действий при вычислениях ДОЗ на фиг. 5 и бреальная ДХ преднамеренно изображена с большими нелинейностями.

Кроме того, одновременно вычисляются значения тангенсов tgα_p наклона измеряемой ДХ путем вычисления в вычитателе 9 разностей значений dT_0p-dT_0p-1 и деления их на шаг квантования 1/Р в делителе 10. Полученные таким образом значения ДОЗ и тангенсов наклона ДХ заносятся в ЗУ. Их значения являются константами, так как зависят только от основных неизменяемых параметров зондирующего сигнала и могут использоваться, например, при контроле устройства. Рядом с этими данными хранятся рассчитанные заранее точные (эталонные) значения dt_ЭT ДХ. Эталонная ДХ представляет собой прямую, идущую под углом 45° от значения -0,5 к значению +0,5. Записанные в ЗУ 8реальнаяДХ, тангенсы наклона и эталонная ДХ хранятся постоянно. После записи их в ЗУ 8 подготовительный этап заканчивается.

В штатном режиме сопровождения для обеспечения успешного функционирования временного дискриминатора фильтр сжатия должен работать по целеуказаниям, которые выражаются значениями ожидаемых задержек Nз отраженных импульсов в тактах дискретизации Т_д, с которыми от фильтра сжатия на вход дискриминаторапоступают амплитуда ЛМ Am, и соседние амплитуды Am-1, Аm+1. По формуле (1) на выходе делителя 5 сформируетсязначение ДОЗdT_n. А в квантователе 6 сформируется с шагом 1/Р значение ДОЗdT_nкв=pdT_n/P. Это значение преобразуется шифратором 7 в адрес, по которому из таблицы эталонных значений dt_ЭТ, хранящейся в ЗУ 8, выбирается близкое кистинному значение dt ДОЗ, неимеющее нелинейных ошибок. На фиг. 5 величины dT_n и dT_nкв отмечены толстыми вертикальными стрелками с координатами по оси абсцисс ~0,87 и p/Р=0,8 соответственно. Для повышения точности проводится линейная экстраполяция измеренного значения dt ДОЗ с использованием дробной части Δ_с, вычисляемой по формуле (4) и поступающей вместе с целеуказанием в блок 11, после чего в нем вычисляется dt_уточн по формуле (3).

Примечание. Выражение (3) имеет место в силу равнобедренности прямоугольного треугольника с катетами [dt_уточн-dt] и [0,87-0,8].

Затем в блоке 13, куда от блока 11 передается значение dt_уточн, вычисляется измеренное абсолютное значение задержки Nз изм по формуле (8), при этом использование уточненного значения ДОЗ, не содержащего нелинейных ошибок, обеспечивает вычисление полного значения задержки ОИ с повышенной точностью. Значение Nз изм поступает на выход устройства.

В режиме обнаружения временной дискриминатор функционирует без предварительного целеуказания. Используются только текущие измеренные значения амплитуд Am, Аm-1, Аm+1 и код m временного положения ЛМ. Как и в сопровождении, вычисляются разность и сумма амплитуд, их отношение dTn и dTnкв. По шифрированному коду dTnкв в блоке 7 вырабатывается адрес, по которому считываются значение ДОЗ dt и тангенс tgα, которые поступают соответственно в блоки 11 и 12, В блоке 12 производится вычисление величины Δ по формуле (7), а в блоке 11 - вычисление dt_уточн. Затем в блоке 13, куда от блока 11 передается значение dt_уточн, вычисляется измеренное полное значение задержки Nз изм по формуле (8), при этом использование уточненного значения ДОЗ, не содержащего нелинейных ошибок, обеспечивает вычисление полного значения задержки ОИ с повышенной методической точностью.

Значение Nз изм поступает на выход устройства.

Пиковое значение амплитуды сжатого ОИ достаточно просто рассчитывается по измеренному значению ДОЗ dt_уточнн и измеренной амплитуде ЛМ, так как форма ГЛ сжатого сигнала не зависит ни от длительности ОИ, ни от значения доплеровского смещения, а также независит от случайныхначальных частот и фаз. На фиг. 3 видна линейная зависимость амплитуды сигнала от временного смещения в пределах одного такта дискретизации от -0,5 до 0 и от 0 до +0,5 Нормированный треугольник на фиг. 3 имеет пиковое значение в вершине А_пик=1, а в основании - Аосн=0,8.

Из подобия треугольников (А_пик_0_Аосн) и (А_лм_(0,5-dt_уточнн,)_Аосн) имеем: А_лм/(0,5-dt_уточнн)=А_ПИК/0,5, где А_лм - измеренная амплитуда ЛМ. Из этого выражения вытекает формула (9) для определения амплитуды пика.

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

При моделировании были взяты следующие исходные параметры:

- длительность составного ЗИ и интервала приема отраженного импульса Т_ЗИ=Тпр=40 мс, период повторения ЗИ Т_п=2Т_ЗИ=80 мс;

- значение полной девиации частоты ЗИ F_дев=0,5 МГц;

- частота дискретизации ОИ и частота следования комплексных сигналов на выходе приемника F_д=1 МГц (такт Тд=1 мкс);

- число элементарных ЛЧМ импульсов (ЭИ) в ЗИN=2000;

- длительность ЭИ T_ЭИ=T_ЗИ/N, Т_эи=20 мкс;

- значение девиации частоты ЭИ F_девЭИ=300 кГц;

- число каналов дальности N_ЗИ=T_зиF_д=40000.

За счет девиации частоты коэффициент сжатия сигнала (база сигнала) составила значение Б=Т_зиF_дев=40000.

В согласованном фильтре отраженный импульс с длительностью Т_ОИ=40 мс сжимается до длительности Т_сж=2Т_зИ/Б=2 мкс, что может являться оценкой разрешающей способности по дальности. При моделировании задавались разные значения задержки и доплеровской частоты сигнала на входе фильтра сжатия и проводились измерения задержки.

Результаты моделирования представлены графиками дискриминаторных характеристик при полной длительности отраженного импульса на фиг. 1 и при укороченной длительности на фиг. 2, графиком формы вершины главного лепестка сжатого составного ЛЧМ импульса на фиг. 3 и результатами сравнения заданных и измеренных задержек, приведенные в таблице 1

Ниже дана оценка технической реализации этого способа.

Оценка выполнена при указанных выше параметрах сигнала и обработки.

Число действительных арифметических операций Q при вычислении дискриминаторной оценки задержки и абсолютного значения задержки одного отраженного импульса на выходе фильтра сжатия составит ориентировочно не более 20 арифметических операций. Если предположить, что среди 160000 каналов дальности окажется одновременно 2000 локальных максимумов целей, что маловероятно, то на все вычисления уйдет 40000 операций.

Эта обработка по сравнению с объемом вычислений в фильтре сжатия несравнимо мала. Если фильтр сжатия выполняется на программируемых интегральных микросхемах (ПЛИС) разработки фирмы Xilinx (США)типа Virtax-7 XC7VSX1140T, где суммирование двух действительных 16-разрядных чисел реализуется с тактом 5 нс, то для всех дискриминаторных вычислений потребуется время Т=40000*5 нс=200 мкс. При допустимой длительности всей обработки, равной длительности сигнала 40 мс, эта обработка составляет 0,5%. То есть, данная обработка не имеет проблем в ее практической реализации.

Таким образом, результаты моделирования и оценки технической реализации показывают, что предлагаемый способ, заключающийся в том, что предварительно осуществляется измерение реальной нелинейной временной дискриминаторной характеристики (ДХ) составного ЛЧМ фильтра сжатия, ее квантование с заданным шагом в пределах одного такта дискретизации входного сигнала, вычисление тангенсов наклона реальной ДХ в точках квантования, формирование линейной эталонной ДХ и запоминание этих характеристик и значений тангенсов; затем в процессе штатной работы осуществляется измерение по каждой цели реальной дискриминаторной оценки задержки (ДОЗ), вычисление квантованного значения этой ДОЗ, проводится ассоциативный выбор по реальной ДОЗ соответствующего эталонного значения ДОЗ, егоуточнение с использованием значения тангенса наклона ДХ ирасчет по уточненному значению ДОЗ искомого абсолютного значения задержки; при этом обеспечивается вычисление ДОЗ, не содержащей нелинейных ошибок, ДХ становится линейной и обеспечивается вычислениеполного значения измеряемой задержки, а также пикового значения амплитуды сигнала с повышенной методической точностью.

Литература

1. Радиоавтоматика: Учеб. пособие для студ. вузов спец. «Радиотехника» / В.А. Бесекерский, А.А. Елисеев, А.В. Небылов и др.; Под ред. В.А. Бесекерского. - М.: Высш. шк., 1985. - 271 с.

2. СПОСОБ СЛЕЖЕНИЯ ЗА ВРЕМЕННОЙ ЗАДЕРЖКОЙ СИГНАЛА И УСТРОЙСТВО ДЛЯ ЕГО РЕАЛИЗАЦИИ, Патент РФ №2157593, H04L 7/02, Н04В 1/10, опубл. 10.10.2000, авторы Гармонов А.В.; Манелис В.Б.; Савинков А.Ю.

3. "Применение цифровой обработки сигналов" под ред. Э. Оппенгейма, М, 1980 г., стр. 325.

4. "Цифровой дискриминатор задержки", патент РФ №2059254, опубл. 27.04.1996, автор Наумов С.В.

Способ дискриминаторного измерения задержки отраженного импульса (ОИ) на выходе фильтра сжатия составного ЛЧМ импульса, предназначенный для оценки положения пика сжатого отраженного импульса (ОИ) по дальности и заключающийся в вычислении разности модулей двух амплитуд и делении ее на сумму этих модулей, отличающийся тем, что предварительно вычисляется реальная временная дискриминаторная характеристика (ДХ) составного ЛЧМ фильтра сжатия путем первоначального задания целочисленной задержки N_з ОИ от начала приема на входе фильтра сжатия и последующего изменения задержки с малым шагом квантования р=1/Р от нуля до единицы интервала дискретизации сигнала Т_Д, с вычислением безразмерных дискриминаторных оценок задержки (ДОЗ) dT_0p по формуле:

где p - номер точки (шага) квантования на интервале дискретизации Т_Д, m - номер отсчета локального максимума (ЛМ) в амплитудном рельефе сжатого сигнала на выходе фильтра сжатия от начала приема, причем знак плюс перед dT_0p ставится тогда, когда соседняя амплитуда слева (A_m-1) больше соседней амплитуды справа (A_m+1), иначе - ставится знак минус; кроме того вычисляются тангенсы углов наклона ДХ tgα_p в точках квантования интервала дискретизации по формуле:

измеренные заранее значения реальной временной ДХ преобразуются в порядковые номера (адреса), по которым в постоянное или оперативное запоминающее устройство (ЗУ) заносятся и там хранятся реальные значения dT_0p, тангенсы углов наклона ДХ tgα_p и точные (эталонные) значения dt_ЭТ ДХ; эталонная ДХ представляет собой прямую, идущую под углом 45° от значения минус 0,5Т_Д к значению плюс 0,5Т_Д; при проведении штатной работы в режиме сопровождения по формуле (1) по вновь измеренным значениям амплитуд A_m, A_m±1, A_m±1 вычисляется ДОЗ dT_n, где n - номер шага квантования; квантованное значение этой ДОЗ dT_nкв является адресом, по которому из ЗУ выбирается соответствующее истинное значение ДОЗ dt; при этом в режиме сопровождения проводится вычисление уточненного значения ДОЗ dt_уточн с использованием дробной части Δ_С заданной задержки N_з по формулам:

где dN_з - остаток заданной задержки, dN_з.кв - квантованное значение остатка заданной задержки, выраженное в долях интервала дискретизации Т_Д; значение Δ_C может поступать в составе целеуказания; в режиме обнаружения проводится уточнение значения ДОЗ с использованием остатка измеряемой задержки Δс учетом тангенса угла наклона ДХ по формулам:

при этом измерение полного значения задержки N_Зuзм в обоих режимах обнаружения и сопровождения осуществляется путем экстраполяции полного измеренного значения задержки m уточненным значением ДОЗ dt_уточн по формуле:

где m - номер отсчета, в котором находится амплитуда ЛМ сжатого ОИ; исходя из линейной зависимости амплитуды сжатого импульса в пределах одного такта дискретизации вычисляется пиковое значение амплитуды по формуле:

где А_ЛМ - амплитуда ЛМ, А_ПИК - искомое значение амплитуды пика;

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