logo
vstyp_mpdo

540. 3Ведення розв'язання одноетапної статичної задачі стохастичного програмування до детермінованої задачі лінійного програмування.

Одноетапна задача стохастичного програмування використовується в тому разі, коли рішення приймаються на підставі відомих характеристик розподілу ймовірностей випадкових параметрів умови задачі до спостережень за їхніми реалізаціями. У такому разі має прийматися найкраще в середньостатистичному розумінні рішення. Тобто випадкові параметри задачі замінюють їх середніми величинами і початкову задачу стохастичного програмування зводять до детермінованої.

Розглянемо лінійну одноетапну задачу стохастичного програмування в такій постановці: визначити план Х, для якого max M {Сума(j=1,n) cj(ω)xj}, P{Сума(j=1,n)aij(ω)xi<bi(ω)}>p (i=1,m), xj>0, ωє Ω (j=1,n),

де вектор коефіцієнтів при змінних у цільовій функції, матриця коефіцієнтів при змінних у системі обмежень, а також вектор є випадковими величинами; ω — випадковий параметр, Ω — множина значень ω, що з’являються з певною ймовірністю. Нехай A(ω) — нормально розподілена випадкова величина з математичним сподіванням aij і дисперсією сігма, а B(ω) i C(ω) — нормально розподілені випадкові величини з математичними сподіваннями і дисперсіями.

Оскільки в обмеженнях задачі матриця A(ω) та вектор B(ω) є нормально розподіленими випадковими величинами, то їх різниці також є випадковими величинами з нормальним розподілом, математичним сподіванням і дисперсією.

Обмеження P{Сума(j=1,n)aij(ω)xi<bi(ω)}>p(i=1,m) еквівалентні нерів­ностям P{дельта(X)<0}>p(i=1,m). Враховуючи, що дельта(Х) нормально розподілена випадкова величина, використаємо функцію нормального закону розподілу, внаслідок чого наведену нерівність можна записати так: 1/Корень(2ПИ)сигма(X)*Интеграл exp{(Е-Дельта)2/2сигма2(Х)}dЕ>p, (i=1,m)

Позначимо: Ф(t)=1/ Корень(2ПИ)* Интеграл е(Е2/2) dЕ. Тоді останню нерівність зведемо до вигляду: Ф(-Дельта(X)/сигма(X))>p, звідки Дельта(X) +Ф-1(p)сигма(X)<0.

Підставивши в цю нерівність значення Дельта(X) і сигма (Х), отримаємо: Ф-1(p)Корень (Сума(j=1,n)Сигма2x2+O2)<b-Сума(j=1,n)aijxj

Отже, початкову стохастичну задачу зведено до детермінованого аналогу з лінійною цільовою функцією та нелінійними обмеженнями: max F=сума(j=1,n)cjxj за умов: Ф-1(p)Корень (Сума(j=1,n)Сигма2x2+O2)<b-Сума(j=1,n)aijxj

Таку задачу можна розв’язати одним з відомих методів розв’язування задач нелінійного програмування, наприклад, методом множників Лагранжа.