Способ управления режимами дозирования и смешивания с применением вейвлет-преобразований Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

УДК 621.929.2/.9:511.336.2

К.А. Дацук, Е.И. Князьков, Д.Б. Федосенков, А.А. Ямпольский, Б.А. Федосенков

СПОСОБ УПРАВЛЕНИЯ РЕЖИМАМИ ДОЗИРОВАНИЯ И СМЕШИВАНИЯ С ПРИМЕНЕНИЕМ ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЙ

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

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

Введение

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

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

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

Научная новизна состоит в преобразовании Ю-материалопотоковых сигналов в 2В/3Б-формат в виде времячастотных распределений, полученных средствами вейвлет-преобразований. На практике такое преобразование позволяет повысить эффек-

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

Методы исследования

При использовании алгоритма вейвлет-поиска соответствия [1] на основе определенной базисной функции g(t, т, s, £) генерируется информационная среда некоторого вейвлет-тезауруса (словаря) в виде семейства вейвлет-функций путем масштабирования (s), перевода (т) и модуляции (£) последних. В качестве базисной функции была выбрана функция Габора (рис. 1а), наилучшим образом соответствующая локальным структурам анализируемых сигналов расхода:

gy (t) = S -U 5g

't-т'

s

exp( j&)

(1)

где

g

- функция Гаусса; у =

V ^

- индекс параметров вейвлета Г абора.

На рис. 1б представлен такой сигнал расхода (на выходе порционного дозатора) и аппроксимирующий его вейвлет Габора (см. рис. 1а).

.т. -4 т-

а)

...... ......

_

■:

Рис. 1. Осциллограммы сигналов: а) адаптированный к процессам дозирования вейвлет; б) сигнал расхода порционного дозатора, спроецированный на вейвлет-функцию Габора

Отметим, что на рис. 1 изображен центрированный сигнал расхода /0(г) :

/0(г) = /(г) - /(Г), (2)

где /(г) - среднее значение постоянной состав-

ляющей сигнала расхода,

‘рег

рег 0

где ‘ - интервал регистрации сигнала.

На рис. 2 представлен реальный материалопотоковый сигнал спирального дозатора (а) и сигнал расхода в центрированной форме (б).

Опишем процедуру реконструкции (восстановления) исходного сигнала на основе его адаптивной аппроксимации рядом вейвлет-функций.

Когда словарь чрезмерно избыточен, поиск векторов наилучшего соответствия сигналу или его остаткам может быть лимитирован [2] в основном подсловарем Ва = ) с В . Предположим

что Га - конечный индекс множества, входящего в состав Г так, что для любого сигнала / е Н, где Н - пространство сигнальных функций с финитной энергией:

™р|(/, Яг)Ыа ь’пр\(/, gr}\ (3)

уеГа уеГ

Здесь (^/,g^ - скалярное произведение сигнала

и вейвлета с индексом У .

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

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

^ для итерации с номером п > 0.

Здесь под Кп/ понимается остаточный сигнал после п-й итерации в рамках выполняемой аппроксимации. В подсловаре Ва ищется такой элемент gr

(словарная волновая форма с индексом уп ), что

Рис. 2. Центрирование сигнала расхода

|{«"/,&")| = иир\{Кп/,£у)\. (4)

Для того чтобы найти словарный элемент, соответствующий / лучше, чем , по методу Ньютона [3] ищется индекс уп по соседству с уп в Г,

где

/,%"

достигает локального максимума. Оче-

видно, что

I («■/ ■ gyn) I г I (*'/■ gy. > |г а™р| (*'/■ gy) |. (5)

Отметим, что функция выбора определяется косвенным образом по вышеуказанной концепции двойного поиска. После выбора вектора g^ вычисляется скалярное произведение нового остатка Яп+1 / с любым ^ е В , поэтому имеем

ЯП+7,g^ = ^Яп/^у)-(Яя/,gЛg■tn^у) . (6)

Так как (яп/,^^ и (яп/,^^ предварительно были сохранены, то обновление состоит лишь в вычислении . Как правило, словари построе-

ны таким образом, что это скалярное произведение получается при малом числе операций. Количество времени на разложение остатков данного сигнала / зависит от желаемой степени точности аппроксимации £. Количество итераций - это минимум р , такой, что остаток вектора после р-й итерации:

р-1 . .

/-К К-"/,Я")Я" "=0

/II.

Из

условия

консервации

(7)

энергии

А\2 =К| (Я"/’^") +|№ЯА следует выра-

жение

р-1,

"=0

/II2-К (*"/-8")'

(8)

Поскольку остаток Яп/ не вычисляется, то для остановки разложения на каждой итерации проверяется справедливость выражения (7). Иными словами, итеративный процесс прекращается при выполнении условия ||^п/|| — £/(г)||, то есть когда остаточный вектор Яп/ после п-й итерации становится меньше заданной ошибки аппроксимации £/(*А сигнала. Количество итераций р зависит

от скорости уменьшения 1йп/ . Оно может изме-

няться в широких пределах в зависимости от вида сигналов.

Уменьшение остаточного вектора ||^"/|| зависит от корреляции между остатками и элементами словаря. Определим корреляционное отношение функции / є Н относительно Б как

Л{/ ) = sup

уєГ

/,ё у

І=1ІЇ

2 \0,5

К(/2)

-1,Ы

-. (9)

Для любого / е Н , Л{/) - больше, чем некоторая положительная постоянная величина. Кроме того, доказано [2], что ||^п/1| уменьшается экспоненциально и что чем ниже корреляционные отношения конкретного сигнала / и его остатков, тем медленнее затухает норма последних. Если сигнал / представляет собой совокупность высокоэнергетических компонент, принадлежащих словарю, то корреляционные отношения / и его остатков высоки, поэтому их норма быстро уменьшается. Эти высокоэнергетические компоненты можно рассматривать как «когерентные структуры» относительно словаря. Если же остатки сигнала имеют низкие корреляционные отношения, то их норма уменьшается медленно, а сам сигнал / необходимо разложить по многим словарным векторам для получения хорошей аппроксимации. Это означает, что информация сигнала / «размыта» по словарю.

Чтобы осуществлять на каждой итерации быстрое вычисление скалярного произведения, воспользуемся аналитическим выражением для определения скалярного произведения двух дискретных (решетчатых) габоровских функций. Для этого рассмотрим сигналы в виде решетчатых функций / и ка с периодом Ы:

^ 4 (10)

/а (") = К1 (" + )

д=—<я ,+да

К(" )= К н(" + яЫ ).

(11)

я=-да ,+да

Сформируем скалярное произведение двух решетчатых сигналов:

^’К) = ^ {п)К {п) =

п=1,М

К

п=1,Ы

К / (п+ч^) К к ( п+шЫ )

д=-¥,+¥ ш=-¥,+¥

= К ї (р) К к (р+шЫ)•

р=-<х>,+<х> ш=-ю,+ю

"=0

Таким образом, перепишем так:

Л. К) =

К К-\Л(іУ1 (і + шЫ )£(і - р ) йі

(13)

Ш=-ГО, +ГО р=-да, +ГО -го

По формуле Пуассона

К5(і -Р)= КехР(І2жЧР). (14)

р=-ГО ,+ГО

^=-ГО,+ГО

Подставив (14) в (13), получим выражение для быстрого расчета скалярного произведения двух решетчатых функций сигналов:

(Л>К) = х Л {п)К {п) =

К К ІЛ Iі)к Iі+шЫ) єхр (і2лчр)йі

(15)

ш=-го,+го р=-го,+го -го

Следовательно, для двух определенных вейвлет-функций g{t) словаря Габора с индексами

"1 =\ Яі>Рі

2жк, і ( 2жк7 і

11 и "2 = 1 ^,р2, --2 I мож-

N

N

но записать скалярное произведение двух решетчатых габоровских сигналов:

' 2^1^2

£"2) =

ехр

2ж(к2 - к )^

< ^ К.

ш=-го,+го р=-го,+го

-№2

ехр

N

-ж

(р2 - рі + ШЫ)2

2 2 *1 + ^2

(16)

ехр

-ж

х ехр

]■

(к2- кі + дЫ)2

N2 (^2 + *22)

(к2 - к1 + ЧЫ)(р2 - рі + ШЫ)

Если вейвлеты В и Ву представляют собой

либо функции Дирака (т.е. )), либо комплексные

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

При задании точности вычислений е для каждого вейвлета Габора gy^ существуют о{ы^1о££\) других таких векторов g , для которых скалярным

произведением , gу^| пренебречь нельзя. Мож-

но показать при этом [2], что выражение (16) потребует о{ы\log^^2) операций для расчета скалярного

произведения некоторого вейвлет-атома ^ со всеми другими атомами ^ из пространства Га.

Полная вычислительная сложность расчетов на одной итерации в алгоритме вейвлет-поиска соответствия имеет порядок O{NlogN). Таким образом, при табуляции гауссовых и гармонических функций на каждую итерацию требуется примерно столько же времени работы центрального процессорного устройства, сколько уходит на реализацию быстрого преобразования Фурье [3] при анализе решетчатого сигнала с N отсчетами.

Рассмотрим подробнее технологию итеративной

процедуры. На первом шаге выбирается вектор ,

дающий наибольшее скалярное произведение с сигналом /(г):

' ' ’ (17)

1 = (/,Е")Е"0 + К11.

Затем остаточный вектор Я1, полученный после аппроксимации / в направлении В , раскладывается подобным же образом. Итеративная процедура повторяется по последующим получаемым остаточным векторам Я"/ , " = 0,Я . Таким образом, на каждой итерации (фазе итеративного процесса) выбирается только одна вейвлет-функция Вг. , где

(і+1) - номер итерации для индекса " ; при этом отбираемый вейвлет Вг. вводится в аппроксимативное выражение (17) по критерию максимума скалярного произведения выбираемой из словаря вейвлет-функции £ и остаточного вектора Я1+1 / на

(і +1) -й итерации. Следовательно,

Я"/ =< Я"/,£у> 8- Я"+1/,

(18)

где Я/ = Я0/ = /(і) - остаточный вектор на нулевой итерации, равный исходному анализируемому сигналу / (і).

Таким способом сигнал раскладывается в сумму времячастотных атомов, выбранных оптимально соответствующими остаткам сигнала (см. схему на рис. 3):

/ =К< Я"/,£" > В" + Я"'1/. (19)

"=0,я

я-1 1-я *7 2-я К-/ К / К/ ІГ4/

жиг/ ► ► г-я т-я

итерация итерация итерация итерация

уровень остаточных векторов і?"/, п = 0,т

І.

1.

ку —►

► (ДМ/.&м>

уровень коэффициентов аппроксимации (вейвлет-коэффициентов)

Рис. 3. Двухуровневая схема адаптивной аппроксимации сигнала

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

В [2] доказано, что данная процедура сходится к /(г) при числе волновых форм избыточного словаря т ^ да, то есть

їія\Яг‘

= 0.

(20)

где ЯЯ/ - остаточный вектор на я-й итерации; ||ЯЯ/|| - норма остаточного вектора в виде решетчатой функции

/ [(кГ )0 ] = К/(кТ \8(і - кТ),

(кТ)»

к=0Л

где {К! )0 - дискретное время для центрированного решетчатого сигнала / [(кТ)0 ]; Ь - размерность

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

выбираемых в направлении индексного вектора уп , определяется как

/(і) = К< Я"/,В" > 8"

(21)

"=0,я—го

где вектор номера вейвлет-функции п также решетчатая функция с единичным шагом.

При этом квадрат нормы сигнала равен

ІІ/ІІ2 = К< Я"/, В у

>

(22)

"=0

Это значит, что сумма квадратов дискрет решетчатого сигнала /(г) равна сумме квадратов скалярных

произведений остаточных векторов Яп/, разложенных в направлении вейвлет-функций g и бау п

зисных функций g .

у п

Результаты и их обсуждение

Пример технологии алгоритма вейвлет-поиска соответствия представлен на рис. 4.

Здесь приведена последовательная процедура обработки исходного сигнала путем итеративного отбора из словаря Б вейвлет-функций Габора и получения остаточных сигналов после 1-й, 2-й и 3-й фаз итеративного процесса. Римскими цифрами I, II, III, IV обозначены фрагменты сигнала /(і), которые аппроксимируются тремя времячастотными атомами Габора, причем последовательность аппроксимации с учетом энергетической плотности этих фрагментов и их остатков такова: II ——^Ш. Фрагмент IV не аппроксимируется, поскольку остаточный сигнал на рис. 4г меньше по норме заданной

ошибки аппроксимации «|/(і)||.

На рис. 5 отображена процедура адаптивной аппроксимации в виде алгоритмической блок-схемы.

В блоке 1 (см. рис. 5) выполняется нормирование вектора сигнала /(і), иными словами, производится операция приведения нормы сигнала /(і) к единичному значению. Следовательно, для решетчатого центрированного сигшла/[(кТ)0] = К/(кТ)о-3({ - кТ),

к=0^

норма которого ||у[(кТ)0]|| Ф 1, нормализация ведет к операции вида --------1------її(кТ)п1. В остальных

\\/[(кТ)0]\\

блоках операции соответствуют выражениям (17)...(21). Заметим, что в блоках 6, 7, 8 ошибка аппроксимации задается также в виде нормированного

вектора «|/(і)||.

"

2

Рис. 4. Технология получения остаточных векторов по формуле (18): а) исходный сигнал; б), в), г) - остаточный сигнал после 1-й, 2-й и 3-й итерации

В блоке 8 происходит сравнение нормы остаточного сигнального вектора ||^п/|| с нормой ошибки

£/(г)\ . Если логическое условие не выполняется

(выход «нет»), то итеративный процесс продолжается. Если условие выполняется (выход «да»), то итеративный процесс прекращается, последний вейвлет занимает свое место в векторе-аппроксиманте, который в блоке 9 фиксируется в виде восстановленного сигнала /(г) = X а ' 8((), где а1 - вейвлет-

1=Ты

коэффициенты.

Опишем численное применение метода вейвлет-поиска соответствия (ВПС) при использовании в процессе анализа материалопотоковых сигналов избыточного вещественного словаря Габора в качестве функционального тезауруса. Для некоторого

. 2жк

индекса у = | т,Б,----

п

и начальной фазы ф , ле-

[0, 2ж], вещественные дискретные времячастотные атомы связаны с комплексными атомарными функциями соотношением

жащей в диапазоне

В(" ф) = К("ф) ^хрііф^" + ехР(- іф\" ).

Рис. 5. Алгоритмическая блок-схема процедуры вейвлет-поиска соответствия

Можно показать, что константа нормирования К {у,фф имеет вид:

(24)

где

ЄХ]

К(",ф) = ^21 + Яе{єхрВ" ,В^)) °'5, Яеіехр(і2фУ\ В" в-)) - вещественная часть

(Р(І2ф\8- &-).

Тогда для остаточного многочлена сигнала Я"/ модуль скалярного произведения равен

(Я"/,В (",ф))| = К (у,ф) Яе(е-“{ Я"/,В")]. (25)

Задавая ф , равную комплексной фазе фу ска-

(23)

лярного произведения

(Я"/,В^

получим

Яе(ехР (- іФ"Xяnf, В—)]| = \{Я"f, В")

(26)

Таким образом осуществляется поиск такого индекса —и , при котором соответствующим вейвлетом

В максимизируется модуль скалярного произведе-'Я"/,В—

ния {К / ,%„) для индекса у в подпространстве Га .

При этом с помощью алгоритма Ньютона [3] исследуется окрестность словарного индекса уи в простран-

2лк,„

стве Г для индекса " = I Тп,$п,

N

при

котором модуль скалярного произведения (К1 / ,ВУ

достигает локального максимума (формула (5)).

Поскольку остаточный многочлен анализируемого сигнала на (п+1)-й итерации равен

Я"

/ = Я"/-(Я"/,В—"ф„)В)

(27)

где g{Уп,ФУ„ ) - словарная функция с индексом уп и начальной фазой ф , то при следующей итерации следует рассчитать для некоторого словарного индекса у е Га скалярное произведение

(яи+1/, Ву) = {Яя/, Ву)-

-(Яя/, В {Уп А))( 8уп ,фу)

Следовательно, в итоге оценивается и подсчиты вается скалярное произведение следующего вида:

(28)

ё— =

К ("п А)

2 .(29)

ехР { А )• (В^ , Ву ) + ехр {-]фуп ) • (Вуп, ёу)

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

(т) к нижнему (ф) - диаграммой (б).

Для получения времячастотного отображения сигнала, реконструированного в ходе ВПС, воспользуемся распределением Вигнера.

Перекрестное распределение Вигнера [1, 2] двух

функций / {г) и определяется как

Ж [ /, Н](і,а) = і-ТО л Сі+г)к Сі-г] ехр (-]ат^йт

Распределение Вигнера функции / {г):

ж/ {г,

)=ж [/,/ ]г,Ф). Так как оно квадратичное, то в итоге получаем

ш/{—2|(«п/,Вуп)Г^ (г,-®) + 2 2 («п/,Вуп){кшГ,Ву.) (31)

п=0 п=0 т=0;т*п

Ш [Ву,Ву, ]М

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

2

Е/(г,ф) = 2|<>| WgIJt,ф). (32)

п=0

На основании известных свойств распределения Вигнера, таких как растяжение и смещение, а также уравнения времячастотного атома получаем для

У = {?,*,£)

Wgу(t,©) = (33)

и, следовательно,

V п

. (34)

У

Распределение Вигнера также удовлетворяет условию

+го +го

І І Жё (і,т)dtdm = \\g\f = 1, (35)

поэтому в соответствии с уравнением локализации энергии имеем

+ТО+ТО

11Е/(1,т)ё1ёт = ||/||2. (36)

-го -го

ТО-ТО

б)

Рис. 6. Задание значений параметров вейвлетов Габора

Функцию Е/ (г, ф) можно интерпретировать как

плотность энергии / на времячастотной плоскости

(г, Ф). В отличие от распределений класса Коэна [4] данный вид распределения не содержит перекрестных членов. Оно также остается положительным,

если (г ,ф) - положительно, т.е. когда g(t) -гауссова функция. Чтобы получить разложение с вещественными коэффициентами, когда сигнал

/ (г) также вещественен, необходимо использовать словари вещественных времячастотных атомов. Для любого у = (т^,£) с ^ ^ 0 и любой фазы ф е [0,2ж] определим

ё(",ф) = К(—ф ё

і - г

+ ф). (37)

Постоянную К(у,ф) устанавливаем такой, чтобы задать ||в(у, ф)| = 1. Начальная фаза ф , скрытая

ранее в комплексных числах, теперь отчетливо прослеживается как один из четырех параметров вещественных атомов. Словарь вещественных времяча-стотных атомов определяется выражением

° = ($(у ф)\у,фег40.2х]. Поиск соответствия, выполненный с данным словарем, разлагает любой вещественный сигнал / (') в ряд:

Л () = 2 ( Яп/, В (Уп, фп)) В (Уп, фпХг) ’ (38)

п=0

где индексы (Уп , фп )=(Jn, 8п , £п , фп ) выбираются для наилучшего соответствия остаткам / .

Для любого у = (г,£,^) вещественные атомы связаны с комплексными атомами следующим выражением:

В (У, фф() = К (у,ф) (е‘фВу(г) + е~‘фВу-(г^ (39)

где у~ = (т,Ь',-£).

Времячастотное распределение энергии вещественного сигнала /(г) получается из ВПС-разложения при суммировании распределений Вигнера комплексных атомов:

Е/ (г,ф) =

+“.,2 и \ .(40)

= Ц(Яп/,В(Упф, })| т(Шв,М'Ш-М

п=0 2

При подстановке в это выражение уравнения (33) получим

V V п

2 п=0

Л (т-Еп )

і - г

V *п

■(т + Еп )

.(41)

Это распределение также удовлетворяет характеристике Е/ (г, ф) как плотности энергии (34).

Данное энергетическое распределение используется для визуализации итоговой сигнальной информации [5] в автоматизированной системе управления смесеприготовительным агрегатом.

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

*

1. Федосенков, Б.А. Управление смесеприготовительным агрегатом на базе вейвлет-преобразований / Б.А. Федосенков, А.С. Назимов, А.В. Шебуков // Автоматизация и современные технологии. Автоматизация научно-исследовательских и производственных процессов. - 2004. - № 8. - C. 7-13.

2. Mallat, S. Matching pursuit with time-frequency dictionaries / S. Mallat and Z. Zhang // IEEE Transactions on Signal Processing. - 1993. - Vol. 41, № 12. - Р. 3397-3415.

3. Корн, Г. Справочник по математике (для научных работников и инженеров) / Г. Корн, Т. Корн. - М.: Наука, 1977. -

832 с.

4. Cohen, L. Time-frequency distributions - A review / L. Cohen // Proc. IEEE. - 1989. - Vol. 77, № 7. - P. 941-981.

5. Использование технологии вейвлет-мониторирования как средства управления динамикой стационарных и нестационарных процессов / К.А. Дацук, Е.Н. Карнадуд, Д.Б. Федосенков и др. // Исследование, разработка и применение высоких технологий в промышленности: Сб. тр. 9-й междунар. науч.-практ. конф., Санкт-Петербург, 22-23.04.2010. - Т. 3. - С. 256-258.

ГОУ ВПО «Кемеровский технологический институт пищевой промышленности», 650056, Россия, г. Кемерово, б-р Строителей, 47.

Тел./факс: (3842) 73-40-40 е-mail: [email protected]

SUMMARY

K.A. Datsuk, E.I. Kniazkov, D.B. Fedosenkov, A.A. Yampolsky, B.A. Fedosenkov The method of controlling the dosing and mixing processes on the basis of wavelet transforms

This article presents an approach which can be used for the representation of material flow processes while dealing with the problem of dosing and mixing regimes control. In this technique the matching pursuit algorithm has been adapted for approximating the flow signals registered at the output of a dosing devices block. As a result of this procedure, the wavelet approximated signals are afterwards transformed into the so-called Wigner distributions (maps) which are the material energy flow signals density distributions over time-frequency space. The advantages, limitations, and possible extensions of the approach considered are outlined. An example on the reconstruction of a flow signal has been reviewed.

Dry material dosing, wavelet function, Wigner distribution, Gabor wavelet parameters, wavelet matching pursuit algorithm.

Kemerovo Institute of Food Science and Technology 47, Boulevard Stroiteley, Kemerovo, 650056, Russia Phone/Fax: +7(3842) 73-40-40 e-mail: [email protected]