автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.01, диссертация на тему:Модели и алгоритмы управления объектами со стохастическим механизмом формирования показателей качества готовой продукции

кандидата технических наук
Белоус, Ирина Геннадьевна
город
Волгоград
год
2010
специальность ВАК РФ
05.13.01
цена
450 рублей
Диссертация по информатике, вычислительной технике и управлению на тему «Модели и алгоритмы управления объектами со стохастическим механизмом формирования показателей качества готовой продукции»

Автореферат диссертации по теме "Модели и алгоритмы управления объектами со стохастическим механизмом формирования показателей качества готовой продукции"

004610154 На правах рукописи

БЕЛОУС ИРИНА ГЕННАДЬЕВНА

МОДЕЛИ И АЛГОРИТМЫ УПРАВЛЕНИЯ ОБЪЕКТАМИ СО СТОХАСТИЧЕСКИМ МЕХАНИЗМОМ ФОРМИРОВАНИЯ ПОКАЗАТЕЛЕЙ КАЧЕСТВА ГОТОВОЙ ПРОДУКЦИИ (НА ПРИМЕРЕ ШАРОВОЙ МЕЛЬНИЦЫ СУХОГО ПОМОЛА ЦЕМЕНТА)

05.13.01 - Системный анализ, управление и обработка информации

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

- 7 ОПТ 2010

Волгоград - 2010

004610154

Работа выполнена в Камышинском технологическом институте (филиал) Государственного образовательного учреждения высшего профессионального образования «Волгоградский государственный технический университет»

Научный руководитель кандидат технических наук, профессор

Крушель Елена Георгиевна.

Официальные оппоненты: доктор технических наук, профессор

Большаков Александр Афанасьевич, доктор технических наук, профессор Шевчук Валерий Петрович.

Ведущая организация Астраханский государственный

технический университет.

Защита диссертации состоится октября 2010г. в часов на заседании диссертационного совета Д 212.028.04 при Волгоградском государственном техническом университете по адресу: 400131 г. Волгоград, пр. Ленина, 28, ауд. 209.

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

Автореферат разослан » сентября 2010 г.

Ученый секретарь

диссертационного совета У^Ю^ 7 Водопьянов В.И.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. Методы решения задач автоматизации объектов в детерминированной постановке хорошо изучены. Однако существует класс промышленных объектов массового распространения, для которых детерминированные модели и методы управления недостаточно пригодны. Объекты этого класса обладают следующими особенностями: во-первых, непрерывный контроль прямых показателей качества готового продукта в реальном времени не всегда возможен (осуществляется лишь эпизодический лабораторный контроль); во-вторых, не удается получить постоянные во времени показатели качества конечного продукта из-за особенностей технологии и, в-третьих, информационные запаздывания в каналах «управляющие воздействия - показатели качества продукта» настолько велики, что непосредственное использование этих показателей для определения управляющих воздействий проблематично.

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

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

В предшествующих работах по автоматизации процесса помола (A.C. Бо-ронихин, Я.Е. Гельфанд, И.Б. Гинзбург, Ю.С. Гризак, Ю.И. Дубинин, ЮЛ. Крахтанов, A.M. Шейнин и др.) в системе управления использовались сигналы, косвенно связанные с тонкостью помола и доступные для непрерывного измерения вблизи внешней оболочки мельницы (например, сигнал микрофона, измеряющего уровень шума). В ряде работ (В.П. Живоглядов, Е.Г\ Крушель* Б.М. Миркин) исследовались задачи управления процессом помола с использованием сигнала распределенного контроля косвенных показателей. Во всех перечисленных работах структура системы управления постулировалась.

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

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

Для решения этих задач потребовалось развитие методов стохастического управления, предложенных ранее в классических работах (Р. Беллмана, А. Брайсона, Ю.И. Дегтярева, Р. Калмана, A.A. Красовского, Б.М. Миллера, П.В. Пакшина, Хо Ю-Ши, А.Н. Jazwinski и др.) и получивших дальнейшее развитие в исследованиях (Б.Ц. Бахшияна, В.Б. Колмановского, М.Н. Красилыцикова, АЖ Матасова, В.С.-Пугачева, Ю.П. Пытьева, И.Н. Синицина, Ф.Л. Черноусько и др.). Известно, что законченные результаты получены для задачи управления линейными динамическими объектами с квадратическим критерием качества управления в условиях, когда на управляющие воздействия и переменные состояния объекта не накладываются ограничения-неравенства, параметры объекта являются детерминированными, входные воздействия и измерительные помехи - аддитивные и гауссовские (условное название данной группы результатов - линейно-квадратично-гауссовская, ЛКГ-теория).

Однако существующие предположения ЛКГ-теории на практике не всегда выполняются. В частности, информация о стохастических характеристиках обычно является неполной; для преодоления фактора стохастической неопределенности предлагаются методы адаптивного управления (Я.3. Цыпкин, Л.С. Понтрягин, В.Ю. Рутковский, A.A. Фельдбаум, A.JI. Фрадков, В.А. Якубович, L. Ljung и др.). Кроме того, управляющие воздействия и переменные состояния практически всегда ограничены, но теоретические основы учета этих ограничений в настоящее время разработаны недостаточно.

Актуальность темы работы в решении теоретических вопросов состоит в выработке подхода к расширению области использования результатов ЛКГ-теории стохастических оптимальных систем управления, позволяющего учесть естественные технологические ограничения на диапазоны изменений управляющих воздействий и показателей качества конечного продукта в условиях, когда вероятностные характеристики внешних воздействий априори известны неточно и подлежат уточнению в ходе работы системы в реальном времени.

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

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

Дня достижения цели в работе решены следующие задачи: 1. Создание имитационной модели примера объектов со стохастическим механизмом формирования показателей качества продукции и со значительными

запаздываниями в каналах «управляющее воздействие - выход» - шаровой мельницы сухого помола цемента.

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

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

4. Разработка методики оценки эффективности предлагаемых модификаций алгоритмов управления и фильтрации. Проведение вычислительных экспериментов для исследования эффективности предлагаемых решений.

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

Достоверность результатов. Обоснованность и достоверность полученных результатов доказывается результатами вычислительных экспериментов, в которых проводится сопоставление показателей качества и характеристик системы, алгоритмическое обеспечение которой синтезировано на базе предложенных методов, с одноименными показателями, достижимыми в системе управления, синтезированной на модели, полученной в работах-предшественниках (В.П. Живоглядов, Е.Г. Крушель, Б.М. Миркин).

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

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

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

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

Практическая значимость результатов.

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

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

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

1. При использовании системы оценки показателя качества продукции по косвенным измерениям переменных состояний (в контуре разомкнутого управления, т.е. без автоматической обратной связи по каналу «показатель качества - производительность») - для раздельного складирования выходного продукта с выделением объемов, соответствующих различным маркам.

2. При наличии контура обратной связи и системы субоптимального управления - для стабилизации выхода объекта (для снижения дисперсии тонкости помола не менее чем на 7% по сравнению с уровнем, достижимом при ручном управлении без нарушений показателей качества цемента).

На защиту выносятся:

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

2. Нелинейный алгоритм управления с прямым учетом ограниченности диапазона изменения управляющих воздействий, позволяющий обобщить результаты ЛКГ-теории в области синтеза структур оптимальных систем на класс объектов со стохастическим механизмом формирования качества продукта.

3. Алгоритм оценивания переменных состояния и неконтролируемых возмущений, построенный на отклонении оценки выхода объекта от задания.

Внедрение. Результаты разработки алгоритмического и программного обеспечения внедрены на ОАО «Себряковский цементный завод» для использования в системе управления производительностью цементной шаровой мельницы с учетом технологических требований к качеству готовой продукции.

Учебный вариант математической модели объекта, алгоритмы управления и средства компьютерной поддержки используются в Камышинском технологическом институте (филиал) Государственного образовательного учреждения высшего профессионального образования «Волгоградский государственный технический университет» для подготовки инженеров по специальности «Автоматизированные системы обработки информации и управления».

Апробация работы. Основные положения и отдельные результаты диссертационной работы докладывались и обсуждались на III и VI всероссийских конференциях «Инновационные технологии в обучении и производстве» (Камышин, 2005; Камышин, 2009), VI международной научно-методической конференции «Информатика: проблемы, методологи, технологии» (Воронеж, 2005), VIII международной научно-практической конференции «Экономико-

организационные проблемы проектирования и применения информационных технологий» (Кисловодск, 2005), XIII международной конференции «Математика. Компьютер. Образование» (Дубна, 2006), X международной конференции «Информатика: проблемы, методологии, технологии» (Воронеж, 2010).

Публикации. Основные результаты работы опубликованы в девяти статьях и материалах конференций, три из которых опубликованы в изданиях из Перечня ВАК. Всего по теме диссертации опубликовано 10 работ.

Объем и структура работы. Работа состоит из введения, четырех глав, заключения, списка литературы и 11 приложений. Основной текст изложен на 163 страницах, содержит 21 рисунок, 14 таблиц, список литературы из 145 наименований публикаций отечественных и зарубежных авторов.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ Во введении обоснована актуальность исследуемой проблемы, сформулированы цель и задачи работы, описана структура диссертации.

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

Показано, что ряд задач автоматизации объектов с недетерминированным механизмом формирования качества готовой продукции может быть эффективно формализован в терминах стохастической теории оптимальных систем управления, наиболее полные результаты которой получены применительно к задаче управления объектом со стохастическими входными воздействиями и с измерениями переменных состояния при наличии помех в ЛКГ-постановке.

Для того, чтобы обеспечить возможность переноса ряда результатов ЛКГ-теории на класс объектов со стохастическим механизмом формирования качества готовой продукции, потребовалось создание математической модели процессов рассматриваемого класса, которая, с одной стороны, с достаточной полнотой отражает свойства реального объекта (как объекта управления) и, с другой стороны, не противоречит допущениям о моделях объекта, для которых применимы теоретические результаты.

Были учтены следующие особенности объектов рассматриваемого класса:

- Наличие запаздываний в каналах передачи входных воздействий на выход объекта.

- Технологические ограничения на измерения показателя качества продукции.

- Распределенный характер обработки материала по длине объекта.

- Неполнота информации о возмущениях и переменных состояния.

Применительно к модели шаровой мельницы сухого помола цемента были конкретизированы общие особенности объектов со стохастическим механизмом формирования качества готовой продукции:

1. Оперативное управление производительностью мельницы в ручном режиме производится по экспресс-показателю качества готового продукта- тонкости помола. В системе автоматического управления тонкость помола не используется, т.к. время на отбор и обработку пробы превышает время пребывания материала в мельнице.

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

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

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

Управляющее воздействие и[$] предполагается скалярным; в качестве управляющего воздействия рассматривается расход сырья, поступающего в горловину мельницы.

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

Динамика значений каждого компонента вектора состояния при изменениях управляющего и/или возмущающего воздействия имитируется последовательным соединением звена запаздывания и апериодического звена 1-го порядка. Распределенный характер процесса приближенно учитывается зависимостью параметров динамических звеньев от расстояния между точкой измерения соответствующего компонента вектора состояния и горловиной мельницы. Предполагается, что величины запаздываний т[1] и постоянных времени Т[1] линейно увеличиваются в зависимости от этого расстояния (1).

Т[1] = а + с-1-Ы, в[1) = Ь + й-1-Ы, (1)

Ю.где дI - дискрета отсчета координаты по длине мельницы, м; Ьм - длина мельницы, м; 1=0.. Ьм - количество дискрет д/, укладывающееся в общую длину мельницы; Д/ - дискрета отсчета времени, мин; а, Ъ,с,<1- параметры линейных зависимостей постоянных времени и запаздываний от I, мин.

Модель рассматривается в дискретном времени, причем интервал дискретизации времени выбирается настолько малым (Д* « Т[/]), что эффект дискретизации можно не учитывать.

Связь между косвенными показателями, измеряемыми непрерывно, и экспресс-показателем качества готового продукта рассматривается как стохастическая, На рис. 1а и 16 приведены результаты синхронизированных во вре-

мени измерений косвенных показателей и тонкости помола в статике (при неизменном положении ножа питателя).

0123456789 10 И 1213 14 ^ 0 1 2 3 4 5 6 7 8 9 10111213 14 Л | Расстояние от горловины меяышцы Расстояние от горловины мельшшы до ее выхода, м Д° ее выхода, м

а) Косвенный показатель б) Тонкость помола

Рисунок 1 - Результаты синхронизированных во времени измерений Для введения стохастического механизма связи между тонкостью помола ^тонк.пом Ь» 0 и значениями косвенных показателей ^КОсв.пок[5>(1 текущее значение каждого косвенного показателя (2) формируется при моделировании аддитивно, с использованием генератора высокочастотного случайного процесса хр[5, ¿], не зависящего от тонкости помола цемента (3). ^косв.пок!^'

I] = л,м

■ -^тонк.пом [3,1]+1Р{8,11 (2)

tttyts.il = 0, вф[1] = Оф-е (3)

где q - параметр, учитывающий среднее уменьшение сигнала датчика на единицу длины мельницы, [с/]=мВ/м.

Среднеквадратичное отклонение аг[1] этого процесса (имитируемого чисто случайной последовательностью) зависит от координаты расположения датчика и убывает по длине мельницы (3).

Предполагается, что математическое ожидание ^хкосв.пак [5>'] косвенного показателя линейно зависит от математического ожидания тх [я, Л тонко-

^гокк.пс.ч -

ста помола в точке размещения датчика I (параметры линейной зависимости для каждого косвенного показателя различны) (4).

г"х^Л]=кА1]-тХптШ]+щрШ], т = (4)

где к,[1\ - коэффициент, показывающий, какая величина тонкости помола будет достигнута после завершения переходного процесса при изменении задающего воздействия на выход объекта, зависит от / экспоненциально, различается для каждой из /камер мельницы, /=1 .Атщ, /„„(-числокамер мельницы; V, а у - параметры, учитывающие значение тонкости помола загружаемого в мельницу материала, среднее уменьшение тонкости помола на единицу длины и величину уменьшения тонкости помола после окончания переходного процесса (при изменении задающего воздействия на выход объекта) соответственно, [г]=% остатка на сите после просева цементной пробы, И=1/м, [у]=безразмерная;

тгц, [5,1] - математическое ожидание помехи формирования косвенного показателя 1р[з, /], мВ;

5 = 0-" 5тах - диапазон изменения тактов дискретного времени.

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

Тонкость помола определяется соотношениями (5-7).

!-^тонк.пом 1> - 1 Л], если 5 < т[I] ¿Ш •*то„к.пом[5 -1,1]+В[1]-и[5-т[1]]+ , (5) '/клинкер [5,4, если 5 > т[/]

^тонк.пом[0,/] = г-е^4, (6)

А[1] = В[1] = к[1]-^, (7)

где АИ, В[/], СЩ - параметры дискретного звена, соответствующего апериодическому звену 1-го порядка;

УИ - возмущающее воздействие (физическая природа которого - случайные изменения силы удара шаров, случайные колебания твердости сырья -клинкера, поступающего в мельницу), нормально распределенная чисто случайная последовательность с заданными ненулевыми математическим ожиданием и средним квадратическим отклонением, мВ; г, } - параметры, учитывающие значение тонкости помола загружаемого в мельницу материала и среднее уменьшение тонкости помола на единицу длины мельницы соответственно, [г]=% остатка на сите после просева цементной пробы, [/]=1/м;

1/клинкер[5Л) - возмущения, вызывающие стохастичность зависимости тонкости помола от расхода клинкера в виде нормально распределенной чисто случайной последовательности с заданным средним квадратическим отклонением и нулевым математическим ожиданием, мВ.

Законы распределения случайных компонентов модели предполагаются гауссовскими с известными статистиками.

Моделирование динамики изменения косвенных показателей происходит из случайных начальных условий, задаваемых для каждой камеры; законы распределения начальных условий предполагались гауссовскими.

В общем случае модель в состояниях рассматривается в следующем виде: *[5 + 1] = ФИ • х|>] + В • ф] + Г[в] • 5 = О,...^^-!, (8) где х - расширенный вектор состояния (п+т+1)х1, учитывающий значения управляющих воздействий в тактах запаздывания т, значения косвенных показателей и значение выхода объекта в данный такт времени, п - число косвенных показателей; ФЭД - (п+т+1)*(п+т+1)-матрица; В - (п+т+1)хт-матрица; Г[У) - (п+т+1)хг-матрица; эти матрицы неслучайны; т - число управляющих воздействий, г - число возмущающих воздействий; «И - га-мерный детерминированный вектор управляющих воздействий; - гаус-совская чисто случайная последовательность г-мерных векторов с заданными статистиками (математическим ожиданием и матрицей ковариаций).

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

представлена на языке, принятом в ЛКГ-теории, с точностью, достаточной для решения задач синтеза и моделирования системы управления.

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

Известно, что задача оптимального стохастического управления при выполнении допущений, принятых в ЛКГ-теории, распадается на независимые подзадачи стохастической фильтрации и детерминированного оптимального управления (рис. 2). Контур 1 рис. 2 представляет собой контур оценки неизме-ряемых переменных состояния и неконтролируемых возмущений и выполняет задачу стохастической оптимальной фильтрации.

Для оценивания неизвестных переменных состояния используется фильтр Калмана, формирующего оценку +1] вектора состояния по результатам расчета предшествующей оценки измерений (с помехами) и доступных компонентов выхода объекта>{^+1]. Р, 9 - параметры фильтра - матрицы ко-вариации ошибки измерения, ковариации измерительных помех и ковариаций входного возмущения соответственно. -^ыхИ - оценка выхода объекта.

w[s]

Объект

x[s]

Формирование выхода

4M Hx[s]

■УМ

ф]|

Фильтр

yW

*«М

"1

L.

Модель объекта ад

Модель формирования выхода

Формирование оптимального управления

---1 Контур 1

Фильтрация

Контур 2 Управление

Рисунок 2 - Структура стохастической оптимальной системы управления

Контур 2 рис. 2 представляет собой контур управления и выполняет задачу детерминированного оптимального управления, алгоритм для которого синтезируется методом аналитического конструирования оптимальных регуляторов (АКОР), базирующегося на оценках, полученных в контуре 1. В системе обеспечивается минимизация квадратического критерия

У = (Фтах] - Xz{Smaxïf " ' Otw] ~ *z[wD + ~

xz[s]7 • Q[s] ■ (x[s] - xz[s]) + (u[s] - uz[s])r ■ R[s] ■ (u[s] - u,M)], (9)

где лДл] -n-мерные векторы задающих воздействий (детерминированные); m2[s] -r-мерные векторы желательных значений управляющих воздействий;

KfSmax], Q[s], s=0,..., smax-l - неотрицательно определенные nxn симметричные матрицы штрафов за отклонение вектора состояния от задания в конце и по ходу процесса управления;

R[s], s=0,..., Snuw-1 - положительно определенная rxr симметричная матрица штрафа за отклонение управляющих воздействий от желательных значений.

Методы, развитые в ЛКГ-теории, обеспечивают возможность формализованного синтеза структуры и алгоритмов управления объектами, которым свойственно наличие факторов неопределенности. Однако практическое применение некоторых результатов ЛКГ-теории встречает следующие трудности:

1. Не всегда удается проследить связь между квадратическим критерием ЛКГ-задачи и инженерными показателями качества.

2. Составляющая ЛКГ-критерия, соответствующая штрафам за отклонения управляющих воздействий от номинальных значений и введенная для косвенного учета ограниченности диапазона управлений, может привести к нежелательному (а иногда и недопустимому) смещению среднего значения выхода объекта относительно задающего воздействия.

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

В качестве иллюстрации недостаточной точности воспроизведения задающего воздействия при использовании управлений, рассчитанных методами ЛКГ-теории, на рис. 3 приведены результаты отработки ступенчатого изменения задания по тонкости помола; параметры второй составляющей критерия (9) выбраны так, чтобы гарантировать отсутствие нарушений допустимого диапазона изменений управляющих воздействий. При этом оказалось, что среднее значение тонкости после достижения стационарного режима отклоняется от задания на 7%.

Тонкость 1«г

помола,

„м

Х,М

10'

60

S0

S. нив

80 s, мив

б) управляющее воздействие

а) выход объекта

Рисунок 3 - Реакция объекта на ступенчатое изменение задания по тонкости помола, используя ЛКГ-теорию * В представляемой работе результаты ЛКГ-теории были использованы только для синтеза структуры системы управления объектами рассматриваемого класса. В алгоритмы расчета управляющих воздействий и оценивания не-измеряемых компонентов вектора состояния были внесены модификации, расширяющие область практического применения теории стохастического управления применительно к рассматриваемому классу объектов

Предлагаемая модификация критерия ЛКГ-задачи состоит в резком уменьшении вклада второй составляющей квадратического критерия (9) по

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

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

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

Для подтверждения эффективности введенной модификации ЛКГ-алгоритма проведено моделирование процесса отработки ступенчатого изменения задающего воздействия. Результаты моделирования (рис. 4), во-первых, показывают, что достигнуто существенно более высокое качество отслеживания задающего воздействия в стационарном режиме (по сравнению с графиком на рис. 3) и, во-вторых, иллюстрируют кратковременность режима, при котором действует нелинейное ограничение управляющего воздействия.

ОД, 1«г

1Ч»С

И-

20

40

60 >0 5, НИН

б) управляющее воздействие

40 60 и

5, МИН

а) выход объекта

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

Результаты улучшения точности управления в системе с модифицированным алгоритмом в системе управления тонкостью помола показаны на рис.5.

Для подтверждения этого предположения были рассчитаны значения среднего квадратического отклонения выхода объекта от задания при использовании двух вариантов алгоритмов

управления процессом помола: первый следовал методам ЛКГ-теории, а во втором использовалась предложенная нелинейная модификация.

Рисунок 5

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

В представляемой работе рассматриваются и сравниваются три схемы фильтров, предназначенных для оценки вектора состояния объекта по результатам измерения части компонентов с помехами:

Схема 1: Классический фильтр Калмана-Бьюси (ФКБ), в котором используются априорные оценки вероятностных характеристик (схема используется как база для сравнения с модифицированными вариантами, предлагаемыми в работе). Схема 2: Фильтр с оценкой вероятностных характеристик в реальном времени по укороченным выборкам предшествующих зашумленных измерений (известный вариант реализации ФКБ для случая, когда входные возмущения и помехи нестационарны (A.A. Дегтярев, Ш. Тайль)).

Схема 5: Новый вариант фильтра, построенный на отклонениях оценки выхода объекта от его желаемого значения. При разработке данной схемы учтено назначение фильтра: использование в замкнутой системе управления выхода объектов рассматриваемого класса.

Суть предлагаемого метода заключается в том, что пересчет оценок вероятностных характеристик и параметров фильтра следует производить в том случае, если ошибка 5[s] оценки выхода объекта (отклонение оценки выхода ■^выхМ от задающего воздействия xz[s]) превышает установленную величину f.

¿M = *zM-W4 №И>?. (Ю)

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

Шаг 1. До начала работы алгоритма производится расчет оценок необходимых статистических характеристик на основе некоторого числа априорных *" данных, либо данных, полученных в результате моделирования системы. В случае, когда количество доступных данных меньше выбранного значения N, статистики рассчитываются на основе тех измерений, которые

¡1, а а я

w с i

sl:

Я Н Я и

Я 2 «

üa

0,8

0,7

0,6

а 0,5 <о

|0,4 о 0,3

|о,1

ЛКГ-алгоритм

, Модифицированный ,'^АКОР с прямым учетом

ограниченности управляющего воздействия

10 15 20 25 30 Нижняя граница ограничения управления,!

доступны на данный момент. В противном случае статистики рассчитываются для N измерений. Шаг 2. Производится расчет параметров ФКБ, в том числе и матрицы ковариа-

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

Шаг 4. Расчет ошибки оценки измерения и сравнение ее с заданной величиной погрешности Если |£[$]| > переход к шагу 1, на котором производится пересчет статистик для выборки из N прошлых измерений. Если 15Ы1 < <Г, то переход к шагу 3.

Приведены результаты сравнения работы рассматриваемых алгоритмов фильтрации при отработке разных типов помех измерению (структура системы показана на рис. 2, рассматриваемые варианты алгоритмов оценки состояния предназначены для контура «Фильтрация»).

На рис. 6 приведен график измерительной помехи с изменяющимися вероятностными характеристиками.

|| 'о

300 мин

Рисунок 6 - Форма помехи измерению Графики соответствующего изменения выхода объекта для рассматриваемых алгоритмов фильтрации показаны на рис. 7.

« « 0.05

а * я

£

и к I 3 5 I

Я 0.05

60

. Хтошлом!5!. фнлыр по схеме 1

120

фильтр по схеме 2

180

240

-тонхлом^» фильтр по схеме 3

X.

300 5, мин

Рисунок 7 - Сравнение реакции фильтров 1,2,3 на изменения помехи измерению Графики сравнения среднего квадратического отклонения (с.к.о.) выхода объекта от задающего воздействия для всех рассматриваемых типов фильтрации, приведены на рис. 8. кинт - интенсивность помех относительно первоначальной при ее изменении для всего исследуемого промежутка времени.

с.к.о.тонкости помола от задающего воздействия, % остатка на сите 0,02

после просева цементной пробы

0,01

-- фильтр по схеме 1

— фнгатр по схеме 2

— фильтр по схеме 3 0.

1 1,5 2 2,5 3

*вэт

а) с.к.о. выхода объекта от задания

с.г.о.

производительности мельницы от номинального 0,16 значения, т/час

0,15

---фильтр по схеме 1

фильтр по схеме 2 - - фильтр по схеме 3 0,14

1 1,5 2 2,5 3

б) с.к.о. управления от номинального значения Рисунок 8 - Сравнение с.к.о. выхода объекта от задания для фильтрации по схемам 1,2,3

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

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

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

Задача 1. «Сортировка цемента по тонкости помола. Моделирование процесса помола в режиме ручного управления». 12,5Т

Максимально допустимая тонкость

Низкокачественный цемент

Цемент среднего качества

Высококачественный цемент

250 300 350 время, мин

Рисунок 9 - Колебания тонкости помола при постоянном расходе сырья без автоматизации

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

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

* 12

г

Максимально допустимая тонкость

Средняя тонкость помола

so введения автоматизации §

S 12

Максимально лопустимая тонкость

I

! П

Источник повышения производительности за счет введения автоматизации

Средни тонкость помола до введения автоматизации

50

150

104-

200 о 50 150 200

время, мин время, мин

а) б)

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

Задача 3. «Определение мест расположения датчиков косвенных показателей, используемых в системе автоматического управления». Метод основан на оценке тонкости помола по значениям датчиков косвенных показателей, зависящим от мест установки датчиков как от параметра. Места установки выбираются на основе компромисса между стремлениями уменьшить дисперсию косвенного показателя и уменьшить влияние временного запаздывания в канале «управляющее воздействие - косвенный показатель».

Разработана методика оценки эффективности предложенных субоптимальных алгоритмов управления (глава 2) и оценки вектора состояния (глава 3).

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

Первый алгоритм реализуется как классический ЛКГ-алгоритм, в котором параметры штрафа за отклонения управляющих воздействий от номинальных значений таковы, что гарантируют соответствие рассчитанных управлений допустимому диапазону.

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

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

1. Оценка математического ожидания отклонения тонкости помола от задания.

2. Оценка среднего квадратического отклонения (с.к.о.) выхода объекта и других компонентов вектора состояния от заданных значений.

3. Относительная (по отношению к задающему воздействию) погрешность среднего квадратического отклонения тонкости помола от задания.

4. Относительная (в % к задающему воздействию) погрешность математического ожидания отклонения тонкости помола от задающего воздействия.

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

1. Оценка математического ожидания управляющего воздействия - производительности мельницы.

2. Оценка среднего квадратического отклонения управления от номинального значения.

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

- В системе управления, синтезированной на базе ЛКГ-теории, не всегда удается достичь требуемой точности отслеживания задающих воздействий по тонкости помола, поскольку соблюдение требования технической реализуемости управляющего воздействия приводит к смещенности оценки математического ожидания тонкости помола относительно задающего воздействия.

- Использование предлагаемого алгоритма ведет к улучшению качества достижения цели за счет устранения смещенности оценки математического ожидания тонкости помола относительно задающего воздействия. Отсутствие смещенности наблюдается и при применении алгоритма управления, основанного на использовании ПИ-регулятора, но показатели качества достижения цели у данного варианта алгоритмического обеспечения, хуже, чем с использованием предлагаемого модифицированного ЛКГ-алгоритма.

- Применение модифицированного ПИ-закона управления (с прямым учетом ограниченности управляющего воздействия) и различных вариантов алгоритмов оценки вектора состояния (схемы 1-3, приведенные в главе 3) в контуре фильтрации позволяет достичь лучших показателей качества по сравнению с вариантом применения классического ПИ-закона. Однако показатели качества управления заметно уступают достижимым при использовании модифицированного ЛКГ-алгоритма с фильтрацией по схемам 1-3.

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

Результаты моделирования подтверждают предположения об эффективности предлагаемого варианта алгоритмического обеспечения.

Наиболее заметное улучшение качества управления достигнуто при использовании следующего варианта алгоритмического обеспечения: модифицированный ЛКГ-алгоритм, базирующийся на оценках вектора состояния объекта с помощью фильтра с оценкой вероятностных характеристик в реальном времени по укороченным выборкам предшествующих зашумленных измерений (схема 2, глава 3) либо фильтра, построенного на отклонениях оценки выхода объекта от его желаемого значения (схема 3, глава 3).

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

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

стохастической теории оптимальных систем. Разработанный подход применен для разработки математической модели одного из объектов данного класса (шаровой мельницы сухого помола цемента).

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

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

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

5. Рекомендации по использованию результатов исследования для усовершенствования системы управления процессом помола цемента в шаровых мельницах переданы на ОАО «Себряковцемент».

Список опубликованных работ по теме диссертации:

Публикации в изданиях, рекомендованных ВАК РФ:

1. Белоус, И.Г. О точности воспроизведения задания в стохастических системах с управляющими воздействиями, ограниченными по диапазону / И.Г. Белоус, Е.Г. Крушель // Известия ВолгГТУ. - 2009. №12(60). - С. 44-48.

2. Белоус, И.Г. Оценка неизмеряемых показателей качества технологического процесса и входных возмущений по результатам автоматического контроля косвенных показателей / И.Г. Белоус, Е.Г. Крушель // Известия ВолгГТУ. -

2009.№12(60).-С. 71-74.

3. Белоус, И.Г. Синтез структуры и алгоритмического обеспечения управления стохастическим объектом / И.Г. Белоус, Е.Г. Крушель // Вестник АГТУ. -

2010.№2.-С. 78-86. Публикации в других изданиях:

4. Крушель, Е.Г. Учебная задача проектирования системы автоматического управления для студентов-системотехников (на примере цементной шаровой мельницы) / Е.Г. Крушель, И.Г. Шляхтина (Белоус) // Информатика: проблемы, методология, технологии: матер. VI междунар. научно-методической конф., г. Воронеж / ВоронежГУ/ - Воронеж, 2006. - С. 217222.

5. Шляхтина (Белоус), И.Г. Постановка учебной задачи проектирования системы стохастического оптимального управления процессом помола цемента

(для студентов-системотехников) / И.Г. Шляхтина // Экономико-организационные проблемы проектирования и применения информационных технологий: сб. науч. тр. VIII международной научно-практической конф еренции РГЭУ «РИНХ» - Ростов-на-Дону, 2006. - С. 90-92.

6. Крушель, ЕГ. Выбор оптимального места установки датчика / Е.Г. Кру-шель, И.Г. Шляхтина (Белоус) // Инновационные технологии в обучении и производстве : матер. III всероссийской кон})., г.Камышин, 1.2 I - Камышин, 2005.-С. 130-134.

7. Шляхтина, И.Г. Компьютерное управление объектом с неполным измерением вектора состояния /И.Г. Шляхтина (Белоус), Е.Г. Крушель //Математика. Компьютер. Образование. Сборник тезисов 13 международной конф. под ред. Г, Ю. Резниченко/-Ижевск, 2006, выпуск 13.-С. 187.

8. Крушель, ЕГ. Влияние ограничения диапазона управляющих воздействий на качество управления / ЕГ. Крушель, И.Г. Белоус // Инновационные технологии в обучении и производстве : матер. VI всероссийской научно-практической конф.,г.Камышин, Т. 4 /-Волгоград, 2010.-С. 97-100.

9. Крушель, ЕГ. Стохастическое управление работой шаровой мельницы сухого помола цемента с оценкой статистик шумов / ЕГ. Крушель, И.Г. Белоус // Инновационные технологии в обучении и производстве : матер. VI всеросс. научно-пракг. конф „г. Камышин, Т. 4 /-Волгоград, 2010. -С. 101-104.

10. Белоус, И.Г. Постановка задачи и структура математической модели цементной шаровой мельницы как примера объекта со стохастическим механизмом формирования показателей качества готовой продукции / И.Г. Белоус, ЕГ. Крушель // Информатика: проблемы, методология, технологии: матер. X международной научно-методической конф., г. Воронеж Т. 1 / -Воронеж, 2010. -С. 87-90.

Белоус Ирина Геннадьевна

Модели и алгоритмы управления объектам и со стохастическим механизмом формирования показателей качества готовой продукции (на примере шаровой мельницы сухого помола цемента)

Автореферат

Подписано в печать 06.09.2010. Формат 60x84 1/16. Усл.печ.л. 1,0. Тираж 130 экз. Отпечатано _ © МАУ «Медиа-холдинг»,

403882, г. Камышин, ул. Октябрьская, 48.

Оглавление автор диссертации — кандидата технических наук Белоус, Ирина Геннадьевна

Введение.

1. Особенности объектов управления с недетерминированным механизмом функционирования.

1.1. Обзор постановок задач автоматизации объектов с недетерминированным механизмом функционирования и обоснование целесообразности совершенствования методов автоматизации процессов управления.

1.1.1. Примеры объектов с недетерминированным механизмом функционирования.

1.1.2. Обзор результатов, достигнутых в области автоматизации объектов с недетерминированным механизмом функционирования.

1.1.3. Сведения о теоретических методах, использованных при синтезе систем автоматизации процесса управления объектами с недетерминированным механизмом функционирования.

1.1.4. Обзор результатов теории оптимальных стохастических систем и обоснование целесообразности развития теоретических методов синтеза методов алгоритмизации задач управления.

1.2. Разработка имитационной модели объекта со стохастическим механизмом формирования показателя качества готового продукта.

1.2.1. Разработка имитационной модели объекта рассматриваемого класса.

1.2.2. Допущения, принятые при создании модели.

1.3. Математическая модель цементной шаровой мельницы как пример объекта со стохастическим механизмом формирования показателей качества готовой продукции „

1.3.1. Работы-предшественники по автоматизации процесса помола, га недостатки и предложения, вносимые в представляемой работе.-.

1.3.2. Особенности технологии помола цемента и контроля качества готового продукта как стохастической системы.

1.3.3. Допущения, принятые при создании модели мельницы.

1.3.4. Математическая модель динамики прогресса помола, прямых и косвенных показателей качества, возмущающих воздействий и измерительных помех.

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

1.4. Оценка адекватности языка теории стохастических оптимальных систем задачам автоматизации контроля процессом помола.

1.4.1. Сравнение описания модели цементной мельницы со стохастической теорией.

1.4.2. Анализ: сходство и отличия, недостатки.

1.4.3. Постановка задач развития методов теории оптимальных стохастических систем для решения практических задач автоматизации объектов со стохастическим механизмом формирования качества.

1.5. Выводы.

2. Синтез алгоритмов управления объектами со стохастическим механизмом формирования показателей качества.

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

2.1.1. Особенности контроля и управления в условиях стохастической неопределенности.

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

2.1.3. Сведения об использовании теоремы разделения.

2.1.4. Структура управления объектами со стохастическим механизмом формирования показателей качества продукции.

2.1.5. Алгоритмы стохастического оптимального управления.

2.2. Критика положений JIKT-теории с позиций практического применения.

2.2.1. Критика критерия.

2.2.2. Замечания о точности воспроизведения задающих воздействий.

2.3. Предлагаемые в работе усовершенствования системы управления. Синтез субоптимальной системы управления, учитывающей ограниченность допустимого диапазона управляющих воздействий.

2.3.1. Описание предлагаемого в работе подхода учета ограниченности управляющих воздействий.

2.3.2. Постановка задачи обоснования эффективности применения субоптимальной системы.

2.3.3. Иллюстрации эффективности введенных модификаций алгоритмов управления.

2.4. Выводы.

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

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

3.1.1. Обоснование г{елесообразности постановки задачи синтеза алгоритмов оценивания вектора состояния и неконтролируемых возмущающих воздействий.

3.1.2. Обзор результатов, достигнутых в области теории фильтрации и прогнозирования.

3.1.3. Направления развития теории фильтрации.

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

3.2.1. Сводка основных результатов теории оптимальной фильтрации. Особенности и свойства фильтра Калмана.

3.2.2. Обзор результатов адаптивной фильтрации. Область применения адаптивных фильтров.

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

3.3.1. Проверка условий применимости фильтра Калмана для фильтрации и прогнозирования вектора состояния в прогрессе помола в шаровой мельнице сухого помола цемента.

3.3.2. Использование алгоритмов фильтрации при автоматизации процесса помола цемента в шаровой мельнице.

3.3.3. Фильтр с использованием априорных оценок вероятностных характеристик.

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

3.3.5. Предлагаемый фильтр, построенный на отклонении оценки вектора состояния от его эюелаемого значения.

3.3.6. Постановка задачи выработки методики обоснования эффективности предлагаемых алгоритмов фильтрации.

3.4. Сравнение эффективности вариантов алгоритмов оценки векторов состояния и возмущающих воздействии.

3.4.1. Схемы вычислительных экспериментов.

3.4.2. Анализ и сопоставление результатов моделирования для рассматриваемых вариантов оценки вектора состояния и возмущающего воздействия.

3.5. Выводы.

4. Исследование вариантов алгоритмов управления объектами со стохастическим механизмом формирования показателей качества. Предложения по использованию разработанных алгоритмов на производстве и в обучении.

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

4.1.1. Краткие сведения об объекте и особенностях управления технологическим процессом.

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

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

главах 2,3.

4.1.4. Исследования, выполненные для обоснования предложений, перечисленных в подразделе 4.1.3.

4.2. Методика выбора варианта алгоритмического обеспечения.

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

главах.

4.2.2. Методика обоснования эффективности предложений по совершенствованию алгоритмического обеспечения.

4.2.3. Классификация и описание показателей качества управления процессом помола, используемых для сопоставления алгоритмов управления согласно методике п.4.2.2.

4.2.4. Перечень сопоставляемых вариантов алгоритмического обеспечения системы управления процессом помола цемента в шаровой мельнице.

4.2.5. Схема проведения экспериментов.

4.3. Результаты вычислительных экспериментов и предложения по составу алгоритмического обеспечения системы автоматизации процесса помола.

4.3.1. Перечень сценариев вычислительных экспериментов.'.

4.3.2. Обобщенный для всех вариантов алгоритмического обеспечения алгоритм управления.

4.3.3. Описание используемых показателей качества.

4.3.4. Сводка показателей качества, полученных при проведении основных экспериментов.

4.3.5 Анализ полученных результатов.

4.4. Методика выбора числа и мест установки датчиков измерения косвенных показателей

4.5. О возможности использования разработанных средств моделирования и алгоритмизации управления процессом помола в учебном процессе.

4.6. Выводы.

Введение 2010 год, диссертация по информатике, вычислительной технике и управлению, Белоус, Ирина Геннадьевна

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

Методы решения задач автоматизации объектов в детерминированной постановке хорошо изучены. Однако существует класс промышленных объектов массового распространения, для которых детерминированные модели и методы управления недостаточно пригодны. Объекты этого класса обладают следующими особенностями: во-первых, непрерывный контроль прямых показателей качества готового продукта в реальном времени не всегда возможен (осуществляется лишь эпизодический лабораторный контроль); во-вторых, не удается получить постоянные во времени показатели качества конечного продукта из-за особенностей технологии и, в-третьих, информационные запаздывания в каналах «управляющие воздействия - показатели качества продукта» настолько велики, что непосредственное использование этих показателей для определения управляющих воздействий проблематично.

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

Можно отметить два противоположных подхода к проектированию систем автоматизации промышленных объектов. С одной стороны, накоплен большой инженерный опыт создания и эксплуатации систем управления, который может быть использован как база для проектирования новых систем. При несомненных достоинствах инженерного опыта можно указать следующее:

- Выбор проекта-аналога в большой степени зависит от квалификации проектировщика.

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

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

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

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

- Часть допущений, принятых на этапе постановки задачи синтеза, может быть связана не столько с особенностями прикладной задачи, сколько с возможностями предлагаемого метода решения. Это замечание относится как к выбору формы критерия оптимальности, так и к составу учитываемых ограничений (в том числе и к выбору формы модели объекта). Из-за возможного несоответствия принятых допущений конкретной прикладной задаче отсутствует гарантия соответствия инженерных показателей качества управления требованиям технического задания на проектирование системы.

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

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

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

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

В предшествующих работах по автоматизации процесса помола (Я.Е. Гельфанд, И.Б. Гинзбург, Ю.И. Дубинин, Ю.Я. Крахтанов, С.Б. Непомнящий, А.Б. Смолянский, A.M. Шейнин, В. Дуда и др.) в системе управления использовались сигналы, косвенно связанные с тонкостью помола и доступные для непрерывного измерения вблизи внешней оболочки мельницы (например, сигнал микрофона, измеряющего уровень шума). В ряде работ (В.П. Живоглядов, Е.Г. Крушель, Б.М. Миркин) исследовались задачи управления процессом помола с использованием сигнала распределенного контроля косвенных показателей. Во всех перечисленных работах структура системы управления постулировалась.

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

Представляет интерес дальнейшее развитие работ в данном направлении, состоящее в отказе от предопределенности структуры системы управления, т.е. в попытке получить решение задачи максимизации производительности объекта в формальной постановке. Это решение может оказаться полезным для формализации процесса проектирования структуры и алгоритмического обеспечения не только применительно к рассматриваемому объекту, но и для класса объектов со сходными особенностями формирования показателей качества готовой продукции. Попутно могут быть решены также вопросы определения мест расположения датчиков косвенных показателей по длине объекта и опре-' деления настроечных параметров алгоритмов.

Для решения этих задач потребовалось развитие методов стохастического управления, предложенных ранее в классических работах (Р. Беллмана, А. Брайсона, Ю.И. Дегтярева, Р. Калмана, A.A. Красовского, H.H. Красовского, Б.М. Миллера, П.В. Пакшина, A.B. Пантелеева, Ю.С. Попкова, Хо Ю-Ши, Jazwinski А.Н., Kashyap R.L., Rao A.R.) и получивших дальнейшее развитие в исследованиях (Б.И. Ананьева, Б.Ц. Бахшияна, И.Я. Каца, В.Б. Колмановского, М.Н. Красилыцикова, А.П. Крищенко, В.В. Малышева, А.И. Матасова, B.C. Пугачева, Ю.П. Пытьева, H.H. Синицина, Г.А. Тимофеевой, В.М. Хаметова, Ф.Л. Черноусько, Cramer Н., Leadbetter М и др.). Известно, что законченные результаты получены для задачи управления линейными динамическими объектами с квадратическим критерием качества управления в условиях, когда на управляющие воздействия и переменные состояния объекта не накладываются ограничения-неравенства, параметры объекта являются детерминированными, входные воздействия и измерительные помехи - аддитивные и гауссовские (условное название данной группы результатов - линейно-квадратично-гауссовская, ЛКГ-теория).

На практике предположения ЛКГ-теории не всегда выполняются. В частности, информация о стохастических характеристиках обычно является неполной; для преодоления фактора стохастической неопределенности предлагаются методы адаптивного управления (Я.З. Цышсин, Б.Т. Поляк, В.Н. Фомин, СЛ. Урясьев, A.B. Назин, L. Ljung и др.). Кроме того, управляющие воздействия и переменные состояния практически всегда ограничены, но теоретические основы учета этих ограничений в настоящее время разработаны недостаточно.

Актуальность темы работы в решении теоретических вопросов состоит в выработке подхода к расширению области использования результатов ЛКГ-теории стохастических оптимальных систем управления, позволяющему учесть естественные технологические ограничения на диапазоны изменений управляющих воздействий и показателей качества конечного продукта в условиях, когда вероятностные характеристики внешних воздействий априори известны неточно и подлежат уточнению в ходе работы системы в реальном времени.

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

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

Для достижения цели в работе решены следующие задачи:

1. Анализ особенностей объектов управления со стохастическим механизмом формирования показателей качества продукции и со значительными запаздываниями в каналах «управляющие воздействия — выход».

2. На основе данного анализа - разработка математической модели одного из объектов рассматриваемого класса — шаровой цементной мельницы.

3. Обзор и анализ применимости методов теории стохастических оптимальных систем и оптимального оценивания векторов состояния и неизмеряемых возмущающих воздействий для алгоритмизации управления объектами рассматриваемого класса.

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

5. Сравнительный анализ алгоритмов стохастического оценивания векторов состояния и возмущающих воздействий в условиях неполной информации о вероятностных характеристиках объектов анализа. Выбор методов и алгоритмов оценивания для систем автоматизации объектов рассматриваемого класса.

6. Разработка методики моделирования систем управления объектами рассматриваемого класса, проведение вычислительных экспериментов для исследования эффективности предложенных алгоритмов.

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

Достоверность результатов. Обоснованность и достоверность полученных результатов доказывается результатами вычислительных экспериментов, в которых проводится сопоставление показателей качества и характеристик системы, алгоритмическое обеспечение которой синтезировано на базе предложенных методов, с одноименными показателями, во-первых, системы управления, синтезированной методами ЛКГ- теории, и, во-вторых, системы управления, в алгоритмическом обеспечении которой использован популярный пропорционально-интегральный закон управления. Достоверность также подтверждается экспертными заключениями работников внедряющей организации — предприятия Себряковский цементный завод.

Научная новизна.

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

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

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

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

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

Практическая значимость результатов.

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

1. При использовании системы оценки показателя качества продукции по косвенным измерениям переменных состояний (в контуре разомкнутого управления, т.е. без автоматической обратной связи по каналу «производительность - показатель качества») - для раздельного складирования выходного продукта — цемента с выделением объемов, соответствующих различным маркам.

2. При наличии контура обратной связи и системы субоптимального управления - для стабилизации выхода объекта - тонкости помола цемента не менее чем на 7% по сравнению с уровнем, достижимом при ручном управлении без нарушений'показателей качества цемента.

На защиту выносятся:

1. Модель описания транспортирования измельчаемого материала по длине цементной шаровой мельницы.

2. Методика управления процессом помола с учетом ограничений на управляющие воздействия.

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

4. Методика обоснования эффективности предложенных субоптимальных алгоритмов оценивания вероятностных характеристик и управления.

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

Учебный вариант математической модели объекта, алгоритмы управления и средства компьютерной поддержки используются в Камышинском технологическом институте (филиал) Государственного образовательного учреждения высшего профессионального образования «Волгоградский государственный технический университет» для подготовки инженеров по специальности «Автоматизированные системы обработки информации и управления».

Апробация работы.

Основные положения и отдельные результаты диссертационной работы докладывались и обсуждались на III и VI всероссийских конференциях «Инновационные технологии в обучении и производстве» (Камышин, 2005; Камышин, 2009), VI международной научно-методической конференции «Информатика: проблемы, методологи, технологии» (Воронеж, 2005), VIII международной научно-практической конференции «Экономико-организационные проблемы проектирования и применения информационных технологий» (Кисловодск, 2005), XIII международной конференции «Математика. Компьютер. Образование» (Дубна, 2006), X международной конференции «Информатика: проблемы, методологии, технологии» (Воронеж, 2010).

Публикации.

Основные результаты работы опубликованы в девяти статьях и материалах конференций [14, 15, 16, 55, 56, 62, 64, 68, 130], три из которых [16, 62, 64] опубликованы в изданиях из Перечня ВАК. Всего по теме диссертации опубликовано 10 работ.

Структура и объем работы.

Диссертационная работа состоит из введения, четырех глав, заключения, списка литературы и приложений. Работа изложена на 163 страницах основного текста, содержит 21 рисунок, 14 таблиц, список литературы из 145 наименований публикаций отечественных и зарубежных авторов.

Заключение диссертация на тему "Модели и алгоритмы управления объектами со стохастическим механизмом формирования показателей качества готовой продукции"

4.6. Выводы

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

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

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

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

5. Разработана методика выбора оптимального места установки датчика измерения косвенного показателя.

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

ЗАКЛЮЧЕНИЕ

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

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

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

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

5. Предложенная методика использована для оценки эффективности предлагаемых алгоритмов управления и фильтрации в системе управления процессом помола в шаровой мельнице сухого помола цемента. Рекомендации по использованию результатов исследования для усовершенствования системы управления процессом помола цемента в шаровых мельницах переданы на ОАО «Себряковцемент».

Список терминов, условных обозначений и сокращений

АКОР — аналитическое конструирование оптимальных регуляторов.

АСОИУ - автоматизированные системы обработки информации и управления.

ЛКГ - линейно-квадратично-гауссовский. м.о. - математическое ожидание.

П - пропорциональный.

ПИ - пропорционально-интегральный.

ПИД - пропорционально-интегрально-дифференциальный.

САУ - система автоматического управления. с.к.о. - среднее квадратическое отклонение.

Библиография Белоус, Ирина Геннадьевна, диссертация по теме Системный анализ, управление и обработка информации (по отраслям)

1. Адаптивные фильтры: пер. с англ. / Под ред. К. Ф. Н. Коуэна и П.М. Гранта. -М.: Мир, 1988.

2. Александров, А.Г. Оптимальные и адаптивные системы / А.Г. Александров. М.: Высшая школа, 1989. - 262 с.

3. Алексеев, Б.В. Производство цемента / Б.В. Алексеев, Г.К. Барбашев. -М., 1985.-250 с.

4. Алексеев, Б.В. Технология производства цемента: учебник для сред, проф.-техн. училищ / Б.В. Алексеев. М. : Высш. школа, 1980.— 266 с.

5. Алтунин, А.Е. Модели и алгоритмы принятия решений в нечетких условиях: монография / А.Е. Алтунин, М.В. Семухин. Тюмень. Изд-во Тюменского государственного университета, 2000. - 325 с.

6. Андреев, Н.И. Теория статистически оптимальных систем управления / Н.И. Андреев. М. : Наука, 1980.

7. Аоки, М. Введение в методы оптимизации / М. Аоки : пер. с англ., М. : Мир, 1977.

8. Афанасьев, В.Н. Математическая теория конструирования систем управления: учеб. для вузов / В.Н. Афанасьев, В.Б. Колмановский, В.Р. Носов. Изд. 3-е, испр. и доп. М. : Высшая школа, 2003. - 614с.: ил.

9. Афанасьев, В.Н. Управление стохастическими системами / В.Н. Афанасьев, В.Б. Колмановский М. : Изд-во МИЭМ, 1989.

10. Балакришнан, А. Теория фильтрации Калмана/ А. Балакришнан : Пер. с англ. М. : Мир, 1988. - 168 е., ил.

11. Бауман, В.А. Механическое оборудование предприятий строительных материалов, изделий и конструкций: учебник для строительных вузов / В.А. Бауман, Б.В. Клушанцев, В.Д.Мартынов. М. : Машиностроение, 1981. -324 с.

12. Бахилина, И.М. Синтез грубых линейных квадратичных гауссовских регуляторов / И.М. Бахилина, С.А. Степанов // Автоматика и телемеханика, 1998.-№7.-С. 96-106.

13. Беллман, Р. Динамическое программирование / Р. Беллман. М. : 1960 400с

14. Белоус, И.Г. Синтез структуры и алгоритмического обеспечения управления стохастическим объектом / И.Г. Белоус, Е.Г. Крушель // Вестник АГ-ТУ.-2010. №2.-С. 78-86.

15. Беседин, П.В. Исследование и оптимизация процессов в технологии цементного клинкера: монография / П.В. Беседин, П.А. Трубаев : .Под общ. ред. П. В. Беседина. Белгород: Изд-во БелГТАСМ, БИЭИ, 2004. - 420 с.

16. Бесекерский, В.А. Теория систем автоматического регулирования / В.А. Бесекерский, Е.П. Попов. М. : Наука, 1975. - 768 с.

17. Богданов, B.C. Основные процессы в производстве строительных материалов / B.C. Богданов, И.А. Семикопенко, A.C. Ильин. Белгород, 2008. - 550с.

18. Богданов, B.C. Шаровые барабанные мельницы / B.C. Богданов. Белгород, 2002 г.

19. Бойчук, JI.M. Оптимальные системы автоматического регулирования / Л.М. Бойчук. Киев : Наукова думка, 1965. — 84 с.

20. Брайсон, А. Прикладная теория оптимального управления / А. Брайсон, Хо Ю-Ши. М. : Мир, 1972. - 544 с.

21. Браммер, К. Фильтр Калмана-Бьюси / К. Браммер, Г. Зиффлинг, М. : Наука, 1982.

22. Востриков, A.C. Оптимальные и адаптивные системы: учеб. пособие /. A.C. Востриков. Новосиб. электротехн. ин-т. Новосибирск, 1977. - 65 с.

23. Востриков, Л.С. Управление динамическими объектами: учеб. пособие / Л.С. Востриков. Новосиб. электротехн. ин-т. Новосибирск, 1979. - 112с.

24. Гилл, Ф. Практическая оптимизация / Ф. Гилл, У. Мюррей, М. Райт : Пер/ с англ. М. : Мир, 1985.

25. Гельфанд, Я.Е. Автоматическое регулирование процессов дробления и помола в промышленности строительных материалов / Я.Е. Гельфанд, И.Е. Гинзбург. Л. : Стройиздат, 1969. - 176 с.

26. Гинзбург, И.Б. Автоматизация цементного производства: справ, пособие / И.Б. Гинзбург, А. Б. Смолянский. Л. : Стройиздат Ленингр. отд-ние, 1986.-С. 190

27. Городецкий, А .Я. Информационные системы. Вероятностные модели и статистические решения: учеб. пособие / А.Я. Городецкий. СПб. : Изд-во СПбГПУ, 2003. - 326 с.

28. Гольцов, A.C. Применение методов искусственного интеллекта в управлении проектами/ A.C. Гольцов, Л.А. Растригин / Под ред. Соколова А.Ю. X.: НАУ им. Н.Е. Жуковского ХАИ, 2002. - 474 с.

29. Граничин, О.Н. Введение в методы стохастической оптимизации и оценивания: Учеб. пособие / О.Н. Граничин. СПб. : Изд-во С.-Петербургского ун-та, 2003. - 131с.

30. Гулько, Ф.Б. Решение нестационарных задач фильтрации и упреждения при произвольной помехе методами моделирования / Гулько Ф.Б., Новосельцева Ж.А. // Автоматика и телемеханика. 1966. № 10. - С 153-168.

31. Дамдинова, Д.Р. Машины и оборудование для измельчения и сортировки строительных материалов: учеб. пособие / Д.Р. Дамдинова, В.Г. Донду-ков. Улан-Удэ : Изд-во ВСГТУ, 2004.

32. Дахнович, A.A. Дискретные системы и цифровая обработка сигналов: учеб. пособие / A.A. Дахнович. — Тамбов : Изд-во Тамб. гос. техн. ун-та, 2007.- 100 с.

33. Дегтярев, A.A. Элементы теории адаптивного расширенного фильтра Калмана / A.A. Дегтярев, Ш. Тайль. М. : Препринт ИПМ им. М.В. Келдыша РАН, 2003.

34. Деревицкий, Д.П. Прикладная теория дискретных адаптивных систем управления / Д.П. Деревицкий, A.JI. Фрадков. М. : Наука, 1981.-216с.

35. Драбкин, Г.С. Автоматизация цементных заводов / Г.С. Драбкин, И.П. Бровар, Я.Е. Гельфанд, Э.Л. Ицкович. -М. : Госстройиздат, 1961.

36. Дуда В. Цемент / В. Дуда : Пер. с нем. Е.Ш.Фельдмана. М. : Стройиздат, 1981.-464с.

37. Дуда, В. Цемент, электрооборудование, автоматизация, хранение, транспортирование: справ, пособие. -М. : Стройиздат, 1987. -373 с.

38. Еремин, Е.Л. Гиперустойчивость системы управления нелинейным объектом с запаздыванием / E.JI. Еремин // Автоматизация технологических процессов : сб. ст. Фрунзе : Фрунзенск. политех, ин-т, 1987.

39. Заде, JI. Теория линейных систем. (Метод пространства состоянии) / JI. Заде, Ч. Дезоер : Пер. с англ. М. : Наука, 1970. - 703 с.

40. Казакевич, В.В. Системы автоматической оптимизации / В.В. Казакевич, А.Б. Родов. М.: Энергия, 1977. - 288 с.

41. Калман, Р. Новые результаты в линейной фильтрации и теории предсказания / Р. Калман, Р. Бьюси. // Техн. механика. Сер. Д. 1961. Т. 83, № 1.

42. Калман, Р. Очерки по математической теории систем / Р. Калман, П. Фалб, М. Арбиб : пер. с англ. / Под ред. Я.З. Цыпкина. Предисл. Э.Л. Наппельбаума. Изд. 2-е, стереот. М.: Едиториал УРСС, 2004. - 400с

43. Кваркернаак, X. Линейные оптимальные системы управления / X. Квар-кернаак, Р. Сиван -М. : Мир, 1977. 650 с.

44. Козлов, В.Н. Управление энергетическими системами. Теория автоматического управления / В.Н. Козлов, В.Е. Куприянов, В.Н. Шашихин : под ред. В.Н. Козлова. СПб. : Изд-во Политехи, ун-та, 2008. - 255с.

45. Колмановский, В. Б. Задачи оптимального управления / В.Б. Колманов-ский // Соросовский Образовательный Журнал 1997. - № 6. - С. 121127.

46. Колмановский, В.Б. Задачи оптимального оценивания / В.Б. Колмановский // Соросовский образовательный журнал, 1999. №11 С. 122-127.

47. Колмановский, В.Б. Приближенный метод для решения задачи оптимального управления в системах с последействием / В.Б. Колмановский, А.И. Матасов. Доклады Академии наук, 1997, Т. 354, №4, С. 465-468.

48. Красовский, A.A. Универсальные алгоритмы оптимального управления непрерывными прочесами / A.A. Красовский, В.Н. Буков, B.C. Шендрик. -М.: Наука, 1977.-272с.

49. Красовский H.H. Задачи управления с гарантированным результатом / H.H. Красовский, В.Е. Третьяков Свердловск : Средне-Уральское книжное изд-во, 1986. - с.

50. Красовский, H.H. Некоторые задачи теории устойчивости движения / H.H. Красовский. М. : Гос. изд-во физ.-мат. лит-ры, 1959. - 211 с.

51. Кругликов, C.B. О принципе дуальности для задач гарантированного априорного управления и оценивания / C.B. Кругликов //Докл. РАН. — 1994. -Т. 335, №5.-С. 570-572.

52. Крушель, Е.Г. Влияние ограниченности диапазона управляющих воздействий на качество управления / Е.Г. Крушель, И.Г. Белоус // Инновационные технологии в обучении и производстве : тез. докл. VII всерос. конф. Камышин. Камышин, 2010. Т. 4. - С. 97-100.

53. Крушель, Е.Г. Выбор оптимального места установки датчика / Е.Г. Кру-шель, И.Г, Шляхтина (Белоус) // Инновационные технологии в обучении и производстве : тез. докл. III всерос. конф. Камышин. Камышин, 2005. -Т. 2. - С. 130-134.

54. Крушель, Е.Г. Информационное запаздывание в цифровых системах управления : Монография / Е.Г. Крушель, И.В. Степанченко. Волгоград : Изд-во ВолгГТУ РПК "Политехник", 2004. - 120 с.

55. Крушель, Е.Г. О сходимости процесса настройки коэффициентов распределенного контроля / Е.Г. Крушель // Идентификация и управление в системах с неполной информацией : сб. ст. — Фрунзе : Изд-во «Илим» Академии наук Киргизской ССР, 1968. С. 23-28.

56. Крушель, Е.Г. О точности воспроизведения задания в стохастических системах с управляющими воздействиями, ограниченными по диапазону / Е.Г. Крушель, И.Г. Белоус // Известия ВолгГТУ. 2009. №12(60). - С. 4448.

57. Крушель, Е.Г. Оценка неизмеряемых показателей качества технологического процесса и входных возмущений по результатам автоматического контроля косвенных показателей / Е.Г. Крушель, И.Г. Белоус // Известия ВолгГТУ. 2009. №12(60). - С. 71-74.

58. Крушель, Е.Г. Способ формирования сигнала распределенного контроля для автоматизации помола цемента в шаровых мельницах / Е.Г. Крушель, В.П. Живоглядов, Б.М. Миркин, Н.П. Шарапов // Журн. Внедренные изобретения. 1975. — №1. - М.

59. Крушель, Е.Г. Способ формирования сигнала распределенного контроля для-автоматизации помола цемента в шаровых мельницах / Е.Г. Крушель, В.П. Живоглядов, Б.М. Миркин, Н.П. Шарапов Авторское свид. №374098

60. Кузин, JI.T. Расчет и проектирование дискретных систем управления / Л.Т. Кузин. М. : ГНТИ Маш. лит., 1962. - 683 с.

61. Кузнецов, Б.Ф. Стохастические модели и методы анализа информационно-измерительных систем АСУ ТП: Монография. / Б.Ф. Кузнецов. Ангарск : Изд-во Ангарской государственной технической академии, 2007. — 180 е.: ил.

62. Куликовский, Р. Оптимальные и адаптивные процессы в системах автоматического регулирования / Р. Куликовский : Пер. с польск. М. : Hay- < ка, 1967. - 380 с.

63. Куо, Б. Теория и проектирование цифровых систем управления/ Б. Куо : Пер. с англ. М.: Машиностроение, 1986. - 448 с.

64. Куржанский, А.Б. Управление и наблюдение в условиях неопределенности / А.Б. Куржанский. -М. : Наука, 1977. 365с.

65. Летов, A.M. Аналитическое конструирование регуляторов / A.M. Летов// Автоматика и телемеханика. 1960. №¡4. - С. 436-441; №5. - С. 561-568; №6. - С. 661-665; 1961, №4. - С. 425-435.

66. Липцер, Р.Ш. Статистика случайных процессов / Р.Ш. Липцер, А.Н. Ширяев.-М. :, 1974

67. Макаров, И.М. Линейные автоматические системы / И.М. Макаров, В.М. Менский. М. : Машиностроение, 1982. - 504 с.

68. Мееров, М.В. Синтез структур систем автоматического регулирования высокой точности / М.В. Мееров. М. : Гос. изд. физ. мат. Лит, 1959. — 284 с.

69. Месарович, М. Общая теория систем: математические основы / М. Меса-рович, Я. Такахара : Пер. с англ. М. : Мир, 1978. - 312 с.

70. Методы классической и современной теории автоматического управления: учебник в 3-х т. Т.2: Синтез регуляторов и теория оптимизации систем автоматического управления / Под ред. Н.Д. Егупова. М. : Изд-во МГТУ им. Н.Э. Баумана, 2000. - 736 с.

71. Методы классической и современной теории автоматического управления: учебник в 3-х т. Т.З: Методы современной теории автоматического управления / Под ред. Н.Д. Егупова. М. : Изд-во МГТУ им. Н.Э. Баумана, 2000.-748 с.

72. Миллер, Б.М. Задача оптимального стохастического управления потоком данных по неполной информации / Б.М. Миллер, К.Е. Авраченков, К.В. Степанян, Г.Б. Миллер // Проблемы передачи информации. 2005. - № 2. -С 89-110.

73. Миллер, Б.М. Теория случайных процессов в примерах и задачах / Б.М. Миллер, А.Р. Панков. М. : ФИЗМАТЛИТ, 2002. - 320 с.

74. Мирошник, И. В. Нелинейное и адаптивное управление сложными динамическими системами / И.В. Мирошник, В.О. Никифоров, А.Л. Фрадков., СПб. : Наука, 2000. - 550с.

75. Митрахович, М.М. Интеграция методов при синтезе сложных систем в условиях априорной неопределенности / М.М. Митрахович // Автоматика. Автоматизация. Элекротехнические комплексы и системы. Херсон: Херсон НТУ, 2008. - №2(22). - С. 46-53.

76. Мишкин, Э. Приспосабливающиеся автоматические системы / Э. Мишкин, Л. Браун. : Пер. с англ. М. : ИИЛ, 1963. - 672 с.

77. Новиков Д. Механизмы функционирования многоуровневых организационных систем. М.: Фонд "Проблемы управления", 1999.- 161 с.

78. Окунев, Л.Я. Высшая алгебра / Л.Я. Окунев. — М. : Учпедгиз, 1958.

79. Основы теории управления: учебн. пособие / Ю.Ю. Громов и др.. Тамбов : Изд-во Тамб. гос. техн. ун-та, 2008. - 240 с.

80. Основы управления технологическими процессами / С.А. Анисимов и др.: под ред. Н.С. Райбмана. М. : Наука, 1978. - 440 с.

81. Острем, К. Системы управления с ЭВМ / К. Острем, Б. Виттенмарк : Пер. с англ. М.: Мир, 1987. - 480 с.

82. Пантелеев, A.B. Методы оптимизации в примерах и задачах / A.B. Пантелеев, Т.А. Летова. Изд. 2-е, испр. - М. : Высш. шк., 2005. - с.

83. Первозванский, A.A. Курс теории автоматического управления: учеб. по-соб. / A.A. Первозванский М. : Наука, 1986. - 616 с.

84. Первозванский, A.A. Математические модели в управлении производством / A.A. Первозванский. М. : Наука, 1975. - 616 с.

85. Понтрягин, Л.С. Математическая теория оптимальных процессов / Л.С. Понтрягин, В.Г. Болтянский, Р.В. Гамкрелидзе, Е.Ф. Мищенко. М. : Наука, 1976.-392 с.

86. Проектирование цементных заводов: учеб. пособие. Под редакцией Зозули П.В., Никифорова Ю.В. СПб. : Синтез, 1995 г. - 445 с.

87. Проскурня, Ю.С. Прогноз динамики продаж на основе адаптивного фильтра Калмана / Ю.С. Проскурня, Б.С. Гривко // Науковий вюник КУЕ1ТУ. 2009. № 2 (24). - С. 76-81.

88. Пугачев, B.C. Стохастические дифференциальные системы. Анализ и фильтрация / B.C. Пугачев, И.Н. Синицын. М. : Наука, 1990. Изд.2, доп. - 632с.

89. Пугачев, B.C. Теория случайных функций и ее применение к задачам автоматического управления / B.C. Пугачев. М. : Гос. изд. физ.-мат. литры., 1962.-883с.

90. Пугачев, B.C. Условно оптимальная фильтрация и экстраполяция непрерывных процессов / B.C. Пугачев // АиТ. 1984. №. - С.

91. Пытьев, Ю.П. Методы математического моделирования измерительно-вычислительных систем / Ю.П. Пытьев. М. : Наука, 2002. - 384 с.

92. Рей, У. Методы управления технологическими процессами / У. Рей. М. : Мир, 1983.

93. Ротач, В.Я. Об уточнении основных положений теории автоматического управления недетерминированными объектами / В.Я. Ротач // Теория и практика построения и функционирования АСУ : сб. науч. тр. М. : МЭИ, 1998.

94. Ротач, В.Я. Теория автоматического управления: учебник для вузов / В.Я. Ротач. 4-е изд., стереот. - М. : Издательский дом МЭИ, 2007. - 400 е., ил.

95. Сапожников, М.Я. Механическое оборудование предприятий строительных материалов, изделий и конструкций / М.Я. Сапожников. М. : Высшая школа, 1971. - 382 с.

96. Саридис, Дж. Самоорганизующиеся стохастические системы управления / Дж. Саридис. М. : Наука, 1980. - 400 с.

97. Системы автоматического управления с запаздыванием: учеб. пособие / Ю.Ю. Громов и др.. Тамбов. : Изд-во Тамб. гос. техн. ун-та, 2007. -76с.

98. Соколов, Н.И. Аналитический метод синтеза линеаризованных систем автоматического регулирования / Н.И. Соколов. М. : Машиностроение, 1966.-328 с.

99. Справочник по теории автоматического управления / Под ред. A.A. Кра-совского. -М. : Наука, 1987. - 712 с.

100. Срагович, В.Г. Адаптивное управление / В.Г. Срагович. М. : Наука, 1981.-384 с.

101. Сысоев, Л.П. Оценка параметров, обнаружение и различение сигналов / Л.П. Сысоев. М. : Наука, 1969. - 230с.114115,116,117,118119120121122123124125126

102. Теория автоматического управления: учеб. пособие для вузов: В 2-х ч.:

103. Под ред. Воронова A.A. М. : Высшая школа, 1977. - 303 с.

104. Теория вероятностей: учеб. для вузов. 3-е изд., испр. / A.B. Печинкин,

105. О.И. Тескин, Г.М. Цветкова и др.; Под ред. B.C. Зарубина, А.П. Крищенко. М. : Изд-во МГТУ им. Н.Э.Баумана, 2004. - 456 с.

106. Тюкин, И.Ю. Адаптация в нелинейных динамических системах / И.Ю.

107. Тюкин, В.А. Терехов. М. : ЛКИ, 2008. - 384 с.

108. Управление и наведение беспилотных маневренных летательных аппаратов на основе современных информационных технологий / К.К. Вереме-енко и др. / Под ред. М.Н. Красилыцикова, Г.Г. Себрякова. М. : ФИЗМАТЛИТ, 2003. - 280 с.

109. Фельдбаум, A.A. Основы теории оптимальных автоматических систем /

110. A.A. Фельдбаум. М. : Физматгиз, 1963. - 552 с.

111. B.Н. Фомин. Л.: Изд-во ЛГУ, 1985. - 336 с.

112. Фомин, В.Н. Рекуррентное оценивание и адаптивная фильтрация / В.Н. Фомин. М.: Наука, 1984г. - 288 с.

113. Черноусько, Ф.Л. Вариационные задачи механики и управления / Ф.Л. Черноусько, Н.В. Баничук. М.: Наука, 1973. - 237с.

114. Черноусько, Ф.Л. Оптимальное управление при случайных возмущениях / Ф.Л. Черноусько, В.Б. Колмановский. М. : Наука, 1978. 352 с.

115. Шейнин, А М Цементобетон для дорожных и аэродромных покрытий / A.M. Шейнин. -М. : Транспорт, 1991. 151с.

116. Athans, М. The Role and Use of the Stochastic Linear-Quadratic-Gaussian Problem in Control System Design / M Athans. IEEE Trans. Automat. Contr., 1971. AC-16. P. 529—552.

117. Birge, J.R. Introduction to Stochastic, Programming, Springer / J.R. Birge, F. Louveaux-NY, 1997. Vol. 1,3.

118. Cramer, H. Stationary and Related Stochastic Processes. Sample Function Properties and Their Application / H. Cramer, M. Leadbetter. New York, London, Sidney : John Wiley, 1967.

119. Fradkov, A. L. Nonlinear and Adaptive Control of Complex Systems / A.L. Fradkov, I.V. Miroshnik, V.O. Nikiforov. Series: Mathematics and Its Applications. Vol. 491. - Kluwer, Dordrecht, 1999. - P 528.

120. Grewal, M.S. Kalman Filtering Theory and Practice / M.S. Grewal, A.P. Andrews : John Wiley Sons, Inc., 2002.

121. Gusev, M.I. On the stability of solution of the inverse problems in control system dinamics / M.I. Gusev // Probl. Control and Inform. Theory. 1988. -Vol.17, №.5.-P.297-310.

122. Jazwinski, A.H. Stochastic processes and filtering theory / A.H. Jazwinski. -NY. : Academic Press, 1970.

123. Josef, P.D. Optimum design of linear multivariate digital control systems / P.D. Josef, J.T. Tou // On linear control theory. 1961. AIEE Trans, on Appl. and Ind., pt. II, v.80. - P. 193-198.

124. Kruglikov, S.V. On the separation principle in guaranteed control / S.V. Krug-likov // Modeling, Estimation and Control of Systems with Uncertainty. Boston etc.: Birkhauser, 1991. - Vol. 335, №.5. - P. 570-572.

125. Landau, I.D. Adaptive control systems: the model reference approach / I.D. Landau. N.Y.: Marsel Dekker, 1979.

126. Luenberger, D.G. Canonical forms for linear multivariable systems / D.G. Lu-enberger // IEEE Trans. Automatic Control, Vol. AC- 12, №3, 1967. P. 290293.

127. Sastry, S.S. Adaptive Control of Linearizable Systems / S.S. S as try, A. Isidori. // IEEE Trans. Automatic Control, Vol. AC- 34, №11, 1989. P. 1123-1131.

128. Simon, D. Kaiman Filtering / D. Simon. Embedded Systems Programming, 2001.-P. 72-79.

129. Utkin, V.l. Sliding mode control in discrete-time and difference systems / V.l. Utkin // Variable Structure and Lyapunov Control / Ed. by A.S.I Zinober. Springer-Verlag, 1994.