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

t, шагов ВКА

Рис. 5. Зависимость изменения массы образца от времени

Примечание: 1, 2, 3 - номера графиков, соответствующие номерам гипотез о структуре модели.

На следующем этапе оценим влияние химических процессов на параметры системы. Для этого примем гипотезу № 2 о том, что кислород может диффундировать в твердом теле и окислять атомы металла. В данном случае в структуру добавляются модули, обозначенные на рисунке 4 пунктирной линией. В результате моделирования такой системы с помощью разработанного программного продукта получаем определенную зависимость (рис. 5, график 2). Из результатов моделирования видно, что в данных конкретных условиях процессы окисления рассматриваемого металла играют существенную роль.

Для проверки гипотезы № 3 о влиянии процессов вторичного окисления металла с образованием второго оксида в структурную модель системы включим процесс вторичного окисления и второй оксид (затонированные блоки на рис. 4). Проведя моделирование указанной системы, получаем результат, приведенный на рисунке 5

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

Таким образом, в результате проведенного анализа разработчик может остановиться на структуре модели, предложенной в гипотезе № 2.

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

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

Литература

1. Агафонов А.Н. [и др.]. Разработка физических принципов и алгоритмов компьютерного моделирования базовых процессов формирования микроструктур методами вероятностного клеточного автомата // Вестн. СамГТУ. Сер.: Физ.-мат. науки, 2007.№ 1. С. 99-107.

2. Коныгин С.Б., Лесухин С.П. Стохастическая модель окисления металлов в газовой среде // Изв. Самарского науч. центра РАН. 2010. Т. 12. № 1 (2). С. 377-380.

3. Коныгин С.Б. Применение метода вероятностного клеточного автомата к моделированию гетерогенных систем // Строительство нефтяных и газовых скважин на суше и на море. 2011. № 4. С. 30-31.

УДК 519.85:519.676

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

Е.С. Семёнкин, д.т.н.

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

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

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

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

Задача поиска структурной модели

Рассмотрим основные этапы структурного анализа. После установления формулы соединения и регистрации дифрактограммы происходят определение (индицирование) и уточнение параметров решетки, а затем пространственной группы симметрии. Далее выполняется полнопрофильная декомпозиция интенсивности перекрытых рефлексов дифрактограммы с использованием известных методов ЛеБэйля или Поули. Следующий этап - поиск модели неизвестной кристаллической структуры. Третьим этапом структурного анализа является уточнение найденной модели.

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

Эффективным средством решения многопараметрических оптимизационных задач в различных областях науки и техники зарекомендовали себя генетические алгоритмы (ГА). Известны работы,

например, Харриса и ряда других авторов по применению классического ГА для определения кристаллических структур в прямом пространстве. Так, для решения задачи определения кристаллической структуры поликристаллического вещества по данным порошковой дифракции предложено использовать двухуровневый гибридный ГА, комбинированный с эффективными математическими методами оптимизации и подходящими дифракционными, кристаллохимическими и физическими критериями оценки структурных моделей [2]. Данный алгоритм при правильной настройке может обеспечивать эволюцию моделей с устойчивой сходимостью в глобальный оптимум, отвечающий истинной структуре, при существенном увеличении размерности структурных задач [3]. Однако определение правильной настройки является очень серьезной проблемой в прикладных задачах.

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

Задача 1-го уровня гибридного алгоритма -найти удовлетворительные (в смысле величины ^-фактора) начальные приближения параметров модели. В процессе выполнения 2-го уровня происходит поиск пути спуска на гиперповерхности ^-фактора из точек, полученных с 1-го уровня ГА. Эффективность сходимости к решению обеспечивает метод минимизации производной разности (МПР) расчетной и экспериментальной ди-фрактограмм. Однако для достижения глобального оптимума необходимо, чтобы начальные приближения, поступившие на 2-й уровень, были достаточно близки к решению. Таким образом, успех в достижении решения обеспечивает ГА, применяемый на 1-м уровне.

Коэволюционный асимптотический ГА

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

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

Существует несколько подходов к решению этой проблемы. Например, прямое кодирование параметров ГА нижнего уровня, решающего задачу оптимизации, в хромосому мета-ГА более высокого уровня, решающего задачу настройки параметров. Другой подход состоит в гибридизации ГА со специальным алгоритмом локального поиска в бинарном пространстве. В некоторых приложениях это приводит к низкой чувствительности эффективности ГА к выбору параметров, что позволяет избежать их настройки. Следующий подход - адаптация стратегии ГА за счет конкурирующих подпопуляций. Каждая подпопуляция имеет собственную стратегию (комбинацию параметров алгоритмов). Перераспределение ресурсов обеспечивает доминирование подпопуляции с наиболее подходящей для решаемой задачи стратегией поиска. Этот подход может рассматриваться как один из примеров коэволюционных ГА, отличающихся идеями организации взаимодействия стандартных ГА. В частности, ранее был предложен подход, использующий как конкуренцию, так и сотрудничество индивидуальных ГА, каждый из которых имеет свою комбинацию параметров [4]. Сотрудничество индивидуальных ГА достигается благодаря миграции лучших решений во все под-популяции, что обеспечивает рост эффективности коэволюционного алгоритма за счет положительного эффекта взаимодействия подпопуляций. Проблемой применения коэволюционных алгоритмов является неопределенность выбора кандидатов для включения в коэволюцию. Число таких кандидатов не может превышать полутора-двух десятков алгоритмов, в то время как выбирать их надо из сотен.

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

Еще меньше настраиваемых параметров содержит асимптотический вероятностный ГА [5], в котором применяется адаптивная мутация, не требующая предварительной настройки. Дополнительным преимуществом асимптотического алгоритма является способ выполнения генетических операторов, который сводится к изменению распределения вероятностей компонент, что значительно сокращает время на проработку алгоритма. Фактически асимптотический вероятностный ГА имеет всего несколько комбинаций настроек параметров.

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

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

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

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

Результаты решения задач поиска структуры с помощью КАГА

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

В первом случае было необходимо по реальной дифрактограмме определить структуру соединения Pd(NH3)2(NO2)2. Единственный на ячейку атом Pd, как и при решении в [2], зафиксирован в

центре инверсии, в общих позициях размещались 4 атома N и О (12 искомых координат), атомы водорода исключены из поиска. Независимая часть ячейки для каждого атома определяет следующую область поиска: хе [0; 0,5], у, zе [0; 1].

Результаты тестирования. Сравнение КАГА с индивидуальными алгоритмами

Использовалась точность кодирования 5 бит на переменную. Выраженная преимущественная ориентация не моделировалась для усложнения задачи поиска, что привело к величине целевой функции (^-фактора) в точке минимума, равной 19,08.

Был осуществлен ряд пусков алгоритма, в которых варьировались размеры популяций и количество поколений. Во всех полученное решение было достаточно близким к известному [2], так что на втором уровне оно сводилось алгоритмом уточнения структуры к точному. Типичный пример сходимости процесса коэволюционного поиска продемонстрирован на рисунке 1. Каждая из линий соответствует какому-либо алгоритму, входящему в коэволюцию. Отдельная точка показывает значение пригодности лучшего индивида популяции в конкретный момент.

Вторая задача поиска кристаллических структур химических соединений заключается в отыскании структуры соединения К^пОф Количество переменных равно 27 - по 3 координаты на девять атомов. Точность кодирования - 5 бит

на переменную. Поиск проводился по всей независимой части ячейки: х, уе [0; 1], zе [0; 0,5].

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

R-фактор

^^^^^ пропорциональная ^ • «ранговая

»малый турнир •••••• средний турнир

ее* большой турнир

Рис. 1. График сходимости КАГА Примечание: ось абсцисс - номер поколения коэволюции, ось ординат - пригодность лучшего индивида для каждого из 5 коэволюционных алгоритмов

R-фактор

^^^^^ пропорциональная • "ранговая

«малый турнир •••••• средний турнир

««« большой турнир

Рис. 2. Работа алгоритма коэволюции для второй задачи

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

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

Функция Надежность лучшего индивидуального ГА Скорость ГА Надеж-ность КАГА Скорость КАГА

AdditivePotential 1 2160 1 1062

DeJong2 1 3040 1 2226

GaussianQuartic 1 2990 1 1861

Griewank 0,98 1428 0,6 3813

Ackley 1 2745 1 2014

Ackley2 1 2775 1 1850

Griewank2 1 2130 1 2146

Himmelblau 0,72 6764 0,4 2237

HyperEllipsoid 1 2595 1 2180

HyperEllipsoid2 1 3875 1 2073

Rastrigin 1 3490 1 3881

Katkovnik 1 3310 1 2278

Multiextremal3 1 3415 1 1747

Multiextremal4 1 4330 1 1838

MultiplicativePotential 1 1845 1 1335

Rastrigin2 1 2415 1 1964

Rastrigin3 1 3425 0,98 4259

Rastrigin WithTurning 1 2970 1 1948

Rosenbrock 0,36 6666 0,44 6502

SphereModel 1 2880 1 1861

Среднее значение 0,953 3262 0,921 2454

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

Литература

1. David W.I.F., Shankland K. Structure determination from powder diffraction data // Acta Cryst. 2008. A64, pp. 52-64.

2. Yakimov Y.I., Semenkin E.S., Yakimov I.S. Two-level genetic algorithm for a fullprofile fitting of X-ray powder patterns //

Z. Kristallogr. Suppl. 2009. Vol. 30, pp. 21-26.

3. Якимов Я.И., Семёнкин Е.С. Автоматизированная система рентгеноструктурного анализа поликристаллов // Программные продукты и системы. 2009. № 4 (88). С. 82-84.

4. Sergienko R., Semenkin E. Competitive Cooperation for Strategy Adaptation in Coevolutionary Genetic Algorithm for Constrained Optimization // IEEE World Congress on Computational Intelligence (WCCI'2010). Barcelona, Spain, 2010, pp. 1626-1631.

5. Галушин П.В., Семёнкин Е.С. Асимптотический вероятностный генетический алгоритм // Вестн. СибГАУ им. акад. М.Ф. Решетнева. 2009. Вып. 4 (25). С. 37-42.

УДК 519.254

ПРОГРАММНЫЙ КОМПЛЕКС ОБНАРУЖЕНИЯ АНОМАЛИЙ ФОРМЫ РЕЛЬСОВЫХ ПУТЕЙ

Г.М. Соломаха, д.ф.-м.н.

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

Ключевые слова: программный комплекс, критерий, аномалия, рельсовая нить, моделирование.

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

Типовыми АФРН, обнаруживаемыми в первую очередь, являются аномалии следующих видов:

- односторонняя просадка (локальное вертикальное искривление) одной рельсовой нити относительно другой;

- одновременная и одинаковая по форме просадка (локальное вертикальное искривление) обеих рельсовых нитей;

- горизонтальное (в плане) одновременное и одинаковое искривление обеих рельсовых нитей.

Возможны комбинации этих простейших АФРН.

В результате предварительных исследований в созданной системе регистрации центральным элементом определен миниатюрный прибор ADIS 16364 фирмы Analog Device, включающий в себя по одному трехосевому датчику линейных ускорений (акселерометр) и угловых скоростей (гироскоп). Поэтому в данной работе учитываются только зависимости линейных ускорений движения вдоль каждой из трех осей декартовой системы координат и угловых скоростей поворота вокруг этих же осей.

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

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

Кроме того, в комплексе предусмотрен режим полунатурного моделирования. В этом режиме