Методы интеллектуального управления процессами обжарки кофе
На правах рукописи
БОЛДИНА ЕЛЕНА АЛЕКСАНДРОВНА
МЕТОДЫ ИНТЕЛЛЕКТУАЛЬНОГО УПРАВЛЕНИЯ
ПРОЦЕССАМИ ОБЖАРКИ КОФЕ
Специальность 05.13.06 – «Автоматизация и управление технологическими
процессами и производствами (в пищевой промышленности)»
АВТОРЕФЕРАТ
диссертации на соискание ученой степени
кандидата технических наук
Москва – 2010
Работа выполнена на кафедре электрооборудования и автоматики
Российского государственного аграрного заочного университета (РГАЗУ)
Научный руководитель: доктор технических наук, профессор
Солдатов Виктор Владимирович
Официальные оппоненты: доктор технических наук, профессор
Судник Юрий Александрович;
кандидат технических наук, доцент
Татаринов Александр Владимирович
Ведущая организация: НИИ «Агроприбор»
Защита состоится «____»_________ 2010 г. в ______ часов на заседании Совета по защите докторских и кандидатских диссертаций Д 212.148.02 при ГОУ ВПО «Московский государственный университет пищевых производств» по адресу: 125080, Москва, Волоколамское шоссе, 11, в ауд. _____
С диссертацией можно ознакомиться в библиотеке Московского государственного университета пищевых производств.
Автореферат разослан «______» ____________ 2010 г.
Ученый секретарь Совета,
к.т.н., доцент Н.О. Воронина
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность проблемы. В России кофе – второй по популярности после чая напиток. Его пьют 83% россиян, в основном – в крупных городах (на долю Москвы и Санкт-Петербурга приходится до 30% всех продаж кофе в России). Так, в Москве потребление кофе в год составляет 8 – 12 тыс. тонн (т.е. около 1 кг в год), а в среднем по России на душу населения приходится около 0,6 кг. По данным социологических исследований, Россия занимает четвертое место в мире по потреблению кофе на душу населения – 0,6 кг в год после Финляндии (11,8 кг), США (3,9 кг) и Англии (2,4 кг).
При переработке кофе важное значение имеет его обжарка, представляющая собой сложнейший процесс, в котором отсутствует определенная стандартная технология производства и, следовательно, требуются глубокие практические знания и умения специалистов.
Интенсивное тепло в процессе обжарки активизирует сложную химическую реакцию. Кофейное зерно теряет около 20% влаги, а на поверхности зерна начинают выделяться углеводы и эфирные масла, содержащиеся внутри, что и придает кофе яркий аромат и вкус. Если же процесс обжарки продолжить, сахар начинает обугливаться, из-за чего постепенно теряются специфические свойства зерна и проявляется горечь обугленного сахара. При повышении температуры начинается разрушение целого ряда химических веществ, влияющих на вкус кофе, появляется горький привкус, пропадает свойственный каждому сорту кофе ароматический букет. Вкус разных сортов кофе начинает сближаться. Чтобы этого не происходило, необходимы строгий контроль температуры и длительности обжарки.
Поэтому целесообразно разработать методы интеллектуального управления температурой обжарки кофе, эффективные при отсутствии достоверной информации о действующих на объект возмущениях, а также уточнить в процессе обжарки каждой партии кофе ее оптимальные температуру и длительность, с целью достижения наилучших показателей качества готового продукта.
Объект исследования – технологические процессы обжарки кофе.
Предмет исследования – методы управления процессами обжарки кофе, эффективные в условиях информационной неопределенности.
Цель и задачи исследования. Цель настоящей работы заключается в разработке методов интеллектуального управления обжаркой кофе, эффективных в условиях информационной неопределенности относительно динамических характеристик объекта и статистических характеристик возмущающих воздействий.
Для достижения этой цели поставлены следующие задачи:
1. Разработать метод робастного управления теплообменными процессами в ростере, эффективный в условиях информационной неопределенности относительно статистических характеристик возмущающих воздействий и наличии погрешностей параметров передаточных функций объекта.
2. Разработать методы коррекции управления по отношению к задающим и контролируемым возмущающим воздействиям.
3. Разработать метод расчета параметров настройки регуляторов в системах с неточно заданными параметрами объекта.
4. Разработать метод идентификации систем управления обжаркой кофе на основе нечеткого логического вывода с использованием модели Мамдани.
5. Разработать метод идентификации систем управления обжаркой кофе на основе нечеткого логического вывода с использованием модели Сугэно.
6. Разработать метод синтеза нечетких регуляторов, с учетом существующих технических и технологических ограничений.
Методическая база и методы исследования. Решение поставленных задач проведено на основе применения теорий автоматического управления, вероятностей и математической статистики, линейной алгебры, оптимизации, интегральных преобразований, дифференциальных уравнений, а также методов системного и функционального анализов.
Достоверность теоретических положений подтвердилась проверкой полученных результатов на компьютерных имитационных моделях, а также испытаниями разработанных методов интеллектуального управления на действующем оборудовании предприятия пищевой промышленности ООО "КАФФА ИНДАСТРИЗ" (г. Москва).
Научная новизна исследований заключается в следующем:
1. Разработан метод робастного управления теплообменными процессами в ростере с использованием интегрального квадратичного критерия, эффективный в условиях информационной неопределенности относительно статистических характеристик возмущающих воздействий.
2. Разработаны методы коррекции управления по отношению к задающим и контролируемым возмущающим воздействиям.
3. Разработан метод расчета параметров настройки регуляторов в системах с неточно заданными параметрами объекта.
4. Разработан метод идентификации систем управления обжаркой кофе на основе нечеткого логического вывода с использованием моделей Мамдани и Сугэно.
5. Разработан метод синтеза нечетких регуляторов с учетом существующих технических и технологических ограничений.
Внедрение результатов исследований. Результаты исследований использованы:
1. При создании АСУТП обжаркой кофе на предприятии пищевой промышленности ООО "КАФФА ИНДАСТРИС" (г. Москва).
2. В учебном процессе, т.е. в методических пособиях, лабораторных работах и лекциях для студентов Российского государственного аграрного заочного университета (РГАЗУ).
Апробация. Основные положения и результаты исследований по теме диссертации доложены, обсуждены и одобрены на:
Всероссийской научно-практической конференции «Интеллект. Инновации. Информация. Инвестиции. Инфраструктура», г. Москва, 2010 г.;
II Ежегодной Всероссийской научно-практической конференции «Перспективы развития информационных технологий», г. Новосибирск, 2010 г.
Публикации. Теме диссертации посвящены 8 научных работ, в которых изложены ее основные положения.
Структура и объем работы. Диссертация состоит из введения, 4 глав, общих выводов, списка литературы и приложения. Она изложена на 164 страницах, включая 7 таблиц, 52 рисунка, список литературы из 165 наименований и два приложения.
На защиту выносятся следующие основные положения:
1. Метод робастного управления температурой внутренней среды машины для обжарки кофе (ростера).
2. Методы коррекции управления по отношению к задающим и контролируемым возмущающим воздействиям.
3. Метод синтеза нечетких ПИД-регуляторов с учетом существующих технических и технологических ограничений.
4. Метод идентификации систем управления обжаркой кофе на основе нечеткого логического вывода.
Совокупность сформулированных и обоснованных научных положений, а также результаты их практической реализации и внедрения в производство представляют собой решение актуальной задачи по разработке методов интеллектуального управления обжаркой кофе в условиях информационной неопределенности относительно динамических характеристик объекта и статистических характеристик возмущающих воздействий.
СОДЕРЖАНИЕ РАБОТЫ
Во введении обоснована актуальность решаемой проблемы, изложено краткое содержание глав диссертации, приведены основные положения, выносимые на защиту.
В первой главе анализируются литературные источники, рассматриваются особенности управления обжаркой кофе в обжарочных машинах (ростерах) фирмы "Пробат".
Ростер работает по принципу пропускания через зерна кофе горячего воздуха, уменьшая, таким образом, риск подгорания внешней оболочки. Воздух нагревается с помощью газовых или масляных горелок. Для удаления пыли из отводимого воздуха используются циклоны.
На рис. 1 представлен ростер типа Sirocco (Пробат UG 15), состоявший из котла или вращающихся цилиндров (барабанов), через которые продувают нагретый воздух, вызывающий турбулентность и, как следствие, – непрерывный нагрев зерен кофе. После обжарки зерна кофе выгружают на металлический перфорированный лист, обдуваемый холодным воздухом.
Рис. 1. Ростер периодического действия типа Sirocco.
Во время обжаривания необходимо воздействовать на зерна кофе большим количеством тепла по возможности равномерно и за относительно короткий промежуток времени. Каждое отдельное зерно должно иметь равномерный коричневый цвет, интенсивность которого заранее определяется программой, и быть без единого пятнышка, являющегося показателем того, что зерно пережарено. Поэтому во время обжарки важно обеспечить весьма точное поддержание в ростере заданной температуры нагретого воздуха, чтобы сохранить ароматический букет у обжаренного кофе. С этой целью анализируется эффективность различных методов управления технологическими процессами. Рассматриваются также критерии управления и ограничения при их оптимизации. Представлен обзор программного обеспечения SCADA, предназначенного для создания программируемых систем отображения информации о технологических процессах в реальном масштабе времени. В заключение главы даны выводы, определяющие конкретные задачи диссертационной работы.
Во второй главе рассматриваются методы робастного управления температурой внутренней среды ростера.
Передаточную функцию ростера, как объекта управления температурой его внутренней среды, представим в виде
, (1)
где – комплексная переменная.
Для определения значений параметров, входящих в выражение (1) для ростера известной немецкой фирмы «Пробат» на основе экспериментальных данных была получена усредненная кривая разгона, представленная на рис. 2.
Рис. 2. Усредненная кривая разгона ростера по каналу
"расход газа – температура воздуха на выходе из ростера".
Воспользовавшись полученной кривой разгона исследуемого объекта для параметров передаточной функции (1) определим следующие численные значения:
; ; . (2)
Для управления температурой внутренней среды ростера был разработан метод, основанный на оптимизации интегрального квадратичного критерия в условиях информационной неопределенности относительно характеристик действующих на объект возмущений.
Интегральный квадратичный критерий , получил широкое распространение ввиду своей простоты в математическом отношении, а также благодаря тому, что во многих практических задачах управления он является удовлетворительной мерой успешности их решения.
Обозначим ошибку управления системы как , тогда интегральный квадратичный критерий определяется следующим выражением:
. (3)
Поскольку в случае линейных стационарных систем удобно использовать частотные методы анализа и синтеза, то воспользуемся известным соотношением
, (4)
где – Фурье-изображение ошибки управления ; – циклическая частота; – мнимая единица.
Для односвязной системы, функциональная схема которой представлена на рис. 3, выполняется равенство
, (5)
где – Фурье-изображение приведенного к выходу объекта возмущающего воздействия с неизвестными характеристиками; и – комплексные частотные характеристики (КЧХ) регулятора и объекта.
Рис. 3. Односвязная система управления: – управляющее воздействие; – сигнал задания; – сигнал ошибки управления; – возмущение, приведенное к выходу объекта; – возмущение действующее на объект; – выходной сигнал объекта в отсутствии возмущений; – результирующий выходной сигнал объекта в присутствии возмущений.
C учетом равенства
, (6)
а также соотношений (4) и (5)) выражение (3) представим в виде
. (7)
Ввиду выполнения неравенства
,
для выражения (7) получим мажорирующую оценку
, (8)
где – амплитудно-частотная характеристика (АЧХ) замкнутой системы, соответствующая каналу передачи сигналов .
Поскольку , то из неравенства (8) следует, что
, (9)
где – показатель колебательности замкнутой системы,
На основании оценок (8) и (9) установим, что при действии на рассматриваемую систему возмущения с неизвестными характеристиками для минимизации критерия целесообразно выполнить требование
(), (10)
где – резонансная частота системы, при которой выполняется равенство .
Примем во внимание, что из равенства (6) можно получить неравенство
, (11)
где – АЧХ замкнутой системы, соответствующая каналу передачи сигналов .
Поскольку из неравенства (11) вытекает, что
,
то требование (8) можно заменить следующим:
(). (12)
Исследования показали, что значения функции минимизируются в интервале частот ( – граничная частота данного интервала), если доминирующие, т.е. оказывающие наибольшее влияние на качество управления, корни характеристического уравнения замкнутой системы расположены на комплексной плоскости следующим образом.
, ; (13)
, ; (14)
,, (15)
где – показатель относительного демпфирования свободного движения замкнутой системы, а – собственная частота ее свободных колебаний.
Отметим, что расположения корней (13), (14) и (15) соответствуют системам с ПИ, ПИД и многопараметрическими регуляторами, причем передаточные функции последних задаются выражением
,
где – коэффициент передачи регулятора; и – постоянные времени дифференцирования и интегрирования соответственно; и – дополнительные постоянные времени.
Интервал частот , в котором при выполнении требований (13) – (15) минимизируются значения функции , определяется при решении уравнения
, (16)
где , а принимает значения 1, 2 или 3 в зависимости от того, какое из требований (17), (18) или (19) выполняется.
Величина в равенствах (18) – (20) и в уравнении (16) выбирается таким образом, чтобы выполнить условие , т.к. в этом случае значения функции минимизируются в заданном интервале частот , а значит, выполняется требование (12).
Следовательно, в случае отсутствия информации о возмущающем воздействии для минимизации мажорирующих оценок (8) и (9) критерия необходимо выбрать значения параметров настройки используемого регулятора таким образом, чтобы доминирующие корни характеристического уравнения замкнутой системы располагались на комплексной плоскости в соответствии с требованиями (13), (14) или (15), в зависимости от типа используемого регулятора, причем значение в равенствах (13) – (15) задаются так, чтобы решение уравнения (21) удовлетворяло условию .
Для оценки эффективности разработанного метода оптимизации применим его к системам, в которых объект (1), (2) управляется ПИ, ПИД и многопараметрическим регулятором.
Графики функций и представлены на рис. 4 и 5.
Рис. 4. АЧХ замкнутой системы по каналу при использовании
следующих регуляторов: 1 – ПИ; 2 – ПИД; 3 – многопараметрический.
Рис. 5. АЧХ замкнутой системы по каналу при использовании
следующих регуляторов: 1 – ПИ; 2 – ПИД; 3 – многопараметрический.
На основании графиков на рис. 4 установим, что при использовании ПИ, ПИД и многопараметрического регуляторов показатель колебательности системы принимает соответственно следующие значения: 1,027; 1,24; 1,567.
Согласно рис. 5 наилучшими фильтрующими свойствами в интервале низких частот, на котором выполняется неравенство , обладает система с многопараметрическим регулятором.
Для оценки эффективности парирования системой возмущающих воздействий построены графики переходных процессов (рис. 6), возникающих при скачкообразном изменении сигнала задания .
Рис. 6. Переходные процессы на выходе системы при использовании
следующих регуляторов: 1 – ПИ; 2 – ПИД; 3 – многопараметрический.
В соответствии с графиками на рис. 6, приходим к выводу, что наиболее быстрое затухание переходных процессов обеспечивается в системе с многопараметрическим регулятором, но при этом перерегулирование примерно на 20 % больше, чем для системы с ПИД-регулятором.
Следовательно, при оптимизации на основе разработанного метода систем управления с неизвестными характеристиками возмущений, целесообразно использовать многопараметрические регуляторы, если не существует достаточно жестких ограничений на величину перерегулирования. В противном случае необходимо применять ПИД или ПИ-регуляторы.
Отметим, что разработанный метод оптимизации робастных систем позволяет выбрать такие компромиссные значения показателей и , при которых достигается правильный баланс между фильтрующими свойствами системы в области низких частот и усилением высокочастотных возмущений в области резонансных частот.
В целях повышения эффективности метода робастного управления разработан также метод расчета параметров настройки регуляторов в системах с неточно заданными параметрами объекта.
Третья глава посвящена разработке методов идентификации систем управления обжаркой кофе на основе нечеткого логического вывода, применение которых позволяет в процессе обжарки оптимизировать ее режимы. При этом на основании анализа экспериментальных исследований установлено существенное влияние продолжительности процесса обжарки кофе и температуры, при которой она производится, на физико-химические показатели готового продукта, такие как, влажность, значение рН, потеря массы кофейных зерен. Причем зависимость этих показателей от продолжительности t и температуры T обжарки пока недостаточно изучена. Для оптимизации режимов обжарки с целью повышения качества выпускаемой продукции использован метод идентификации зависимости показателя pH от продолжительности t и температуры T обжарки кофе с помощью нечеткого логического вывода.
Идентификация, т. е. построение математической модели по результатам наблюдений, является важной и сложной задачей, возникающей при управлении различными объектами. Причем в современной теории идентификации все более важную роль начинают играть методы, привлекающие лингвистическую информацию при построении моделей сложных зависимостей. Одним из наиболее разработанных в инженерном отношении инструментов учета лингвистической информации является теория нечетких множеств и нечеткая логика.
В диссертационной работе использовался метод двухэтапной идентификации нелинейных зависимостей с помощью нечетких баз знаний.
Первый этап – структурная идентификация – представляет собой формирование нечеткой базы знаний, которая грубо отражает взаимосвязь между входами и выходом с помощью лингвистических правил «ЕСЛИ – ТО». Лингвистические правила получаются в результате экстракции нечетких знаний из экспериментальных данных. На втором этапе проводится параметрическая идентификация исследуемой зависимости путем нахождения таких параметров нечеткой базы знаний, которые минимизируют отклонение модельных и экспериментальных результатов.
Типовая структура модели на основе нечеткого логического вывода показана на рис. 7.
Рис. 7.. Типовая структура модели нечеткого логического вывода.
Нечеткая модель содержит следующие блоки:
- фаззификатор, преобразующий фиксированный вектор влияющих факторов в вектор нечетких множеств , необходимых для выполнения нечеткого логического вывода;
- нечеткая база знаний, содержащая информацию о зависимости
в виде лингвистических правил типа «ЕСЛИ – ТО»;
- машина нечеткого логического вывода, которая на основе правил базы знаний определяет значение выходной переменной в виде нечеткого множества , соответствующего нечетким значениям входных переменных ;
- дефаззификатор, преобразующий выходное нечеткое множество в четкое число .
При определении зависимости показателя pH от режимных параметров t и T использовались два вида нечетких моделей – типа Мамдани и типа Сугэно, реализованные в популярном пакете прикладных программ MATLAB. Эти модели отличаются форматом базы знаний и процедурой дефаззификации.
Идентифицируемая зависимость была представлена выборкой экспериментальных данных вида «входы-выход»:
(Tr, t r, pH r), , (17)
где Tr, t r и pH r – соответственно температура, длительность обжарки кофе и pH в -том эксперименте, причем N – объем выборки, т.е. количество выполненных экспериментов.
Задача идентификации состоит в нахождении нечеткой модели , обеспечивающей минимальное значение среднеквадратической невязки:
, (18)
где – значение выхода нечеткой модели при значениях ее входов, заданных режимными параметрами Tr и t r.
Выход нечеткой модели зависит от ее структуры – базы знаний и параметров: функций принадлежностей, реализаций логических операций, метода дефаззификации, а также коэффициентов линейных функций в заключениях правил для модели типа Сугэно. Нахождение структуры и параметров нечеткой модели, обеспечивающих выполнение требования (18) и является задачей идентификации на основе нечеткого логического вывода..
Пакеты Fuzzy Logic Toolbox и Optimization Toolbox, входящие в состав системы MATLAB, использовались для идентификации с помощью нечеткого логического вывода зависимости pH от температуры T и продолжительности обжарки кофе t, причем в качестве исходных данных для идентификации применялась обучающая выборка в форме (17).
Полученные при этом графики идентифицируемой зависимости представлены на рис. 8.
Отметим, что при малых обучающих выборках качество идентификации оказалось существенно выше для модели типа Мамдани. Это объясняется тем, что исходная, основанная на экспертных высказываниях, нечеткая модель уже отражает основные особенности идентифицируемой зависимости. С увеличением объема обучающей выборки лучшее качество идентификации обеспечивает модель типа Сугэно. При больших выборках точность идентификации модели типа Сугэно выше, чем для модели типа Мамдани. Однако после обучения модель типа Мамдани остается прозрачной: ее параметры – функции принадлежности – легко интерпретируются лингвистическими термами. Для моделей типа Сугэно типовое явление - сложность содержательной интерпретации ее параметров.
Рис. 8. Поверхности «входы-выход»:
a) нечеткой модели типа Мамдани; b) нечеткой модели типа Сугэно.
Отметим, что наряду с нечеткими моделями для идентификации зависимости pH от режимных параметров обжарки использовались также традиционные полиномиальные модели. Исследования показали, что при одинаковых объемах обучающей выборки (17) ошибка идентификации для нечетких моделей значительно меньше, чем для полиномиальных.
В целях повышения эффективности управления при отсутствии информации о динамике объекта разработан метод синтеза нечетких регуляторов, с учетом существующих технических и технологических ограничений.
В четвертой главе решаются задачи коррекции управления по отношению к задающим и контролируемым возмущающим воздействиям, применение которой способствует повышению его качества.
Рассмотрим коррекцию сигнала задания. На рис. 9 приводится функциональная схема системы управления с коррекцией сигнала задания.
Рис. 9. Система управления с коррекцией сигнала задания.
Как видно из рис. 9, сигнал поступает на вход специального устройства, называемого корректирующим, которое преобразует его в сигнал , фактически исполняющий роль задающего воздействия.
Рассмотрим метод синтеза передаточной функции корректирующего устройства с учетом расположения доминирующих полюсов системы управления.
Поскольку передаточная функция от входа до выхода, представленной на рис. 9 системы, определяется выражением
, (19)
где и – изображения по Лапласу сигналов и соответственно, то передаточную функцию можно выбрать так, чтобы ее числитель и доминирующие полюсы замкнутой системы взаимно сократились. Тогда в знаменателе передаточной функции (19) доминировать станут полюсы более удаленные от мнимой оси, чем те, которые удалось сократить. Вследствие этого качество переходных процессов в системе при изменениях сигнала резко повышается.
Пусть , – совокупность всех доминирующих полюсов замкнутой системы, тогда синтезируемую передаточную функцию можно представить в виде
, (20)
где и – соответственно коэффициент передачи и постоянная времени корректирующего устройства.
Отметим, что устройство с передаточной функцией (20) технически реализуемо, т.к. степень полинома ее числителя не превышает степени полинома знаменателя.
Выберем значения постоянных и так, чтобы выполнялись следующие требования:
; (21)
. (22)
При выполнении требований (21) и (22) корректирующее устройство мало влияет как на медленно, так и на весьма быстро изменяющиеся сигналы , пропуская их практически без изменений. Зато, все остальные сигналы изменяются им таким образом, чтобы переходный процесс завершился наиболее быстро.
Исходя из выражения (20) и требования (21) получим
, (23)
а с учетом требования (22) имеем
. (24)
Используя выражения (23) и (24) можно определить неизвестные параметры передаточной функции корректирующего устройства (20).
Чтобы убедиться в эффективности данного устройства, построим графики АЧХ каналов передачи задающего воздействия с коррекцией и без нее, которые представлены на рис. 10 и 11.
Как видно из рис. 10, применение синтезированного корректирующего устройства позволило полностью устранить резонансный пик АЧХ канала передачи сигналов , что способствует резкому уменьшению колебательности системы, без снижения ее быстродействия, т.к. значение показателя при этом остается неизменным.
Согласно рис. 11, с помощью корректирующего устройства можно резко улучшить фильтрующие свойства системы по каналу передачи сигналов , т.к. у АЧХ данного канала исчезает резонансный пик, и при этом не возрастают ее значения в области низких частот.
Рис. 10. АЧХ канала передачи сигналов :
1 – с коррекцией; 2 – без коррекции.
Рис. 11. АЧХ канала передачи сигналов :
1 – с коррекцией; 2 – без коррекции.
Повышение качества переходных процессов в системе при использовании корректирующего устройства иллюстрируется рис. 12.
Рис. 12. Переходные процессы при действии единичного ступенчатого возмущения
по каналу передачи сигналов :
1 – процесс с коррекцией; 2 – процесс без коррекции.
Анализируя графики, представленные на рис. 12, приходим к выводу, что применение корректирующего устройства позволяет полностью устранить перерегулирование и колебательность.
Отметим, что устранение колебательности особенно полезно при управлении температурой, т.к. в таких системах обычно отсутствуют устройства, обеспечивающие искусственное охлаждение объекта при превышении температурой заданного значения, а при его естественном охлаждении требуется довольно длительное время для устранения, возникшего перерегулирования.
Рассмотрим теперь коррекцию управления по отношению к контролируемому возмущающему воздействию . В этом случае функциональная схема системы управления представлена на рис. 13.
При решении данной задачи обычно стараются выбрать передаточную функцию корректирующего устройства так, чтобы выполнялось тождество
, (25)
где – передаточная функция объекта по каналу действия возмущения на выходную величину .
Рис. 13. Система с коррекцией управления по отношению к возмущению .
Чтобы определить условия, при которых выполняется тождество (25), получим выражение для передаточной функции , воспользовавшись следующими равенствами:
; (26)
; (27)
; (28)
; (29)
, (30)
где и – изображения по Лапласу выходных сигналов объекта и , обусловленных воздействиями и соответственно; и – передаточные функции объекта по каналам передачи и .
Поскольку
,
то с учетом равенств (26) – (30) имеем
. (31)
Следовательно, тождество (25) выполняется, если
. (32)
Чтобы корректирующее звено, синтезируемое согласно выражению (32), было технически реализуемо, канал передачи объекта должен быть более инерционным, чем его канал передачи . В таком случае степень полинома числителя передаточной функции (32) не превышает степени полинома ее знаменателя.
В тех случаях, когда корректирующее устройство с передаточной функцией (32) технически реализуемо, то его синтез осуществляется в два этапа.
На первом этапе выбирается закон управления, т.е. определяется из тех или иных соображений тип регулятора и значения его параметров настройки.
На втором этапе, когда передаточная функция регулятора уже известна, в соответствии с выражением (32) определяется передаточная функция корректирующего устройства.
Однако на практике синтез корректирующих устройств на основе выражения (32) зачастую оказывается невозможным, т.к. канал объекта, по которому действует возмущение , является менее инерционным, чем канал передачи управляющих воздействий . Поэтому в таких случаях приходится искать другие способы решения задачи коррекции управления по отношению к контролируемым возмущениям.
Рассмотрим один из подходов к решению данной проблемы, основанный на взаимном сокращении нулей числителя и доминирующих полюсов знаменателя передаточной функции , заданной выражением (31).
В этом случае на первом этапе, когда выбирается закон управления, и определяются значения параметров настройки регулятора необходимо использовать методы расчета, базирующиеся на каком-либо из требований: (13), (14) или (15), к расположению доминирующих полюсов замкнутой системы.
Целесообразность использования указанных методов объясняется тем, что они позволяют определить не только значения параметров настройки выбранного регулятора, но и расположение доминирующих полюсов , замкнутой системы.
Поскольку для знаменателя передаточной функции (31) выполняются равенства
, ,
то для взаимного сокращения нулей и доминирующих полюсов системы числитель этой передаточной функции при , также должен обращаться в нуль, т.е.
, . (33)
Поэтому на втором этапе синтеза корректирующего устройства его передаточная функция ищется в виде
, (34)
где , – параметры настройки корректирующего устройства; – заданная постоянная времени корректирующего устройства.
Ввиду того, что степень полинома числителя передаточной функции (34) не превышает степени полинома ее знаменателя, то синтезируемое корректирующее устройство технически реализуемо.
Так как после выполнения первого этапа синтеза расположение доминирующих полюсов системы становится известным, то согласно требованиям (13), (14) и (15) получают определенные значения входящие в них величины и .
В таком случае выберем значение постоянной времени , исходя из условия
, (35)
выполнение которого гарантирует, что корректирующее устройство не будет существенно повышать инерционность канала передачи сигналов , т.к. это привело бы к нежелательному затягиванию возникающих при изменениях возмущающего воздействия переходных процессов на выходе системы.
Пусть
,
тогда требования (33) можно представить в виде
, . (36)
В результате решения системы уравнений (36) при заданном значении определяются значения параметров настройки корректирующего устройства.
Для обоснованного выбора постоянной времени целесообразно исследовать, возникающие на выходе системы управления при действии по каналу единичного ступенчатого возмущения , переходные процессы , графики которых представлены на рис. 14.
Рис. 14. Переходные процессы на выходе системы управления:
1 – без коррекции; 2 – ; 3 – .
Как видно из рис. 14, минимальные значения выброса и длительности переходного процесса достигаются при .
Следовательно, применение разработанного метода коррекции позволяет значительно улучшить показатели качества переходных процессов в системе управления.
Отметим, что разработанные методы интеллектуального управления нашли применение при создании АСУТП обжаркой кофе на предприятии пищевой промышленности ООО «КАФФА ИНДАСТРИЗ" (г. Москва), причем испытания показали, что после внедрения АСУТП процесс обжарки кофе происходит при оптимальных режимах, позволяющих экономить топливно-энергети-ческие ресурсы и сократить потери выпускаемой продукции.
ОБЩИЕ ВЫВОДЫ
1. Разработан метод робастного управления теплообменными процессами в ростере с использованием интегрального квадратичного критерия, эффективный в условиях информационной неопределенности относительно статистических характеристик возмущающих воздействий.
2. Разработаны методы коррекции управления по отношению к задающим и контролируемым возмущающим воздействиям.
3. Разработан метод расчета параметров настройки регуляторов в системах с неточно заданными параметрами объекта.
4. Разработан метод идентификации систем управления обжаркой кофе на основе нечеткого логического вывода с использованием модели Мамдани.
5. Разработан метод идентификации систем управления обжаркой кофе на основе нечеткого логического вывода с использованием модели Сугэно.
6. Разработан метод синтеза нечетких регуляторов, с учетом существующих технических и технологических ограничений.
7. Разработанные методы приняты к использованию в учебном процессе Российского государственного аграрного заочного университета (РГАЗУ) и предприятием пищевой промышленности ООО "КАФФА ИНДАСТРИЗ" г. Москва при внедрении АСУТП обжарки кофе.
ПУБЛИКАЦИИ АВТОРА
Статьи, опубликованные в периодических изданиях, рекомендованных ВАК:
1. Болдина Е.А. Коррекция управления по отношению к задающим воздействиям // Международный научный журнал. – 2010. – № 1. – С. 58–61.
2. Болдина Е.А., Солдатов В.В. Расчет параметров настройки регуляторов в системах с неточно заданными параметрами объекта // Международный научный журнал. – 2010. – № 1. – С. 61–63.
3. Болдина Е.А., Солдатов В.В., Жиров М.С. Коррекция возмущений связанных с изменениями задающих воздействий в замкнутых робастных системах управления технологическими процессами // Автоматизация в промышленности. – 2010. – № 4. – С. 15–16.
Статьи, опубликованные в других изданиях:
4. Болдина Е.А., Солдатов В.В. Идентификация нелинейных функциональных зависимостей с использованием нечеткого логического вывода // Вестник Российского государственного аграрного заочного университета. –2010. – № 1. – С. 1–4.
5. Болдина Е.А., Солдатов В.В. Математическое моделирование зависимости влажности зерен кофе от продолжительности их обжарки // Вестник Российского государственного аграрного заочного университета. – 2010. – № 1.
– С. 1–3.
6. Болдина Е.А. Разработка нечетких ПИД-регуляторов для управления технологическими процессами / II Ежегодная Всероссийская научно-практическая конференция с международным участием. Перспективы развития информационных технологий, г. Новосибирск, 2010. Сб. науч. тр. – С. 181–187.
7. Болдина, Е. А. Особенности управления объектами с распределенными параметрами / I Всероссийская научно-практическая конференция. Интеллект. Инновации. Информация. Инвестиции. Институты. Инфраструктура, г. Москва, 2010. Сб. науч. тр. – С. 115-118.
8. Болдина Е.А., Солдатов В.В. Методы интеллектуального управления технологическими процессами АПК. – М.: МГУТУ, 2010. – 127 с.