Maximization with an objective variable in the function
显示 更早的评论
Suppose I have a profit (z) equation
z = PY - WL - z^(-iota)*L^(nu)
where P, W, iota and nu are given and iota is some decimal.
I want to maximize z.
Would I be able to find Y, L that maximizes z using fmincon?
If not, which function command is used to solve this problem?
采纳的回答
You can use fmincon,
fmincon(@(x) -x(3), [Y0,L0,z0], [],[],[],[],lb,ub, @(x) nonlcon(x,P, W, iota , nu));
function [c,ceq] = nonlcon(x,P, W, iota , nu)
[Y,L,z]=deal(x(1),x(2),x(3));
ceq = z - P*Y - W*L - z^(-iota)*L^(nu); c=[];
end
but with such a nonlinear constraint, you will need good initial guesses of the unknown parameters Y0,L0,z0 in order to avoid local minima. Because you only have 3 unknowns, it shouldn't be hard to vectorize a discrete grid search for an approximate solution, given appropriate lower and upper bounds lb and ub on the unknowns.
[YY,LL,zz]=ndgrid( linspace(lb(1):ub(1),300) ,...
linspace(lb(2):ub(2),300) , ...
linspace(lb(3):ub(3),300));
con=abs( zz - P.*YY - W.*LL - zz.^(-iota).*LL.^(nu) )<=tolerance;
i=find( zz(con)==max(zz(con)) ,1);
[Y0,L0,z0]=deal(YY(i), LL(i),zz(i)); %approximate solution - intial guess
更多回答（1 个）
Maximilian Schönau
2021-1-4
编辑：Maximilian Schönau
2021-1-4
0 个投票
It seems to me, that fmincon can only find local minima. It also has much math I dont understand, an easier (to understand) function which goes in your direction would be fminsearch.
The problem is, that you probably need a global minimum, which as far as I know is not as easy to find with fminsearch.
I would use fminsearch to find a local minumum, and plot in 2D your function with all the relevant Y and L. That enables you to verify, if the found minimum is the global minimum (for the reasonable values of Y,L). If fminsearch did not find the "global minumum" of your plot, change the starting vector of fminsearch to a value near of the minimum you got out of your plot.
That is the way I have solved this kind of problems, maybe a smart mathematician has a better solution. Maybe there is a smart way you can solve your Problem with symbolic math?
类别
在 帮助中心 和 File Exchange 中查找有关 Calculus 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
