Компьютерное моделирование и оптимизация технологических процессов и оборудования. Дворецкий С.И - 142 стр.

UptoLike

где
при 0>θ
k
,
θ
θ
θ
θ
=
+
θ
;
~
;
~
k
k
k
k
C
C
C
k
при 0
<
θ
k
,
k
kkk
k
zdCzdC
C
θ
θθ±θ
θ
±
),,(),
~
,(
~
~
**
.
В этом случае коррекция найденных значений
),
~
,(
~
*
θzdC
может быть выполнена по формулам
()
=
θ=
K
k
k
CzdCC
1
*
)(),
~
,(
~
ˆ
,
(4.68)
где
∫∫
+
θ
θθθ
θ
+θθθ
θ
=
0
0
)(
ˆ
)(
ˆ
kkk
k
kkk
k
dP
C
dP
C
C
k
.
(4.69)
Если ошибка в определении
k
θ подчиняется нормальному закону распределения, то интегралы в
последнем выражении можно переписать в виде:
{}
{
}
(
)
{}
θ
θ
θθ
θπ
θ
0
2
exp
2
1
k
k
kk
k
k
d
D
M
D
;
{}
{
}
(
)
{}
θ
θ
θθ
θπ
θ
0
2
exp
2
1
k
k
kk
k
k
d
D
M
D
.
Блок-схема вычисления
CCO
i
I
с использованием (4.67) – (4.68) представлена на рис. 4.8.
В блоке 1 генерируются случайные величины
k
θ
θ
θ
,...,,
21
с использованием соответствующих плот-
ностей распределения. В блоке 2 для каждой реализации
)(ν
θ рассчитывается оптимальный режим
функционирования ХТП
*
z
и оптимальные значения векторной целевой функции.
В блоке 3 последовательно вычисляются квазиоптимальные режимы
*
~
z
для неопределенных пара-
метров
kk
θ±θ . При этом блок 3 работает последовательно
K
2 раз, каждый раз для отдельной незави-
симой группы неопределенных параметров (здесь
K
размерность вектора θ ).
В блоке 4 вычисляются значения компонент вектора
),
~
,(
~
*
θzdC и
k
C
θ
по формулам (4.68), (4.69).
Блоки 5 – 9 предназначены для вычисления вероятностей и проверки выполнения технологических
условий с заданной вероятностью
[
]
ν+ν=
ννν
/
ˆ
)1()()(
)()1(CCO)(CCO
iii
CdIdI ,