⎧⎪⎪⎨
⎪⎪⎩
Y
R,i
(
ξ
)
при
V
ξ,i
+
V
ξ,i
+1
>
0
,
Y
R,i
(
ξ
) +
Y
L,i
+1
(
ξ
)
2
при
V
ξ,i
+
V
ξ,i
+1
= 0
,
Y
L,i
+1
(
ξ
)
при
V
ξ,i
+
V
ξ,i
+1
<
0
и
⎧⎪⎪⎨
⎪⎪⎩
Y
R,j
(
η
)
при
V
η,j
+
V
η,j
+1
>
0
,
Y
R,j
(
η
) +
Y
L,j
+1
(
η
)
2
при
V
η,j
+
V
η,j
+1
= 0
,
Y
L,j
+1
(
μ
)
при
V
η,j
+
V
η,j
+1
<
0
.
Далее по тексту будем обозначать расширенный массив так же, как
и раньше —
Y
i,j
.
В каждой ячейке с номером
i
для каждой восстанавливаемой ве-
личины
Y
i,j
осуществляется расчет индекса немонотонности
Ind
(
Y
)
:
Ind
(
Y
) =
1
12
|−
Y
i
+2
,j
+ 16
Y
i
+1
,j
−
30
Y
i,j
+ 16
Y
i
−
1
,j
−
Y
i
−
2
,j
|
1
2
|−
Y
i
+2
,j
+ 4
Y
i
+1
,j
−
3
Y
i,j
|
+
1
2
|
3
Y
i,j
−
4
Y
i
−
1
,j
+
Y
i
−
2
,j
|
+
θ
,
где
θ
— малый параметр.
Далее для центра и границ
i
-й ячейки уточняется индекс немоно-
тонности:
Ind
C
=
(
Ind
(
Y
i
+1
,j
) +
Ind
(
Y
i,j
) +
Ind
(
Y
i
−
1
,j
))
3
;
Ind
L
=
(
Ind
(
Y
i
−
1
,j
) +
Ind
(
Y
i,j
))
2
, Ind
R
=
(
Ind
(
Y
i
+1
,j
) +
Ind
(
Y
i,j
))
2
.
Регулируется нижняя граница ограничителей:
R
min
= (1
−
Ind
C
)
W
max
+
Ind
C
W
min
.
Определяются значения ограничителей
R
(
ξ, Ind
)
, R
(
η, Ind
)
на ле-
вой и правой границах ячейки с номером
i
:
R ξ
=
ξ
i
−
1/2
,j
, Ind
= (1
−
Ind
L
) +
Ind
L
R
min
,
R ξ
=
ξ
i
+1/2
,j
, Ind
= (1
−
Ind
R
) +
Ind
R
R
min
.
Далее проводится минимизация значений ограничителей
R
(
Ind
)
:
R
Lim
(
Ind
) = min
R ξ
=
ξ
i
−
1/2
,j
, Ind , R ξ
=
ξ
i
+1/2
,j
, Ind .
Рассчитывается среднее по
i
-му распределению значение ограни-
чителя
R
Lim
(
Ind
)
R
cp
=
N
i
R
Lim,i
N
.
58 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2009. № 4
