Моделирование тенденций поведения элементов систем уникальных объектов Текст научной статьи по специальности «Физика»

Вычислительные технологии

Том 10, № 5, 2005

МОДЕЛИРОВАНИЕ ТЕНДЕНЦИЙ ПОВЕДЕНИЯ ЭЛЕМЕНТОВ СИСТЕМ УНИКАЛЬНЫХ

ОБЪЕКТОВ

А. Б. Логов, Р.Ю. Замараев, А. А. Логов Институт угля и углехимии СО РАН, Кемерово, Россия e-mail: {a.logov,r.zamaraev}@kemsc.ru, logoff@rambler.ru

We present a further development of the entropy analysis method for the complex multifactor systems of unique objects. The concept of classes of behavior of an elements of the system in a phase space is introduced. The principles of the construction of behavior models along with trends revealing are presented. Some features of behavior modeling applied to the analysis of Mendeleev periodic system of chemical elements are demonstrated

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

Здесь обсуждается основная и простейшая модель — тенденция в поведении элемента линейно зависит от начальных координат Ui.o и U2.o на плоскости состояний и обусловлена величиной отклонения от эллиптической границы между видами состояния [2].

Поведение на границе в параметрической форме (U1(t), U2(t)) в точности может быть описано следующей системой:

AB (1 + r2) ^ulw + (A2 - B2) r ад) = (a2 + B2r2) ад)

I U2W,

(1)

АВ (1 + г2) ^Ц2^ - (а2 - в2) г ад) = - (А2г2 + в2) ад),

где А(т, г) и В(т, г) — полуоси граничного эллипса, изображающего гармоническое решение в системе мощностью т с коэффициентом корреляции г = М[Ц • Ц2].

В произвольной точке окрестности эллипса можно оценить параметр отклонения

= Ц2'° = , (2)

Ui .ell (to) U .ell (to) с учетом которого модель поведения представляется в виде

AB(1 + r2) Udtt) + [(A2 - B2)r - ABZ(1 + r2)] ад) = (A2 + в2r) ад), AB(1 + r2) - [(A2 - B2)r + ABZ(1 + r2)] ад) = -(A2r + B2) Us(t).

© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2005.

МОДЕЛИРОВАНИЕ ПОВЕДЕНИЯ ЭЛЕМЕНТОВ СИСТЕМ

55

При тестировании метода на примере анализа периодической системы химических элементов [2] нетривиальный результат выделения видов состояния был получен для комбинации фундаментальных характеристик: Е(2) — величина энергии связи, Е(3) — длина связи в двухатомной молекуле (полный перечень из 56 рассмотренных характеристик приведен в [2]).

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

Р(,72-3)- ^1П(Д(3))] К(*' 2) (4)

Р(г/2'3)= а[1п(Е(2))] ■ Ж*)' (4)

где

т _

К(г/з)- Е(г/зШ Я(1/з) _ ; ¿=1 1 т

а2[1п(Д(з))]--- Е1п2(К(г/з)).

т — 1 ¿=1

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

1

и2(г/2;3)-1п(Р(г/2; 3))^мтал, (5)

где

т

а2[1п(Р(2; 3))] - -г$>2(Р(г/2; 3)).

т — 1

г=1

Как показано в [2], производная от модели потенциала связи имеет вид

М[Е-1(2)] — Е-1(2) М[Е-1(3) — Е-1 (3)]

с(г/2'3)--^[Е^---¿Ечвд-, (6)

56 А. Б. Логов, Р. Ю. Замараев, А. А. Логов

и, соответственно, получаем отображение комбинации на ось ординат

= G§§, (7)

где

т

а2[С(2;3)] = -С2(г/2;3).

т - 1

г=1

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

(Ш

0.2119 + Цу (') = 0.6637 Ц2 ('),

/ (8) 0.1799—^ - Ц(') = -0.4927 Ц(')

с решением

Цу(0) = 2.80 • в°'420 • ео8(0 - 0.093), Ц(0) = 2.62 • в°'420 • ео8(0 + 0.255).

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

Здесь достаточно характерно сопоставление свойств углерода С(6) и фтора Е(9), чьи изображающие точки расположены достаточно близко, но разделены прямой линейной регрессии. Начальная тенденция углерода — сближение по комбинированной характеристике с остальной системой; фтор обнаруживает тенденцию к выделению из совокупности. Такое же свойство проявляется у азота

Портрет группы инертных газов не только радикально отличается от остальной системы, но и расположен почти ортогонально линии связи. Смысл соответствующей независимости иллюстрируется моделью поведения криптона Кг(36).

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

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

[1] Логов А.Б., Кочетков В.Н., Рожков А.А. Энтропийный подход к моделированию процесса реструктуризации угольной отрасли. М.: Недра, 2001. 324 с.

[2] Логов А.Б., Замараев Р.Ю., Логов А.А. Анализ функционального состояния промышленных объектов в фазовом пространстве. Кемерово: Институт угля и углехимии СО РАН, 2004. 168 с.

Поступила в редакцию 11 февраля 2005 г., в переработанном виде —12 мая 2005 г.