Графические модели (курс лекций)/2013/Задание 6

Материал из MachineLearning.

Перейти к: навигация, поиск


Начало выполнения задания: 29 апреля 2013 г.
Срок сдачи: 26 мая 2013 г. (воскресенье), 23:59.

Программная среда для выполнения задания — MATLAB. Неэффективная реализация кода может негативно отразиться на оценке.

Модель Изинга

Прямоугольная система соседства
Прямоугольная система соседства
Треугольная система соседства
Треугольная система соседства

Модель Изинга — математическая модель статистической физики, предназначенная для описания магнитных свойств вещества. Каждой вершине кристаллической решётки (рассматривается двухмерный случай) сопоставляется число, называемое спином и равное +1 или −1 («поле вверх»/«поле вниз»). Каждому из 2^N возможных вариантов расположения спинов (где N — число атомов решётки) приписывается энергия, состоящая из взаимодействия спинов соседних атомов J и действия внешнего магнитного поля H:

E(X) = -\left( J \sum_{(i,j) \in\mathcal{E}} x_i x_j + \sum_{i=1}^N H_i x_i \right),

где x_i — переменные, соответствующие спинам, \mathcal{E} — система соседства (в данном задании рассматриваются две системы соседства: прямоугольная и треугольная), каждое ребро в \mathcal{E} считается один раз. Вероятность нахождения в каждом конкретном состоянии X задается распределением Гиббса:

P(X) = \frac{1}{Z} \exp{\left( -\beta E(X) \right)}, \qquad \beta = \frac{1}{kT},

где Z — нормировочная константа, T — температура, k — параметр.

Если J = 1 , то вещество называется ферромагнетиком. Если J = -1, то вещество называется антиферромагнетиком.

Формулировка задания

Пример иллюстрации состояния модели Изинга размера 20 на 20.
Пример иллюстрации состояния модели Изинга размера 20 на 20.
  1. Провести исследование модели Изинга с помощью схемы Гиббса:
    • Вывести формулы для одномерных условных распределений вида p(x_i | X_{\backslash{i}});
    • Реализовать с помощью схемы Гиббса процедуру оценки математического ожидания энергии \frac{1}{N}\mathbb{E}E, стандартного отклонения энергии \frac{1}{N}\sqrt{\mathbb{D}E}, математического ожидания квадрата общей намагниченности модели \sqrt{\mathbb{E}(\mu^2)}, где \mu = \left(\sum_{i=1}^N x_i \right) / N для заданных параметров \beta и внешнего магнитного поля H. Требования по эффективности реализации: 1000 итераций метода Гиббса для прямоугольной решетки размера 20 на 20 и ста значений параметра \beta должны выполняться не более 20 секунд. В отчете привести соответствующие замеры времени работы кода;
    • Построить графики зависимости \frac{1}{N}\mathbb{E}E,\ \frac{1}{N}\sqrt{\mathbb{D}E},\ \sqrt{\mathbb{E}(\mu^2)} от температуры для треугольной и четырехугольной систем соседства, ферромагнетика и антиферромагнетика (всего 4 модели). Проинтерпретировать полученные результаты (в частности идентифицировать, локализовать и изучить фазовый переход). Для построения графиков использовать следующие значения параметров: размер решетки 20 на 20 (N = 400), k = 1, 10000 итераций метода Гиббса для оценки статистик, внешнее магнитное поле H_i = 0, температуры T = 0.5 : 0.1 : 10;
    • Для ферромагнетика с четырехугольной системой соседства отобразить характерные отдельные конфигурации X (см. рис.) в зависимости от температуры (низкая температура, окрестность фазового перехода, высокая температура). Проинтерпретировать результаты;
    • Исследовать влияние константного внешнего магнитного поля на фазовый переход в ферромагнетике с прямоугольной системой соседства. Рассмотреть также внешнее магнитное поле со следующей структурой: на половине решетки H = 1, на другой половине H = -1;
  2. Провести исследование модели Изинга с помощью вариационного подхода, где в качестве факторизованного семейства для q(X) выступает независимое распределение по всем компонентам q(X)=\prod_{i=1}^Nq_i(x_i):
    • Вывести формулу для пересчета отдельного фактора q_i(x_i), формулу для функционала \mathcal{L}(q), формулы для оценки \frac{1}{N}\mathbb{E}E,\ \frac{1}{N}\sqrt{\mathbb{D}E},\ \sqrt{\mathbb{E}(\mu^2)}, где мат.ожидания вычисляются по распределению q(X);
    • Реализовать с помощью вариационного подхода процедуру оценки величин \frac{1}{N}\mathbb{E}E,\ \frac{1}{N}\sqrt{\mathbb{D}E},\ \sqrt{\mathbb{E}(\mu^2)}, а также нижней оценки \mathcal{L}(q) для логарифма нормировочной константы \log Z для заданных параметров \beta и внешнего магнитного поля H. Требования по эффективности реализации: для решетки размера 20 x 20 сто итераций поиска вариационного приближения для ста значений параметра \beta должны выполняться не более одной секунды, включая этап оценки всех статистик. В отчете привести соответствующие замеры времени работы кода;
    • Построить графики зависимости \frac{1}{N}\mathbb{E}E,\ \frac{1}{N}\sqrt{\mathbb{D}E},\ \sqrt{\mathbb{E}(\mu^2)},\ \mathcal{L}(q) от температуры для треугольной и четырехугольной систем соседства, ферромагнетика и антиферромагнетика (всего 4 модели). Проинтерпретировать полученные результаты (в частности идентифицировать, локализовать и изучить фазовый переход). Для построения графиков использовать следующие значения параметров: размер решетки 20 на 20 (N = 400), k = 1, внешнее магнитное поле H_i = 0, температуры T = 0.5 : 0.1 : 10;
    • Для ферромагнетика с четырехугольной системой соседства отобразить найденное распределение q(X) (в шкале серого) в зависимости от температуры (низкая температура, окрестность фазового перехода, высокая температура). Проинтерпретировать результаты;
    • Исследовать влияние константного внешнего магнитного поля на фазовый переход в ферромагнетике с прямоугольной системой соседства. Рассмотреть также внешнее магнитное поле со следующей структурой: на половине решетки H = 1, на другой половине H = -1;
  3. Сравнить результаты схемы Гиббса с результатами вариационного подхода. Рассмотреть ферромагнетик с прямоугольной системой соседства. Привести графики математического ожидания и дисперсии энергии, корня из математического ожидания квадрата намагниченности в одних осях для двух подходов.
  4. Написать отчет в формате PDF с описанием всех проведенных исследований.

Рекомендации по выполнению задания

  1. Для схемы Гиббса и вариационного подхода:
    • Рекомендуется реализовывать алгоритмы вывода с векторными операциями по параметру \beta, т.е. проводить вычисления для всех температур сразу;
    • На этапе тестирования процедур рекомендуется проводить эксперименты на неквадратной решетке и со случайным внешним магнитным полем H;
  2. Для схемы Гиббса:
    • Для оценки статистик распределения с помощью схемы Гиббса следует отбрасывать конфигурации X, полученные на первой трети итераций;
    • В качестве примеров конфигураций X лучше брать ситуации с последних итераций схемы Гиббса;
    • Одной из возможных проверок на корректность оценки статистик по сгенерированному набору конфигураций X является следующая: сделать большое число итераций по схеме Гиббса, рассмотреть несколько различных подмножеств сгенерированных конфигураций, по каждому из подмножеств оценить статистики. При правильной реализации статистики, оцененные по разным подмножествам конфигураций, должны быть близки между собой;
  3. Для вариационного подхода:
    • Рекомендуется запускать вариационную оптимизацию из нескольких начальных приближений для q: случайные, все вероятности одинаковые и близки к единице, шахматная доска и т.д. При этом для текущей температуры наилучшим вариационным приближением признается такое, при котором значение нижней границы \mathcal{L}(q) является максимальным;
    • Одной из проверок на корректность оптимизации в рамках вариационного подхода является монотонное возрастание в итерациях значения нижней границы \mathcal{L}(q);
    • Одной из возможных проверок на корректность вычисления различных статистик по приближению q является использование метода Монте Карло, при котором из распределения q генерируется набор конфигураций X^j, по которым затем оцениваются необходимые статистики.

Оформление задания

Выполненное задание следует отправить письмом по адресу bayesml@gmail.com с заголовком письма «[ГМ13] Задание 6, Фамилия». Убедительная просьба присылать выполненное задание только один раз с окончательным вариантом. Новые версии будут рассматриваться только в самом крайнем случае. Также большая просьба строго следовать указанным ниже прототипам реализуемых функций.

Присланный вариант задания должен содержать в себе:

  • Текстовый файл в формате PDF, содержащий описание проведенных исследований;
  • Все исходные коды с необходимыми комментариями.

Для проверки корректности прототипов реализованных функций рекомендуется воспользоваться специальной процедурой. Задания, для которых данная процедура будет выходить с ошибкой, проверяться не будут!

Метод Гиббса для оценки статистик распределений
[E, D, M, S] = gibbsIsing4(H, J, betaAll, num_iter) — прямоугольная система соседства
[E, D, M, S] = gibbsIsing6(H, J, betaAll, num_iter) — треугольная система соседства
ВХОД
H — внешнее магнитное поле, матрица размера vS x hS;
J — параметр модели, равен 1 или -1;
betaAll — вектор значений параметра \beta (вектор-строка длины \beta_0);
num_iter — количество итераций схемы Гиббса;
ВЫХОД
E — значения мат.ожиданий энергии на один спин \frac{1}{N}\mathbb{E}E для каждой температуры, вектор-строка длины \beta_0;
D — значения стандартных отклонений энергии на один спин \frac{1}{N}\sqrt{\mathbb{D}E} для каждой температуры, вектор-строка длины\beta_0;
M — значения средней магнетизации на один спин \sqrt{\mathbb{E}\mu^2} для каждой температуры, вектор-строка длины \beta_0;
S — примеры конфигураций X для всех температур, массив размера vS x hS x \beta_0.

 

Вариационный подход
[E, D, M, L] = varIsing4(H, J, betaAll, opt_params) — прямоугольная система соседства
[E, D, M, L] = varIsing6(H, J, betaAll, opt_params) — треугольная система соседства
ВХОД
H — внешнее магнитное поле, матрица размера vS x hS;
J — параметр модели, равен 1 или -1;
betaAll — вектор значений параметра \beta (вектор-строка длины \beta_0);
opt_params — (необязательный параметр) параметры оптимизационного процесса, структура со следующими полями:
'max_iter' — максимальное количество итераций, по умолчанию = 300;
'tol_crit' — необходимая точность по значению нижней границы, по умолчанию = 10^{-4};
'num_start' — количество различных начальных приближений, по умолчанию = 1;
ВЫХОД
E — мат.ожидание энергии на один спин \frac{1}{N}\mathbb{E}E для каждой температуры, вектор-строка длины \beta_0;
D — стандартное отклонение энергии на один спин \frac{1}{N}\sqrt{\mathbb{D}E} для каждой температуры, вектор-строка длины \beta_0;
M — средняя магнетизация \sqrt{\mathbb{E}\mu^2} для каждой температуры, вектор-строка длины \beta_0;
L — нижние границы для логарифмов нормировочных констант для каждой температуры, вектор-строка длины \beta_0.
Личные инструменты