One-Variable Optimization Problem - MathOverflow most recent 30 from http://mathoverflow.net2013-05-21T16:15:47Zhttp://mathoverflow.net/feeds/question/45296http://www.creativecommons.org/licenses/by-nc/2.5/rdfhttp://mathoverflow.net/questions/45296/one-variable-optimization-problemOne-Variable Optimization ProblemVenous0072010-11-08T12:07:06Z2010-12-18T09:36:17Z
<p>$W_{opt}=\arg {\max(\pi_0 F_{L_0}(W)-\frac{\pi_1}{W}\int_0^W F_{L_1}( \alpha )d \alpha )}$ </p>
<p>subject to $\quad \int_0^W F_{L_0} (\alpha)d\alpha <\xi$ </p>
<p>We should find analytically the optimal $W >0$ which maximize the first equation subject to the second equation, where $F( \cdot )$ is comulative distribution function (CDF), and $L_0$ and $L_1$ are positive random variables. $\xi$, $\pi_0$, $\pi_1$ are constant. Also, $0<\pi_0, \pi_1<1$ and $\pi_0 + \pi_1 =1$. All variables are real. Further, if needed, we can assume that, for example, $L_0$ and $L_1$ may have Erlang or exponential distribution.</p>
http://mathoverflow.net/questions/45296/one-variable-optimization-problem/45313#45313Answer by Gilead for One-Variable Optimization ProblemGilead2010-11-08T15:20:20Z2010-11-08T16:19:49Z<p>Hmm, you could rewrite it this way: (I'm going to assume that $\pi_0, \pi_1, \xi'', \mu, \lambda$ are pre-defined constants, and $\alpha,W$ are variables)</p>
<p><code>$$
\begin{align}
&\max_{W}\; \pi_{0} (1 - e^{-\mu W}) - \frac{\pi_{1}}{W}z_{1}(W)\\
s.t.\;& \frac{dz_{0}(\alpha)}{d\alpha} = 1 - e^{-\mu \alpha},\quad z_{0}(0) = 0\\
& \frac{dz_{1}(\alpha)}{d\alpha} = 1 - e^{-\lambda \alpha},\quad z_{1}(0) = 0\\
& z_{0}(W) + \epsilon \leq \xi''\\
& W \geq 0
\end{align}
$$</code>
where $z_{0},z_{1}$ are auxiliary variables and $\epsilon$ is a numerical tolerance value. This then becomes a DAE (differential-algebraic equation) optimization problem, which can be solved numerically. (though given that your decision variable is also the independent variable in the differential equations, some further bilinear transformations may be required. See <a href="http://dx.doi.org/10.1016/j.na.2005.03.066" rel="nofollow">http://dx.doi.org/10.1016/j.na.2005.03.066</a>. Some software packages do this automatically.).</p>
<p>Edit: I just realized, if indeed $\xi''$ is a constant as I have assumed, the inequality constraint can easily be converted into a bound.
<code>$$
\begin{align}
\int_{0}^{W^U} 1-e^{-\mu \alpha}\,d\alpha = \xi'' - \epsilon\\
\frac{e^{-\mu W^{U}}}{\mu} + W^{U} - \frac{1}{\mu} = \xi'' - \epsilon\\
\end{align}
$$</code>
Solve for $W^{U}$, and replace the above inequality constraints with:
$$ 0 \leq W \leq W^{U}
$$
You may be able to solve this using optimal control methods.</p>