О гладкой аппроксимации вероятностных критериев в задачах стохастического программирования
Ключевые слова:
стохастическое программирование, вероятностный критерий, квантильный критерий, аппроксимация, численные методы, сигмоидальная функцияАннотация
В данной работе исследуется один из возможных вариантов гладкой аппроксимации вероятностных критериев в задачах стохастического программирования. Исследование проведено в приложении к задачам оптимизации функции вероятности и функции квантили для функционала потерь, зависящего от вектора управления и одномерной абсолютно непрерывной случайной величины. В данной работе исследуется один из возможных вариантов гладкой аппроксимации вероятностных критериев в задачах стохастического программирования. Исследование проведено в приложении к задачам оптимизации функции вероятности и функции квантили для функционала потерь, зависящего от вектора управления и одномерной абсолютно непрерывной случайной величины. Основная идея аппроксимации – замена разрывной функции Хевисайда в интегральном представлении функции вероятности на гладкую функцию, обладающую такими свойствами как непрерывность, гладкость, а также имеющую легко вычислимые производные. Примером такой функции является функция распределения случайной величины, распределенной по логистическому закону с нулевым средним и конечной дисперсией – сигмоида. Величина, обратно пропорциональная корню из дисперсии, при этом является параметром, обеспечивающим близость исходной функции и ее аппроксимации. Такая замена позволяет получить гладкое приближение функции вероятности, для которого легко могут быть найдены производные по вектору управления и иным параметрам задачи. В статье доказана сходимость аппроксимации функции вероятности, полученной при замене функции Хевисайда на сигмоидальную функцию, к исходной функции вероятности, и получена оценка погрешности такой аппроксимации. Далее получены приближенные выражения для производных функции вероятности по вектору управления и параметру функции, доказана их сходимость к истинным производным при выполнении ряда условий на функционал потерь. С помощью известных соотношений между производными функции вероятности и функции квантили получены приближенные выражения для производных функции квантили по вектору управления и уровню вероятности. Рассмотрены примеры, демонстрирующие возможность применения предложенных оценок к решению задач стохастического программирования с критериями в форме функции вероятности и функции квантили, в том числе в случае многомерной случайной величины.
Литература
Лысенко И.В. Нечеткая оптимизация: новый подход к постановке и решению задач // Труды СПИИРАН. 2004. 1. № 2 (1). C. 90--118.
Бабуров В.И., Иванцевич Н.В., Олянюк П.В., Саута О.И. Нечеткие критерии и понятия в задаче выбора приемника ГНСС для навигационно-пилотажного комплекса // Труды СПИИРАН. 2012. 1. № 23 (4). C. 357--368.
Фильченков А.А. Меры истинности и вероятностные графические модели для представления знаний с неопределенностью // Труды СПИИРАН. 2012. 1. № 23 (4). C. 254--295.
Кибзун А.И., Кан Ю.С. Задачи стохастического программирования с вероятностными критериями.М.: Физматлит, 2009. 372 с.
D. Dentcheva, S. Penev, A. Ruszczyński. Statistical estimation of composite risk functionals and risk optimization problems // Annals of the Institute of Statistical Mathematics. 2017, Vol. 69, Issue 4, pp. 737--760.
D. Dentcheva, A. Ruszczyński. Risk forms: representation, disintegration, and application to partially observable two-stage systems // Mathematical Programming. 2019, pp. 1--21.
Кибзун А.И., Наумов А.В., Норкин В.И. О сведении задачи квантильной оптимизации с дискретным распределением к задаче смешанного целочисленного программирования // Автоматика и телемеханика, 2013, № 6, С. 66--86.
Кибзун А.И., Игнатов А.Н. Сведение двухшаговой задачи стохастического оптимального управления с билинейной моделью к задаче смешанного целочисленного линейного программирования // Автоматика и телемеханика, 2016, № 12, С. 89--111.
Иванов С.В., Кибзун А.И. Выборочная аппроксимация двухэтапной задачи стохастического линейного программирования с квантильным критерием // Тр. ИММ УрО РАН, 23, № 3, 2017, С. 134--143.
Кибзун А.И., Соболь В.Р. Двухшаговая задача хеджирования европейского колл-опциона при случайной длительности транзакций // Тр. ИММ УрО РАН, 21, № 3, 2015, С. 164--174.
Райк Э. Дифференцируемость по параметру функции вероятности и стохастический псевдоградиентный метод для ее оптимизации // Изв. АН ЭССР, физ.-мат., 1975, 24, № 1, С. 3--8.
Кибзун А.И., Третьяков Г.Л. Дифференцируемость функции вероятности. // Докл. РАН, 1997, 354, № 2, C. 159--161.
Uryas'ev S.. Derivatives of probability functions and some applications. // Annals of Operations Research. 1995, Vol. 56 Issue 1-4, pp. 287--311.
Uryas'ev S.. Derivatives of probability functions and integrals over sets given by inequalities. // Journal of Computational and Applied Mathematics. 1994, Vol. 56, Issue 1-2, pp. 197--223.
Rene Henrion. Gradient estimates for Gaussian distribution functions: application to probabilistically constrained optimization problems // Numerical Algebra, Control and Optimization, 2012, Vol. 2, 4, pp. 655--668.
Wim van Ackooij, ReneHenrion. (Sub-)Gradient Formulae for Probability Functions of Random Inequality Systems under Gaussian Distribution // SIAM/ASA J. Uncertainty Quantification, 2017, Vol. 5, pp. 63--87.
Georg Ch. Pflug, Heinz Weisshaupt. Probability gradient estimation by set-valued calculus and applications in network design // SIAM Journal on Optimization, 2005, Vol. 15, pp. 898--914.
Yu C, Zelterman D. A general approximation to quantiles. Commun Stat Theory Methods. 2017, 46(19), pp. 9834--9841.
H. Okagbue, M. Adamu, T. Anake. Ordinary Differential Equations of the Probability Functions of the Weibull Distribution and their Application in Ecology. International Journal of Engineering and Future Technology. 2018, 15, pp. 57--78.
J. Garniera, A. Omraneb, Y. Rouchdyc. Asymptotic formulas for the derivatives of probability functions and their Monte Carlo estimations // European Journal of Operational Research, 2009, Vol. 198, Issue 3, pp. 848--858.
A. Mafusalov, S. Uryasev. Buffered Probability of Exceedance: Mathematical Properties and Optimization // SIAM Journal on Optimization, 2018, Vol. 28, 2, pp. 1077--1103.
T. Zhang, S. Uryasev, Y. Guan. Derivatives and subderivatives of buffered probability of exceedance // Operations Research Letters, 2019, Vol. 47, Issue 2, pp 130--132
Ширяев А.Н. Вероятность. Кн. 1. 4-е изд., перераб. и доп. - М., 2007, 552 с.
On smooth approximation of probabilistic criteria in stochastic programming problems
Abstract
In this paper we study one of the possible variants of smooth approximation of probability criteria in stochastic programming problems. The research is applied to the optimization problems of the probability function and the quantile function for the loss functional depending on the control vector and one-dimensional absolutely continuous random variable. In this paper we study one of the possible variants of smooth approximation of probability criteria in stochastic programming problems. The research is applied to the optimization problems of the probability function and the quantile function for the loss functional depending on the control vector and one-dimensional absolutely continuous random variable. The main idea of the approximation is to replace the discontinuous Heaviside function in the integral representation of the probability function with a smooth function having such properties as continuity, smoothness, and easily computable derivatives. An example of such a function is the distribution function of a random variable distributed according to the logistic law with zero mean and finite dispersion, which is a sigmoid. The value inversely proportional to the root of the variance is a parameter that provides the proximity of the original function and its approximation. This replacement allows us to obtain a smooth approximation of the probability function, and for this approximation derivatives by the control vector and by other parameters of the problem can be easily found. The article proves the convergence of the probability function approximation obtained by replacing the Heaviside function with the sigmoidal function to the original probability function, and the error estimate of such approximation is obtained. Next, approximate expressions for the derivatives of the probability function by the control vector and the parameter of the function are obtained, their convergence to the true derivatives is proved under a number of conditions for the loss functional. Using known relations between derivatives of probability functions and quantile functions, approximate expressions for derivatives of quantile function by control vector and by the level of probability are obtained. Examples are considered to demonstrate the possibility of applying the proposed estimates to the solution of stochastic programming problems with criteria in the form of a probability function and a quantile function, including in the case of a multidimensional random variable.
Опубликован
Как цитировать
Раздел
Copyright (c) 2020 Виталий Романович Соболь, Роман Олегович Торишный

Это произведение доступно по лицензии Creative Commons «Attribution» («Атрибуция») 4.0 Всемирная.
Авторы, которые публикуются в данном журнале, соглашаются со следующими условиями: Авторы сохраняют за собой авторские права на работу и передают журналу право первой публикации вместе с работой, одновременно лицензируя ее на условиях Creative Commons Attribution License, которая позволяет другим распространять данную работу с обязательным указанием авторства данной работы и ссылкой на оригинальную публикацию в этом журнале. Авторы сохраняют право заключать отдельные, дополнительные контрактные соглашения на неэксклюзивное распространение версии работы, опубликованной этим журналом (например, разместить ее в университетском хранилище или опубликовать ее в книге), со ссылкой на оригинальную публикацию в этом журнале. Авторам разрешается размещать их работу в сети Интернет (например, в университетском хранилище или на их персональном веб-сайте) до и во время процесса рассмотрения ее данным журналом, так как это может привести к продуктивному обсуждению, а также к большему количеству ссылок на данную опубликованную работу (Смотри The Effect of Open Access).