Оптимизации методом покоординатного спуска

clear all; clc;
g = 0.05; % постоянная шага
d = 0.01; % дельта
% Начальная точка
x1 = -3; 
x2 = 3; 
k = 1;  % Счетчик шагов
kmax = 100; % Предельное число шагов, 
% задается для предотвращения зацикливания
% Массивы для хранения промежуточных координат
x1trace = [x1]; 
x2trace = [x2]; 
i = 2; 
while k < kmax 
% Спуск по первой координате
gr1 = exp(-5*x2^2)*(6*cos(3*x1)*sin(3*x1) - 1); 
x1 = x1 + g*gr1
% Сохранение координат
x1trace(i) = x1; 
x2trace(i) = x2; 
i = i + 1; 
% Спуск по второй координате
gr2 = 10*x2*exp(-5*x2^2)*(x1 - sin(3*x1)^2); 
x2 = x2 + g*gr2
% Сохранение координат
x1trace(i) = x1; 
x2trace(i) = x2; 
i = i + 1; 
% Проверка условия останова
if sqrt(gr1^2 + gr2^2) <= d; 
break; % Выход из цикла в случае выполнения условия
end
k = k + 1; 
end
% Построение графика
x = -5:0.1:5; 
y = -5:0.1:5; 
[X, Y] = meshgrid(x, y); 
Z = ((sin(3.*X)).^2.-X)./exp(4.*Y.^2.+Y.^2); 
[C, h] = contour(X, Y, Z); 
clabel(C, h); 
% Отображение меток на линиях уровня
hold on;
plot(x1trace, x2trace, '-'); 
% Вывод начальной точки на график
text(x1trace(1) + 0.2, x2trace(1) + 0.5, 'M0'); 
% Вывод решения на график
text(x1 + 2, x2, ...
strvcat(['x1 = ' (num2str(x1))], ...
        ['x2 = ' (num2str(x2))], ...
        ['k = '  (num2str(k))])); 
hold off;

clear all; clc;
g = 0.05; % постоянная шага
d = 0.01; % дельта
% Начальная точка
x1 = -2; 
x2 = 2; 
k = 1;  % Счетчик шагов
kmax = 100; % Предельное число шагов, 
% задается для предотвращения зацикливания
% Массивы для хранения промежуточных координат
x1trace = [x1]; 
x2trace = [x2]; 
i = 2; 
while k < kmax 
% Спуск по первой координате
gr1 = exp(-x1^2-4*x2^2)*(2*x1^2-x1+3*sin(6*x1)+x1*cos(6*x1) - 1); 
x1 = x1 + g*gr1
% Сохранение координат
x1trace(i) = x1; 
x2trace(i) = x2; 
i = i + 1; 
% Спуск по второй координате
gr2 = -8*x2*exp(-x1^2-4*x2^2)*(sin(3*x1)^2-x1); 
x2 = x2 + g*gr2
% Сохранение координат
x1trace(i) = x1; 
x2trace(i) = x2; 
i = i + 1; 
% Проверка условия останова
if sqrt(gr1^2 + gr2^2) <= d; 
break; % Выход из цикла в случае выполнения условия
end
k = k + 1; 
end
% Построение графика
x = -2:0.1:2; 
y = -2:0.1:2; 
[X, Y] = meshgrid(x, y); 
Z = ((sin(3.*X)).^2.-X)./exp(4.*Y.^2.+X.^2); 
[C, h] = contour(X, Y, Z); 
clabel(C, h); 
% Отображение меток на линиях уровня
hold on;
plot(x1trace, x2trace, '-'); 
% Вывод начальной точки на график
text(x1trace(1) + 0.2, x2trace(1) + 0.5, 'M0'); 
% Вывод решения на график
text(x1 + 2, x2, ...
strvcat(['x1 = ' (num2str(x1))], ...
        ['x2 = ' (num2str(x2))], ...
        ['k = '  (num2str(k))])); 
hold off;

x1 = -1; 
x2 = 1;