Стохатические модели
Теория стохатических моделей
Теория перколяции
Модели дендритов
Теория клеточных автоматов
Модель изинга
Генетические алгоритмы
Случайная функция –случайная величина, зависящая не только от случайного события ω, но и от какого-либо параметра. Если этот параметр – время, то случайная функция называется случайным процессом и обозначается ξ (t,ω).
Величина ξ может быть как скалярной (скалярный случайный процесс), так и векторной (векторный или многомерный случайный процесс). Она может принимать как дискретные значения (процесс с дискретными состояниями), так и пробегать непрерывный ряд значений
Методы Монте-Карло
Метод Монте-Карло - общее название группы численных методов, основанных на получении большого числа реализаций стохастического (случайного) процесса, который формируется таким образом, чтобы его вероятностные характеристики совпадали с аналогичными величинами решаемой задачи.
Алгоритм Буффона
Суть метода была в бросании иглы длиной L на плоскость, расчерченную параллельными прямыми, расположенными на расстоянии r друг от друга.
Интегрирование методом Монте Карло
Ограничим функцию прямоугольником (n-мерным параллелепипедом в случае многих измерений), площадь которого S0 можно легко вычислить;
«набросаем» в этот прямоугольник (параллелепипед) некоторое количество точек (N штук), координаты которых будем выбирать случайным образом;
определим число точек (K штук), которые попадут под график функции;
площадь области, ограниченной функцией и осями координат, S даётся следующим выражением S = S0 K / N
Прохождение нейтронов сквозь пластинку
Взаимодействие нейтронов с веществом характеризуется в рассматриваемом случае двумя постоянными Σc и Σs , которые называются сечением поглощения и сечением рассеяния (индексы с и s—это первые буквы английских слов capture -- захват и scattering -- рассеяние). Сумма этих сечений Σ = Σc + Σs называется полным сечением.
Задача о протекании
Пусть есть прямоугольная решетка, состоящая из N*N свободных и занятых узлов.
Пусть вероятность того, что узел занят – p, и, соответственно, вероятность того, что узел свободен (1-p).
Будем называть кластером соседние занятые узлы и стягивающим кластером (перколяционным кластером) такой кластер, который начинается на одной границе и заканчивается на противоположной границе решетки.
Образование стягивающего кластера называется перколяцией.
Метод Монте-Карло
- Задать значение вероятности занятости узла p.
- Разыграть состояние каждого узла при заданном значении p, тем самым построив конкретную реализацию решетки.
- Определить, есть ли в этой реализации хотя бы один стягивающий кластер. Если есть, то увеличить счетчик числа стягивающих кластеров Mc на единицу.
- Повторить пункты 2-3 M раз.
- Получить оценку вероятности образования стягивающего кластера
- Повторить 1-5 для ряда значений вероятности p в интервале от 0 до 1.
Источники задач
Физика (исторически первый источник) – задачи просачивания жидкостей в пористой среде; некоторые задачи, связанные с проводимостью.
Математика – случайные графы.
Химия.
Другие области - для описания возникновения связных структур в случайных средах (кластеров), состоящих из отдельных элементов.
Порог перколяции
Порог pᶜ перколяции – доля занятых узлов, при которой возникает перколяционный кластер.
Бесконечная квадратная решетка
pᶜ =0,5 для задачи связей;
pᶜ ≈ 0,59275 для задачи узлов.
Методы перколяции для сетевых алгоритмов
Лавинная маршрутизация на основе перколяции для сенсорных сетей
Перколяционная маршрутизация для оптической кластерной сети
Перколяционный поиск в сетях со степенным распределением связности
мобильных сетях, уровень изменений топологии в которых настолько высок, что лавинная доставка данных является единственным способом маршрутизации.
Применение теории перколяции для оценки общей устойчивости и уязвимостей транспортных сетей
Явления макромира - математические модели (бесконечные и непрерывные).
Основная особенность: возможность одновременного (параллельного) изменения состояния всей системы, в то время как каждый участок системы взаимодействует только со своими непосредственными соседями. Это свойство позволяет при моделировании связать события, происходящие на микроуровне, с изменениями макроуровневого моделируемого объекта.
Автоматы
Если состояние системы в произвольный момент времени характеризуется лишь ее предыдущим состоянием и набором правил, регламентирующих ее переход, то она называется автоматом.
Клеточные автоматы широко применяются для моделирования систем, в которых важную роль играет пространственное взаимодействие между элементами.
В физике, например, клеточные автоматы применяются для анализа явлений переноса (теплопроводности, диффузии и вязкости) и моделирования твердого тела.
Классической системой с мелкозернистым параллелизмом является клеточный автомат.
Классификация по типам поведения
Класс 1: Результатом эволюции начальных условий является быстрый переход к гомогенной стабильности. Любые негомогенные конструкции быстро исчезают.
Класс 2: Результатом эволюции начальных условий является быстрый переход в неизменяемое негомогенное состояние либо возникновение циклической последовательности. Большинство структур начальных условий быстро исчезает, но некоторые остаются. Локальные изменения в начальных условиях оказывают локальный характер на дальнейший ход эволюции системы.
Класс 3: Результатом эволюции почти всех начальных условий являются псевдо-случайные, хаотические последовательности. Любые стабильные структуры, которые возникают почти сразу же уничтожаются окружающим их шумом. Локальные изменения в начальных условиях оказывают неопределяемое влияние на ход эволюции системы.
Класс 4: Результатом эволюции являются структуры, которые взаимодействуют сложным образом с формированием локальных, устойчивых структур. В результате эволюции могут получаться некоторые последовательности Класса 2, описанного выше. Локальные изменения в начальных условиях оказывают неопределяемое влияние на ход эволюции системы. Некоторые клеточные автоматы этого класса обладают свойством универсальности по Тьюрингу.
Вероятностные правила
Такие клеточные автоматы называются стохастическими. В таких системах задается вероятность, что на следующем шаге клетка сменит свой цвет на другой.
Криптография
Одно из правил - возможный блочный шифр.
Двумерные клеточные автоматы используются для генерации случайных чисел.
Клеточные автоматы предложены для использования в криптосистемах с открытым ключом.
В этом случае односторонняя функция является результатом эволюции конечного клеточного автомата, первоначальное состояние которого сложно найти.
По заданному правилу легко найти результат эволюции клеточного автомата, но очень сложно вычислить его предыдущие состояния.
Тьюрмит - это машина Тьюринга, которая имеет ориентацию в пространстве, текущее состояние и «ленту», состоящую из бесконечного двухмерного массива ячеек.
В настоящее время для сборки устройств из отдельных молекул используют атомно-силовой микроскоп, но это ручная и очень непроизводительная работа. Если же удастся создать молекулярный аппарат типа тьюрмита, он сможет по командам своего "мозга" быстро строить нужные молекулярные структуры.
Законы живого мира
Приспособляемость
Естественный отбор
Наследование
среди генетически различающихся особей одной и той же популяции выживают и оставляют потомство только наиболее приспособленные к окружающей среде.
Применение принципов эволюции к формированию и представлению знаний
Эволюционируют не сами живые организмы, а их хромосомы (концентрированно содержат в себе все свойства организма). Эволюция совершенствует хромосомы
Выживаемость организма определяется свойствами хромосом
Требуется оценка каждой хромосомы, отражающая и свойство приспособляемости
Терминология
Ген (свойство, знак, детектор) – атомарный элемент генотипа (набор хромосом, характеризует особь)
Хромосома – упорядоченная последовательность генов (цепочка генов)
Локус (лок) – позиция данного гена в хромосоме
Аллель – значение конкретного гена, определяющее свойство организма
Генотип - набор хромосом, характеризует особь
Фенотип – набор внешних свойств, соответствующих данному генотипу. В формализме ГА фенотип – декодированная из хромосом структура, описывающая реальные параметры задачи.
Популяция - конечное множество особей с одинаковым генотипом
Базовые принципы
А. Имеется некоторая задача, ее решение – новое знание, которое является логическим следствием значений параметров, описывающих данную задачу.
Б. В ГА обрабатываются не сами параметры задачи, а их закодированная форма в виде хромосом.
В. На множестве хромосом определяется целевая функция приспособляемости, оценивающая уровень приспособленности д/каждой хромосомы.
Г. Решение задачи ищется в закодированном виде, обладающих наиболее высокими значениями целевой функции приспособляемости.
Д. При поиске решения используется только целевая функция без ее производных и какой-либо другой дополнительной информации.
Е. Поиск решения происходит не из одной точки (хромосомы), а из некого их множества.
Ё. Порождение нового поколения хромосом осуществляется на основе вероятностных механизмов выбора, а не детерминированности.
Стандартный ГА
Операторы
Селекция – метод рулетки
Оператор скрещивания – одноточечное
Оператор мутации - для каждой хромосомы определяется вероятность её мутации
Сравнительная эффективность
переборные методы
локальные методы – метод Ньютона, градиентного спуска, алгоритм Гаусса
генетические алгоритмы
Применение
Оптимизация функций
Оптимизация запросов в базах данных
Разнообразные задачи на графах (задача коммивояжера, раскраска, нахождение паросочетаний)
Настройка и обучение искусственной нейронной сети
Задачи компоновки
Составление расписаний
Игровые стратегии
Теория приближений
Искусственная жизнь
Биоинформатика (фолдинг белков)
Синтез конечных автоматов
Настройка ПИД регуляторов
Модель Изинга - математическая модель статистической физики, предназначенная для описания намагничивания материала.
Пусть рассматриваются N спинов σi, i=1, 2,... на решетке, каждый может быть направлен либо вверх, либо вниз:
Каждый спин взаимодействует с внешним полем H:
Спины, ближайшие друг к другу также взаимодействуют
Параллельно выстраиваться спинам "выгодно" - этому соответствует энергия -J, J>0(b)
Антипараллельные спины энергетически невыгодны - на образование каждой такой пары требуется положительная энергия +J
Энергия
Статсумма и термодинамические переменные
Статсумма - сумма по всем микроскопическим состояниям системы:
Фазовый переход
Пусть H=0 (нет выделенного направления в пространстве - симметричный случай). Энтропия пытается разупорядочить спины, но параллельные состояния более энергетически выгодны.
При больших T побеждает энтропия: M=0
При малых T перебарывает энергия:
В пограничной точке - точке Кюри Tc - имеет место фазовый переход.
Кабельная теория Ролла - множество предположений и результатов, которые относятся к распространению и взаимодействия электрических сигналов в дендритных деревьях.
Основы кабельной теории - математический анализ затухания сигнала в трансатлантическом телефонном кабеле между Америкой и Англией
два морфологически идентичных дендритных дерева, но с различными электрическими свойствами, могут иметь совершенно разные вычислительные характеристики.