О МНОГОМЕРНОМ ПРЕДСТАВЛЕНИИ СИГНАЛОВ АВТОМАТИЧЕСКИХ ДОЗИРУЮЩИХ УСТРОЙСТВ
Авторы:
1.
Сибирская генерирующая компания
2.
Кемеровский государственный университет
3.
Кемеровский государственный университет
4.
Сибирский государственный индустриальный университет
(д-р техн. наук, профессор; профессор кафедры автоматизации и информационных систем)
Тип:
Статья
Страницы:
с 26
по 32
Статус:
Опубликован
Получено:
05.10.2018
Одобрено:
25.10.2018
Опубликовано:
25.10.2018
Классификаторы:
ВАК 05.12.2013 Системы, сети и устройства телекоммуникаций
ВАК 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
ВАК 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
ВАК 05.13.10 Управление в социальных и экономических системах
ВАК 05.13.18 Математическое моделирование, численные методы и комплексы программ
ВАК 05.13.19 Методы и системы защиты информации, информационная безопасность
УДК 681.5:517.443
ГРНТИ 20.01 Общие вопросы информатики
ГРНТИ 28.01 Общие вопросы кибернетики
ГРНТИ 49.01 Общие вопросы связи
ГРНТИ 50.01 Общие вопросы автоматики и вычислительной техники
ГРНТИ 82.01 Общие вопросы организации и управления
ВАК 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
ВАК 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
ВАК 05.13.10 Управление в социальных и экономических системах
ВАК 05.13.18 Математическое моделирование, численные методы и комплексы программ
ВАК 05.13.19 Методы и системы защиты информации, информационная безопасность
УДК 681.5:517.443
ГРНТИ 20.01 Общие вопросы информатики
ГРНТИ 28.01 Общие вопросы кибернетики
ГРНТИ 49.01 Общие вопросы связи
ГРНТИ 50.01 Общие вопросы автоматики и вычислительной техники
ГРНТИ 82.01 Общие вопросы организации и управления
Язык материала:
русский
Ключевые слова:
процесс автоматического дозирования, 1D-сигнал, 2D/3D-пространство, мгновенная частота, время-частотное отображение, многомерное распределение, мультикомпонентный сигнал
Аннотация (русский):
Предметом исследования является разработка специфического подхода, основанного на применении многомерных вейвлет-распределений для реализации функций мониторингового управления процессами дозирования в составе смесеприготовительного агрегата. В качестве ключевого компонента этого подхода используется многомерное время-частотное распределение Вигнера - Вилле, входящее в состав распределений класса Коэна. Исследование фокусируется на сигналах, характеризующих процессы массопереноса в виде материалопотоковых измерительных сигналов расхода в соответствующих узлах агрегата. Показано, что распределение Вигнера - Вилле в терминах времени представляет собой преобразование Фурье произведений частей рассматриваемого сигнала для прошедшего и будущего моментов времени, а соответствующее распределение для частотного спектра - преобразование Фурье произведений частей сигнала для высокочастотного и низкочастотного фрагментов спектра. Отмечено, что при использовании комплексной модели сигнала дозирования дискретные отсчеты рассматриваются как его вещественные значения. Описание параметров сигнала (амплитуды, фазы, частоты) проведено с привлечением преобразования Гильберта. В распределения класса Коэна, которыми отображаются одномерные нестационарные сигналы расхода, введено понятие «мгновенной частоты». Представлено графическое пояснение преобразования технологического сигнала расхода из одномерной временной области во время-частотное 2 D/ 3 D -пространство. Подробно разобрана технология формирования многомерного отображения в виде распределения Вигнера для одномерных сигналов расхода дозаторов непрерывного действия спирального или шнекового типа. Рассмотрены свойства поддержки распределения Вигнера, позволяющие судить о наличии или отсутствии время-частотных элементов распределения на интервале регистрации сигнала. Продемонстрировано получение распределения Вигнера для сигнала дозирования непрерывно-прерывистого типа. Сделан вывод о том, что при определенном характере сигнала для нулевых фрагментов последнего на распределении возникают ненулевые, т. е. виртуальные (паразитные), время-частотные элементы. Помимо распределения Вигнера, рассмотрены два других распределения - Рихачека и Пейджа, - отображающие тот же сигнал и также содержащие виртуальные элементы, но в разных частях время-частотного пространства. Аналитически представлено обобщенное многомерное распределение многокомпонентного сигнала с наличием ядра распределения, в состав которого входит параметр коррекции, позволяющий управлять интенсивностью виртуальной энергии сигнала.
Предметом исследования является разработка специфического подхода, основанного на применении многомерных вейвлет-распределений для реализации функций мониторингового управления процессами дозирования в составе смесеприготовительного агрегата. В качестве ключевого компонента этого подхода используется многомерное время-частотное распределение Вигнера - Вилле, входящее в состав распределений класса Коэна. Исследование фокусируется на сигналах, характеризующих процессы массопереноса в виде материалопотоковых измерительных сигналов расхода в соответствующих узлах агрегата. Показано, что распределение Вигнера - Вилле в терминах времени представляет собой преобразование Фурье произведений частей рассматриваемого сигнала для прошедшего и будущего моментов времени, а соответствующее распределение для частотного спектра - преобразование Фурье произведений частей сигнала для высокочастотного и низкочастотного фрагментов спектра. Отмечено, что при использовании комплексной модели сигнала дозирования дискретные отсчеты рассматриваются как его вещественные значения. Описание параметров сигнала (амплитуды, фазы, частоты) проведено с привлечением преобразования Гильберта. В распределения класса Коэна, которыми отображаются одномерные нестационарные сигналы расхода, введено понятие «мгновенной частоты». Представлено графическое пояснение преобразования технологического сигнала расхода из одномерной временной области во время-частотное 2 D/ 3 D -пространство. Подробно разобрана технология формирования многомерного отображения в виде распределения Вигнера для одномерных сигналов расхода дозаторов непрерывного действия спирального или шнекового типа. Рассмотрены свойства поддержки распределения Вигнера, позволяющие судить о наличии или отсутствии время-частотных элементов распределения на интервале регистрации сигнала. Продемонстрировано получение распределения Вигнера для сигнала дозирования непрерывно-прерывистого типа. Сделан вывод о том, что при определенном характере сигнала для нулевых фрагментов последнего на распределении возникают ненулевые, т. е. виртуальные (паразитные), время-частотные элементы. Помимо распределения Вигнера, рассмотрены два других распределения - Рихачека и Пейджа, - отображающие тот же сигнал и также содержащие виртуальные элементы, но в разных частях время-частотного пространства. Аналитически представлено обобщенное многомерное распределение многокомпонентного сигнала с наличием ядра распределения, в состав которого входит параметр коррекции, позволяющий управлять интенсивностью виртуальной энергии сигнала.
Ключевые слова:
процесс автоматического дозирования, 1D-сигнал, 2D/3D-пространство, мгновенная частота, время-частотное отображение, многомерное распределение, мультикомпонентный сигнал
процесс автоматического дозирования, 1D-сигнал, 2D/3D-пространство, мгновенная частота, время-частотное отображение, многомерное распределение, мультикомпонентный сигнал
Объектом предлагаемого исследования являются процессы автоматического дозирования сухих мелкодисперсных или зернистых материалов, которые востребованы в пищевой и пищеперерабатывающей промышленности, стройиндустрии, фармацевтической промышленности, в ряде секторов аграрно-промышленного комплекса. Предметная область исследования представляет собой разработку подхода, основанного на применении многомерных вейвлет-распределений для реализации семантически прозрачных функций мониторинга и синхронного автоматического управления режимами работы дозирующих устройств (дозаторов) в составе смесеприготовительных агрегатов. В качестве основного компонента разрабатываемого подхода используется многомерное время-частотное распределение Вигнера - Вилле, относящееся к так называемому классу Коэна [1]. При этом внимание фокусируется на сигналах, характеризующих режимные процессы массопереноса в виде материалопотоковых измерительных сигналов расхода в соответствующих технологических точках (узлах) смесеприготовительного агрегата. Время-частотное представление измерительного сигнала Исходя из функционального структурного представления исследуемого сигнала, можно разработать некоторые простые эмпирические правила, чтобы выяснить характер распределения Вигнера [2-4], отображающего сигнал во время-частотном пространстве. Распределения Вигнера - Вилле в терминах времени и частотного спектра соответственно выглядят следующим образом: W(t, ω) = (2π)-1(t -τ/2)exp(-jωτ)s(t + τ/2)dt; (1) W(t, ω) = (2π)-1 (ω + θ/2)exp(-jtθ)s(ω - θ/2)dθ, (2) где интегралы представляют собой преобразования Фурье автокорреляционных функций сигнала, а в отсутствие экспонент интегралы являются свертками сигнала s(t) для прошедшего и будущего интервалов времени, а также для высокочастотного и низкочастотного диапазонов спектра сигнала относительно текущей мгновенной частоты; при этом первый интеграл - суть преобразование Фурье произведения частей сигнала s(t) для прошедшего и будущего моментов времени, а второй - то же для высокочастотной и низкочастотной частей спектра; символом * отмечена операция сопряжения комплексной сигнальной функции s(t); t, τ - текущие моменты времени; ω, θ - мгновенные угловые частоты в спектре сигнала. При использовании комплексной (аналитической) модели сигнала h(t) = x(t) + jy(t), j = (-1)0,5, значения дискретных отсчетов сигнала x(t) рассматриваются как его вещественные части. Значения мнимой части y(t) аналитического сигнала однозначно определяются из x(t) с помощью преобразования Гильберта [5, 6]. Преобразование Гильберта функции f(t) и его обратное преобразование определяются в виде FHi(t) = (1/π) f(τ)dt/(τ - t) = (-1/πt)*f(t); f(t) = -(1/π) FHi(ξ)/(ξ - t) = -(1/πt)*FHi. С использованием преобразования Гильберта аналитический сигнал, связанный с вещественной функцией f(t), интерпретируется комплексной функцией s(t) = f(t) - jFHi(t). Тогда мгновенные амплитуда, фаза и частота сигнала s(t) будут определяться следующим образом: Ai(t) = [f 2(t) + FHi2(t)]0,5; Ф(t) = arctg[-FHi(t)/f(t)]; ωi = 2πfi = dФ[ω(t)]/dt = Im[dlns(t)/dt]= Im[(1/s)(ds(t)/dt)]. Понятие «мгновенной частоты» здесь основано на обобщенном представлении технологических сигналов в виде нестационарных по частоте процессов, т. е. на присутствии в составе сигнала в каждый текущий момент времени гармонической составляющей определенной варьируемой частоты. Заметим, что в случае использования сигнала в виде вещественной функции вместо сопряженной берется собственно сигнальная функция. Для идентификации мультикомпонентных нестационарных сигналов, к которым относятся сигналы расхода на выходе блока дозаторов, и состоящих из ряда 1D-нестационарных сигналов воспользуемся методом их анализа на время-частотной основе. Время-частотный анализ выполняется с использованием специфических время-частотных преобразований. Рис. 1 поясняет способ представления сигнала во время-частотном пространстве. Рис. 1. Отображение 1D-сигнала на 2D-время-частотную область (при учете амплитудного фактора сигнала - на 3D-пространство) Время-частотное отображение технологического сигнала представляет собой преобразование сигнала из одномерной (1D-) временной области во время-частотную, двумерную/трехмерную (2D-/3D-). При этом преобразование производится в интегральной форме. В частности, на рисунке условно показано преобразование одномерного импульсного сигнала, заданного во временной области (сигнал с теневым фоном), в 3D-отображение время-частотного типа с двумя компонентами разной частоты и амплитуды в трехкоординатной системе (с осями «время», «частота», «амплитуда»). Особенности формирования 2D/3D-распределения Вигнера на примере сигналов дозирования Из уравнения (1) видно, что в некоторый текущий момент времени складываются части, сформированные из произведений сигнала в прошедший момент, умноженные на сигнал в будущий момент, причем прошедший интервал времени будет равен интервалу времени в будущем. Поэтому, чтобы увидеть, является ли распределение Вигнера нулевым в какой-либо момент времени, следует мысленно свернуть часть сигнала слева от данного момента с правой частью и посмотреть, есть ли какое-либо перекрытие. Если это так, распределение Вигнера не будет равно нулю, в противном случае оно будет нулевым. Рассмотрим ненулевой непрерывный сигнал конечной длительности (например, сигнал расхода дозатора непрерывного действия: спирального или шнекового) на интервале (t1, t2). При любых моментах времени, меньших t1 (для которых сигнал отсутствует), свертка такого сигнала с сигналом в моменты времени, большие t1, будет равна нулю. Следовательно, для сигналов с конечной длительностью распределение Вигнера равно нулю вплоть до момента возникновения сигнала. Это позитивное свойство распределения, поскольку на карте Вигнера не должно отображаться никаких ненулевых время-частотных элементов распределения, если сигнал равен нулю. В любой точке справа от t1, но со значениями, меньшими t2, свертка сигналов дает ненулевые значения. Аналогичные рассуждения справедливы для всех моментов времени, больших t2. Поэтому для ограниченного по времени сигнала распределение Вигнера всегда равно нулю до начала сигнала и после окончания сигнала, т. е. W(t, ω) = 0 для t ≤ t1 или t ≥ t2, если сигнал s(t) отличен от нуля только в диапазоне (t1, t2). Вследствие аналогичности структур уравнений (1) и (2) те же соображения относятся и к частотной области. Если анализируемый сигнал (например, сигнал материалопотока в той или иной технологической точке смесеприготовительного агрегата) имеет ограниченную полосу частотного спектра, распределение Вигнера будет нулевым для всех частот, которые не включены в эту полосу, т. е. W(t, ω) = 0 для ω ≤ ω1 или ω ≥ ω2, если частотный спектр сигнала s(ω) отличен от нуля только в диапазоне (ω1, ω2). Эти аспекты - свойства поддержки распределения Вигнера. При этом диапазон поддержки равен интервалу регистрации сигнала материалопотока дозирования, на котором производится цифровая обработка с целью расчета соответствующего время-частотного распределения. Сравнительное представление разных время-частотных отображений 1D-сигнала процесса дозирования Далее рассмотрим сигнал расхода для дозирования непрерывно-прерывистого типа, соответствующий сигналу материалопотока на выходе дозатора непрерывного действия, работающего в импульсно-периодическом режиме, для которого, как в вышеприведенном случае, его значения - ненулевые в диапазоне (0, t1), а также при (t2, t3) (рис. 2). а б в Рис. 2. Распределения Вигнера (а), Рихачека (б), Пейджа (в) В диапазоне же от t1 до t2 в сигнале возникает пауза, сигнал s(t) = 0. Рассмотрим интервал паузы и некоторую точку (момент времени) внутри этого интервала. Свертывая левую часть сигнала, лежащую слева от произвольной точки внутри интервала, с правой частью сигнала s(t) (здесь в обеих частях сигнала содержатся его ненулевые значения), получаем эффект перекрытия. Следовательно, распределение Вигнера для интервала паузы - ненулевое, несмотря на то, что сигнал здесь отсутствует. Отсюда вывод: время-частотное распределение Вигнера не равно нулю, когда сигнал равен нулю; это обстоятельство вызывает значительные трудности при интерпретации время-частотных отображений, к которым относится и распределение Вигнера. На рис. 2 помимо распределения Вигнера показаны для сравнения два других распределения (Рихачека и Пейджа [1, 7]), которым свойственны те же эффекты - возникновение виртуальных элементов для тех интервалов исследуемого сигнала, на которых соответствующее распределение должно быть равно нулю. Здесь сигнал дозирования (в виде переменной составляющей расхода) на выходе спирального дозатора возникает в нулевой момент с частотой ω1 и блокируется в момент t1, затем процесс дозирования снова начинается в момент t2 уже с частотой ω2 ˃ ω1 и прекращается при t = t3. Паразитные элементы распределений отображаются на разных распределениях по-разному: на распределении Вигнера - в диапазоне паузы на полусуммарной частоте [1] двух интервалов дозирования, на распределении Рихачека - в диапазонах работы дозатора, но проявляют себя на обратных частотах по отношению к реальным интервалам дозирования. На распределении Пейджа виртуальная энергия сигнала дозирования сосредоточена во временном диапазоне более высокочастотного процесса дозирования (с частотой ω2), однако проявляет себя в этом диапазоне на меньшей частоте ω1. Коррекция распределения Вигнера Подставив выражение для мультикомпонентного сигнала , где k - номер сигнальной компоненты, N - число компонент в составе сигнала, в запись обобщенного распределения [7, 8] класса квадратичных время-частотных распределений (3) где Ф(θ, τ) - ядро распределения, получим (4) где и - соответственно автономные и перекрестные (интерференционные, паразитные) члены, формирующие полную энергию сигнала. В соответствии с [9] ядро распределения Вигнера, способное минимизировать паразитную энергию сигнала, имеет вид (5) где σ - коэффициент коррекции, управляющий (способный минимизировать) интенсивностью виртуальной энергии сигнала. Подставив запись ядра (5) в обобщенное выражение распределения и проинтегрировав по θ, получим в итоге: (6) Данное выражение представляет собой модифицированную форму корректирующего распределения Вигнера - Вилле (распределения Чуи - Уилльямса) с использованием коррекции в виде экспоненциального ядра (5). Выводы 1. Таким образом, негативной особенностью время-частотного распределения Вигнера является возникновение в нем (при определенной структуре технологического сигнала) виртуального отображения несуществующих сигнальных компонент. 2. На практике, при выполнении расчетов распределений Вигнера - в рамках алгоритма оценки текущих режимов функционирования технологического оборудования (здесь - дозирующих устройств непрерывного или дискретного действия) - паразитные элементы распределения вырезаются непосредственно из формул (3) и (4) и в дальнейших расчетах не участвуют. 3. В качестве альтернативного способа получения отображений информационных сигналов в автоматизированных смесеприготовительных комплексах рекомендуется пользоваться время-частотным распределением Чуи - Уилльямса (6), варьируя при этом входящий в его состав коэффициент коррекции (позволяющий менять величину виртуальной энергии сигнала) в желаемом диапазоне.
1. Cohen L. Time-frequency distributions - A Review // Proceedings of the IEEE. July 1989. Vol. 77. N. 7. P. 941-981.
2. Cohen L. On a fundamental property of the Wigner distribution // IEEE Trans. Acoust., Speech, Signal Processing. 1987. Vol. ASSP-35. P. 559-561.
3. Debnath L. Recent development in the Wigner-Ville distribution and time-frequency signal analysis // PINSA. January 2002. 68, A. N. 1. P. 35-56.
4. Федосенков Б. А. Научно-технические основы создания и моделирования автоматизированных систем управления непрерывными смесеприготовительными процессами: автореф. дис. … д-ра техн. наук. М.: МГУПП, 2005. 55 с.
5. Huang N., Shen S. Hilbert-Huang transform and its applications. Singapore: World Scientific, 2005. 324 p.
6. Мартуганова Е. Р. Модель веб-сервиса по специализированной обработке данных на основе жадных алгоритмов. М.: МГУ им. М. В. Ломоносова, 2014. 86 с.
7. Малла С. Вейвлеты в обработке сигналов. М.: Мир, 2005. 673 с.
8. Debnath L. Wavelet transforms and their applications. Boston: Birkhauser, 2002. 565 p.
9. Swindlehurst A. L., Kailath T. Near-field source parameter estimation using a spatial Wigner distribution approach // In Advanced Algorithms and Architectures for Signal Processing III. 1989. Vol. 975. P. 86-92.
