close

Вход

Забыли?

вход по аккаунту

?

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

код для вставкиСкачать
Кибернетика. Информационные системы и технологии
Библиографический список
1. Кухарев Г.А. Биометрические системы: Методы и сред5. Шапиро Л., Стокман Дж. Компьютерное зрение. М.: Бином,
ства идентификации личности человека. СПб.: Политехника,
2006. 752 с.
2001. 240 с.
6. Вуде Р., Гонзалес Р. Цифровая обработка изображений.
2. Даминов О.А. и др. Выделение признаков изображений
М.: Техносфера, 2005. 752 с.
отпечатков пальцев при идентификации личности // При7. Яне Б. Цифровая обработка изображений. М.: Техносфекладная информатика и компьютерное моделирование: мат.
ра, 2007. 583 с.
Всероссийской научно-практ. конф. Уфа, 2012. Т.2. С.46–48.
8. Эсонтурдиев М.Н., Аширханов А.К. Алго-ритмы предвари3. Журавлев Ю.И. Избранные научные труды. -М.: Магистр,
тельной обработки изображений отпечатков пальцев // Ин1998. 420 с.
теллектуальные системы принятия решений и проблемы
4. Абдугаффаров И.А. Системы обработки, анализа и иденвычислительного интеллекта: мат. междунар. научн. конф.
тификации отпечатков пальцев: автореф. дис. … канд. техн.
Херсон: ХНТУ, 2012. С.443–444.
наук. Ташкент: НТЦ “Современные информационные технологии” АН РУз, 2002. 18 с.
УДК 676
МОДЕЛИРОВАНИЕ АЛГОРИТМОВ УПРАВЛЕНИЯ ПРИ ИСПОЛЬЗОВАНИИ СЛУЧАЙНЫХ
ВЕЛИЧИН
© В.К. Соломина1
Иркутский государственный технический университет,
664007, Россия, г. Иркутск, ул. Лермонтова, 83.
Рассмотрена задача моделирования алгоритмов адаптивного управления с использованием случайных величин.
Для моделирования разработан алгоритм адаптивного управления процессом обжига извести с использованием
разностного динамического уравнения прогнозирования выходной координаты объекта. Параметры модели корректируются по мере поступления информации о входах и выходе объекта. Полученные результаты позволяют
оценивать работу алгоритма в режиме управления.
Ил. 2. Табл. 1. Библиогр. 1 назв.
Ключевые слова: алгоритм; моделирование; случайные величины; обжиг извести.
CONTROL ALGORITHMS MODELING WHEN USING STOCHASTIC VARIABLES
V.K. Solomina
Irkutsk State Technical University,
83 Lermontov St., Irkutsk, Russia, 664074.
The paper deals with the problem of adaptive control algorithm modeling with the use of stochastic variables. The alg orithm of adaptive control of lime burning has been developed with the use of dynamic difference equation for predicting
the object output coordinate. Model parameters are adjusted with the information on the object inputs and outputs. The
obtained results allow to assess the algorithm performance in the control mode.
2 figures. 1 source.
Key words: algorithm; modeling; stochastic variables; lime burning.
Создание автоматизированных систем в отраслях
народного хозяйства потребовало разработки сравнительно простых, но высокоэффективных (по точности,
быстродействию, приспособлению к изменяющимся
условиям) алгоритмов управления. Для достижения
этих целей создана и развивается теория и практика
адаптивного управления. Все работы этого направления можно условно разделить на две группы. Для
первой группы характерно использование регуляторов
фиксированной структуры с адаптивной (по мере поступления новой информации) коррекцией их параметров (рис. 1). Особенностью второй группы [63] является использование адаптивной модели (с переменными непрерывно корректируемыми параметрами)
для непрерывного синтеза (на основе модели и локальных критериев оптимальности) управляющих воз-

действий (рис. 2). Через x обозначена желаемая
траектория движения объекта. Преимущества второго
подхода наиболее сильно проявляются при управлении распределенными в пространстве технологическими процессами. Именно этот класс объектов рассматривается в работе. Так как датчики для основных
переменных процесса находятся в определенных точках распределенного технологического процесса, то
динамику его удобно описывать разностными уравнениями с чистыми запаздываниями по возмущающим и
управляющим входным воздействиям. На каждом такте управления во времени по текущим (и, естественно, предшествующим) измерениям входов и выходов
объекта производится коррекция параметров (из-за
этого параметры модели становятся переменными).
___________________________
1
Соломина Вера Константиновна, кандидат технических наук, доцент кафедры информатики, тел.: (3952) 405183,
89641268859.
Solomina Vera, Candidate of technical sciences, Associate Professor of the Department of Computer Science, tel.: (3952) 40518 3,
89641268859.
34
ВЕСТНИК ИрГТУ №12 (71) 2012
Кибернетика. Информационные системы и технологии
z
h
u
h
Объект
управления
x
Регулятор
фиксированной
структуры
u
z
Объект
управления
x
Синтезируемый
регулятор
x*


x*
Блок
перестройки
параметров
модели
Блок
перестройки
параметров
регулятора
Устройство управления
Рис. Рис.
В.11
Затем при расчете управлений производится (на основе той же модели) прогнозирование выходных переменных объекта с учетом чистых запаздываний по
управляющим входным переменным. А при первом
подходе для учета запаздываний надо в регулятор
включать упредители. Сложность их синтеза и реализации намного превышает сложность всех алгоритмов
адаптивного управления (с идентификацией), получаемых при реализации второго подхода. Здесь подчеркнуто только одно из преимуществ второго подхода, хотя имеются и другие. Реализация адаптивных
систем с идентификатором (АСИ) позволяет осуществлять достаточно эффективное управление малоизученными объектами, для которых не были известны ни структура, ни параметры их моделей. В
процессе реализации АСИ выяснилось также, что одни и те же алгоритмы работы АСИ можно применять
для широкого спектра объектов различной природы.
Это обстоятельство позволяет решать задачу автоматизированного
проектирования
математического
обеспечения АСИ. Включение идентификатора в цепь
обратной связи приводит к тому, что идентификация
происходит в замкнутом контуре управления в реальном масштабе времени и главная цель идентификации направлена только на достижение цели управления (например, осуществление движения выходных
координат системы (объект–устройство управления)

по заданной траектории x с минимальными уклонениями от нее).
В отличие от работ, основанных на использовании
разностных аналогов соответствующих уравнений в
частных производных, в работе промышленные распределенные динамические объекты описываются
существенно более простыми разностными уравнени-
Устройство управления
Рис. В.2Рис. 2
ями с различными чистыми запаздываниями по входным воздействиям (управляющим и наблюдаемым
неуправляющим), и для этого подхода необходимо
строить новые алгоритмы адаптивного управления.
Они должны быть просты в реализации и достаточно
эффективны.
Обжиговые вращающиеся печи являются объектом с распределенными параметрами. В процессе
обжига извести во вращающихся печах изменяются
условия протекания технологического процесса. Из-за
этого возникают пережог материала и значительные
колебания качества извести. Необходимо так управлять процессом, чтобы, несмотря на изменяющиеся
условия, обеспечивать заданное качество продукта.
Указанные выше сложные задачи управления потребовали включения в алгоритмическое и программное обеспечение АСУ ТП новых алгоритмов адаптивного управления с идентификацией.
Решена задача управления обжигом извести во
вращающихся печах. Показателем качества извести
является x(t ) [CaO%] . Собраны экспериментальные данные в процессе эксплуатации печи и методом
статистического анализа построена адекватная дискретно-временная динамическая модель с чистыми
запаздываниями во входных (управляемом и наблюдаемом) воздействиях:
y(t )   0  1 x(t  1)   2 z2 (t  1) +
+  3u1 (t  1)   4u2 (t  1) +  5u3 (t  1) + (1)
+  6 z3 (t  2)   7 u4 (t  3)   8 z5 (t  3) .
Здесь x ( CaO% ) – качество конечного продукта; u1 – расход мазута; u 2 – расход воздуха I; u 3 –
ВЕСТНИК ИрГТУ №12 (71) 2012
35
Кибернетика. Информационные системы и технологии
расход воздуха II; u 4 – расход шлама; z1 – температура в зоне обжига, измеряемая термопарой; z2 –
температура в зоне обжига, измеряемая пирометром;
z3 – температура (измеряемая с помощью термопары) в зоне сушки; z 4 – разрежение в зоне сушки;
z5 –
плотность шлама.
Модель (1) положена в основу динамической модели с перестраиваемыми параметрами

y(k |  (t ))   0 (t )  1 (t ) x(k  1)   2 (t ) z2 (k 1) +
+  3 (t )u1 (k  1)   4 (t )u2 (k  1) +  5 (t )u3 (k  1) +
+  6 (t ) z3 (k  2)   7 (t )u4 (k  3)   8 (t ) z5 (k  3) ,(2)
используемой при синтезе алгоритма адаптивного
управления (с текущей идентификацией (на t -ом так
те) вектора параметров  (t ) ). Аргумент k этой модели принимает значения k  t , t  1,  при иден-

 (t ) . При этом за счет

наилучшего подбора параметров  (t ) минимизирутификации параметров
ются невязки выходов объекта

x(k ) и модели
y(k |  (t )) на одном или нескольких тактах работы
системы. При синтезе управлений u1 (t )  u4 (t ) производится прогноз выхода модели (обычно на несколько тактов вперед). Индекс k модели принимает
значения t  k . Сама последовательность операций
при расчете управляющих воздействий представляет
собой основную новизну результатов работы (для
вышеуказанного процесса обжига извести, а также для
всех рассматриваемых в дальнейшем технологических процессов).
Цель управления состоит в наилучшем приближе*
нии выхода объекта x к желаемому значению x в
каждый текущий момент времени. В работе приводятся результаты численных исследований алгоритма
адаптивного управления с идентификацией.
Осуществлен анализ использования проекционного алгоритма идентификации и рекуррентного алгоритма наименьших квадратов в контуре управления.
По результатам испытаний делается вывод, что при
наличии плохих начальных приближений параметров
лучше
использовать
рекуррентный
алгоритм
наименьших квадратов. Этот алгоритм также более
устойчив к помехам. При работе алгоритмов в окрестности "оптимального" режима (при решении типичной
задачи стабилизации) и при сравнительно малом
уровне помех оба алгоритма обеспечивают получение
примерно одинаковых итоговых результатов, но проекционный алгоритм более прост в вычислительном
отношении.
Любая исследовательская, проектная деятельность связана с моделированием. При описании
сложных объектов, на работу которых влияют случайно действующие факторы, применяют процедуру
идентификации, когда модель строится по результа-
36
там наблюдений за входными и выходными переменными объекта. Идентификация связана с методами
статистической обработки результатов наблюдений.
Разработаны лабораторные работы, в которых
предлагается исследовать методы первичной статистической обработки информации. Работа состоит из
двух частей, в первой рассматриваются методы имитационного моделирования случайных величин, процессов, во второй – соответствующие процедуры оценивания.
Случайные величины, их распределение и
числовые характеристики. Случайной величиной
называется такая переменная, которая в результате
испытаний принимает то или иное заранее непредсказуемое значение. Для количественного описания случайной величины необходимо задать область ее существования и способ количественного определения
вероятности попадания величины в произвольную
часть области существования, то есть установить
связь между значением случайной величины и вероятностью появления этого значения. Последнее задается с помощью закона распределения случайной
величины.
Законы распределения могут быть дискретными
(для дискретных случайных величин) и непрерывными
(для непрерывных случайных величин). Случайная
величина имеет дискретное распределение, если может принимать конечное число значений x1 , x2 ,...xn с
px1 , px2 ,.... pxn  и дискретный
закон полностью определяется совокупностью xi и
pxi  .
вероятностями
Для непрерывных случайных величин вопрос: какова вероятность того, что данная величина примет
точно значение x – лишен смысла, и закон характеризуется плотностью вероятности:
f x   lim
x0
px  x  x  x 
.
x
Функция F x  , характеризующая вероятность того, что случайная величина не превзойдет значения
x , называется интегральной функцией распределения (интегральным законом):
F x  
x
 f x dx .
x
Тогда f x  называется дифференциальной
функцией распределения (дифференциальным законом). Кроме формул (1) и (2) связь между f x  и
F x  определяется так:
f x  
d
F X  .
dx
Законы распределения характеризуются следующими свойствами:
ВЕСТНИК ИрГТУ №12 (71) 2012
Кибернетика. Информационные системы и технологии
f x   0;
px1  x  x2  
x2

x1

  f x dx;
 f x  dx  1
F x2   F x1 ,
,
DCx   C 2 Dx  ;
Da  bx   b2 Dx  .
x2  x1 ;
если
Px1  x  x2   F x2   F x1 ;
F    1;
m1 a  bx   a  bm1 x  ;
DC   0 ;
F   
.
Закон распределения в интегральной или дифференциальной формах задается с помощью таблиц,
аналитически или графически.
Закон распределения является исчерпывающей
характеристикой одномерной случайной величины.
Однако в большинстве случаев достаточно указать
лишь некоторые числовые характеристики этого распределения, среди которых наиболее часто используются моменты, определяемые выражениями:
для непрерывного распределения

M g x  
 g x  f x dx ;

для дискретного
n
M g x    g xi pxi  .
i 1
Для заданного вида распределения f x  момент
M g x  является неслучайным числом и для
g x   x k называется начальным моментом k  го
mk x  , а для
порядка
и
обозначается
g x   x  m1  – центральным моментом k  го
порядка M k x  . Наиболее важными числовыми хаk
рактеристиками являются следующие (для непрерывных величин):
математическое ожидание (среднее значение)
 xf x dx ;

дисперсия

D x   M 2 
 x  m  f x dx ;
2
1

асимметрия
1 
M 5 x 
M2
эксцесс
2 
3
2
x 
;
M 4 x 
.
M 22 x 
Математическое ожидание и дисперсия обладают
следующими свойствами:
m1 c   c ;
m1 Cx   Cm1 x  ;
M 2 x   Dx   m2 x   m12 x  .
Моменты связаны с параметрами законов распределения.
Моменты имеют физический смысл. Математическое ожидание m1 x  определяет среднее значение
случайной величины и является абсциссой центра
тяжести фигуры, ограниченной кривой f x  . Центральный момент второго порядка (дисперсия) является характеристикой разбросов (вариации, оси) случайной величины относительно среднего значения.
Показатель  1 характеризует симметрию кривой
f x  относительно среднего значения. Если  1  0 ,
то асимметрия правосторонняя, если  1  0 , – то
левосторонняя. Коэффициент аксцесса (крутость)
 2  0 у нормального закона распределения дает
возможность сравнивать вершину симметричной кривой f x  с нормальным законом. Если  2  0 , то
вершина данного закона более острая и высокая, чем
у нормального с тем же средним значением и дисперсией.
Среди немоментных характеристик наибольшее
распространение получили:
медиана M e , определяемая как такое значение
x , для которого
ME

f x dx 


m1 
Центральные моменты могут вычисляться через
начальные, например:
мода

 f x dx  0,5 ;
ME
M O , равная такому x , при котором f x  до-
стигает максимального значения.
Для нормального закона распределения справедливо соотношение m1  M O  M e .
При моделировании работы различных объектов
управления необходимо искусственно получать случайные величины, распределенные по заданным законам распределения. Это достигается с помощью
применения таблиц случайных чисел или построения
специальных генераторов, которые реализуются либо
физически, либо программно на ЭВМ. В последнем
случае величины, генерируемые с помощью рекуррентных алгоритмов, называются псевдослучайными,
так как в отличие от случайных получаются по определенным правилам. Псевдослучайные числа можно
получать различными способами, например путем
деления одного числа на другое и использования
остатков возведения двух n - значных чисел в квадрат и выбора n средних цифр (метод Неймана), ис-
ВЕСТНИК ИрГТУ №12 (71) 2012
37
Кибернетика. Информационные системы и технологии
пользования разложения иррациональных чисел и т.п.
Полученные псевдослучайные числа проверяются на
случайность с помощью специальных критериев.
Чаще всего на ЭВМ реализуется генератор
(датчик) равномерно распределенных на интервале
(0,1) случайных чисел xoi . В состав стандартного математического обеспечения входит процедура генерации xoi с именем RND(x).
Моделирование псевдослучайных чисел, распределенных по нормальному закону с нулевым математическим ожиданием и единичной дисперсией, осуществляется по алгоритму
12
xN   xoi  6 .
i 1
Формула следует из центральной предельной
теоремы Ляпунова, согласно которой случайная величина, состоящая из суммы однородных случайных
величин, распределена по нормальному закону.
Для получения (из равномерной xoi ) случайной
величины
2

1 N 

S  D x  
 xi  x  ,

N  1 i1 

где N – объем выборки; xi – значения случайной ве2

личины. Для определения закона распределения
строится эмпирическая гистограмма, которая представляет собой графическое изображение интегрального вариационного ряда в виде прямоугольников.
Способы получения случайных величин
Наименование закона
Способ получения
случайной величины
Нормальный
x  M  sx
Н
x  M  S  M xo
Равномерный
Релея
x  S  2 ln xo
Односторонний
нормальный
Коши
x  S xН
x  Stg xo  0,5
x  S sin  xo  0,5
Арксинуса
xi с заданной плотностью распределения
f x  обычно используют обратную функцию
F 1 x  к интегральному закону F x  , задаваемому
формулой, и xi вычисляют в соответствии с соотношением
xi  F 1 xoi  .
F 1 xoi  для некоторых законов
Вычисление
осуществляется достаточно легко, однако во многих
случаях формула не дает аналитического решения,
поэтому используют механизм появления того или
иного закона распределения. Например, величина xi ,
равная сумме двух независимых равномерно распределенных случайных чисел, распределена по закону
Симпсона, следовательно, для получения xi может
быть использована формула
xi  xoi1xoi .
Формулы для получения псевдослучайных чисел,
распределенных по ряду законов, приведены в таблице.
Обработка и анализ выборки случайной величины,
полученной экспериментально или путем моделирования работы объектов, сводятся к вычислению оценок моментных характеристик, определению вида закона распределения и числовых значений его параметров. Вычисление моментов производится с учетом
дискретного характера выборки, например:
После построения гистограммы вычисляют моменты по сгруппированным данным: по виду гистограммы или из теоретических соображений (иногда
используют специальные диаграммы Ринда) выбирают тип распределения. По формулам, связывающим
моменты с параметрами, вычисляют оценки параметров закона распределения.
Для проверки гипотезы о близости теоретического
закона распределения к эмпирической гистограмме
применяют специальные критерии согласия, которые
отвечают на вопрос: вызвано ли расхождение распределений случайными причинами, связанными с недостаточным объемом выборки, или оно существенно и
выбранный теоретический закон не соответствует реальному распределению. Наиболее распространенными являются следующие критерии.
Критерий согласия Пирсона (  -квадрат). Вычисляется оценка
2
  
 nj  nj 
2
 ,
 
n
j 1
j

где
n j – эмпирические частоты (из гистограммы); n j –
теоретические,
которые
определяются
как
n j  f x j x ; x j – середина интервала.
Затем по таблице вероятностей для
2
при дан-
 
ном числе степеней свободы r находится P  , и
если она значительно отлична от нуля, то расхождение считается случайным и гипотеза о соответствии
теоретического закона выполняется. Число степеней
свободы
r  k  l 1 ,
где l – число параметров закона распределения.
2
для оценки среднего значения


x  m1 x  
для дисперсии
38
1 N
 XI ;
N I 1
ВЕСТНИК ИрГТУ №12 (71) 2012
Кибернетика. Информационные системы и технологии
Критерий Пирсона применяют, если
N  50
и
к  5.
Критерий согласия Колмогорова  оценивает
близость теоретического и эмпирического законов.
Недостатком рассмотренных критериев является
отсутствие прямого ответа на вопрос о применимости
к эмпирической гистограмме выбранного закона распределения. Этот недостаток устраняется при исполь-
зовании критерия Б.С. Ястремского, который базируется на
 2 - критерии Пирсона:
I
2 k
2k  4Q
,
где Q при k  20 принимается равным 0.6.
Если I<3, то гипотеза о близости законов распределения принимается.
Библиографический список
1. Рубан А.И., Соломина В.К. Адаптивное управление обжигом извести // Бумажная промышленность.1991. №12. С.21 –23.
ВЕСТНИК ИрГТУ №12 (71) 2012
39
Документ
Категория
Без категории
Просмотров
15
Размер файла
2 635 Кб
Теги
величины, моделирование, алгоритм, случайных, использование, pdf, управления
1/--страниц
Пожаловаться на содержимое документа