Constraint function in optimization toolbox

2020 6 30
1 个回答
回答已采纳
更新时间：2020 6 30
4 次查看（30 天）

0 个投票

One of the variables in my model is Measured Power which is being calculated by an external black box model connected with MATLAB. The way its calculated is below:

4 个评论
显示 2更早的评论 隐藏 2更早的评论

Rik
Rik 2020-6-30
Why are you using so many globals? They should be avoided as much as possible anyway, and the function you posted is ignoring most of them, as well as the input.
Rik
Rik 2020-6-30
Comment sent as email:
Can I send you an email with more details as the program is proprietary?
Rik
Rik 2020-6-30
Since this isn't really my field I don't know if I can solve the actual problem. I just noticed the globals, which tend to be a bad sign.
Also, this forum is a hobby, private consultation is a job, and it isn't my job. There are plenty of people you could hire.
Hussain Ja
Hussain Ja 2020-6-30
I have skipped all the lines of code to show my constraint which is "Powermeasm". I am using the same code for the main program and the constraint function. The constraint function program has just 2 additional lines declaring the two constraints; c and ceq. I would like to know if this procedure is correct? Can I caclulate "Powermeasm" in the main program and provide a constraint on it through another constraint program? If so, why isn't it working?

请先登录，再进行评论。

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

 采纳的回答

Alan Weiss
Alan Weiss 2020-6-30

0 个投票

I think that you missed the step where you call the external function to calculate the nonlinear constraint.
function [c,ceq] = con_EbsOptimize4(para)
global asmInfo app oc model objects Pr Hu Pmeasm Hmeasm Power Powermeasm;
model.Simulate(errors);
Power = oc.CastToComp11(objects.Item('Generator_1'));
Powermeasm = Power.QREAL.Value;
c = Powermeasm-192000;
ceq = [];
end
There are two other things that you might want to keep in mind for this:
  1. Often, the function that calculates the objective function is the same as the function that calculates the nonlinear constraint. If this is true for your function, consider using the technique in Objective and Nonlinear Constraints in the Same Function.
  2. If you are using fmincon as your solver, keep in mind the problems in finite differences mentioned in Optimizing a Simulation or Ordinary Differential Equation.
Alan Weiss
MATLAB mathematical toolbox documentation

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

Hussain Ja
Hussain Ja 2020-6-30
编辑：Hussain Ja 2020-6-30
I am using the same code for the main program and the constraint function. The constraint function program has just 2 additional lines declaring the two constraints; c and ceq. I would like to know if this procedure is correct?
When I run the code you have shown above I get the following error:
Optimization running.
Error running optimization.
Output argument "c" (and maybe others) not assigned during call to "con_EbsOptimize4".
I also get the following warning:
Warning: The value of local variables may have been changed to match the globals. Future versions of MATLAB will require
that you declare a variable to be global before you use that variable.
Alan Weiss
Alan Weiss 2020-6-30
I'm sorry, I canot tell what you did. Please include your code. Feel free to not include all the code that calculates the values, but I do need to see the fmincon call and the constraint function definition, at least the first few lines of that function.
Alan Weiss
MATLAB mathematical toolbox documentation
Hussain Ja
Hussain Ja 2020-6-30
The main program is running fine. The constraint function doesn't appear to be working as the output I get is violating the constraint.
Alan Weiss
Alan Weiss 2020-6-30
I see the objective function calculation and the constraint function calculation. I also need to see the fmincon call, and any options settings.
Alan Weiss
MATLAB mathematical toolbox documentation
Hussain Ja
Hussain Ja 2020-6-30
I call it as follows:
code11(2,[1.8,0.1],[3,0.9],15,3)
Alan Weiss
Alan Weiss 2020-6-30
Things look OK to me. Try setting MaxGenerations_Data to 1e5 and see if that helps. The nonlinear constraint solver takes many, many generations to work.
I suggest that you might want to try patternsearch on this problem, and maybe surrogateopt if you have a recent MATLAB version, as ga is not very efficient when handling nonlinear constraints.
Alan Weiss
MATLAB mathematical toolbox documentation
Hussain Ja
Hussain Ja 2020-6-30
My concern is not the optimality of the solution. It's just that the solution that it returns clearly violates the supplied constraint by a significant margin. So I have a hunch that somehow it's not calculating the constraint violation correctly and giving it as zero when clearly the violation is much much larger
Alan Weiss
Alan Weiss 2020-6-30
You didn't say before that "...it's not calculating the constraint violation correctly and giving it as zero when clearly the violation is much much larger." Can you give iterative display data showing what you mean, or something else to help me understand what might be happening? And did you try my previous suggestion, to set MaxGenerations_Data to 1e5?
Alan Weiss
MATLAB mathematical toolbox documentation
Hussain Ja
Hussain Ja 2020-6-30
Thanks Mr. Weiss. I figured out the problem. The code is working fine now. I had written the wrong generator name.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Get Started with Optimization Toolbox 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by