problem running fmincon, 'Input arguments to function include colon operator.'

Question

Ekin Ceyda Cetin 2017-1-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/318706-problem-running-fmincon-input-arguments-to-function-include-colon-operator

评论： Walter Roberson 2019-4-30

采纳的回答： Walter Roberson

在 MATLAB Online 中打开

I am trying to run fmincon and i get this error:

'Input arguments to function include colon operator.

To input the colon character, use ':' instead.'

I cannot understand why I get this error.

My code is below.

Anyone any ideas? Thank you.

 % code
  low_b = zeros(size(param));
options = optimoptions('fmincon', ...
    'Algorithm','interior-point', ...   
    'DerivativeCheck', 'on',...   
    'GradObj', 'on',...
    'HessFcn', @hessfh, ...
    'Hessian', 'user-supplied',...
    'Display','final');
[xfinal,fval,exitflag,output] = fmincon(@fh,[0.1;0.1],[],[],[],[],low_b,[],[],options);
 % code 
function [f,gradfh]=fh(param)
syms param1 param2 real
param=[param1;param2];
f3 =((191352748222723/(7128637342969177/1125899906842624)^((param1 - 1)/param1)/(7128637342969177/1125899906842624)^(1/param1))/9511602413006487552 + ((2733703199318633/70368744177664 - param2)^(1 - param1)*(2733703199318633/70368744177664 - param2)^param1)/(561*(param2 - 2733703199318633/70368744177664)))^2 + (log(2733703199318633/70368744177664 - param2)/561 + 2074186557694359/(1125899906842624*param1) - log(-(param2 - 2733703199318633/70368744177664)*(2733703199318633/70368744177664 - param2)^param1)/(2733703199318633/70368744177664 - param2)^param1)^2 
f = matlabFunction(f3,'vars',{param});
gradf3 = jacobian(f3,param).';
gradfh = matlabFunction(gradf3,'vars',{param});
end
    % code
function hessf=hessfh(param)
syms param1 param2 real
param=[param1;param2]
f3 =((191352748222723/(7128637342969177/1125899906842624)^((param1 - 1)/param1)/(7128637342969177/1125899906842624)^(1/param1))/9511602413006487552 + ((2733703199318633/70368744177664 - param2)^(1 - param1)*(2733703199318633/70368744177664 - param2)^param1)/(561*(param2 - 2733703199318633/70368744177664)))^2 + (log(2733703199318633/70368744177664 - param2)/561 + 2074186557694359/(1125899906842624*param1) - log(-(param2 - 2733703199318633/70368744177664)*(2733703199318633/70368744177664 - param2)^param1)/(2733703199318633/70368744177664 - param2)^param1)^2 
gradf3 = jacobian(f3,param).';
hessf3 = jacobian(gradf3,param);
hessf = matlabFunction(hessf3,'vars',{param});
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2017-1-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/318706-problem-running-fmincon-input-arguments-to-function-include-colon-operator#answer_249075

The second output of the objective function should be the gradient evaluated at the input point, not a function handle. Likewise the hessian should be the value of the hessian, not a function handle.

11 个评论
显示 9更早的评论隐藏 9更早的评论

Ekin Ceyda Cetin 2017-1-1

在 MATLAB Online 中打开

If I do not use the gradient and hessian as function handle,

then i get the error from 'deal'.

'Error using deal

The number of outputs should match the number of inputs.'

Any suggestions in this situation?

    % code
syms param1 param2
param=[param1;param2];
f3 =((191352748222723/(7128637342969177/1125899906842624)^((param1 - 1)/param1)/(7128637342969177/1125899906842624)^(1/param1))/9511602413006487552 + ((2733703199318633/70368744177664 - param2)^(1 - param1)*(2733703199318633/70368744177664 - param2)^param1)/(561*(param2 - 2733703199318633/70368744177664)))^2 + (log(2733703199318633/70368744177664 - param2)/561 + 2074186557694359/(1125899906842624*param1) - log(-(param2 - 2733703199318633/70368744177664)*(2733703199318633/70368744177664 - param2)^param1)/(2733703199318633/70368744177664 - param2)^param1)^2;
gradf3 = jacobian(f3,param).';
fh = matlabFunction(f3,gradf3,'vars',{param});
hessf3 = jacobian(gradf3,param);
hessfh = matlabFunction(hessf3,'vars',{param});
low_b = zeros(size(param));
options = optimoptions('fmincon', ...
    'Algorithm','interior-point', ...   
    'DerivativeCheck', 'on',...   
    'GradObj', 'on',...
    'HessFcn', hessfh, ...
    'Hessian', 'user-supplied',...
    'Display','final');
[xfinal,fval,exitflag,output] = fmincon(fh,[0.1;0.1],[],[],[],[],low_b,[],[],options);

Ekin Ceyda Cetin 2017-1-2

在 MATLAB Online 中打开

I tried to do as you say, I am getting colon operator error as I get in the beginning.

'Input arguments to function include colon operator. To input the colon character, use ':' instead.'

What am I doing wrong? Obviously I am missing sth but don't know what.

% code
syms param1 param2
param=[param1;param2];
f3 =((191352748222723/(7128637342969177/1125899906842624)^((param1 - 1)/param1)/(7128637342969177/1125899906842624)^(1/param1))/9511602413006487552 + ((2733703199318633/70368744177664 - param2)^(1 - param1)*(2733703199318633/70368744177664 - param2)^param1)/(561*(param2 - 2733703199318633/70368744177664)))^2 + (log(2733703199318633/70368744177664 - param2)/561 + 2074186557694359/(1125899906842624*param1) - log(-(param2 - 2733703199318633/70368744177664)*(2733703199318633/70368744177664 - param2)^param1)/(2733703199318633/70368744177664 - param2)^param1)^2;
gradf3 = jacobian(f3,param).';
hessf3 = jacobian(gradf3,param);
hessfh = matlabFunction(hessf3,'vars',{param});
low_b = zeros(size(param));
options = optimoptions('fmincon', ...
  'Algorithm','interior-point', ...   
  'DerivativeCheck', 'on',...   
  'GradObj', 'on',...
  'HessFcn', hessfh, ...
  'Hessian', 'user-supplied',...
  'Display','final');
[xfinal,fval,exitflag,output] = fmincon(@(param) fh(param, f, gradfh),[0.1;0.1],[],[],[],[],low_b,[],[],options);
    % code
function [f,gradfh]=fh(param, f, gradfh)
syms param1 param2
param=[param1;param2];
f3 =((191352748222723/(7128637342969177/1125899906842624)^((param1 - 1)/param1)/(7128637342969177/1125899906842624)^(1/param1))/9511602413006487552 + ((2733703199318633/70368744177664 - param2)^(1 - param1)*(2733703199318633/70368744177664 - param2)^param1)/(561*(param2 - 2733703199318633/70368744177664)))^2 + (log(2733703199318633/70368744177664 - param2)/561 + 2074186557694359/(1125899906842624*param1) - log(-(param2 - 2733703199318633/70368744177664)*(2733703199318633/70368744177664 - param2)^param1)/(2733703199318633/70368744177664 - param2)^param1)^2 
f = matlabFunction(f3,'vars',{param})
if nargout > 1 % gradient required
    gradf3 = jacobian(f3,param).';
    gradfh = matlabFunction(gradf3,'vars',{param});
end
end

Ekin Ceyda Cetin 2017-1-2

在 MATLAB Online 中打开

Do you mean, i need to use the function handle just to make the choice of calculating only objective function or calculating gradient as well if necessary? I tried that too, not enough input arguments. Is this what you mean?

So sorry for being such a trouble :'

% code
syms param1 param2
param=[param1;param2];
f3 =((191352748222723/(7128637342969177/1125899906842624)^((param1 - 1)/param1)/(7128637342969177/1125899906842624)^(1/param1))/9511602413006487552 + ((2733703199318633/70368744177664 - param2)^(1 - param1)*(2733703199318633/70368744177664 - param2)^param1)/(561*(param2 - 2733703199318633/70368744177664)))^2 + (log(2733703199318633/70368744177664 - param2)/561 + 2074186557694359/(1125899906842624*param1) - log(-(param2 - 2733703199318633/70368744177664)*(2733703199318633/70368744177664 - param2)^param1)/(2733703199318633/70368744177664 - param2)^param1)^2;
f = matlabFunction(f3,'vars',{param});
gradf3 = jacobian(f3,param).';
gradfh = matlabFunction(gradf3,'vars',{param});
hessf3 = jacobian(gradf3,param);
hessfh = matlabFunction(hessf3,'vars',{param});
low_b = zeros(size(param));
param0=[3.6;2.6];
options = optimoptions('fmincon', ...
  'Algorithm','interior-point', ...   
  'DerivativeCheck', 'on',...   
  'GradObj', 'on',...
  'HessFcn', hessfh, ...
[xfinal,fval,exitflag,output] = fmincon(@(param1, param2) oneOrTwo(param1, param2,f,gradfh),param0,[],[],[],[],low_b,[],[],options);
% code
function [a, b] = oneOrTwo(c, d, h, g) 
syms c d
a = h(c,d);
if nargout > 1
 b = g(c,d) ;
end
end

Walter Roberson 2017-1-2

在 MATLAB Online 中打开

% code
syms param1 param2
param=[param1;param2];
f3 =((191352748222723/(7128637342969177/1125899906842624)^((param1 - 1)/param1)/(7128637342969177/1125899906842624)^(1/param1))/9511602413006487552 + ((2733703199318633/70368744177664 - param2)^(1 - param1)*(2733703199318633/70368744177664 - param2)^param1)/(561*(param2 - 2733703199318633/70368744177664)))^2 + (log(2733703199318633/70368744177664 - param2)/561 + 2074186557694359/(1125899906842624*param1) - log(-(param2 - 2733703199318633/70368744177664)*(2733703199318633/70368744177664 - param2)^param1)/(2733703199318633/70368744177664 - param2)^param1)^2;
f = matlabFunction(f3,'vars',{param});
gradf3 = jacobian(f3,param).';
gradfh = matlabFunction(gradf3,'vars',{param});
hessf3 = jacobian(gradf3,param);
hessfh = matlabFunction(hessf3,'vars',{param});
low_b = zeros(size(param));
param0=[3.6;2.6];
options = optimoptions('fmincon', ...
  'Algorithm','interior-point', ...   
  'DerivativeCheck', 'on',...   
  'GradObj', 'on',...
  'HessFcn', hessfh, ...
[xfinal,fval,exitflag,output] = fmincon(@(params) oneOrTwo(params, f,gradfh), param0, [], [], [], [], low_b, [], [], options);
% code
function [fun_val, grad_val] = oneOrTwo(params, fun_h, grad_h) 
  fun_val = fun_h(params);
  if nargout > 1
   grad_val = grad_h(params) ;
  end

Rustem Devletov 2019-4-26

Hello, friend. I'm really really feeling despondent. I don't know how to implement the following idea. My supervisor told me to supple the gradient to a function from another file. How can I implement this? Can you suggest anything?