Идентификация структуры и параметров модели химической реакции методом сетевого оператора Текст научной статьи по специальности «Математика»

НАУКА и ОБРАЗОВАН И Е

Эл № ФС77 - 48211. Государственная регистрация №0421200025. КБМ 1994-0408

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

# 12, декабрь 2012

Б01: 10.7463/0113.0521994

Дивеев А. И., Семиков М. В., Софронова Е. А.

УДК 681.5.015

Вычислительный центр им. А.А. Дородницына Российской академии наук

Российский университет дружбы народов

aidiveev@mail.ru semikovmv@yandex.ru sofronova ea@mail.ru

Введение

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

В частности метод генетического программирования^] и его развитие, метод грамматической эволюции [2] и метод аналитического программирования [3] могут быть использованы для получения математических моделей различных систем и процессов. К этому же классу методов относится метод сетевого оператора [4-6]. Различие между этими методами заключается в том, что с помощью метода сетевого оператора можно получить математические выражения, включающие только унарные и бинарные операции. Методы генетического программирования, грамматической эволюции и аналитического программирования предназначены для синтеза текста программ, поэтому они могут включать любые функции и операторы. С другой стороны, метод сетевого оператора при поиске структуры и параметров математического выражения использует базисную модель для данного выражения, которая может задаваться на основании опыта. Такой подход позволяет реализовать целенаправленный поиск оптимального решения в окрестности базисного, что существенно сокращает время поиска и пространство возможных решений.

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

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

1. Постановка задачи

Рассмотрим задачу идентификации математической модели химической реакции.

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

X = (to, у0 , у н)),

>у j,...,Цн,у )Ь (1)

где ук - наблюдаемый в момент ^ вектор показателей концентрации химической реакции,

Ук ••• Укш*к-1 < гк, к = 1уN;

Ук1 =Ф 1 (х), к = 1, N, у = 1, т; (2)

ук =

х = Х ... Хп^, X! - концентрация вещества I химической реакции, фу(х) - известные

функции, описывающие значения показателя концентраций.

Известно, что в химической реакции первоначально при детальном рассмотрении механизма реакции по стадиям I < п . Остальные п — I веществ получаются в результате реакции.

Для описания изменения концентраций веществ во времени определена система дифференциальных уравнения

х = f ^ (ХЬ..., XI)), (3)

где Х1 )=[м^1 {х^ Хн ) ... , Хуг )]Г, 1 < 1]ь..., 1Г , уг < I.

Функции Ы’р (х;- , Ху ) описывают закон изменения концентраций Хг- , Ху при

взаимодействии веществ и ]' :

Р ^ Р

р ■'р

Здесь qp - неизвестный по величине постоянный параметр, р = 1, Г.

(4)

Предполагаем, что структура функции г(qp, Хг- , Ху ) и значения параметров qp,

р ■'р р = 1, Г неизвестны.

Заданы ограничения в виде алгебраических уравнений баланса

а у(х)= 0, у = 1,1. (5)

Необходимо найти структуру функций (4) и параметры qp, р = 1, Г, чтобы в результате решения системы дифференциальных уравнений (3) для заданных начальных значений

а)} - весовые коэффициенты

2. Метод сетевого оператора

Для решения задачи используем метод сетевого оператора [4-6]. Метод позволяет находить математические выражения, построенные на следующих множествах:

- множество переменных

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

Для представления сетевого оператора в памяти компьютера используем матрицу сетевого оператора

где Ь - число узлов в сетевом операторе.

Матрица сетевого оператора совпадает по структуре с матрицей смежности графа сетевого оператора. Вместо единиц матрицы смежности, указывающих на дуги графа,

*1 (°)= хЬ -, х1 (°)= X0, X+1 (°)= 0,•••,хп(°)= 0

(6)

минимизировать значение критерия

(7)

где х(^ ) - решение системы (3) в момент ^ и удовлетворять уравнениям баланса (5)

(8)

(9)

- множество параметров

(1°)

- множество унарных операций

(11)

- множество бинарных операций

(12)

(13)

установим номера связанных с соответствующими дугами унарных операций, у г- у = к, если дуга ] ) связана с операцией Р к ( г). На главной диагонали ставим номера бинарных

операций, у і і = т, если узел І связан с бинарной операцией X т (г' , г").

В сетевом операторе узлы нумеруют так, чтобы, номер узла, откуда дуга выходит, был меньше номера узла, куда дуга входит, V(i, ]) І < ] , Тогда матрица сетевого оператора будет иметь верхний треугольный вид. По матрице сетевого оператора можно вычислить значение математического выражения, для этой цели вводят дополнительные три вектора для указания номеров узлов для переменных, параметров и выходов.

Приведем наиболее часто употребляемые унарные и бинарные операции.

1_

^ Р3 (г) = -г , Р4 (г) = ^п(г)лЩ ,

Pl(z )= z, P2 (z ) =

-1 I I

S , если z >

z2, иначе,

P5 (z

sgn(z)

, если z < S,

S

, P6(z

—, иначе, z

S , если z > -ln(s).

e , иначе,

Pl(z)

ln(s), если - ln|z > —, lnl z|, иначе,

S

P8 (z )

sgn(z), если |z| > -ln(s),

1 - e

- z

1 + e

-z

иначе,

11, если г > 0, , ч , ч

Р9 (г , Р10 (г )= ^п(г),

0, иначе,

Pll(z )= cos(z), Pl2 (z ) = sin(z), Pl3 (z) arCtan(z ), P14 (z)

sgn(z) I I 1

——, если z >

z 3 , иначе

Pl5(z)

Чє, если |z| < S, ¥~z, иначе,

P16(z)4z ТГ ^<1j Pll(z)=sgn(z)ln(zl +1), sgn(z ), иначе,

Pl8(z)

sgn(z)

если z >

-ln(s),

S і I I \ P19(z)

sgn(z)\e'z' - ll иначе,

0, если г| > -Іи(^), sgn(г)е г, иначе,

1

S

1, если Z >8

0, если Z < 0 , р21 ^):

иначе

3z 2 2 z 3

„2 „3

88

1, если z > —

2

1 8 -1, если z < —

2

3

3z 4 z'

23

88

-, иначе

Р22\ z

(z ):

0, если 21 > -1п8),

— 2

е 1 ', иначе,

Р23(z)

2 =«

sgn(z) I К 1

------—, если 2 >

8

z — z , иначе,

38

Р24 (2 )=<

8

1 + 8 1

, если z > —1п(8).

1 + e

— z

, иначе,

(15)

Х0(z",2") = 2' + z ", Х1 (2',z ") = 22 , X2(2,z") = шах^Х}, X3(2',z")= ШШ^"^"},

X4(г',2")= г' + г" — 22 , Х5(г',г")= sgn(z' + + (г") ,

X 6 (2 ', 2 " )= Sgn(z ' + 2 " )(| 2 " + 2 "|), X 7 (2 ', 2 " )= Sgn (2 ' + 2 " )шах^ \\х "} . (16)

Для бинарных операций (16) имеем следующие единичные элементы:

е0 = 0, е1 = 1, е2 = —, е3 =-----, е4 = 0, е5 = 0 , еб = 0, е7 = 0 .

8 8

(17)

Матрица ^ сетевого оператора размерностью Ь х Ь построена на множествах из N переменных, Р параметров, Ж унарных операций и V бинарных операций. Матрица позволяет описывать

Ь ]—1

К0 = VЬ—N—Р п Е Ф )жк

]=N+Р+1 к=1

(18)

где

О (к ):

() — 1)

(] — 1 — к )к!

математических выражений.

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

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

б(гк, V1 )=£ ]Г эйку -У,' ] )

І = 1 ] = 1

(19)

Vк, V1 є 5(1, Ь,с, и), Vк = ^ ] , V1 = [р/у ] , І, ] = 1,1.

Если Vк, V1 є Н(£, Ь,С,и) и б(тк, V1 )= А, то говорим, что расстояние между

сетевыми операторами V и V равно А .

А-окрестностью сетевого оператора VІ є е(ь, Ь , С, и) называем подмножество

І, а)с 5(Ь,Ь, с, и), сетевых операторов, которые находятся на расстоянии не более А от

сетевого оператора

Vі , VТ] є «(VІ,а), 5^І, V] )< А.

Пусть мы имеем функцию /0 (V) на множестве сетевых операторов. Тогда считаем, что функция /0 (V) достигает локального минимума на множестве сетевых операторов, если

такой, что V* є п(ї, 1), /0 (V )> /0 (V).

Пусть w оператор вариации или просто вариация сетевого оператора

О V Ф V . (20)

Элементарной вариацией сетевого оператора называем вариацию, удовлетворяющую условию

5^ о V, V)= 1. (21)

Условие локального минимума функции /3 (V) имеет вид

/0^ О V)-/0 (V )> 0. (22)

Элементарные вариации сетевого оператора могут отличаться друг от друга. Для сетевого оператора достаточно использовать следующие вариации: 0 - замена унарной операции, 1 -замена бинарной операции, 2 - добавление унарной операции, 3 - удаление унарной операции.

При выполнении элементарной вариации 3 следует учитывать, что столбцы и строки сетевого

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

w = И И 2 Из И 4 ]^, (23)

где И1 - номер вариации, и2 - номер строки, Из - номер унарной или бинарной операции.

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

Wi = ^ Z,1,^, w1,1), (24)

Щ о Т = w i,/ о...о w ^ о Т (25)

Алгоритм поиска оптимального математического выражения включает следующие

операции:

1. Задаем сетевой оператор Т , который, по мнению исследователя, наиболее близко описывает искомое решение. Данный сетевой оператор называем базисным.

2. Генерируем наборы вариаций Wi = ^^1,...,w7,1), i = 1,Н , и вычисляем значения

целевой функции для каждого нового сетевого оператора /^ (Щi о Т0 ).

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

4. Если /0 (Щ о Т0 )< /0 (Т0 ), то Т0 = Щ о Т0 .

Вероятность Ра нахождения оптимального сетевого оператора с помощью

рассмотренного алгоритма зависит от количества Н сгенерированных наборов вариаций:

Пт РА = 1. (26)

Н ——^

3. Пример

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

*1 = — И — И — ™3, х2 = — "И^ + И 2 ,

•*3 = И — И2 — И3 ,

*4 = И2,

•*5 = И2 ,

•*6 = И3,

*7 = И3 ,

где И = g(^1,*1,*2), и2 = g(#2,*1,*з), и3 = g(^з,*1,*3), g^,а,Ь) - неизвестная функция трех аргументов, qi - неизвестные значения постоянных параметров, i = 1,2,3 .

Для системы дифференциальных уравнений заданы начальные условия:

*1 (0)= 0,8818, *2 (0)= 0,1182, * (0) = 0, i = 3,7.

Уравнения баланса (5) имеют вид

3*1 + 3*3 + 3*5 + 3*6 + 3*7 = 3*1 (0),

6 *1 + 6 *3 + 6 *5 + 6 *6 + 6 *7 = 6 *1 (0),

2 *1 + 2 *3 + 2 *4 + *5 + 2 *6 + 2 *7 = 2 *1 (0),

*2 + *3 + *7 = *2 (0).

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

*1*2 У *1 (0)*2(0) .

Экспериментальные значения интенсивностей и моменты времени, в которые она определена, приведены в табл.1.

Для решения задачи использовали метод сетевого оператора. Для поиска решения использовали вариационный генетический алгоритм с множественным базисом. Алгоритм имел следующие значения параметров: количество возможных решений в начальной популяции -512; число поколений - 128; число возможных скрещиваемых пар в поколении - 256; число вариаций в одном решении - 4; вероятность мутации 0.7; число базисов - 5; число элитарных решений 12; число поколений между сменой базисов - 16.

Таблица 1 Экспериментальные данные

t к У к t к У к

0 9.9452с-1 7.4139с-1 4.7645с-1

5.4882e-3 9.9452с-1 7.6284с-1 4.6634с-1

2.1722e-2 9.7851с-1 7.8800с-1 4.5475с-1

4.7246e-2 9.5385с-1 8.1096с-1 4.4443с-1

5.8350e-2 9.4332с-1 8.2000с-1 4.4043с-1

8.0944e-2 9.2225с-1 8.4518с-1 4.2947с-1

1.0172с-1 9.0328с-1 8.5800с-1 4.2401с-1

1.1463с-1 8.9169с-1 8.8418с-1 4.1305с-1

1.3783с- 1 8.7125с-1 9.0740с-1 4.0357с-1

1.4803с- 1 8.6240с-1 9.1841с-1 3.9915с-1

1.7553с-1 8.3901с-1 9.5718с-1 3.8397с-1

1.9994с-1 8.1878с-1 9.5223с-1 3.8588с-1

2.1132с- 1 8.0951с-1 9.7710с-1 3.7640с-1

2.3904с-1 7.8738с-1 9.9004с-1 3.7156с-1

2.7197с-1 7.6188с-1 1.0159с+0 3.6208с-1

2.8449e-1 7.5240e-1 1.0394e+0 3.5366e-1

3.1348e-1 7.3090e-1 1.0556e+0 3.4798e-1

3.2858e-1 7.1995e-1 1.0889e+0 3.3660e-1

3.5197e-1 7.0330e-1 1.1265e+0 3.2417e-1

3.6797e-1 6.9213e-1 1.1415e+0 3.1933e-1

3.8207e-1 6.8244e-1 1.1812e+0 3.0690e-1

4.0297e-1 6.6833e-1 1.2062e+0 2.9932e-1

4.1854e-1 6.5801e-1 1.2218e+0 2.9469e-1

4.2592e-1 6.5317e-1 1.2612e+0 2.8331e-1

4.4218e-1 6.4264e-1 1.2830e+0 2.7721e-1

4.5537e-1 6.3421e-1 1.3092e+0 2.7005e-1

4.6471e-1 6.2831e-1 1.3408e+0 2.6163e-1

4.8162e-1 6.1778e-1 1.3670e+0 2.5489e-1

5.0299e-1 6.0472e-1 1.3819e+0 2.5110e-1

5.1737e-1 5.9609e-1 1.3971e+0 2.4732e-1

5.4788e-1 5.7817e-1 1.4073e+0 2.4480e-1

5.6293e-1 5.6954e-1 1.4229e+0 2.4101e-1

5.8689e-1 5.5605e-1 1.4512e+0 2.3427e-1

6.0408e-1 5.4658e-1 1.4721e+0 2.2943e-1

6.1220e-1 5.4216e-1 1.5057e+0 2.2185e-1

6.2628e-1 5.3458e-1 1.5356e+0 2.1532e-1

6.4578e-1 5.2425e-1 1.5765e+0 2.0669e-1

6.5262e-1 5.2068e-1 1.6224e+0 1.9742e-1

6.7141e-1 5.1099e-1 1.6671e+0 1.8879e-1

6.9223e-1 5.0046e-1 1.6920e+0 1.8416e-1

7.0195e-1 4.9562e-1 1.7631e+0 1.7152e-1

7.1866e-1 4.8740e-1 1.7904e+0 1.6689e-1

1.8185e+0 1.6227e-1 3.3185e+0 3.621e-2

2.0544e+0 1.2817e-1 3.3418e+0 3.537e-2

2.1171e+0 1.2038e-1 3.3718e+0 3.4328e-2

2.1401e+0 1.1764e-1 3.3964e+0 3.3494e-2

2.1821e+0 1.1280e-1 3.4216e+0 3.2661e-2

2.2338e+0 1.0712e-1 3.4880e+0 3.0562e-2

2.2518e+0 1.0523e-1 3.4877e+0 3.0571e-2

2.2986e+0 1.0039e-1 3.5514e+0 2.8683e-2

2.3569e+0 9.4710e-2 3.6515e+0 2.5951e-2

2.3748e+0 9.3033e-2 3.7627e+0 2.3220e-2

2.3953e+0 9.1144e-2 3.8088e+0 2.2175e-2

2.4281e+0 8.8202e-2 3.8179e+0 2.1974e-2

2.4256e+0 8.8422e-2 3.8566e+0 2.1140e-2

2.4619е+0 8.5269е-2 3.8561е+0 2.1150е-2

2.4718е+0 8.4435е-2 3.9068е+0 2.0105е-2

2.4944е+0 8.2547е-2 3.9063е+0 2.0114е-2

2.5175е+0 8.0659е-2 3.9487е+0 1.9280е-2

2.5520е+0 7.7927е-2 3.9481е+0 1.9290е-2

2.5905е+0 7.4985е-2 3.9477е+0 1.9300е-2

2.6572е+0 7.0146е-2 3.9044е+0 2.0153е-2

2.6845е+0 6.8258е-2 3.9467е+0 1.9319е-2

2.7125е+0 6.6370е-2 4.0375е+0 1.7641е-2

2.7414е+0 6.4480е-2 4.0370е+0 1.7650е-2

2.8019е+0 6.0696е-2 4.0365е+0 1.7660е-2

2.8335е+0 5.8808е-2 4.0360е+0 1.7669е-2

2.8333е+0 5.8817е-2 4.0354е+0 1.7678е-2

2.8960е+0 5.5242е-2 4.0349е+0 1.7688е-2

2.8959е+0 5.5251е-2 4.0343е+0 1.7697е-2

2.9875е+0 5.0411е-2 4.1472е+0 1.5809е-2

3.0477е+0 4.7469е-2 4.2896е+0 1.3710е-2

3.0475е+0 4.7478е-2 4.3523е+0 1.2877е-2

3.0473е+0 4.7488е-2 4.4369е+0 1.1832е-2

3.0650е+0 4.6654е-2 4.4361е+0 1.1842е-2

3.0876е+0 4.5610е-2 4.4353е+0 1.1851е-2

3.1299е+0 4.3722е-2 4.4345е+0 1.1861е-2

3.1741е+0 4.1833е-2 4.6079е+0 9.9725е-3

3.1738е+0 4.1843е-2 4.5066е+0 1.1036е-2

3.1736е+0 4.1852е-2 4.5057е+0 1.1046е-2

3.2198е+0 3.9964е-2 4.5049е+0 1.1055е-2

3.2195е+0 3.9974е-2 4.5040е+0 1.1064е-2

3.2624е+0 3.8296е-2 4.5032е+0 1.1074е-2

3.2901е+0 3.725е-2 4.6023е+0 1.0029е-2

3.3185е+0 3.621е-2 4.5805е+0 1.0249е-2

3.3418е+0 3.537е-2 4.6880е+0 9.2046е-3

3.3718е+0 3.4328е-2 4.6867е+0 9.2147е-3

3.3964е+0 3.3494е-2 4.5984е+0 1.0067е-2

В результате был получен сетевой оператор

\

т

0 0 0 23 5 0 0 0

0 0 0 23 14 0 0 0

0 0 0 0 1 0 0 0

0 0 0 1 23 0 0 0

0 0 0 0 1 23 0

0

0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

3

1 23 3 1 15 00

который, соответствует следующему уравнению

g(д,, а, Ь) = ^А - А3 -(а - А3 )3 - 2А + А3 - д1 (а - а3 )ь - Ь3 )+ (а - а3 )3 (ь - Ь3 )3

где

А = д, (ь5 - а2Ь4 - Ь6 + а2Ь6 )- д (аЬ2 - а3Ь2 - аЬ4 + а3Ь4)3 .

а

-, (27)

(28)

Для параметров д, получены следующие значения: ^ = 26,09753, д2 = 0,01590, д3 = 314,67213. График экспериментальных и вычисленных значений интенсивностей, полученных с использованием полученных функций g (д,, а, Ь) приведен на рис. 1.

Рис. 1 Экспериментальные и расчетные значения интенсивностей

Заключение

Из сравнения графиков, представленных на рис. 1 можно сделать вывод о высокой точности полученной математической модели химической реакции. В результате расчетов можно выдвинуть предположение об уточнении закона взаимодействия веществ Гульденберга-Ваага, концентрации веществ при их взаимодействии пропорциональны произведению полиномов от концентраций каждого вещества [7].

Работа выполнена по гранту РФФИ № 12-07-00324-а. «Структурная и параметрическая идентификация кинетических моделей реакций нейтрального металлокомплексного катализа».

3

а

Список литературы

1. Koza J.R. Genetic Programming: On the Programming of Computers by Means of Natural Selection. Cambridge, Massachusetts, London, MA: MIT Press, 1992. 819 p.

2. O’Neill M. Automatic Programming in an Arbitrary Labguage: Evolving Programs in Grammatical Evolution: PhD thesis. University of Limerick, 2001. 163 p.

3. Zelinka I. Analytic programming by Means of Soma Algorithm // Proc. 8th International Conference on Soft Computing Mendel’02. Brno, Czech Republic, 2002. P. 93-101.

4. Дивеев А.И. Метод сетевого оператора. М.: ВЦ РАН, 2010. 178 с.

5. Дивеев А.И., Софронова Е.А. Метод сетевого оператора и его применение в задачах управления. М.: РУДН, 2012. 182 с.

6. Diveev A.I., Sofronova E.A. The Network Operator Method for Search of the Most Suitable Mathematical Equation // In the book “Bio-Inspired Computational Algorithms and Their Applications” / Edited by Shangce Gao. Croatia: InTech, 2012. P. 19-42.

7. Слинько М.Г. Основы и принципы математического моделирования каталитических процессов. Новосибирск: Институт катализа им. Г.К. Борескова СО РАН, 2004. 488 с.

8. Овчинников М.Ю., Юнусов А.А., Хурсан С.Л., Губайдуллин И.М. Моделирование генерации синглетного кислорода при разложении диметилдиоксирана с использованием технологии OpenMP // Параллельные вычислительные технологии (ПаВТ’2012) : труды международной научной конференции (Новосибирск, 26-30 марта 2012 г.). Челябинск: Издательский центр ЮУрГУ, 2012. С. 735-736.

SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU

SCIENCE and EDUCATION

EL № FS77 - 48211. №0421200025. ISSN 1994-0408

electronic scientific and technical journal

Identification of the structure and parameters of a chemical reaction model by a network operator method

# 12, December 2012

DOI: 10.7463/0113.0521994

Diveev A.I., Semikov M.V., Sofronova E.A.

Dorodnicyn Computing Centre of RAS Peoples’ Friendship University of Russia

aidiveev@mail.ru semikovmv@yandex.ru sofronova ea@mail.ru

The authors consider the identification problem for a chemical reaction model. The mathematical model in the form of a system of ordinary differential equations is developed on the basis of the given experimental data on the variation in time index of the concentration of substances that interact in the chemical reaction. To construct this model which clarifies the law of mass action a network operator is used. As a result of application of the network operator method, clarification of the law as a product of polynomials of the concentration of substances that interact in the chemical singlet oxygen reaction is obtained.

Publications with keywords: identification, the method of network operator Publications with words: identification, the method of network operator

References

1. Koza J.R. Genetic Programming: On the Programming of Computers by Means of Natural Selection. Cambridge, Massachusetts, London, MA, MIT Press, 1992. 819 p.

2. O’Neill M. Automatic Programming in an Arbitrary Labguage: Evolving Programs in

Grammatical Evolution: PhD thesis. University of Limerick, 2001. 163 p.

3. Zelinka I. Analytic programming by Means of Soma Algorithm. Proc. 8th International

Conference on Soft ComputingMendel’02. Brno, Czech Republic, 2002. P. 93-101.

4. Diveev A.I. Metod setevogo operator a [Method of network operator]. Moscow,

Computing Center of RAS Publ., 2010. 178 p.

5. Diveev A.I., Sofronova E.A. Metod setevogo operatora i ego primenenie v zadachakh upravleniia [Method of network operator and its application in problems of control]. Moscow, Peoples’ Friendship University of Russia Publ., 2012. 182 p.

6. Diveev A.I., Sofronova E.A. The Network Operator Method for Search of the Most Suitable Mathematical Equation. In book: Shangce Gao, ed. Bio-Inspired Computational Algorithms and Their Applications. Croatia, InTech, 2012, pp. 19-42.

7. Slin'ko M.G. Osnovy i printsipy matematicheskogo modelirovaniia kataliticheskikh protsessov [The basics and principles of mathematical modeling of catalytic processes]. Novosibirsk, Boreskov Institute of Catalysis SB RAS Publ., 2004. 488 p.

8. Ovchinnikov M.Iu., Iunusov A.A., Khursan S.L., Gubaidullin I.M. Modelirovanie generatsii singletnogo kisloroda pri razlozhenii dimetildioksirana s ispol'zovaniem tekhnologii OpenMP [Modeling of generation of singlet oxygen in the decomposition of dimethyldioxirane with the use of OpenMP technology]. Parallel'nye vychislitel'nye tekhnologii (PaVT’2012) : trudy mezhdunarodnoi nauchnoi konferentsii [Proc. of the International scientific conference "Parallel computational technologies (PCT’2012)"], Novosibirsk, 26-30 March, 2012. Cheliabinsk, Publishing center of SUSU, 2012, pp. 735-736.