Maximization with an objective variable in the function
7 views (last 30 days)
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?
Matt J on 4 Jan 2021
Edited: Matt J on 4 Jan 2021
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)
ceq = z - P*Y - W*L - z^(-iota)*L^(nu); c=;
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) , ...
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
More Answers (1)
Maximilian Schönau on 4 Jan 2021
Edited: Maximilian Schönau on 4 Jan 2021
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?