8 / 14
Прогнозирование числа отказов программного обеспечения космических аппаратов
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2016. № 5
31
иска искомых параметров определяется шагом перебора параметров
и .
i
i
Для
повышения точности уменьшается шаг. Алгоритм расчетов и исходные данные
приведены на рис. 5 и 6. Зафиксированные значения
k
M
представлены на рисун-
ках в виде матрицы-столбца
N
. Как следует из рис. 5, для
1, 85
минимум
A
= 6,103 достигается при
0,06.
Соответственно для
0, 06
минимум
A
= 4,109 достигается при
1,6
(см. рис. 6). Таким образом, в рассматриваемом
случае минимум критериальной функции (4) достигается при
1,6,
которому
соответствует
0, 06.
Величина
0,0
M
оценивается по изложенному методу
наименьших квадратов для найденных оптимальных значений
0,06
и
= 1,6.
При этом минимум критериальной функции (4) находится методом перебора по
0,0
M
. Результаты расчета для рассмотренного примера приведены на рис. 8.
0.06
1.55
0
r 1
M n( )
M
0 0
1
a
m
1 e
m 1
(
)
m
m 0 n
for
b
s
1 e
s 1
( )
s
s 0 n
for
M
i 1
j
M
i j
1 b
i j
i 1
j
if
0
i
k
M
i k
b
i k
otherwise
j 0 n
for
i 0 n 1
for
M
i 1
j
M
i j
1 a
i j
i 1
j
if
0
r 1
k
M
i k
b
i k
r
i
k
M
i k
a
i k
otherwise
j r n
for
i r n 1
for
M
i j
M
i j
i j
if
0 otherwise
j 0 n
for
i 0 n
for
M
N
1
1
1
3
3
1
2
1
2
1
0.5
A
1
10
k
M 10( )
(
)
1 1
1
M 10( )
k k
N
k
2
A 4.218
A 1.55 (
) 4.218
A 1.6( ) 4.109
A 1.65 (
) 4.15
Рис. 6.
Алгоритм и результаты расчета критериальной функции
A
для
0,06
и различных значений