Дві часткові стохастичні моделі з розв'язувальними правилами нульового порядку.

Матеріал з Вікі ЦДУ
Перейти до: навігація, пошук

Розглянемо 2 часткові P-моделі з детермінованими обмеженнями і випадковими коефіцієнтами лінійної форми. Побудуємо відповідні детерміновані еквіваленти.

Нехай потрібно

Неможливо розібрати вираз (невідома помилка): P(cx\leq{k})\to{max}

(1)

за умов

Неможливо розібрати вираз (невідома помилка): Ax\geq{b}

           (2)

Неможливо розібрати вираз (невідома помилка): x\geq{0}

            (3)

Елементи матриці А та компоненти вектора b детерміновані, а компоненти вектора c випадкові.

Припускають, що розв’язок задачі (1)-(3) визначається серед детермінованих векторів.

До схеми виду (1)-(3) зводиться задача планування виробництва при випадкових втратах, пов’язаних з реалізацією різних технологічних способів. Оптимальний план повинен максимізувати ймовірність того, що сумарні затрати не перевищать деякої, заданої вищою організацією, величини.

Модель №1.

В цій моделі випадковий вектор Неможливо розібрати вираз (невідома помилка): ~c=c(\omega)

припускається рівним Неможливо розібрати вираз (невідома помилка): ~c=c_{0}+c_{1}\tau(\omega)

, де Неможливо розібрати вираз (невідома помилка): ~c_{0},c_{1}

– детерміновані вектори,Неможливо розібрати вираз (невідома помилка): ~\tau(\omega)
– випадкова величина.

Припускають також, що гіперплощина Неможливо розібрати вираз (невідома помилка): ~c_{1}x=0

не перетинається з многогранною множиною (2)-(3). Нехай для визначеності в точках множини (2)-(3) Неможливо розібрати вираз (невідома помилка): ~c_{1}x>0


Детермінована задача представляє собою наступну задачу дробово-лінійного програмування:

Неможливо розібрати вираз (невідома помилка): {\frac{k-c_{0}x}{c_{1}x}}\to{max}


Неможливо розібрати вираз (невідома помилка): Ax\geq{b} , Неможливо розібрати вираз (невідома помилка): x\geq{0}


Дійсно за прийнятих припущень:

Неможливо розібрати вираз (невідома помилка): P(cx\leq{k})=P\{\omega:c(\omega)x\leq{k}\}=P\{\omega:\tau(\omega)\leq{\frac{k-c_{0}x}{c_{1}x}}\}


Модель №2.

В цій моделі компоненти Неможливо розібрати вираз (невідома помилка): ~c_j

випадкового вектора c припускаються нормально розподіленими з параметрами Неможливо розібрати вираз (невідома помилка): ~\mu_j,\sigma_j

, тобто Неможливо розібрати вираз (невідома помилка): c_j\in N(\mu_j,\sigma_j) .

Припускають також, що точка x=0 не є планом задачі.

Детермінований еквівалент представляє собою задачу нелінійного, точніше неопуклого, програмування:

Неможливо розібрати вираз (невідома помилка): \Big[\Big(k-\sum^{n}_{j=1} {\mu_j x_j}\Big)\Big(\sum^{n}_{j=1} {\sigma_j^2 x_j^2}\Big)^{-1/2}\Big]\to max


Неможливо розібрати вираз (невідома помилка): Ax\geq{b} , Неможливо розібрати вираз (невідома помилка): x\geq{0}


Дійсно за прийнятих припущень:

Неможливо розібрати вираз (невідома помилка): P(cx\leq{k})=P\Bigg(\sum^{n}_{j=1} {c_j x_j}\leq{k}\Bigg)=P\Bigg\{\frac{\sum^{n}_{j=1} {c_j x_j}-\sum^{n}_{j=1} {\mu_j x_j}}{\sqrt{\sum^{n}_{j=1} {\sigma_j^2 x_j^2}}}\leq{\frac{k-\sum^{n}_{j=1} {\mu_j x_j}}{\sqrt{\sum^{n}_{j=1} {\sigma_j^2 x_j^2}}}}\Bigg\}=P(\eta\leq{\eta_0})


Неможливо розібрати вираз (невідома помилка): \eta=\frac{\sum^{n}_{j=1} {c_j x_j}-\sum^{n}_{j=1} {\mu_j x_j}}{\sqrt{\sum^{n}_{j=1} {\sigma_j^2 x_j^2}}}


Неможливо розібрати вираз (невідома помилка): \eta_{0}=\frac{k-\sum^{n}_{j=1} {\mu_j x_j}}{\sqrt{\sum^{n}_{j=1} {\sigma_j^2 x_j^2}}}


В силу прийнятих припущень випадкова величина Неможливо розібрати вираз (невідома помилка): \eta

нормально розподілена: Неможливо розібрати вираз (невідома помилка): \eta\in N(0,1)

. Тому

Неможливо розібрати вираз (невідома помилка): P(cx\leq{k})=P(\eta\leq{\eta_0})=1/\sqrt{2\pi}\int\limits_{-\infty}^{\eta_0}{e^{-t^2/2}}