способ определения периода многочастотного сигнала, представленного в цифровой форме

Формула изобретения

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

где n_min и n_mах – минимальная и максимальная граница ожидаемой области нахождения величины периода;

Тc₀ – ожидаемая величина периода,

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

Описание изобретения к патенту

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

Известен способ определения периода многочастотного сигнала, представленного в цифровой форме (Гольдштейн Е.И., Даниленко Т.Г. Универсальная методика определения периода многочастотного сигнала, представленного в цифровой форме //Измерительная техника. 2001. N 6. с. 47-49.), который заключается в определении периода сигнала Тс по результатам сравнения трёх интегралов:

m>k>1,

где t – период дискретизации сигнала;

u(t₁), u(t₂), u(t_m) – первый, k-й и m-й цифровой отсчёт рассматриваемого сигнала соответственно;

N - размер массива цифровых отсчётов рассматриваемого сигнала;

t_N - время, приходящееся на N–й цифровой отсчёт рассматриваемого сигнала;

k, m - целые числа.

Путём перебора параметров k u m определяют минимальное значение Д от сравнения указанных интегралов:

Минимальному значению А соответствует искомое значение периода сложного многочастотного сигнала Тc=tNmin, где N_min - номер отсчёта, соответствующий минимальному значению .

Недостатками известного способа являются большие вычислительные затраты, связанные с вычислением интегралов.

Известен способ определения периода многочастотного сигнала, представленного в цифровой форме (Шутова И.А. Алгоритмы и программы обработки ваттметрической информации для технического диагностирования скважинных штанговых насосных установок: Автореферат диссертации на соискание учёной степени канд. тех. наук: 05.13.11. Томск, 1992. 18 с.: ил. В надзаг.: Томский политехнический университет. Библиогр.: с. 17-18(13назв.) АР92-11105 МФ-93/кн. ГРНТИ 52.47.19), выбранный в качестве прототипа, заключающийся в том, что определяют размер сравниваемых подмассивов

где N – длина реализации рассматриваемого сигнала.

Далее определяют границы ожидаемой области нахождения искомой величины периода:

где n_min и n_mах – минимальная и максимальная граница ожидаемой области нахождения величины периода;

Тс₀ – ожидаемая величина периода.

Запоминают первые b отсчётов входного сигнала U(t_i), образуя тем самым начальный подмассив U(t_m), где m[0; b] параметр изменения номера отсчёта. Начиная с момента времени, соответствующего n_min, запоминают текущие значения отсчётов входного сигнала U(t_i), образуя тем самым сравниваемый подмассив U(t_m+j) размера b, сдвинутый относительно начального подмассива на интервал времени j[n_min; n_max]. Сравнивают отсчёты начального и сравниваемых подмассивов, суммируют результаты сравнения. Путём перебора всех возможных значений сравниваемых подмассивов в ожидаемой области нахождения величины периода Т_C образуют массив сумм квадратичных оценок D_j:

где U(t_m) – m-й отсчёт начального подмассива рассматриваемого сигнала;

U(t_m+j) – m-й отсчёт сравниваемого подмассива рассматриваемого сигнала;

j[n_min, n_max] – параметр сдвига сравниваемого подмассива.

В соответствии с интервалом времени, который приходится на D_min, выдают искомую величину Тc.

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

Задачей предлагаемого изобретения является упрощение способа определения периода многочастотного сигнала, представленного в цифровой форме.

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

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

n_min=Tc₀-b;

n_max=Tc₀+b,

где n_min и n_max - минимальная и максимальная граница ожидаемой области нахождения величины периода;

Тc₀ - ожидаемая величина периода, между начальным и сравниваемыми подмассивами определяют арифметическое отклонение.

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

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

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

На фиг.1 изображена структурная схема реализации предлагаемого способа определения периода многочастотного сигнала, представленного в цифровой форме.

На фиг.2 изображена зависимость сумм арифметических оценок от времени.

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

В таблице 1 приведены значения сигнала U(t_i), взятые через период Тc.

В таблице 2 приведены рекомендуемые значения размера сравниваемых массивов в зависимости от максимального уровня шумов и максимальной погрешности определения величины Тc.

Способ определения периода многочастотного сигнала, представленного в цифровой форме, реализован с помощью структурной схемы (фиг.1), содержащей запоминающие устройства 1 и 2 (ЗУ1 и ЗУ2 соответственно), на входы которых подаётся сигнал U(t_i) в цифровой форме, устройство вычитания 3 (УВ), входы которого соединены с выходами запоминающих устройств 1 и 2, устройство суммирования 4 (УС), вход которого соединён с выходом устройства вычитания 3, устройство минимизации суммы арифметических оценок 5 (УМ), вход которого соединён с выходом устройства суммирования 4, а выход с входом устройства вывода информации 6 (У Выв).

Запоминающие устройства 1 и 2 (ЗУ1 и ЗУ2) могут быть реализованы на основе оперативного запоминающего устройства микроконтроллера серии AT89S8252, AT89S4D12 или аналогичной серии. Устройство вычитания 3 (УВ), устройство суммирования 4 (УС), устройство минимизации суммы арифметических оценок 5 (УМ) и устройство вывода информации 6 (У Выв) могут быть реализованы на микроконтроллере серии AT89S8252, AT89S4D12 или аналогичной серии.

На входы запоминающих устройств 1 и 2 (ЗУ1 и ЗУ2) был подан следующий тестовый сигнал:

Истинный период данного сигнала составляет Тc=0,2 с, что можно наблюдать по данным из таблицы 1. В данной таблице приведены k отсчётов рассматриваемого сигнала, взятых относительно первого отсчёта (U(t_k)) и через период этого сигнала (U(t_k+Tc))

Сначала выбирают параметр b в зависимости от максимального уровня помех. В таблице 2 приведены рекомендуемые данные, полученные экспериментальным путём, по выбору параметра b в зависимости от максимального уровня шумов R_max. При проведении экспериментов использовался центрированный белый шум. В данной таблице приведена также погрешность _max определения величины Тc в зависимости от параметров R_max и b.

Далее определяют границы ожидаемой Тс₀-области нахождения искомой величины периода:

n_min= Тс₀ - b;

n_mах=Тc₀+b,

где n_min и n_max – минимальная и максимальная граница ожидаемой области нахождения величины периода;

Тс₀ – ожидаемая величина периода.

Для сигнала U(t_i) с Тс₀=0,2 с при b=50 и t=10^-4 секунд получим n_min=1950 и n_max=2050.

С помощью запоминающего устройства 1 (ЗУ1) запоминают первые b отсчётов входного сигнала U(t_i), образуя тем самым начальный подмассив U(t_m), где m[0; b] - параметр изменения номера отсчёта. Начиная с момента времени, соответствующего n_min, с помощью запоминающего устройства 2 (ЗУ2) запоминают текущие значения отсчётов входного сигнала U(t_i), образуя тем самым сравниваемый подмассив U(t_m+j) размера b, сдвинутый относительно начального подмассива на интервал времени j[n_min; n_max]. С помощью устройства вычитания 3 (УВ) сравнивают отсчёты начального и сравниваемых подмассивов без учёта знака разности сравнения, что равносильно вычислению модуля разности. С помощью устройства суммирования 4 (УС) суммируют результаты сравнения. Путём перебора всех возможных значений сравниваемых подмассивов в ожидаемой области нахождения величины периода Тc образуют массив сумм арифметических оценок D_j:

На фиг. 2 приведён график изменения массива D_j в зависимости от текущего времени. С помощью устройства минимизации сумм арифметических оценок 5 (УМ) определяют минимальное значение массива сумм арифметических оценок D_min. На фиг.3 ясно видно, что минимум D_min массива D_j приходится на искомую величину истинного периода рассматриваемого сигнала. Устройство вывода информации 6 (У Выв) в соответствие с интервалом времени, который приходится на D_min, выдаёт искомую величину Тc, которая для рассматриваемого сигнала U(t_i) составила 0,2 секунды.

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