Найти минимум многопеременной функции

Question

Вопрос : Найти минимум f (x, y) = x ^ 2 + y ^ 2-2 * x-6 * y + 14 в окне [0,2] × [2 , 4] с шагом 0,01 для х и у.

Мой подход : Найти первые частные производные fx и fy. Критические точки удовлетворяют уравнениям fx (x, y) = 0 и fy (x, y) = 0 одновременно. найти частные производные второго порядка fxx (x, y), fyy (x, y) и fxy (x, y), чтобы найти D.

clc
clear all
syms x y
fun=x^2+y^2-2*x-6*y+14;
fx=diff(fun,x);
fy=diff(fun,y);
pt=solve(fx==0,fy==0);
sol = struct2array(pt)
fxx=diff(fx,x);
fyy=diff(fy,y);
fxy=diff(fx,y);
D=subs(fxx,[x y],[1 3])*subs(fyy,[x y],[1 3])-(subs(fxy,[x y],[1 3]))^2
fxx_val=subs(fxx,[x y],[1 3])
minimum_value=subs(fun,[x y],[1 3])

Правильно ли я делаю то, что задал вопрос? Кроме того, что насчет окна и инкремента упоминается этот вопрос. Любые советы или решения будут оценены.
Заранее спасибо.

Answer 1

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

Пожалуйста, прочитайте код

f = @(x,y)x.^2+y.^2-2.*x-6.*y+14;

% x range
x_lb = 0;
x_ub = 2;

% y range
y_lb = 2;
y_ub = 4;

step = 0.01;

% lower bound of x, initial guess as xmin
xmin = x_lb;

% lower bound of y, initial guess as ymin
ymin = y_lb;

% f at the lower bounds, initial  fmin 
fmin = f(xmin, ymin);

for x = x_lb:step:x_ub

    for y = y_lb:step:y_ub
        % function evaluation
        fval = f(x, y);

        %replace fmin if the newly evaluated f is less than the actual fmin
        if fval < fmin
            fmin = fval;

            % save current x and y where f is minimum 
            xmin = x;
            ymin = y;

        end
    end

end

Решение

xmin = 1;
ymin = 3;
fmin = 4;

Answer 2

Я бы предложил использовать возможности Matlab для вычисления с помощью матриц.Тогда цикл не требуется.

% your function, look up anonymous functions 
func = @(x,y) x.^2 + y.^2 - 2.*x - 6.*y + 14;

% get matrices for you x- and y-window
[xg, yg] = meshgrid(0:.01:2, 2:0.01:4);

% compute all in one call
result = func(xg,yg);

% find total minimum
minimum = min(result(:));

% find the index of the (first) minimum, for other equations, there might
% be more than one
ind = find(result==minimum, 1);

% Output the result
fprintf('The minimum (%d) is located at x: %d, y: %d.\n', minimum,  xg(ind), yg(ind));