Построение графика во время подбора параметра parfor

В этом примере показано, как выполнить развертку параметра параллельно, и график прогрессируют во время параллельных расчетов. Можно использовать DataQueue контролировать результаты во время расчетов на параллельном пуле. Можно также использовать DataQueue с функциями языка параллельного программирования, такими как parfor, parfeval и spmd.

Пример показывает, как выполнить развертку параметра в классической системе, Осцилляторе Ван дер Поля. Эта система может быть выражена как набор ОДУ, зависящих от двух параметров Осциллятора Ван дер Поля,$\mu$ и$\nu$:

$$\dot{x } =\nu y$$

$$\dot{y } =\mu \left(1-x^2 \right)y -x$$

Можно выполнить параллельную развертку параметра по параметрам$\mu$ и$\nu$ использованию parfor цикл, чтобы узнать средний период$y$ при варьировании их. Следующая анимация показывает выполнение этого примера в локальном кластере.

Настройте значения развертки параметра

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

gridSize = 6;
nu = linspace(100, 150, gridSize);
mu = linspace(0.5, 2, gridSize);
[N,M] = meshgrid(nu,mu);

Подготовьте объемную поверхностную диаграмму визуализировать результаты

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

Z = nan(size(N));
c = surf(N, M, Z);
xlabel('\mu Values','Interpreter','Tex')
ylabel('\nu Values','Interpreter','Tex')
zlabel('Mean Period of  y')
view(137, 30)
axis([100 150 0.5 2 0 500]);

Создайте DataQueue, чтобы выбрать результаты во время развертки параметра

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

D = parallel.pool.DataQueue;
D.afterEach(@(x) updateSurface(c, x));

Выполните развертку параметра и постройте результаты

Используйте parfor выполнять параллельную развертку параметра. Дайте рабочим команду решать систему для каждой комбинации параметров в meshgrid и вычислять средний период. Сразу передайте результат обратно каждой итерации клиенту, когда рабочий закончит расчеты.

parfor ii = 1:numel(N)
    [t, y] = solveVdp(N(ii), M(ii));
    l = islocalmax(y(:, 2));
    send(D, [ii mean(diff(t(l)))]);
end
Starting parallel pool (parpool) using the 'local' profile ...
connected to 6 workers.

Масштабируйте до кластера

Если у вас есть доступ к кластеру, можно увеличить расчет. Для этого удалите предыдущий parpool, и откройте новый с помощью профиля в большем кластере. Код ниже показов кластерный профиль под названием 'MyClusterInTheCloud'. Чтобы запустить этот код самостоятельно, необходимо заменить 'MyClusterInTheCloud' с именем вашего кластерного профиля. Настройте количество рабочих. Пример показывает 4 рабочим. Увеличьте размер полного расчета путем увеличения размера сетки.

gridSize = 25;
delete(gcp('nocreate'));
parpool('MyClusterInTheCloud',4);
Parallel pool using the 'local' profile is shutting down.

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

Функции помощника

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

function [t, y] = solveVdp(mu, nu)
f = @(~,y) [nu*y(2); mu*(1-y(1)^2)*y(2)-y(1)];
[t,y] = ode23s(f,[0 20*mu],[2; 0]);
end

Объявите, что функция для DataQueue обновляет график результатами, которые прибывают от рабочих.

function updateSurface(s, d)
s.ZData(d(1)) = d(2);
drawnow('limitrate');
end
Starting parallel pool (parpool) using the 'MyClusterInTheCloud' profile ...
connected to 4 workers.

Смотрите также

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