СПОСОБ СТАТИСТИЧЕСКОЙ ОБРАБОТКИ СЛУЧАЙНЫХ ПОТОКОВ ИМПУЛЬСОВ С ДИСКРЕТНЫМ ВРЕМЕНЕМ Российский патент 2006 года по МПК G01R23/02 H04B17/00 

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

Известны способы обработки случайных потоков импульсов, выполняемой в интересах получения информации о наличии полезных сигналов, оценки их параметров и выделения потоков импульсных сигналов, принадлежащих отдельным источникам радиоизлучения. Первая группа известных способов основана на методах совпадения импульсных сигналов [1. Финкелыптейн М.Н. Основы радиолокации. - М.: "Сов. радио", 1973, с.254], в том числе реализуемые с применением накопителей на линиях задержки с положительной обратной связью [2. Лезин Ю.С. Оптимальные фильтры и накопители импульсных сигналов. - М.: "Сов. радио", 1969]. Недостатками этих способов являются ограниченные функциональные возможности из-за невыполнения функции фильтрации импульсных потоков, а также потребность в априорной информации о тактовом интервале обнаруженных потоков импульсов.

Другие из известных способов, обладающие возможностью выделения (фильтрации) импульсных потоков, основаны на селекции (отборе) определенных импульсов входящего потока, отличающихся количественными признаками их формы (амплитуды, длительности), а в некоторых случаях и признаками "содержания" радиоимпульсов (например, вида и параметров внутриимпульсной модуляции [3. Справочник по основам радиолокационной техники. - М.: В.И., 1967, с.519]. Недостатком этих способов является невысокие точностные и надежностные характеристики принимаемых решений при совместном обнаружении, оценке определяющих параметров и фильтрации случайных потоков бинарно-квантованных импульсов стандартной формы с неизвестным дискретным временем при наличии мешающих импульсов и ошибок наблюдения.

Периодический поток импульсов является частной моделью случайного потока импульсов с дискретным временем при тактовой вероятности появления импульса, равной единице. С учетом этого наиболее близким по технической сущности к предложенному способу является способ обнаружения периодических импульсных последовательностей и оценки их периода, включающий измерение временных интервалов (ВИ) между последовательными моментами прихода импульсов входящего потока, накопление первичной выборки (массива) из N значений ВИ, вычисление первой, второй,..., (N+1)-й вторичных выборок ВИ между моментом прихода первого, второго,..., (N+1)-го опорного импульса соответственно до моментов прихода всех других импульсов накопленной реализации входящего потока, вычисление в рамках отдельных вторичных выборок ВИ статистических оценок наибольшего общего делителя (НОД) возможных сочетаний из двух и более ВИ, формирование частных рядов распределения оценок НОД для отдельных вторичных выборок ВИ, объединение частных рядов распределения оценок НОД, в которых выявлено группирование их значений, совместное обнаружение случайного потока импульсов с дискретным временем по факту группирования оценок НОД в объединенном ряду распределения их значений и оценку его тактового интервала (ТИ) средним значением выявленной группы оценок НОД [4. Россия, заявка на изобретение №2003115913/28 от 27.05.2003, G 01 R 23/02].

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

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

Поставленная задача решается за счет того, что в известном способе, включающем измерение временных интервалов (ВИ) между последовательными моментами прихода импульсов входящего потока, накопление первичной выборки (массива) из N значений ВИ, вычисление первой, второй,..., (N+1)-й вторичных выборок ВИ между моментом прихода первого, второго,..., (N+1)-го опорного импульса соответственно до моментов прихода всех других импульсов накопленной реализации входящего потока, вычисление в рамках отдельных вторичных выборок ВИ статистических оценок наибольшего общего делителя НОД возможных сочетаний из двух и более ВИ, формирование частных рядов распределения оценок НОД для отдельных вторичных выборок ВИ, объединение частных рядов распределения оценок НОД, в которых выявлено группирование их значений, совместное обнаружение случайного потока импульсов с дискретным временем по факту группирования оценок НОД в объединенном ряду распределения их значений и оценку его ТИ средним значением выявленной группы оценок НОД. Новым является то, что по порядковым номерам вторичных выборок ВИ, при обработке которых выявлено группирование оценок НОД в окрестности найденной оценки ТИ, определяют базовую группу импульсов обнаруженного потока с дискретным временем, вычисляют моменты прихода импульсов базовой группы относительно ее первого импульса и их кратности в единицах оценки ТИ, вычисляют с использованием данных о базовой группе импульсов по формуле несмещенной линии регрессии с известным угловым коэффициентом, равным оценке ТИ, прогнозируемый относительный момент t′_k дискретного времени, превышающий относительное время t′_N+1 прихода последнего импульса первичной реализации потока, формируют с привязкой к относительному моменту времени t′_k симметричный логический строб, управляющий прохождением ожидаемого импульса входящего потока, в случае появления импульса в стробе принимают его за отфильтрованный импульс обнаруженного потока и измеряют относительное время его прихода t′_n+1, которое вместе с его кратностью в единицах оценки ТИ записывают в массив данных о базовой группе импульсов, если общее количество импульсов в базовой группе превысило базовое число n≥2 процесса прогнозирования, то исключают первый импульс базовой группы и корректируют текущее относительное время и данные о частично обновленной базовой группе импульсов, уточняют оценку ТИ по величине углового коэффициента несмещенной линии регрессии данных об обновленной базовой группе импульсов и вычисляют очередной прогнозируемый момент t′_k+1 дискретного времени с уточненной оценкой ТИ, в случае непоявления импульса в стробе вычисляют прогнозируемый момент t′_k+1 дискретного времени при прежнем значении оценки ТИ, если в течение определенного числа сеансов прогнозирования моментов дискретного времени не появился ни один импульс входящего потока в выставляемых стробах, то принимают решение об исчезновении обнаруженного потока и производят накопление новой первичной выборки ВИ, причем значения прогнозируемого момента дискретного времени t′_k и уточняемой оценки ТИ определяют соответственно по формулам:

где - оценка тактового интервала;

t′_i, m_i - относительный момент прихода i-го импульса базовой группы размером из n≥2 импульсов и его кратность в единицах оценки ТИ соответственно (t′₁=0, m₁=0);

(m_n+k) - кратность в единицах оценки ТИ случайного временного интервала от первого импульса базовой группы до прогнозируемого момента времени t′_k, превысившего относительный момент прихода t′_N+1 последнего импульса первичной реализации входящего потока плюс время τ_B<Т, затрачиваемое вычислительным устройством на выполнение необходимых расчетов;

К_г - максимальная "глубина" прогнозирования моментов t′_k в единицах оценки ТИ при отсутствии импульсов в последовательно выставляемых стробах;

За статистическую оценку НОД значений ВИ между импульсами входящего потока принимают максимальный коэффициент линии регрессии, проходящей через начало прямоугольных координат (несмещенная линия регрессии) с наименьшей суммой квадратов уклонений собственных узловых точек полурешетки [τ, m], где τ - ось ординат непрерывных значений ВИ, m=1, 2,... - ось абсцисс целочисленных кратностей ВИ в единицах оценки ТИ обнаруживаемого случайного потока импульсов с дискретным временем [5. Анишин А.С., Батурин Ю.О. Журнал "Радиотехника", 2002, №10, с.73-74].

Факт группирования оценок НОД в ряду распределения их значений выявляют по превышению заданного порога N₀₁ количеством оценок НОД, попадающих в заданный по величине интервал ΔT идентичности (близости) их значений.

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

Более подробно существо предложенного способа рассмотрим с использованием временных диаграмм, приведенных на фиг.2,...,6.

Пусть наблюдается реализация входящего потока стандартных по форме импульсов (фиг.2а), состоящая из пяти импульсов (№2, 3, 4, 6, 7) случайного потока с дискретным временем и трех мешающих импульсов (№1, 5, 8) с непрерывным временем случайных моментов их прихода.

Исходной информацией для предложенного способа являются измеренные значения 7-и ВИ (τ_i,i+1,i=) между последовательными моментами прихода 8-и импульсов входящего потока. Запишем значения ВИ первичной выборки в нулевую строку таблицы (фиг.3).

По значениям ВИ первичной выборки вычислим первую, вторую,..., восьмую вторичные выборки ВИ, отсчитываемых в двух направлениях (в "прошлое" и "будущее") от момента прихода соответственно первого, второго,..., восьмого опорного импульса до моментов прихода каждого импульса накопленной реализации входящего потока. Значения ВИ первой, второй и т.д. вторичных выборок запишем в соответствующие строки таблицы (фиг.3). Случайные ВИ с дискретным, а в общем случае многомодальным распределением (из-за ошибок наблюдения моментов прихода импульсов) в таблице помечены штриховкой.

В рамках каждой вторичной выборки ВИ организуем перебор возможных сочетаний из двух и более ВИ. Каждое сочетание ВИ обработаем с помощью оператора статистической оценки НОД.

По результатам обработки вторичных выборок ВИ сформируем частные ряды распределения (спектры) оценок НОД, основные типовые разновидности которых приведены на фиг.4 и 5. При обработке вторичной выборки ВИ, образованных относительно опорного импульса, находящегося на дискретной шкале времени, спектр оценок НОД преимущественно будет иметь вид, приведенный на фиг.4. Здесь отмечается группирование оценок НОД по факту превышения первого порога N₀₁ количеством "совпадающих" оценок НОД.

Напротив, при обработке вторичной выборки ВИ, образованных относительно опорного мешающего импульса, спектр оценок НОД будет преимущественно иметь вид, приведенный на фиг.5. Здесь группирование оценок НОД не отмечается из-за отсутствия факта превышения порога N₀₁.

Частные спектры оценок НОД, в которых выявлено группирование их значений, объединяем в суммарный спектр оценок НОД (фиг.6), в котором отмечаем группирование оценок НОД по факту превышения второго порога N₀₂ количеством выявленной группы оценок НОД. Следовательно, принимаем решение об обнаружении потока импульсов с дискретным временем. За оценку ТИ обнаруженного потока импульсов с дискретным временем принимаем среднее значение "совпадающих" оценок НОД в суммарном спектре их значений.

По номерам частных спектров оценок НОД, в которых выявлено группирование оценок НОД в окрестности полученной оценки ТИ (предположим, что ими являются частные спектры НОД с номерами 2, 3, 6, 7), определяем базовую группу импульсов №2, 3, 6, 7 обнаруженного потока с дискретным временем. Импульс №4 по результатам обработки ВИ в базовую группу не входит (фиг.2б).

Далее вычисляем моменты прихода импульсов базовой группы относительно ее первого импульса (t′₁=0, t′₂=τ₂₃, t′₃=(τ₂₃+τ₃₄+τ₄₅+τ₅₆), t′₄=(t′₃+τ₆₇)) и их кратности в единицах оценки ТИ (m₁=0, m₂=2, m₃=6, m₄=7), а также - момент прихода последнего импульса первичной выборки ВИ относительно первого импульса базовой группы t′_N+1=(t′₄+τ₇₈) (см. фиг.2б).

С использованием данных о базовой группе импульсов по формуле (1) несмещенной линии регрессии с известным угловым коэффициентом Т, равным оценке ТИ, вычисляем прогнозируемый относительный момент t′_k дискретного времени, который превышает относительный момент t′_N+1 прихода последнего импульса обрабатываемой реализации потока плюс время τ_B, затрачиваемое на выполнение необходимых вычислений (τ_B<Т),

Далее сформируем с привязкой к спрогнозированному моменту t′_k симметричный логический строб, предназначенный для управления прохождением импульса входящего потока (фиг.2в). Импульс, который появляется в стробе, принимаем за отфильтрованный импульс случайного потока с дискретным временем (фиг.2г). Измеряем момент его прихода относительно первого импульса базовой группы t′₅ (фиг.2г) и вместе с его кратностью m₅ в единицах оценки ТИ записываем в массив данных о базовой группе импульсов. Если количество импульсов в базовой группе превышает базовое число процесса прогнозирования, например в рассматриваемом эпизоде n=4(n≥2), то исключаем первый импульс базовой группы и корректируем данные о частично обновленной базовой группе импульсов следующим образом (фиг.2г):

1. t′₁:=t′₂-τ₂₃=0; t′₂:=t′₃-τ₂₃; t′₃:=t′₄-τ₂₃; t′₄:=t′₅-τ₂₃;

2. m₁=0, m₂:=m₃-m₂, m₃:=m₄-m₂, m₄:=m₅-m₂,

где τ₂₃ - интервал (величина) коррекции относительного времени, связанной с исключением первого импульса базовой группы (фиг.2д).

С использованием частично обновленной базовой группы импульсов уточняем по формуле (2) линейной регрессии оценку ТИ

При не появлении импульса в стробе очередной прогнозируемый момент t′_k+1 дискретного времени вычисляем по формуле (2) при прежнем значении оценки ТИ.

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

Таким образом, в заявленном способе совместное обнаружение, оценка ТИ и фильтрация (выделение) импульсов случайного потока с дискретным временем осуществляется с использованием единых физических и математических операций, выполняемых над входящим потоком импульсов и значениями его ВИ соответственно. Способ позволяет определять значения случайного точечного процесса (каким является обнаруженный поток импульсов с дискретным временем) не только на интервале первичной выборки, но и вне этого интервала.

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

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

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

На фиг.1 приведена структурная схема устройства в виде возможного варианта технической реализации способа, на фиг.2 - временные диаграммы, поясняющие существо способа и работу устройства, на фиг.3 - таблица первичной и вторичных выборок ВИ, на фиг.4 приведен экспериментальный частный спектр оценок НОД, в которых выявлено их группирование, на фиг.5 - экспериментальный частный спектр оценок НОД, в котором отсутствует их группирование, на фиг.6 - суммарный спектр оценок НОД, в котором выявлено их группирование.

Устройство, реализующее заявленный способ (фиг.1), содержит блок 1 задания размера N первичной выборки ВИ числом (N+1) импульсов входящего потока, элемент ИЛИ 2, измеритель ВИ 3, микропроцессор 4, таймер (часы) 5 относительного времени, схему сравнения двоичных кодов (компаратор) 6, формирователь строба 7 и элемент И 8, при этом вход блока 1, объединенный с первым входом элемента И 8, является входом устройства, на который поступает входящий поток импульсов, выход блока 1 задания размера первичной выборки ВИ соединен со входом запуска таймера 5 и первым входом элемента ИЛИ 2, выход которого соединен со входом измерителя ВИ 3, разрядный выход которого подключен к разрядному входу микропроцессора 4, второй выход (кода коррекции относительного времени) которого соединен с соответствующим разрядным входом таймера 5, разрядный выход которого соединен с первым входом компаратора 6, второй вход которого подключен к третьему выходу (кода прогнозируемого момента времени) микропроцессора 4, управляющий выход которого соединен со входом запуска блока 1, выход компаратора 6 соединен со входом формирователя 7 логического строба, выход которого соединен со вторым входом элемента И 8, выход которого подключен ко второму входу элемента ИЛИ 2 и является выходом устройства.

Устройство, реализующее данный способ, работает следующим образом. Входящий поток импульсов ϕ(t) (фиг.2а) поступает на вход блока 1 задания размера выборки, который при поступлении сигнала запуска извне либо от микропроцессора 4 пропускает серию из (N+1) последовательно следующих импульсов входящего потока. Первый импульс этой серии запускает таймер 5 и через элемент ИЛИ 2 вместе последовательными импульсами серии воздействуют на вход измерителя ВИ 3. Кодовые эквиваленты N измеренных ВИ записываются в массив исходных данных микропроцессора 4. Одновременно импульсы входящего потока поступают на первый вход закрытого по второму входу элемента И 8 и на выходе устройства не появляются (фиг.2г).

В процессе обработки значений ВИ первичной выборки микропроцессор 4 по факту обнаружения случайного потока импульсов с дискретным временем на основе промежуточных решений, принятых при анализе частных спектров оценок НОД, выявляет базовую группу из 4-х импульсов (фиг.2б), вычисляет величину и формирует код коррекции (интервал τ₁₂ на фиг.2а) относительного времени прихода импульсов базовой группы, предъявляет его (код) на разрядный вход таймера 5. Затем микропроцессор 4 вычисляет с использованием оценки ТИ по формуле линейной регрессии (1) значение и формирует двоичный код прогнозируемого момента дискретного времени t′_k (фиг.2в) и выдает этот код на второй вход компаратора 6.

В момент появления на первом входе компаратора 6 двоичного кода таймера 5, равного коду прогнозируемого момента дискретного времени, на выходе компаратора 6 появляется импульс, который запускает формирователь строба 7. Сигнал строба длительностью τ_ст≪Т открывает элемент И 8 по его второму входу. При появлении во время действия строба импульса входящего потока последний проходит на выход элемента И 8 и на второй вход элемента ИЛИ 2. С выхода элемента ИЛИ 2 импульс проходит на непрерывно работающий измеритель ВИ 3, который выдает кодовый эквивалент ВИ τ₈₉ (фиг.2г) в микропроцессор 4 и частично обновляет массив текущих данных о входящем потоке импульсов.

Импульс входящего потока, который появился на выходе устройства, является отфильтрованным и одновременно базовым импульсом. Его параметры в относительном времени вносят в массив данных о базовой группе импульсов. Если количество базовых импульсов в группе превышает базовое число n (например, n=4) процесса прогнозирования, то исключают первый импульс базовой группы и корректируют данные частично обновленной базовой группы импульсов с учетом изменения на величину ВИ τ₂₃ относительного времени (фиг.2д).

На основе частично обновленной базовой группы импульсов микропроцессор 4 уточняет оценку ТИ по формуле (2) и осуществляет подготовку данных для работы устройства на новом прогнозируемом тактовом интервале обнаруженного случайного потока импульсов с дискретным временем.

Поток импульсов на выходе компаратора 6 представляет собой восстановленную периодическую последовательность импульсов с периодом, равным найденной оценке ТИ. Для повышения точности работы устройства измеритель ВИ 3 и таймер 5 относительного времени целесообразно подключать к общему источнику временных меток (на фиг.1 этот источник опущен).

Блок 1 задания размера первичной выборки ВИ представляет собой известное логическое устройство для формирования пачки из (N+1) импульсов входящего потока [А.с. 1520513 СССР. Устройство для моделирования урны /Бюл. №41, 1989]. Измеритель ВИ 3 может быть построен на основе известного измерителя [Мирский Г.Я. Радиоэлектронные измерения. - М.: "Энергия" 1985, с.175]. Таймер 5 представляет собой двоичный счетчик временных меток со схемой вычитания двоичного кода, корректирующего относительное (внутреннее) время работы элементов устройства. Компаратор 6 двоичных кодов является схемой сравнения значений двух двоичных кодов. Остальные элементы устройства являются типовыми элементами цифровой и импульсной техники.

Таким образом, предложенный способ технически реализуем и по сравнению с прототипом обладает расширенной областью применения за счет возможности совместного выполнения большего числа функций: обнаружения, оценки ТИ и фильтрации (выделения) случайного потока импульсов с дискретным временем. Заявленный способ ориентирован на решение актуальных задач синхронизации приемных и передающих устройств, в том числе конфликтующих радиоэлектронных систем, в условиях априорной неопределенности о параметрах используемых сигналов и при наличии мешающих факторов. Кроме того, новая функция фильтрации импульсных потоков может быть востребована при построении помехозащищенных каналов обработки информации многопозиционных систем определения координат источников радиоизлучений.

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

Способ статистической обработки случайных потоков импульсов с дискретным временем, включающий измерение временных интервалов (ВИ) между последовательными моментами прихода импульсов входящего потока, накопление первичной выборки (массива) из N значений ВИ, вычисление первой, второй,..., (Н+1)-й вторичных выборок ВИ между моментом прихода первого, второго,..., (N+1)-го опорного импульса соответственно до моментов прихода всех других импульсов накопленной реализации входящего потока, вычисление в рамках отдельных вторичных выборок ВИ статистических оценок наибольшего общего делителя (НОД) возможных сочетаний из двух и более ВИ, формирование частных рядов распределения оценок НОД для отдельных вторичных выборок ВИ, объединение частных рядов распределения оценок НОД, в которых выявлено группирование их значений, совместное обнаружение случайного потока импульсов с дискретным временем по факту группирования оценок НОД в объединенном ряду распределения их значений и оценку его тактового интервала (ТИ) средним значением выявленной группы оценок НОД, отличающийся тем, что по порядковым номерам вторичных выборок ВИ, при обработке которых выявлено группирование оценок НОД в окрестности найденной оценки ТИ, определяют базовую группу импульсов обнаруженного потока с дискретным временем, вычисляют моменты прихода импульсов базовой группы относительно ее первого импульса и их кратности в единицах оценки ТИ, вычисляют с использованием данных о базовой группе импульсов по формуле несмещенной линии регрессии с известным угловым коэффициентом, равным оценке ТИ, прогнозируемый относительный момент t'_k дискретного времени, превышающий относительное время t'_N+1 прихода последнего импульса первичной реализации потока, формируют с привязкой к относительному моменту времени t'_k симметричный логический строб, в случае появления импульса в стробе принимают его за отфильтрованный импульс обнаруженного потока и измеряют относительное время его прихода t'_n+1, которое вместе с его кратностью в единицах оценки ТИ записывают в массив данных о базовой группе импульсов, если общее количество импульсов в базовой группе превысило базовое число n≥2 процесса прогнозирования, то исключают первый импульс базовой группы и корректируют текущее относительное время и данные о частично обновленной базовой группе импульсов, уточняют оценку ТИ по величине углового коэффициента несмещенной линии регрессии данных об обновленной базовой группе импульсов и вычисляют очередной прогнозируемый момент t'_k+1 дискретного времени с уточненной оценкой ТИ, в случае непоявления импульса в стробе вычисляют прогнозируемый момент t'_k+1 дискретного времени при прежнем значении оценки ТИ, если в течение определенного числа сеансов прогнозирования моментов дискретного времени не появился ни один импульс входящего потока в выставляемых стробах, то принимают решение об исчезновении обнаруженного потока и производят накопление новой первичной выборки ВИ.