Текст произведения
(PDF):
Читать
Скачать
Введение Поскольку процессы смесеприготовления, в силу системно-технологических причин, являются нестационарными, время-частотно-зависимыми, традиционные методы анализа и выработка на их основе соответствующих способов управления рабочими режимами агрегата не могут быть применимы в данной ситуации. Вейвлет-управление в условиях автоматизированного агрегата, осуществляемое подачей управляющих воздействий со стороны мониторингового комплекса для поддержания начальных (номинальных) режимов, позволяет избежать возникновения неблагоприятных условий при реализации процесса смесеприготовления. Поэтому решение вопросов интенсификации и совершенствования смесеприготовительных процессов в непрерывно действующих агрегатах с центробежными и вибрационными смесителями для переработки сухих дисперсных материалов на базе новых методов управления является актуальной задачей, представляющей значительный практический интерес для ряда ведущих отраслей экономики. Целью настоящей работы является разработка нетрадиционного аппарата для исследования нестационарных время-частотно-зависимых процессов, использующего в своей основе всплесковые преобразования на основе микролокальных волновых составляющих (вейвлет-функций, или всплесков). Объект и методы исследования Объектом исследования являются процессы дозирования в смесеприготовительных агрегатах для производства сухих мелкодисперсных и/или зернистых композиций различного назначения. Предметная область исследований охватывает методы цифрового управления технологическим объектом, в основе которого - формирование так называемых квадратичных обратных связей на базе вейвлет-преобразований и время-частотных распределений. Обобщенная теория время-зависимых (динамических) спектров рассмотрена в [4]. Такие распределения представляют собой зависимость энергии (интенсивности) сигнала (здесь - материалопотока) от времени и частоты одновременно. Частотно-время-зависимые сигналы являются нестационарными. Таким образом, время-частотные распределения являются эффективным средством для анализа нестационарных сигналов. В соответствии с подобным распределением можно выявить, какая доля энергии сигнала лежит в пределах определенного временного и частотного диапазонов, то есть в рамках так называемого допустимого прямоугольника неопределенности (окна Гейзенберга) [5] на время-частотной плоскости. Из теории Фурье-анализа [1] известно, что мгновенная энергия сигнала (в нашем случае - сигнала материалопотока) x(t) интерпретируется как его интенсивность в единицу времени в момент t, то есть, |x(t)|2 а в течение Δt - |x(t)|2 Δt. Интенсивность же сигнала x(t) на единицу угловой частоты ω равна |S(ω)|2, где S(ω) - преобразование Фурье сигнала x(t). Энергия сигнала в интервале частот Δω составляет |S(ω)|2 ∙ Δω. Преобразование Фурье сигнала определяется как: . Совместную время-частотную функцию распределения энергии сигнала в интервалах времени Δt и частоты Δω определим как E(t, ω) ΔtΔω. Следовательно, полная энергия сигнала материалопотока равна: . (1) Термин «функция распределения» (или «распределение») применительно к детерминированным нестационарным сигналам используется для того, чтобы показать, как «распределяется» энергия сигнала в прямоугольнике неопределенности - ячейке Δt, Δω. Эффекты, выражающиеся в появлении интенсивностей сигнала в тех ячейках время-частотного пространства, в которых их не должно быть, являются артефактами, то есть искусственными (виртуальными) проявлениями энергии сигнала. Последние проявляются вследствие билинейной (квадратичной) структуры распределения Вигнера-Вилле [4], которое ведет к появлению так называемых перекрестных (интерференционных) членов, являющихся причиной появления виртуальной интенсивности при анализе многокомпонентных сигналов, к которым принадлежат материалопотоковые сигналы: , где k - номер компоненты; N - количество компонент в структуре сигнала. Подставив данное выражение в запись обобщенного распределения класса квадратичных время-частотных распределений [4] , (2) где Ф(θ,τ) - ядро распределения, получим: , здесь Ekk и Ekl - соответственно автономные и перекрестные члены, формирующие полную энергию сигнала. В интерпретации [4] ядро распределения Вигнера, минимизирующее виртуальную энергию сигнала, имеет вид: Ф(θ,τ) = exp(-θ2τ2/σ), (3) где σ - константа управления долей виртуальной энергии. Подставляя такую запись ядра в общую формулу распределения энергии и интегрируя по θ, получим в итоге: . Данное выражение представляет собой конечную запись распределения Вигнера с использованием коррекции в виде экспоненциального ядра (3). В табл. 1 приведены выражения ряда двумерных время-частотных распределений, которые формируются из обобщенного распределения E(t, ω) с ядром Ф(θ,τ). Здесь: x(t) - материалопотоковый сигнал; ω и τ, t, u -текущие частота и время; S*(ω) - сопряженная спектральная функция сигнала x(t): S*(ω) = Re{S(ω)} - jIm{S(ω)}, где - комплексная спектральная функция сигнала; a - коэффициент формы синуса; σ - коэффициент управления величиной виртуальных членов распределения; h(τ) - импульсная переходная характеристика смесительной системы (или ее фрагмента). Таблица 1 Квадратичные распределения и их ядра в составе (2) Распределение Ядро Ф(θ, τ) Запись распределения Е(t, ω) Вигнера-Вилле 1 Марджено-Хилла Кирквуда-Рихачека Фазонормированный синус Пейджа Чуи-Уилльямса exp(θ2τ2/σ) Спектрограмма Рассмотрим процедуру получения время-частотного распределения материалопотокового сигнала на выходе блока дозирующих устройств, включающего два шнековых дозатора с рабочими частотами ω1 и ω2 > ω1 При этом считаем, что исходные сигналы получены в результате обработки первичных сигналов посредством высокочастотной фильтрации. Это означает, что постоянная составляющая производительности дозаторов исключена из анализируемых сигналов. Итак, двухкомпонентный сигнал состоит из двух гармоник и имеет вид: x(t) = Xm1 ∙ exp(jω1t) + Xm2 ∙ exp(jω2t). Тогда исправленное распределение Вигнера запишется как . (4) Обозначим: . При σ→ ∞ limσ→ ∞ Ω (ω, ω1, ω2 … σ) = δ[ω - ½(ω1 + ω2]. Последнее выражение показывает, что на полусуммарной частоте происходят импульсные всплески в виде дельта-функций, объясняемые наличием перекрестного (третьего) члена в выражении исправленного распределения Вигнера (4). При этом, для σ = ∞ ядро (3) обобщенного распределения (2) Ф(θ,τ) = 1, что соответствует обычному распределению Вигнера с наличием виртуальных компонент. Таким образом, при возрастании/уменьшении σ величина перекрестных составляющих повышается/снижается. В случае сброса/наброса нагрузки на два дозатора шнекового или спирального типов процессы дозирования являются нестационарными, поскольку носят частотно-время-зависимый характер, то есть представляют собой частотно-параметри-ческие процессы. Суммарный сигнал имеет вид: x(t) = Xm1 ∙ exp[(j(β1t2 + ω1t)] + Xm2 ∙ exp[(j(β2t2 + ω2t)] в составе которого присутствуют две мгновенные время-зависимые частоты: ω = β1t + ω1 и ω = β2t + ω2. И в этом случае виртуальные энергетические компоненты путем соответствующего выбора демпфирующего коэффициента σ легко нейтрализуются. Минимизация виртуальных членов возможна при выборе коэффициента управления виртуальной энергией в диапазоне σ = 105 ... 6,5∙105. Для вейвлет-сигналов словаря Габора вида: т.е. сигналов, представляющих собой гармонические сигналы, модулированные функцией Гаусса, распределение Вигнера определяется как: W(t, ω) = exp [-αt2 - (ω - βt -ω0)2 / α]. Отсюда видно, что при такого рода сигналах распределение Вигнера всегда положительно. Данное свойство позитивности распределения соответствует его семантической особенности - способности отображать энергетическую плотность исследуемого сигнала. Если коэффициент α - мал, то энергия сигнала концентрируется эллипсообразно вдоль линии с частотой ω = ω0 + βt, являющейся производной фазы d[(βt/2 + ω0)t]/dt. При α = 0 (при единичной модуляции) получаем частотно-время-зависимый сигнал x(t)α=0 = exp [(jβt/2 +ω0)t], распределение Вигнера которого вырождается в дельта-функцию с наклоном β на время-частотной плоскости (рис. 1а): W(t, ω)α=0 = δ[ω - (βt +ω0)]. (5) Отсюда видно, что вся энергия локализуется вдоль мгновенной частоты (βt +ω0). При β = 0 (и α = 0) энергия концентрируется вдоль горизонтальной линии на частоте ω0 = const (рис. 1б): W(t, ω) = δ(ω - ω0). (6) а) б) Рис. 1. Распределение Вигнера для: а) гармонического сигнала с время-зависимой частотой; б) гармонического сигнала с постоянной частотой а) б) Рис. 2. Расчетные распределения Вигнера реальных сигналов для случаев: а) - формула (5) и б) - формула (6) Фактическая расчетная дискретность и частотно-временной разброс отображения одномерного сигнала в виде его представления в 2D-пространстве (МСМ-карты, то есть карты модифицированного сигнала материалопотока) (рис. 2) объясняется, главным образом, действием следующих факторов: 1. неточностью адаптивной аппроксимации исходного сигнала вследствие неквадратичности скалярного произведения по всем отсчетам аппроксимируемого сегмента осциллограммы, то есть в силу соотношения: где ; k - номер отсчета (дискреты) анализируемого сигнала; Nl - длина дискретного интервала (сегмента осциллограммы); ||x(t)|| - норма сигнала; 2. генетической (присущей по определению) локальной «атомарностью» (ограниченностью величины носителя) вейвлет-функций, используемых для аппроксимации материалопотоковых сигналов; 3. особенностью процедуры вейвлет-поиска соответствия, выражающейся в фрагментарности процесса «проецирования» тех или иных участков осциллограммы исследуемого сигнала на соответствующие словарные поля вейвлет-тезауруса; 4. дискретностью исходного одномерного сигнала, участвующего в численной реализации двумерного распределения; 5. наличием краевых эффектов, вызванных наложением ограничения на ширину расчетного окна (интервала), ведущего к возникновению разрывов первого рода на границах интервала. Таким образом, при двухкомпонентном гармоническом сигнале в распределении Вигнера возникают концентрации энергии на двух частотах (фактическое распределение энергии) и на полусуммарной частоте (виртуальное распределение энергии). При использовании в составе блока дозирующих устройств дозаторов непрерывного типа (шнековых и / или спиральных) в распределении Вигнера возникают виртуальные концентрации энергии на полусуммарных частотах - при определении средней частоты любых двух дозаторов. Следовательно, для N дозаторов непрерывного типа получаем (1/2)N(N-1) локализованных концентраций виртуальной энергии (табл. 2). Таблица 2 Количество дозаторов, N Число виртуальных членов в составе распределения Вигнера 2 3 4 5 1 3 6 10 При расчете время-частотных распределений сложных мультикомпонентных сигналов, последние выражаются в виде дискретных решетчатых функций. В этом случае, в соответствии с теоремой отсчетов [2], распределение Вигнера решетчатого сигнала определяется как: (7) где 1/T - частота отсчетов (частота дискретизации); выбирается из соотношения T ≤ π/2 ωmax, где ωmax - максимальная частота в спектре сигнала. Как видно из выражения (7), дискретное распределение Вигнера - периодично с периодом π (поскольку частота мнимой экспоненты равна 2ω), а не 2π - как в случае непрерывного сигнала. Следовательно, частота дискретизации ωs должна быть в два раза больше частоты Найквиста ωN, то есть: ωs ≥ 2ωN = 4ωmax При решетчатой форме сигнала материалопотока, дискретизированного по времени, и непрерывной частоте ω получаем дискретную версию распределения Вигнера: Данное выражение соответствует условию безутратного - в информационном смысле - восстановлению непрерывного одномерного материалопотокового сигнала из его дискретного распределения. Результаты и их обсуждение В качестве примера картины мониторинга текущих процессов в СМПА для выработки и подачи управляющих воздействий на исполнительный механизм дозатора порционного типа, на рис.3 приведены осциллограмма материалопотока при напряжении U, равном 100 В, на его электроприводе, её вейвлет-аппроксимация и карта Вигнера соответственно. Здесь изображён реальный сигнал дозирования, снятый с пьезоэлектрического преобразователя. В соответствии с теоремой отсчётов первичный недискретизированный сигнал может быть восстановлен из решетчатой функции при выполнении условия: fxd ³ 2fxmax = fN, где fxmax- мгновенная максимальная частота в спектре сигнала материалопотока; fxd - частота временной дискретизации аналогового сигнала в АЦП; fN - частота Найквиста. Анализ осциллограмм первичного сигнала материалопотока показал, что на средний интервал дозы td =1,1 с приходится 70-80 периодов, что соответствует мгновенной частоте в диапазоне 58,5-66,6 Гц. При количестве отсчетов в АЦП, равном 212 = 4096 на интервале Δt = 17,7 с, частота дискретизации составляет 232,6 Гц. В нашем случае имеем: 232,6 Гц>2 (58,5-66,6 Гц), т.е. 236,6 Гц > 117 ...133,2 Гц, что соответствует допустимости дискретизации с такой частотой. При этом количество отсчетов на длительности дозы в 1,1 с составляет 256 значений. Рис. 3. Преобразование одномерного первичного материалопотокового сигнала в 2D-отображение: а) реальный сигнал, б) реконструированный сигнал по 7 вейвлетам, в) карта Вигнера. Напряжение U = 100 В Рис. 12. С целью определения на картах Вигнера (рис. 3в) элементов, характеризующих текущие режимы процесса дозирования, производится привязка время-частотных атомов, соответствующих номинальным режимам, к определенным время-частотным фреймам. Так, на рис. 4 и 5 представлены фреймированные карты Вигнера: для порционного и спирального дозаторов (фреймы показаны в виде пунктирных окон). Создание фреймов необходимо для определения девиации атомов текущего режима относительно их номинальной локализации на карте. По величине девиации и тарировочным характеристикам - fd =Fc(u) (для дозаторов непрерывного типа) или Td= Fb(u) (для дозаторов дискретного типа) - определяется напряжение рассогласования на электроприводе в контуре регулирования соответствующего дозатора, после чего формируется номинальное напряжение, которое через тиристорный управляемый выпрямитель подается на двигатель, приводя последний к номинальному состоянию по частоте вращения. Таким образом реализуется автоматическое управление режимом работы дозирующего устройства путем формирования обратной связи в ждущем режиме в вейвлет-среде. а) б) Рис. 4. Осциллограммы и их время-частотные карты, соответствующие номинальному (а) и текущему (возмущенному - б) режимам работы дозатора порционного типа а) б) Рис. 5. Осциллограммы и их время-частотные карты, соответствующие номинальному (а) и текущему (возмущенному - б) режимам работы спирального дозатора Разработанная векторная автоматизированная система управления [7] смесеприготовительным агрегатом в вейвлет-среде представлена на рис. 6. Она состоит из N скалярных автономных контуров управления дозаторами и смесителем. Функциональное назначение элементов схемы указано внутри блоков. Блоки 1, 2, 3, 4 формируют собственно объект управления - смесеприготовительный агрегат. Автономный аналого-цифровой скалярный контур управления дозаторами включает блоки 5-14. В блоке 5 производится текущая регистрация сигналов материальных потоков с помощью первичных преобразователей (тензометрических и пьезоэлектрических датчиков). В состав цифрового вейвлет-сегмента скалярного контура входят блоки 7-12. В блоке 10 определяется цифровое значение текущего напряжения электропривода дозатора ДУj - при обращении к базе данных тарировочных характеристик двигателей дозирующих устройств, сформированной в компьютере управления агрегатом. В блоках 11, 12 осуществляется определение корректирующих воздействий на двигатели дозаторов по тарировочной базе Tj(fj)=F(uj) и подача откорректированных напряжений на ЦАП, создающих номинальные «фреймированные» МСМ-карты (рис. 4, 5). Рис. 6. Блок-схема векторной системы управления смесеприготовительным агрегатом в вейвлет-среде: - технологические потоки; - скалярные информационные потоки; - векторные информационные (сигнальные) потоки; Tj - период дозирования (для дозаторов дискретного действия); fj - частота дозирования (для дозаторов непрерывного действия); БИИ - блок исходных ингредиентов; БМД - блок мультидозирования; РЕ - рецептурная емкость; ЛРК - локальный рецикл-канал; БУ - блок управления; ВАСУ - векторная автоматизированная система управления: СМУ - смесительное устройство (СНД) В блоке 13 производится комплекс операций по преобразованию цифрового сигнала (кода) потока на выходе модуля «LPT-порт компьютера управления/ ЦАП» в аналоговое напряжение (управляющее воздействие) на зажимах двигателя дозатора. Аналоговое управляющее воздействие формируется в тиристорном управляемом выпрямителе с использованием принципа широтно-импульсной модуляции. Управление смесительным устройством осуществляется реализацией процедур в блоках 15-18. На входе СНД производится мониторирование степени флуктуаций материалопотока путем анализа коэффициента пульсаций КpdΣ, на выходе ведется непрерывный контроль сглаживающей способности S0(ω) смесителя. Для заданных уровней сглаживания и существующем значении коэффициента локального рецикла klr формируется такое управляющее воздействие на клапан отражающего элемента СНД, которому соответствует меньшее значение klr, создающее такую же величину сглаживающей способности смесителя. Последний режим ведет к увеличению интенсивности выходного потока, что, в свою очередь, снижает время приготовления смеси и повышает производительность смесительного устройства при прежних параметрах режима смесеприготовления, характеризующих качество смеси. Заключение 1. Сделан вывод о необходимости применения аппарата вейвлет-преобразований как одного из основных средств, формирующих базовую математическую платформу анализа и моделирования систем управления процессами в смесительных агрегатах и комплексах непрерывного действия. 2. Для выполнения эффективной обработки регистрируемых материалопотоковых сигналов применен алгоритм вейвлет-поиска соответствия, благодаря которому производится аппроксимация сигналов наборами вейвлет-функций из специализированных время-частотных словарей (тезаурусов). С целью возможности визуального семантического представления одномерных переменных (потоков) последние посредством квадратичного распределения модифицируются в двумерные / трехмерные отображения, которые могут интерпретироваться как динамические текущие спектры материальных потоков в агрегате. Описаны процедуры адаптивной аппроксимации технологических сигналов вейвлет-функциями, генерируемыми в рамках время-частотного избыточного словаря; они позволяют точно (во всем частотном диапазоне материалопотоковых сигналов) и оперативно воспроизводить восстановление одномерных расходовых сигналов в вейвлет-среде. 3. В качестве аппроксимирующего базиса при анализе сигналов используется базис вейвлет-функций Габора с тетрадой параметров τ, s, ξ, φ (соответственно: смещения, масштаба, частоты модуляции и начальной фазы), дающий возможность с высокой точностью восстанавливать реальные анализируемые материалопотоковые сигналы. 4. Для целей визуализации и облегчения управления динамикой одномерных материалопотоковых процессов предложен и используется двухступенчатый подход при их обработке: на первом этапе производится аппроксимация сигналов методом вейвлет-поиска соответствия, на втором - их преобразование в 2D-формат.Специфические особенности непрерывных смесеприготовительных процессов требуют применения нетрадиционных аппаратов вейвлет-преобразований и время-частотных распределений - при реализации процессов обработки технологической информации и формирования условий для управления фрагментами смесительного агрегата в рамках цифровых автоматизированных систем с обратными связями по многомерным координатам. 5. Разработана технология адаптации алгоритма вейвлет-поиска соответствия применительно к обработке нестационарных по частоте материалопотоковых сигналов; это позволило эффективно: а) идентифицировать и контролировать специфические режимы работы дозирующих устройств, обусловленные заданной технологией получения смесевых композиций; б) управлять динамикой смесеприготовительного агрегата, используя при этом карту модифицированного сигнала материалопотока (карту Вигнера) в качестве регулируемой двумерной время-частотной координаты, что, в конечном счете, позволяет рационализировать процесс производства высококачественных смесей. 6. Алгоритм вейвлет-поиска соответствия с габоровским словарем обеспечивает наиболее точное описание время-частотных структур среди доступных в настоящее время методов, так как он описывает представленные в сигнале структуры в терминах их времени возникновения, частотного и временного охвата, амплитуды и фазы - с разрешением, которое может быть настроено до теоретических пределов. 7. Поскольку смесительная аппаратура (смеситель, каналы направленной организации потока) аппроксимированы звеньями первого / второго порядков и обладают свойствами низкочастотного фильтра, то определенные высокочастотные составляющие сигнала материалопотока, имеющие место на этапе дозирования, могут отсутствовать в спектре сигнала на выходе СМПА и / или в локальном, глобальном рецикл-каналах, и, следовательно, их невозможно использовать для целей управления, т.е. сигналы материалопотока, регистрируемые на выходе смесительных аппаратов, пригодны для управления дозирующими устройствами, работающими с относительно низкими частотами в спектре материалопотоковых сигналов. 8. На базе вейвлет-анализа разработан способ мониторингового управления производственным процессом смесеприготовления, обеспечивающий непрерывный анализ смесей, управление процессом дозирования и динамикой смесительного узла с воздействием на вектор параметров последнего. Автоматизированное управление соответствующим фрагментом реализуется путем дискретного формирования обратной связи в ждущем режиме. Такие процедуры управления могут использоваться в SCADA-подобных системах. 9. Сформированы концепции создания систем управления с обратными связями на основе многомерных время-частотных переменных состояния. Это позволяет - совместно с функцией управления - реализовывать функцию визуального мониторинга текущих режимов работы агрегата, обогащающую семантическую основу процедуры управления.