考虑以下Matlab函数
function [f, dfx1, dfx2] = optifun(x1,x2)
f = x1(1)^2 + x1(2)^2 + x2^2;
% Gradients
dfx1(1) = 2*x1(1);
dfx1(2) = 2*x1(2);
dfx2 = 2*x2;
我的目标是使用函数`fminunc对x1
和x2
优化上述函数,该函数还包含渐变。如果以上述方式定义,我不知道是否可以优化功能。
我的方法是(但我知道它不会起作用):
options = optimoptions('fmincon', 'SpecifyObjectiveGradient',true);
% Initializing
x10 = [1, 1];
x20 = 1;
[t1, t2] = fminunc(@(x1, x2)optifun(x1,x2), x10, x20, options);
编辑:我已根据用户指出进行了更正:@ m7913d
答案 0 :(得分:0)
主要问题
您的函数签名与fminunc
的期望不符:您只能指定一个(初始)x
向量和一个渐变向量,它应包含所有变量/渐变。如果您不想optifun
(这是首选解决方案),您可以按如下方式定义辅助函数:
function [f, df] = optifun_helper(x)
[f, dfx1, dfx2] = optifun(x(1:2),x(3));
df = [dfx1 dfx2];
end
使用此功能解决您的优化问题:
[x] = fminunc(@(x) optifun_helper(x), [x10 x20], options);
次要问题
您的代码还有其他问题。我得到的第一个错误是:
Error using optimoptions (line 118)
Invalid solver specified. Provide a solver name or handle (such as 'fmincon' or @fminunc).
Type DOC OPTIMOPTIONS for a list of solvers.
如建议的错误消息,您应该将求解器指定为第一个参数:
options = optimoptions('fminunc','SpecifyObjectiveGradient',true);
第二个问题是由于拼写问题(您定义optifun
而非df1
),dfx1
的某些输出变量从未定义过:
function [f, dfx1, dfx2] = optifun(x1,x2)
f = x1(1)^2 + x1(2)^2 + x2^2;
% Gradients
dfx1(1) = 2*x1(1);
dfx1(2) = 2*x1(2);
dfx2 = 2*x2;
end
作为结论,请务必仔细阅读错误消息和文档,并尝试使用预期的语法来填充代码。