Відмінності між версіями «Одноетапні стохастичні задачі з лінійними розв’язувальними правилами. М-модель та V-модель.»

Матеріал з Вікі ЦДУ
Перейти до: навігація, пошук
Рядок 22: Рядок 22:
  
 
де <math> \bar{c} </math> і <math>\bar{b} </math> - математичні очікування векторів ''c'' і ''b''.
 
де <math> \bar{c} </math> і <math>\bar{b} </math> - математичні очікування векторів ''c'' і ''b''.
 +
 +
Зведемо тепер умови (1.2) до еквівалентного детермінованого вигляду. Позначимо через <math>\ a_i=(a_{i1},...,a_{in}) </math> ''i''-ту вектор-строку матриці А і введемо випадкову змінну <math>\ \zeta_i </math>:
 +
 +
<math>\ \zeta_i=\frac{(b_i-\bar{b}_i)-a_iD(b-\bar{b})}{\sqrt{M[(b_i-\bar{b}_i)-a_iD(b-\bar{b})]^2}} </math>.
 +
 +
Легко бачити, що
 +
 +
<math>\ \bar{\zeta}_i=M \zeta_i=0  </math> і <math>\ \sigma^2_{\zeta_i}=M(\zeta_i-\bar{\zeta}_i)^2=1  </math>,
 +
 +
тобто випадкові величини <math>\ \zeta_i </math> мають нульові математичні сподівання і одиничні дисперсії.
 +
 +
Неважко впевнитися безпосередніми обчисленнями, що умови  <math>\ \sum^n_{j=1}a_{ij}x_j \le b_i </math>, або, що те ж саме у припущенні (1.3), <math>\ a_iDb \le b_i </math> еквівалентні співвідношенням
 +
 +
<math>\ \zeta_i \ge \zeta^0_i=\frac{-\bar{b}_i+a_iD\bar{b}}{\sqrt{M[(b_i-\bar{b}_i)-a_iD(b-\bar{b})]^2}} </math>

Версія за 09:45, 18 січня 2013

    Розглянемо наступну М-модель стохастичного програмування з ймовірнісними обмеженнями:

Неможливо розібрати вираз (невідома помилка): \ M(cx) \rightarrow max , (1.1)

Неможливо розібрати вираз (невідома помилка): \ P \left \{ \sum^{n}_{j=1} a_{ij}x_j \right \} \geq p_{i}, i=1,...,m . (1.2)

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

припускаються включеними у систему нерівностей (1.2) з Неможливо розібрати вираз (невідома помилка): \ p_{i}=1 

.

Будемо вважати матрицю Неможливо розібрати вираз (невідома помилка): \ A=||a_{ij}||

детермінованою, а вектори b і c незалежними випадковими векторами. Компоненти кожного із цих векторів можуть бути корелбованими між собою. Наступні обчислення будемо вести, припускаючи, що складові векторо обмежень b розподілені нормально.

Задамо розв'язувальне правило у вигляді

Неможливо розібрати вираз (невідома помилка): \ x=Db , (1/3)

де D - невідома детермінована матриця розміру Неможливо розібрати вираз (невідома помилка): \ n \times m .

Знайти розв'язки задачі (1.1)-(1.3) - означає очислити елементи Неможливо розібрати вираз (невідома помилка): \ d_{ij}

матриці D.

Перетворимо запис задачі (1.1)-(1.3). Підставимо (1.3) у виразі для показника якості (1.1) розв'язку задачі. Враховуючи статистичну незалежність векторів c і b, маємо

Неможливо розібрати вираз (невідома помилка): M(cx)=\overline{cx}=\overline{cDb}=\sum^m_{i=1} \sum^n_{j=1}d_{ij}\bar{b_i}\bar{c_i} ,

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

і Неможливо розібрати вираз (невідома помилка): \bar{b} 
- математичні очікування векторів c і b.

Зведемо тепер умови (1.2) до еквівалентного детермінованого вигляду. Позначимо через Неможливо розібрати вираз (невідома помилка): \ a_i=(a_{i1},...,a_{in})

i-ту вектор-строку матриці А і введемо випадкову змінну Неможливо розібрати вираз (невідома помилка): \ \zeta_i 

Неможливо розібрати вираз (невідома помилка): \ \zeta_i=\frac{(b_i-\bar{b}_i)-a_iD(b-\bar{b})}{\sqrt{M[(b_i-\bar{b}_i)-a_iD(b-\bar{b})]^2}} .

Легко бачити, що

Неможливо розібрати вираз (невідома помилка): \ \bar{\zeta}_i=M \zeta_i=0

і Неможливо розібрати вираз (невідома помилка): \ \sigma^2_{\zeta_i}=M(\zeta_i-\bar{\zeta}_i)^2=1  

,

тобто випадкові величини Неможливо розібрати вираз (невідома помилка): \ \zeta_i

мають нульові математичні сподівання і одиничні дисперсії.

Неважко впевнитися безпосередніми обчисленнями, що умови Неможливо розібрати вираз (невідома помилка): \ \sum^n_{j=1}a_{ij}x_j \le b_i , або, що те ж саме у припущенні (1.3), Неможливо розібрати вираз (невідома помилка): \ a_iDb \le b_i

еквівалентні співвідношенням

Неможливо розібрати вираз (невідома помилка): \ \zeta_i \ge \zeta^0_i=\frac{-\bar{b}_i+a_iD\bar{b}}{\sqrt{M[(b_i-\bar{b}_i)-a_iD(b-\bar{b})]^2}}