Пользовательская функция построения графика

О пользовательских функциях построения графика

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

Создание пользовательской функции построения графика

Чтобы создать функцию построения графика для этого примера, скопируйте и вставьте следующий код в новый файл функции в MATLAB^® Редактор:

function stop = psplotchange(optimvalues, flag)
% PSPLOTCHANGE Plots the change in the best objective function 
% value from the previous iteration.
  
% Best objective function value in the previous iteration
persistent last_best
 
stop = false;
if(strcmp(flag,'init')) 
        set(gca,'Yscale','log'); %Set up the plot
        hold on;
        xlabel('Iteration'); 
        ylabel('Log Change in Values');
        title(['Change in Best Function Value']);
end
 
% Best objective function value in the current iteration
best = min(optimvalues.fval);  
 
 % Set last_best to best
if optimvalues.iteration == 0
last_best = best;
        
else
        %Change in objective function value
			 change = last_best - best; 
        plot(optimvalues.iteration, change, '.r');
end

Сохраните файл как psplotchange.m в папке на пути MATLAB. Код объяснен в Как работы Функции построения графика.

Настройте проблему

Проблема совпадает с Ограниченной Минимизацией Используя patternsearch, и Оптимизируйте Live Editor Тэска. Настраивать проблему:

Введите следующее в командной строке MATLAB:

x0 = [2 1 0 9 1 0];
Aineq = [-8 7 3 -4 9 0];
bineq = 7;
Aeq = [7 1 8 3 3 3; 5 0 -5 1 -5 8; -2 -6 7 1 1 9; 1 -1 2 -2 3 -3];
beq = [84 62 65 1];

Поскольку это - линейно ограниченная проблема, установите PollMethod опция к 'GSSPositiveBasis2N'. Включайте обоих @psplotbestf встроенная функция построения графика и пользовательская функция построения графика @psplotchange в опциях.
```
options = optimoptions('patternsearch',...
    'PlotFcn',{@psplotbestf,@psplotchange},...
    'PollMethod','GSSPositiveBasis2N');
```

Запустите оптимизацию с пользовательской функцией построения графика

Запустите пример путем вызова patternsearch запуск с x0.

[x,fval] = patternsearch(@lincontest7,x0,...
    Aineq,bineq,Aeq,beq,[],[],[],options);

Поскольку шкала y - ось в более низком пользовательском графике является логарифмической, график показывает только изменения, которые больше 0.

Как работы функции построения графика

Функция построения графика использует информацию, содержавшуюся в следующих структурах.

optimvalues — Текущее состояние решателя, структуры
flag — Текущий статус алгоритма, вектора символов

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

Пользовательские операторы функции построения графика

Оператор	Описание
`persistent last_best`	Создает персистентную переменную `last_best`, лучшее значение целевой функции в предыдущем поколении. Персистентные переменные сохраняются по множественным вызовам функции построения графика.
`set(gca,'Yscale','log')`	Настраивает график, прежде чем алгоритм запустится.
`best = min(optimvalues.fval)`	Наборы `best` равняйтесь минимальному значению целевой функции. Поле `optimvalues.fval` содержит значение целевой функции в текущей итерации. Переменная `best` минимальное значение целевой функции. Для полного описания полей структуры `optimvalues`, смотрите Структуру Функций построения графика.
`change = last_best - best`	Устанавливает переменную `change` к лучшему значению целевой функции в предыдущей итерации минус лучшее значение целевой функции в текущей итерации.
`plot(optimvalues.iteration, change, '.r')`	Строит переменную `change` в текущем значении целевой функции, для текущей итерации содержал in`optimvalues.iteration`.

Документация