close

Вход

Забыли?

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

?

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

код для вставкиСкачать
Управление, вычислительная техника и информационные технологии
УДК 681.51.015
ИСПОЛЬЗОВАНИЕ МЕТОДА ГРУППОВОГО УЧЕТА
АРГУМЕНТОВ ДЛЯ ВЫБОРА СТРУКТУРЫ МОДЕЛИ
ДИНАМИЧЕСКОГО ОБЪЕКТА
В.М. Понятский, С.И. Велешки, А.В. Жирнова
Рассмотрены вопросы оценки возможности использования метода группового
учета аргументов для оценки структуры модели динамического объекта. Проведенные исследования показали возможность использования предложенного подхода для
оценки структуры объекта в виде типового динамического звена.
Ключевые слова: модель, оценка, структура.
Введение
При обработке результатов экспериментов важной проблемой является оценка структуры модели исследуемого динамического объекта. Рассмотрим возможность использования для решения указанной проблемы
метода группового учета аргументов.
Пусть объект представляется в виде совокупности присущих ему
состояний S1 , S2 , …, Sk , по которым можно судить о его поведении внутри
системы либо за ее границами.
Каждое
состояние характеризуется
набором
параметров
( x1 , x2 , …, xn , xn 1 ), номенклатура которых устанавливается:
Si  Si (x1 , x2 , …, xn , xn 1 ) .
Среди рассматриваемых параметров выделяются несколько наиболее значимых, которые олицетворяют цель исследования. Именно такой
смысл вкладывают в термин результативный признак и обозначают его
Y0 . Не нарушая общности, можно считать, что таким параметром является
xn1 , поэтому Yo  xn 1 . Исходная информация представляется в виде таблицы.
Исходная информация
Номер
состояния
…
S1
S2
Sk
Параметры системы
x1
x11
x 21
…
…
x2
x12
x 22
…
…
…
…
…
xn
x12
x 2n
…
Yo  xn 1
Y1
Y2
…
x k1
x k2
x kn
Yk
Ставится задача определить математическую зависимость результативного признака Y0 с набором управляемых параметров x1, x2 , …, xn ,
которые предполагаются независимыми.
255
Известия ТулГУ. Технические науки. 2013. Вып. 2
Построение и анализ эффективного уравнения Y   (x o ) , собственно, и являются главным моментом в теории статистических методов,
так как именно результативный признак используется в дальнейшем для
управления объектом, вывода системы на экстремальный уровень, прогнозирования за границами системы и реализации иных приложений.
Для создания исходных данных для метода группового учета аргументов используются разностные уравнения.
Метод группового учета аргумента
При малой величине выборки, когда при построении регрессионных моделей использовать статистические гипотезы о плотности распределения, например, гипотезу о Гауссовском распределении, невозможно,
используется индуктивный подход. В соответствии с этим методом последовательно порождаются модели возрастающей сложности до тех пор, пока не будет найден минимум некоторого критерия качества модели. Этот
критерий качества называется внешним критерием, так как при настройке
моделей и при оценке качества моделей используются разные данные. Достижение глобального минимума внешнего критерия при порождении моделей означает, что модель, доставляющая такой минимум, является искомой.
Описание метода
Индуктивный алгоритм определения модели оптимальной структуры состоит из следующих основных шагов.
Шаг 1. Пусть задана выборка
D   x n , y n  nN1 , x  R P .
Выборка разбивается на тестовую (L) и обучающую (C). Эти множества удовлетворяют условиям разбиения L  C  W , L C   . Матрица
XL состоит из тех векторов–строк xn , для которых индекс n L . Вектор
YL состоит из тех элементов yn , для которых индекс n L . Обычно обучающая выборка составляет 70 % из всей выборки.
Вследствие разбиения получаем
 XL


, Y W   Y L , L  C  W .
X W  
 XC
Y C 
Шаг 2. Назначается базовая модель. Под базовой моделью понимается вид отношения между зависимой переменной y и свободными переменными x с базисным набором опорных функций, используемым в дальнейшем для формирования моделей различной сложности.
Шаг 3. Исходя из поставленной задачи выбирается внешний критерий, описывающий качество модели.
Шаг 4. Порождаются модели–претенденты на основе получения
частных описаний из базовой модели (скрещивание) и отбора лучших из
них (селекция).
256
Управление, вычислительная техника и информационные технологии
Шаг 5. Настраиваются параметры моделей. Для настройки используется внутренний критерий – критерий, вычисляемый с использованием
обучающей выборки. Каждому элементу вектора xn – элементу выборки D
– ставится в соответствие вектор an . Строится матрица A w – набор векторов–столбцов an . Матрица A w разбивается на подматрицы A L и A C .
Предполагается, что наименьший квадрат невязки
|| yC  yCo ||2  ( yC  yCo )Т ( yC  yC 0 ) , где yCo  AC wCo
доставляет значение вектора параметров wCo , которые вычисляются методом наименьших квадратов:


1
wCo  ACT AC ACT yC .
Шаг 6. Вычисляется качество порожденных моделей. При этом используются тестовая выборка и назначенный внешний критерий.
Шаг 7. Модель, доставляющая минимум внешнему критерию, считается оптимальной.
Если значение внешнего критерия не достигает своего минимума
при увеличении сложности модели или значение функции качества неудовлетворительно, то выбирается лучшая модель из моделей заданной сложности. Под сложностью модели подразумевается число настраиваемых параметров модели. Существуют следующие причины, по которым глобальный минимум может не существовать:
1) данные слишком зашумлены,
2) среди данных нет необходимых для отыскания модели переменных,
3) неверно задан критерий выбора,
4) при анализе временных рядов существует значительная временная задержка отыскиваемой причинно–следственной связи.
Внешние критерии
Критерий выбора модели можно назвать внешним, если он получен
с помощью дополнительной информации, не содержащейся в данных, которые использовались при вычислении параметров моделей. Например,
такая информация содержится в дополнительной, тестовой выборке.
Алгоритм МГУА использует и внутренний критерий, и внешний.
Внутренний критерий используется для настройки параметров модели,
внешний критерий используется для выбора модели оптимальной структуры. Возможен выбор модели по нескольким критериям.
Самые распространённые критерии: критерий регулярности и критерий минимального смещения.
Критерий регулярности ∆2(С) включает среднеквадратическую
ошибку на обучающей подвыборке С, полученную при параметрах модели, настроенных на тестовой подвыборке L:
257
Известия ТулГУ. Технические науки. 2013. Вып. 2
2
2 (C )  yC  AC wLo .
Критерий ∆2(С) также обозначается ∆2(С\L), то есть ошибка на подвыборке С при параметрах, получаемых на подвыборке L.
Критерий минимального смещения (критерий непротиворечивости)
модели: модель, которая имеет на обучающей выборке одну невязку, а на
тестовой – другую, называется противоречивой. Этот критерий включает
разность между зависимыми переменными модели, вычисляемыми на двух
различных выборках L и C. Критерий не включает ошибку модели в явной
форме. Он требует, чтобы оценки коэффициентов в оптимальной модели,
вычисленные на множествах L и С, различались минимально.
Критерий непротиворечивости как критерий минимума смещения
имеет вид
2
2  wLo  wCo .
Комбинаторный алгоритм МГУА
В комбинаторных алгоритмах выполняется перебор всевозможных
частных моделей из заданного полиномиального базиса с выбором лучшей
из этих моделей по заданному внешнему критерию. При переборе сложность частных моделей наращивается от 1 до максимального числа n (числа аргументов базисного набора функции). Таким образом, общая схема
комбинаторного алгоритма включает следующие операции: по методу
наименьших квадратов (МНК) определяются коэффициенты всех частных
моделей при сложности s  1, s  2, ..., s  n ; для каждой из них вычисляется значение внешнего индивидуального или комбинированного критерия
селекции; единственная модель оптимальной сложности выбирается по
минимальному критерию.
Структура комбинаторного алгоритма
В структуре алгоритма выделяются три основных блока:
1) преобразования исходных данных в соответствии с выбранной
системой опорных (базисных) функций, в которой ищется модель;
2) генерирования (перебора) полного или неполного множества
усложняющихся частных моделей в выбранном базисе;
3) вычисления значения некоторого критерия селекции, имеющего
свойства внешнего дополнения, и последовательного отбора частных моделей, лучших по этому критерию.
Формирование базиса
Если заданы измерения некоторых входных переменных
z1 , z1 , ..., z v моделируемого объекта и максимальная степень полинома, то
число слагаемых n в полном полиноме степени  max от v переменных
определяется однозначно:
258
Управление, вычислительная техника и информационные технологии
i
n    max .
i
i 1
Полный полином базисной модели записывается в следующем общем виде:
v
n
v
n
i 1
j 1
i 1
y   a i  z j i j   a i x i ,
где каждый обобщенный линейный аргумент x i является нелинейной
функцией исходных переменных z j :
v
xi   z j ij .
j 1
Для формирования степеней аргументов удобно организовать процедуру получения последовательности v –разрядных чисел σ i=(σ1i, σ2i,…,
σvi) с основаниями 1, 2, ..., max  1 и отбирать те из них, которые удовлетворяют указанному ограничению. Например, при v  2 (два аргумента) и
 max  2 получим последовательность n=6 чисел  i , i  1, 2, ..., 6 :
00,10,01,20,11,02, что соответствует полному полиному
y  a1  a 2 z1  a 3z 2  a 4 z12  a 5 z1z 2  a 6 z12 .
Итак, члены xi полного полинома являются базисным набором
опорных функций для комбинаторного алгоритма. Описанная процедура
позволяет сформировать матрицу измерений обобщённых аргументов
X [ N x n] , где N – число точек измерений.
Матрицу X следует разделить на обучающую C длиной N C и тестовую L длиной N L , причем N C  N L  N .
Генерирование (перебор) частных моделей
Основными операциями, выполняемыми в этом блоке, являются:
формирование структуры очередной частной модели; формирование соответствующей нормальной системы уравнений; решение полученной системы (оценка коэффициентов модели).
Формирование структур частных моделей формализуется с помощью двоичного структурного вектора d  (d1 , d1 , ..., d n ) : если элемент d i
этого вектора принимает значение 1, то соответствующий i –й аргумент
включается в частную модель, если – 0, то не включается i  1, 2, ..., n .
Изменение состояний вектора d можно организовать многими способами, но наиболее простым по своей идее является следующий способ:
получать все возможные варианты размещения в векторе d сначала одной
единицы (всего C1n  n вариантов), затем двух (всего C1n  0.5nn  1 вари-
259
Известия ТулГУ. Технические науки. 2013. Вып. 2
антов) и т.д. вплоть до n единиц ( C nn  1 вариант). Таким образом, на k –й
селекции число генерируемых моделей будет равно C kn .
В соответствии с приведенной схемой алгоритм работает так. Сначала определяются все модели при s  1 , т.е. состоящие лишь из одного аргумента:
q1  a1 , q 2  a 2 x1, q 3  a 3 x 2 , ... , q n  a n x n .
Далее рассматриваются всевозможные модели при s  2 , состоящие
из двух аргументов:
q1  a1  a 2 x1, q 2  a1  a 3 x 2 , ...,
q i  a1  a o x 2 x 3 , ..., q j  a 2 x1  a 3 x 3 , ...,
q k 1  a 2 x1  a o x 2 x 3 , ..., q k  a o x1x 3  a o x 2 x 3 .
Аналогично строятся частные модели при s  3 , при s  4 , и т.д. до
n  1 модели при s  n , т.е. до полного полинома.
Cn
Общее число вариантов составит p n = 2 n - 1 различных структур,
т.е. полный перебор. Однако программная реализация этого способа является далеко не быстродействующей, и такой перебор структур целесообразно применять только в алгоритмах неполного перебора.
Очевидно, что p n при увеличении n очень быстро возрастает
(например, p10 = 1024 , а p15 = 32768 ). Поэтому возможности комбинаторного перебора ограничены. Существуют некоторые способы ускорения
счета при полном переборе. Например, оптимальная схема перебора, основанная на вычислении коэффициентов последовательно получаемых частных моделей с помощью метода окаймления, позволяет увеличить n примерно до 23.
При большом n приходится применять целесообразное усечение
перебора. Например, известен способ усечения треугольника перебора:
рассматриваются все модели со сложностью от s  1 до s max  m , а модели
с большим числом аргументов исключаются из перебора. Максимальная
сложность s max задается исходя из возможности ЭВМ.
Наиболее компактной и универсальной является схема изменения
двоичного вектора по принципу работы двоичного счётчика, в последний
разряд которого добавляется единица. При этом, что существенно, соблюдается взаимно однозначное соответствие между порядковым номером
очередной частной модели и состоянием структурного вектора. Этот способ алгоритмически очень прост и весьма удобен именно при полном переборе, несмотря на то, что количество и состав аргументов в частных моделях все время меняются. Более того, он позволяет даже организовать рекуррентную процедуру перебора.
Для формирования нормальной системы уравнений, соответствующей очередному структурному вектору, можно поступить формально: из
260
Управление, вычислительная техника и информационные технологии
столбцов полной матрицы X , указанных единичными элементами d , составляется частная матрица X1 , а затем вычисляются элементы нормальной матрицы X1T X1 , X1T Y . Однако при полном переборе этот формальный
путь является наихудшим, поскольку приводит к многократному вычислению одних и тех же скалярных произведений. Например, для структур 010,
011 и 110 трижды вычисляется одна и та же величина
N
x T2 x 2   x 22 j .
j 1
Поэтому в комбинаторном алгоритме достаточно только один раз вычислить матрицы полной нормальной системы:
 x1T x1 x T2 x1 ... x1T x n 
 x1T y
 T

 T 
x 2 x1 x T2 x 2 ... x T2 x n 
x y
T
T

X X
, X Y 2 .
 ...
 ... 
...
...
... 
 T

 T 
T
T
 x n x1 x n x 2 ... x n x n 
 x n y
Для получения любой частной нормальной системы достаточно
взять элементы матрицы XT X , находящиеся на пересечении строк и
столбцов, указанных единицами вектора d , а также соответствующие элементы вектора X T Y .
Для решения каждой нормальной системы X iT Xi a i  XiT y , т.е. для
вычисления оценок коэффициентов частной модели, можно применять
любые процедуры решения систем алгебраических уравнений (с симметричной матрицей) с хорошими вычислительными свойствами.
Отбор по внешним критериям
Вычисления в этом блоке, как и в предыдущем, организуются с
учетом применяемых внешних критериев.
Селекция (отбор) моделей, лучших по заданному критерию I ,
обычно выполняется не в конце перебора (когда получены все частные модели), а в процессе его. Для этого запоминаются значения критерия для заданного числа F первых моделей, а затем величина критерия каждой последующей модели I i сравнивается с худшим I max из F значений. Если
Ii  I max , то новая модель запоминается вместо худшей (запоминаются
структура, оценки коэффициентов и критерий); если же Ii  I max Ii , то эта
модель пропускается. После окончания перебора оставшиеся F моделей
являются лучшими из всех p n в смысле заданного критерия. Очевидно,
что при таком подходе блоки генерации и отбора алгоритмически объединяются.
Необходимо отметить, что в алгоритмах МГУА обычно выполняется еще один этап вычислений – оценка качества отобранных лучших моде261
Известия ТулГУ. Технические науки. 2013. Вып. 2
лей. При этом вычисляются, например, среднеквадратическая ошибка аппроксимации и ошибка экстраполяции.
На рисунке приведен алгоритм комбинаторного МГУА, реализованный в среде Matlab.
Рассмотрим применение алгоритмов МГУА для оценки структуры
модели динамического объекта первого порядка.
Для уравнения первого порядка разностное уравнение имеет вид
y n  a o x n  a1x n 1  by n 1 .
Для передаточной функции W  3 /(0,01s  1) имеем следующее разностное уравнение:
(1)
y n  2x n 1  0.33y n 1 .
Вычисленная среднеквадратическая ошибка: I = 8,1562e–028.
Проведено тестирование разработанной программы, реализующей
метод МГУА. Для этого запускаем файл «MGUA.m» и вводим следующие
данные: количество селекций (countSelection = 10), количество выводов
лучшей модели (countTopModel = 7), максимальная степень полинома (extendPoly = 2).
Метод МГУА рассматривает все возможные комбинации следующего полинома:
ao  a1 xn  a2 xn1  a3 yn1  a4 xn xn1  a5 xnYn1  a6 xn1 yn1  a7 xn2  a8 xn21  ...
В результате получены следующие модели.
Model 1:
1,4709e-014 +2x n-1-0.333y n-1-1,56e-014x n-1y n-1+1,29e-014x 2n-1- 8,6275e-015y 2n-1;
I = 7,8293e–023 ; Number selection – 6;
Model 2: 2x n -1 - 0,33333yn -1 - 8,0231e - 018x n -1y n -1 ;
I =8,1135e–023; Number selection – 3;
2x n -1 - 0,33333yn -1 - 3,4694e - 016y 2n -1 ;
Model 3:
I =1,1412e–022; Number selection – 3;
Model 4:
2x n-1 - 0,33333yn-1 - 1,5458e - 013x n x n-1 - 1,5684e - 013x 2n - 1,5735e - 013 * X 2n-1 ;
I =1,1912e–022; Number selection – 5;
Model 5: 2x n -1 - 0,33333yn -1 - 1,4311e - 017x n y n -1 ;
I =1,5862e–022; Number selection – 3;
Model 6: 2x n -1 - 0,33333yn -1 ;
I =2,1771e–022; Number selection – 2;
Model 7:
3.69e-014 -7,97e-013xn+2 xn-1-0,33 yn-1+5,80e-015xn-1 yn-1+7,7217e-015 yn-2 1 ;
I =2,6826e–022; Number selection – 6.
262
Управление, вычислительная техника и информационные технологии
1
Начало
v
–
количество
параметров
системы;
σ
–
максимальная
степень
полинома;
F количество вывода лучших
моделей;
m – количество селекции (m≤ n);
X0 - матрица (N x v) содержит
значения параметров;
Y0
–
значения
системы
на
выходе.
2
Ввод:
v,σ, F , m
матрицы
X0,Y0
3
n – количество слагаемых
в полиноме
λ
– матрица (n x v) содержит значение
степеней при z членов. Организованная
процедуру получения последовательности
v-разрядных чисел;
λ = MatrixExtend(n,v,σ )
4
λ=
MatrixExtend
(n,v,σ )
X
–
Матрица
(N
x
n)
измерений
обобщенных
измерений,
т.е.
значения
z1,z2,…,zv;
X = MatrixZ(X0,N,n)
5
X=
CalcMatrixX
(X0,N,n)
Y
–
Матрица
(N+1,
F)
содержит
коэффициенты лучших F моделей, а
последней столбец значение внешнего
критерия;
Y = CalculusFunc(X,Y0,n,N,m)
6
Y=
CalculusFunc
(X,Y0,n,N,m)
7
Осуществляется вывод
моделей
вместе
с
значением
внешнего
критерия.
Вывод
лучших
моделей (Y)
8
а
конец
1
X = MatrixZ(X0,N,n)
Начало
X – Матрица
(N x n) измерений
обобщенных измерений, т.е. значения
z1,z2,…,zv;
2
i=0,N,1
3
7
j=0,n,1
Возврат
матрица X
4
8
конец
X[i,j]=1
5
k=0,v,1
6
X[i,j] *=
X0[I,k]^λ[j,k]
б
Алгоритм МГУА: а – общий алгоритм; б – алгоритм создания
обобщенной матрицы «X» (начало, см. также с. 264 и 265 )
263
Известия ТулГУ. Технические науки. 2013. Вып. 2
Y = CalculusFunc(X,Y0,n,N,m)
Y – Матрица (N+1, F) содержит коэффициенты
лучших F моделей, а последней столбец значение
внешнего критерия;
1
Начало
s – номер селекции. Цикл
выполнения m cелекции.
2
s=1,m,1
9
Возврат
массив Y
3
10
Конец
4
l – количество возможных
моделей на s-ой селекции.
D = Structure
Matrix(n,l,s)
D – структурная матрица (l,n). Содержит
структурные
векторы,
которые
соответствуют
формированных
моделей.
5
i=0,l-1,1
8
SelectionMode
l(I,A,X,s,Y)
6
Отбор F лучших моделей и их
запись в Y
Вычисление
параметров A[i]
модели
Вычисление критерия i-ой
модели.
7
I[i]=
ExternalCrite
rion(X,A[i],Y0
,N)
I[i]= ExternalCriterion(X,A[i],Y0,N)
в
SelectionModel(I,A,X,s,Y)
1
Начало
2
Есть
выбранных
модели
(Y=NULL)
Отбор F лучших моделей и их
запись в Y
Нет
Да
4
3
Сравниваем
значения критерия
старых моделей с
новых моделей
Создаем Y и в ней
записываем F модели
для которых I минимален
Конец
Создаем Y и в ней
записываем F модели
для
которых
I[i]
минимален.
Сравниваем значения критерия старых
моделей с новых моделей.
Величина критерия каждой новой модели Ii
сравнивается с худшим Imax из F значений.
Если Ii<Imax, то новая модель запоминается
вместо худшей (запоминается структура,
оценки коэффициентов и критерий); если
же Ii≥Imax, то эта модель пропускается.
г
Продолжение: в – алгоритм получения F лучших моделей;
г – отбор F лучших моделей
264
Управление, вычислительная техника и информационные технологии
D = StructureMatrix(n,l,s)
1
Начало
D – структурная матрица (l,n). Содержит
структурные векторы, которые соответствуют
формированных моделей.
2
P=2^n - 1
Вычисляем
общее
количество моделей.
Матрица всех
возможных моделей.
3
Создаем матрица
D_Full
4
Создаем матрица
D
Из матрицы D_Full выбираем все
строки в котором количество
единиц равно S. Все выбранные
строки записываются в матрицу D.
5
Возврат
матрицы D
6
д
конец
Окончание: д – алгоритм создания структурной матрицы на s–
ой селекции
Анализ полученных результатов показывает, что разностное уравнение (1) и модель номер 6 совпадают.
Проведем оценку структуры модели динамического объекта второго порядка.
Уравнение второго порядка имеет разностное уравнение вида
y n  a o x n  a1x n  2  a1x n  2  bo y n 1  b1y n  2 .
Для передаточной функции W  3 /(0,05s 2  0.35s  1) получено следующее разностное уравнение:
y n  0,29x n 1  0,29x n  2  1,93y n 1  0,93y n  2 .
(2)
Вычислим среднеквадратическую ошибку: I = 7.2436e–029.
Метод МГУА рассматривает все возможные комбинации следующего полинома:
ao  a1 xn  a2 xn1  a3 xn2  a4 yn1  a5 yn2  a6 xn xn1  a7 xn xn2  a8 xn yn1  ... .
Заданы количество селекций (countSelection = 5), количество выводов лучшей модели (countTopModel = 5), максимальная степень полинома
(extendPoly = 2).
Получены следующие структуры.
Model
1:
0,00293x n -1 + 0,0028638x n - 2 - 1,9305yn -1 + 0,93239yn - 2 ;
I =3,3942e–022; Number selection – 4;
265
Известия ТулГУ. Технические науки. 2013. Вып. 2
Model 2:
0,00293x n -1 + 0,00286x n - 2 - 1,93yn -1 + 0,932yn - 2 - 2,59e - 013x n -1y n - 2 ;
I =1,9816e–020; Number selection – 5;
Model 3:
- 2,037e - 012x n + 0,00293x n -1 + 0,0028638x n - 2 - 1,9305yn -1 + 0,93239yn - 2 ;
I =3,9365e–020; Number selection – 5;
Model 4:
0,00293x n -1 + 0,0028638x n - 2 - 1,9305yn -1 + 0,93239yn - 2 - 8,2107e - 013x n - 2 y n - 2 ;
I =5,8978e–020; Number selection – 5;
Model 5:
0,00293x n-1 + 0,0028638x n-2 - 1,9305yn-1 + 0,93239yn-2 + 1,4483e - 012x n y n-2 ;
I =6,6075e–020; Number selection – 5.
Полученные результаты показывают, что разностное уравнение (2)
и модель номер 1 совпадают.
Проведенные исследования показали, что предложенный подход
обеспечивает оценку структур типовых динамических звеньев.
Необходимо отметить, что c увеличением номера селекции рассматривается все больше и больше комбинаций входных и выходных параметров, поэтому при задании слишком большого количества селекций
может получиться так, что модель разностного уравнения не попадает в
лучшие модели.
Таким образом, рассмотрено использование метода группового учета аргумента для оценки структуры модели динамического объекта. Проведенные исследования показали возможность оценки объекта в виде типового динамического звена.
Список литературы
1. Ивахненко А.Г., Степашко В.С. Помехоустойчивость моделирования. Киев.: Наукова Думка, 1985. 216 c.
2. Себер Дж. Линейный регрессионный анализ. М.: Мир, 1980. 456
с.
3. Интернет ресурс –http:/www.gmdh.net – официальный сайт метода МГУА.
Понятский Валерий Мариафович, канд. техн. наук, доц., pwmru@rambler.ru,
Россия, Тула, Тульский государственный университет,
Велешки Стаян, студент, Россия, Тула, Тульский государственный университет,
Жирнова Алена Владимировна, студентка, Россия, Тула, Тульский государственный университет
266
Управление, вычислительная техника и информационные технологии
USE OF A METHOD OF THE GROUP ACCOUNT OF ARGUMENTS FOR A CHOICE
OF STRUCTURE OF MODEL OF DYNAMIC OBJECT
V.M. Ponyatskiy, S.I. Velechki, A.V. Girnova
Questions of an estimation of possibility of use of a method of the group account of
arguments for an estimation of structure of model of dynamic object are considered. The carried out researches have shown efficiency of the offered approach.
Key words: model, an estimation, structure.
Ponyatskiy Valery Mariafovich, candidate of
pwmru@rambler.ru, Russia, Tula, Tula State University,
technical
sciences, docent,
Velechki Stajan, student, Russia, Tula, Tula State University,
Girnova Alena Vladimirovna, student, Russia, Tula, Tula State University
УДК 66.02:519.771.3
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПИЛОТНОЙ УСТАНОВКИ
ПРОЦЕССА СУСПЕНЗИОННОЙ ПОЛИМЕРИЗАЦИИ СТИРОЛА
М.А. Сафин, А.Г. Лопатин, Д.П. Вент, В.П. Савельянов
Разработана адекватная математическая модель процесса суспензионной полимеризации стирола, которая учитывает влияние основных параметров (концентрации инициатора и температуры) на ход процесса во времени. Модель позволяет подобрать оптимальный температурный режим и создать на её основе систему адаптивного управления процессом.
Ключевые слова: математическая модель, суспензионная полимеризация, пилотная установка, стирол.
Наиболее важную и многообразную группу процессов химической
технологии составляют химические процессы, в которых в результате химической реакции происходит изменение химического состава, внутреннего строения и свойств исходных веществ. Аппараты для их осуществления
– химические реакторы - являются центральным аппаратом в химикотехнологическом процессе, включающем подготовку сырья, химическое
превращение, выделение продуктов. От условий работы в значительной
степени зависит эффективность всего химического производства. Одним
из представителей такого производства является синтез суспензионного
267
Документ
Категория
Без категории
Просмотров
30
Размер файла
813 Кб
Теги
структура, метод, выбор, учет, использование, групповой, pdf, модель, аргументы, динамическое, объекты
1/--страниц
Пожаловаться на содержимое документа