ГОСТ Р 50779.26-2007 Статистические методы. Точечные оценки, доверительные, предикционные и толерантные интервалы для экспоненциального распределения (МЭК 60605-4:2001).

ГОСТ Р 50779.26-2007

Группа Т59

НАЦИОНАЛЬНЫЙ СТАНДАРТ РОССИЙСКОЙ ФЕДЕРАЦИИ

 Статистические методы

 ТОЧЕЧНЫЕ ОЦЕНКИ, ДОВЕРИТЕЛЬНЫЕ, ПРЕДИКЦИОННЫЕ

И ТОЛЕРАНТНЫЕ ИНТЕРВАЛЫ ДЛЯ ЭКСПОНЕНЦИАЛЬНОГО РАСПРЕДЕЛЕНИЯ

 Statistical methods. Point estimates, confidence intervals, prediction intervals

 and tolerance intervals for exponential distribution

ОКС 03.120.30

Дата введения 2008-06-01

Предисловие

Цели и принципы стандартизации в Российской Федерации установлены Федеральным законом от 27 декабря 2002 г. N 184-ФЗ "О техническом регулировании", а правила применения национальных стандартов Российской Федерации - ГОСТ Р 1.0-2004 "Стандартизация в Российской Федерации. Основные положения"

Сведения о стандарте

1 ПОДГОТОВЛЕН Открытым акционерным обществом "Научно-исследовательский центр контроля и диагностики технических систем" (ОАО "НИЦ КД") и Техническим комитетом по стандартизации ТК 125 "Статистические методы в управлении качеством продукции" на основе собственного аутентичного перевода стандарта, указанного в пункте 4

2 ВНЕСЕН Управлением развития, информационного обеспечения и аккредитации Федерального агентства по техническому регулированию и метрологии

3 УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Приказом Федерального агентства по техническому регулированию и метрологии от 27 декабря 2007 г. N 578-ст

4 Настоящий стандарт является модифицированным по отношению к международному стандарту МЭК 60605-4:2001 "Испытания оборудования на надежность. Часть 4. Статистические процедуры для экспоненциального распределения. Точечные оценки, доверительные, предикционные и толерантные интервалы" (IEC 61605-4:2001 "Equipment reliability testing - Part 4: Statistical procedures for exponential distribution - Point estimates, confidence intervals, prediction intervals and tolerance intervals"), путем внесения технических отклонений, объяснение которых приведено во введении к настоящему стандарту.

Наименование настоящего стандарта изменено относительно наименования указанного международного стандарта для приведения в соответствие с ГОСТ Р 1.5-2004 (подраздел 3.5)

5 ВВЕДЕН ВПЕРВЫЕ

Информация об изменениях к настоящему стандарту публикуется в ежегодно издаваемом информационном указателе "Национальные стандарты", а текст изменений и поправок - в ежемесячно издаваемых информационных указателях "Национальные стандарты". В случае пересмотра (замены) или отмены настоящего стандарта соответствующее уведомление будет опубликовано в ежемесячно издаваемом информационном указателе "Национальные стандарты". Соответствующая информация, уведомление и тексты размещаются также в информационной системе общего пользования - на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет

 Введение

В отличие от применяемого международного стандарта в настоящий стандарт не включены ссылки на МЭК 60050-191:1990 "Международный электротехнический словарь. Глава 191. Надежность и качество услуг", который нецелесообразно применять в национальном стандарте из-за отсутствия принятого гармонизированного национального стандарта. В соответствии с этим изменено содержание разделов 1 и 3. Кроме того, изменено содержание раздела 5, добавленные предложения и ссылки на национальные стандарты выделены в тексте стандарта курсивом.

      1 Область применения

Настоящий стандарт устанавливает статистические методы определения точечных оценок, доверительных, предикционных и толерантных интервалов для интенсивности отказов объектов, наработки которых на отказ подчиняются экспоненциальному распределению. Это означает, что интенсивность отказов постоянна во времени. Несмотря на то, что в стандарте рассмотрена интенсивность отказов, приведенные количественные методы применимы также к интенсивности других событий, наработки до появления которых подчиняются экспоненциальному распределению. Таким образом, приведенные в настоящем стандарте методы применимы, например, к постоянным параметрам потока отказов, интенсивности восстановлений. Однако для удобства и ненужных повторений дальнейшее изложение сделано только для отказов и интенсивности отказов.

Использование методов настоящего стандарта должно сопровождаться проверкой предположений о постоянстве интенсивности отказов или параметра потока отказов (см. ГОСТ Р МЭК 60605-6).

Настоящий стандарт применим также в случаях, когда объекты случайной выборки подвергают испытаниям с регистрацией наработок до отказа для оценки показателя надежности.

      2 Нормативные ссылки

В настоящем стандарте использованы нормативные ссылки на следующие стандарты:

ГОСТ Р 50779.10-2000 (ИСО 3534.1:1993) Статистические методы. Вероятность и основы статистики. Термины и определения (ИСО 3534-1:2006 "Статистика. Словарь и условные обозначения. Часть 1. Общие статистические термины и термины, используемые в вероятностных задачах", IDT)

ГОСТ Р 50779.11-2000 (ИСО 3534.2:1993) Статистические методы. Статистическое управление качеством. Термины и определения (ИСО 3534-2:2006 "Статистика. Словарь и условные обозначения. Часть 1. Прикладная статистика", IDТ)

ГОСТ Р МЭК 60605-6-2007 Надежность в технике. Критерии проверки постоянства интенсивности отказов и параметра потока отказов (МЭК 60605-6:1997 "Испытания оборудования на надежность. Часть 6. Критерии проверки постоянства интенсивности отказов или параметра потока отказов", IDТ)

Примечание - При пользовании настоящим стандартом целесообразно проверить действие ссылочных стандартов в информационной системе общего пользования - на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет или по ежегодно издаваемому информационному указателю "Национальные стандарты", который опубликован по состоянию на 1 января текущего года, и по соответствующим ежемесячно издаваемым информационным указателям, опубликованным в текущем году. Если ссылочный стандарт заменен (изменен), то при пользовании настоящим стандартом следует руководствоваться заменяющим (измененным) стандартом. Если ссылочный стандарт отменен без замены, то положение, в котором дана ссылка на него, применяется в части, не затрагивающей эту ссылку.

      3 Термины, определения и обозначения

3.1 Термины и определения

В настоящем стандарте применены термины по ГОСТ Р 50779.10, ГОСТ Р 50779.11 и ГОСТ Р ИСО 16269-8.

3.1.1

двусторонний доверительный интервал (two-sided confidence interval): Если и - две функции наблюдаемых значений таких, что для оценки параметра распределения совокупности вероятность равна (где - константа положительная и меньше 1), интервал между и - это двусторонний доверительный интервал для при доверительной вероятности .   [ГОСТ Р 50779.10, статья 2.57]

3.1.2

односторонний доверительный интервал (one-sided confidence interval): Если (или ) - функция наблюдаемых значений такая, что для оценки параметра распределения совокупности вероятность (или вероятность ) больше или равна (где - константа положительная и меньше 1), интервал от наименьшего возможного значения до (или от до наибольшего возможного значения ) - это односторонний интервал для при доверительной вероятности .   [ГОСТ Р 50779.10, статья 2.58]     Примечание - см. приложение В.

3.1.3

толерантные границы (tolerance limits): Две границы и , относительно которых можно утверждать, что по крайней мере заданная доля совокупности лежит между этими границами с вероятностью .   [ГОСТ Р 50779.11, статья 1.4.3]     Примечание - Интервал от и называют толерантным интервалом.

3.1.4

3.1.5

3.2 Обозначения

В настоящем стандарте приведены следующие обозначения:

	- нижняя или верхняя односторонняя доверительная граница показателя надежности, такого как MTTF (средняя наработка до отказа), интенсивность отказов, вероятность безотказной работы и т.д.;
	- уровень значимости - % это уровень доверия, для которого вычисляют доверительные интервалы и границы;
,	- нижняя и верхняя границы односторонних толерантных интервалов уровня %;
,	- нижняя и верхняя границы двустороннего толерантного интервала уровня %;
	- истинная постоянная интенсивность отказов;
	- оценка постоянной интенсивности отказов;
,	- нижняя и верхняя границы односторонних доверительных интервалов для истинной интенсивности отказов;
,	- нижняя и верхняя границы двустороннего доверительного интервала для истинной интенсивности отказов;
	- истинная средняя наработка до отказа;
	- оценка средней наработки до отказа;
,	- нижняя и верхняя границы односторонних доверительных интервалов для истинной средней наработки до отказа;
,	- нижняя и верхняя границы двустороннего доверительного интервала для истинной средней наработки до отказа;
	- общее количество объектов испытаний;
	- доля будущей совокупности (используется при определении толерантных интервалов);
	- количество обнаруженных отказов;
,	- нижняя и верхняя границы односторонних предикционных интервалов уровня %;
,	- нижняя и верхняя границы двустороннего предикционного интервала уровня %;
	- календарное время окончания испытаний;
	- суммарное время испытаний - общая наработка всех объектов до завершения испытаний (см. приложение С);
	- продолжительность будущего периода (используется при определении предикционных интервалов);
	- продолжительность прошедшего периода (используется при определении предикционных интервалов);
	- квантиль -распределения уровня с и степенями свободы;
	- сумма первых ( ) членов распределения Пуассона со средним равным ;
	- вероятность безотказной работы для экспоненциального распределения и заданного времени ;
	- квантиль -распределения уровня с степенями свободы.

      4 Требования и предположения

Статистические методы, описанные в настоящем стандарте, справедливы только в случае, когда интенсивность отказов рассматриваемых объектов (см. раздел 1) постоянна во времени. Если это требование выполняется, средняя наработка до отказа совокупности, из которой отобраны объекты, является обратной величиной к интенсивности отказов. Методы проверки гипотезы о постоянстве интенсивности отказов приведены в ГОСТ Р МЭК 60605-6.

4.1 Предположения и информация, необходимая для определения точечной оценки и доверительных интервалов

Предполагается, что из генеральной совокупности с постоянной интенсивностью отказов отобрана случайным образом выборка из

объектов, которые поставлены на испытания одновременно или в различные моменты времени. Условия испытаний должны быть одинаковыми для всех объектов, подвергаемых испытаниям, а отказавшие объекты могут заменяться или не заменяться. Испытания могут быть прекращены при достижении заданного времени испытаний

(испытания с ограниченной наработкой) или когда произошло заданное количество отказов

(испытания с ограниченным количеством отказов). В вычислениях используют наработку до отказа каждого отказавшего объекта, а также наработку неотказавших объектов к моменту завершения испытаний.

Примечание - Испытания с ограниченным количеством отказов и с ограниченной наработкой часто упоминаются в литературе как испытания с цензурированием типа I и типа II соответственно.

Используемые данные включают наработки до отказа объектов, поставленных на испытания, и объектов, которыми заменены отказавшие.

Необходимая информация (см. приложения А и В):

- уровень доверия (

);

- суммарная наработка

, наблюдаемая за время испытаний (см. приложение С);

- количество отказов

, обнаруженных за время испытаний;

- признак завершения испытаний (при достижении заданного количества отказов или при достижении заданной наработки).

4.2 Требования и предположения для определения предикционных интервалов

В случае предикционных интервалов справедливы перечисленные выше предположения. Кроме того, предполагается известным количество наработок, попадающих в интервал времени

. Необходимо определить предикционный интервал для количества наработок, которые попадут в будущий период времени

.

Таким образом, необходима следующая информация:

- уровень доверия, соответствующий предикционному интервалу;

- количество отказов

, обнаруженных в течение периода

;

- значения

и

.

4.3 Требования и предположения для определения толерантных интервалов

Предположения для толерантных интервалов аналогичны предположениям для предикционных интервалов. Необходимая информация:

- уровень доверия, для которого определяется толерантный интервал;

- доля

объектов совокупности, для которой прогнозируют отказ;

- суммарная наработка за время испытаний

: общая наработка объектов на момент завершения испытаний (см. приложение С)

;

- количество отказов

, произошедших в течение периода

;

- признак завершения испытаний (количество отказов или значение наработки);

- значения

и

.

________________

В случае экспоненциального распределения наработок до отказа не существует различий, каким образом получена наработка

: в результате длительных испытаний одного объекта или менее продолжительных испытаний нескольких объектов. Однако на практике лучше избегать крайностей, т.к. условия экспоненциальной модели не могут сохраняться слишком долго. Эти утверждения распространяются на весь раздел 4.

      5 Метод вычисления точечных оценок и доверительных интервалов

В данном разделе выделены следующие подразделы:

- испытания с ограниченной наработкой;

- испытания с ограниченным количеством отказов;

- испытания, в которых отказавшие объекты не заменяют;

- испытания, в которых отказавшие объекты заменяют;

- испытания для определения точечных оценок;

- испытания для определения односторонних доверительных границ;

- испытания для определения двусторонних доверительных границ.

Структура раздела приведена на рисунке 1.

Рисунок 1 - Структура раздела 5

5.1 Испытания с ограниченной наработкой

5.1.1 Точечные оценки

Для получения точечной оценки необходимо выполнить следующую процедуру (справедливую для испытаний, как с заменой, так и без замены отказавшего объекта).

Фиксируют количество отказов

и суммарную наработку за время испытаний

. Для получения точечной оценки

(интенсивность отказов) используют следующую формулу:

.                                                                   (1)

Для получения точечной оценки

(MTTF) используют следующую формулу:

.                                                                   (2)

Примечание - Точечная оценка, полученная в соответствии с уравнением (2), является смещенной. Однако для небольших значений

(менее 10) это смещение можно уменьшить, заменяя

на (

). Для больших значений

смещение является приемлемым.

Если в процессе испытаний отказы не произошли, точечная оценка для MTTF не может быть получена. Однако могут быть получены нижняя доверительная граница для MTTF, нижняя граница одностороннего предикционного интервала и верхняя граница одностороннего доверительного интервала для интенсивности отказов.

5.1.2 Доверительные границы

Для получения границ доверительного интервала по результатам испытаний с ограниченной наработкой необходимо знать, проводились ли замены отказавших объектов в процессе испытаний (испытания "с заменой") или они не проводились (испытания "без замены").

5.1.2.1 Испытания с заменой

Ниже приведены расчетные формулы для определения границ доверительного интервала.

а)

Верхняя граница одностороннего доверительного интервала для

и нижняя граница одностороннего доверительного интервала для

Для выбранного уровня доверия

% и значений

и

вычисляют верхнюю границу одностороннего доверительного интервала для интенсивности отказов

:

,                                         (3)

где

- уровень доверия (например 0,90 или 90%). Значения

определяют по таблицам квантилей

-распределения (см. приложение D). Аналогично определяют нижнюю границу одностороннего доверительного интервала средней наработки до отказа

:

,                                         (4)

б)

Нижняя граница одностороннего доверительного интервала для

и верхняя граница одностороннего доверительного интервала для

Нижнюю границу одностороннего доверительного интервала

и верхнюю границу одностороннего доверительного интервала

вычисляют по следующим формулам:

,                                           (5)

.                                           (6)

в)

Границы двустороннего доверительного интервала

и

Нижнюю границу

и верхнюю границу

интенсивности отказов определяют по следующим формулам:

,                                             (7)

.                                           (8)

Если в процессе испытаний отказы не произошли, то можно вычислить только верхнюю границу одностороннего доверительного интервала для

. Нижнюю границу

и верхнюю границу

средней наработки до отказа (MTTF) вычисляют по следующим формулам:

,                                           (9)

.                                                      (10)

Если в процессе испытаний отказы не произошли, то можно определить только нижнюю границу одностороннего доверительного интервала для

.

5.1.2.2 Испытания без замены

Ниже приведены расчетные формулы для определения границ доверительного интервала.

а)

Верхняя граница одностороннего доверительного интервала

Для выбранного уровня доверия

% и значений

и

вычисляют верхнюю границу доверительного интервала для интенсивности отказов

:

,                                           (11)

где

- уровень доверия (например 0,90 или 90%). Значения

вычисляют по таблице D.1 (приложение D).

Аналогично вычисляют нижнюю границу одностороннего доверительного интервала для средней наработки

:

.                                           (12)

б)

Нижняя граница одностороннего доверительного интервала

Вычисляют нижнюю границу

доверительного интервала для интенсивности отказов с уровнем доверия

%:

.                                           (13)

Аналогично определяют верхнюю границу одностороннего доверительного интервала

:

.                                           (14)

в)

Границы двустороннего доверительного интервала -

и

.

Вычисляют нижнюю

и верхнюю

границы доверительного интервала для интенсивности отказов с уровнем доверия

%:

,                                           (15)

.                                           (16)

Аналогично вычисляют нижнюю

и верхнюю

границы для истинной средней наработки до отказа:

,                                           (17)

.                                           (18)

г) Метод для неизвестных наработок

Метод предполагает, что испытания проводились "без замены", и использует следующую исходную информацию:

- количество объектов

, первоначально поставленных на испытания;

- количество отказов

, зафиксированных в процессе испытаний;

- продолжительность испытаний

.

Метод включает следующие шаги:

Шаг 1. Вычисляют верхнюю (

) и нижнюю (

) границы двустороннего доверительного интервала для вероятности безотказной работы:

,                    (19)

.                    (20)

Значения

определяют по таблице F.1.

Шаг 2. Вычисляют верхнюю границу двустороннего доверительного интервала для MTTF (

), подставляя значение

в полученное следующее выражение:

.                                         (21)

Шаг 3. Вычисляют нижнюю границу двустороннего доверительного интервала для MTTF (

), подставляя значение

, в полученное следующее выражение:

.                                         (22)

Примечания

1 Уравнения (19) и (20) следуют непосредственно из выражения

.                                                     (23)

2 Верхнюю и нижнюю границы

одностороннего доверительного интервала определяют, заменяя

в уравнениях (19) и (20) на

. Такие же величины используют для определения соответствующих значений

и

.

5.2 Аналитический метод. Испытания с ограниченным количеством отказов

5.2.1 Точечные оценки

Точечную оценку интенсивности отказов

вычисляют по формуле

.                                                 (24)

Точечную оценку MTTF вычисляют по формуле

.                                                 (25)

Примечание - Точечная оценка, определяемая в соответствии с уравнением (22), является смещенной. Однако для небольших значений

(менее 10) это смещение может быть уменьшено заменой

на

. Для больших значений

смещение является приемлемым.

Если в процессе испытаний отказы не обнаружены, точечная оценка интенсивности отказов не может быть получена, поскольку в этом случае концепция испытаний с заданным количеством отказов теряет смысл.

5.2.2 Доверительные интервалы

При определении доверительных границ по результатам испытаний с ограниченным количеством отказов не нужно знать, заменялись или нет отказавшие объекты в процессе испытаний.

Таким образом, метод применим как для испытаний с заменой, так и для испытаний без замены.

Суммарную наработку за время испытаний

определяют в соответствии с приложением С. Метод определения

не зависит от того, проводились ли замены отказавших объектов в процессе испытаний.

5.2.2.1 Верхняя граница одностороннего доверительного интервала

Для заданного уровня доверия и полученных значений величин

и

вычисляют верхнюю границу интенсивности отказов

:

,                                        (26)

где (

) - уровень доверия (например 0,90 или 90%). Значение

определяют по таблице квантилей

-распределения (см. приложение D). Аналогично определяют нижнюю границу одностороннего доверительного интервала для средней наработки до отказа

:

.                                        (27)

5.2.2.2 Нижняя граница одностороннего доверительного интервала

Для заданного уровня доверия и полученных значений величин

и

вычисляют нижнюю границу интенсивности отказов

,                                        (28)

где (

) - уровень доверия (например 0,90 или 90%). Значение

определяют по таблице D.1 (приложение D). Аналогично вычисляют верхнюю границу одностороннего доверительного интервала для средней наработки до отказа

:

.                                        (29)

5.2.2.3 Границы двустороннего доверительного интервала

Вычисляют нижнюю и верхнюю границы доверительного интервала для

:

,                                        (30)

.                                        (31)

Аналогично вычисляют границы доверительного интервала

и

для средней наработки до отказа:

,                                        (32)

.                                        (33)

      6 Предикционный интервал для количества отказов за будущий период

6.1 Границы предикционного интервала

и

Нижнюю

и верхнюю

границы предикционного интервала, соответствующего уровню доверия

% и

наблюдаемых отказов (в прошлом), определяют в соответствии со следующей процедурой:

Шаг 1. Граница

является наименьшим целым числом, удовлетворяющим неравенству (см. [1]):

.                          (34)

Шаг 2. Граница

является наименьшим целым числом, удовлетворяющим неравенству

.                          (35)

Процедура поиска необходимого целого числа, удовлетворяющего уравнению (34) или (35), может быть компьютеризирована. Таблица квантилей

-распределения приведена в приложении F.

6.2 Границы односторонних предикционных интервалов

Нижнюю и верхнюю границу соответствующих предикционных интервалов получают заменой

на

в уравнениях (34) или (35).

      7 Процедура определения границ толерантного интервала

7.1 Верхняя толерантная граница Пуассона

Верхняя граница вероятности для количества отказов

в

будущих периодов продолжительностью

(или

систем в период продолжительностью

является наименьшим целым числом

(обозначаемым

), удовлетворяющим неравенству

,                        (36)

где

- вероятность того, что количество отказов не будет превышать

;

- функция распределения Пуассона с параметром

.

Примечание - Величина

связана с

-распределением

,

где

,

.

Значения величины

приведены в таблице Е.1 (приложение Е).

Здесь предполагается, что значение

неизвестно, в этом случае

определяют, используя следующую процедуру.

Шаг 1. Определяют верхнюю границу (

) одностороннего доверительного интервала для

с уровнем доверия

%, используя формулы (3), (11) или (26) соответственно.

Шаг 2. Подставляют значение

в неравенство (36) и находят наименьшее целое число

, удовлетворяющее неравенству. Это целое число

является искомой верхней толерантной границей

.

Таким образом, изготовитель может утверждать с уровнем доверия

%, что "по крайней мере,

из

будущих периодов будут содержать не более

отказов системы" или напротив "по меньшей мере

систем будет иметь не более

отказов в период продолжи

тельности

".

7.2 Нижняя толерантная граница Пуассона

Нижняя граница вероятности для количества отказов

в

будущих периодов продолжительности

(или

систем в период продолжительности

) является наибольшим целым числом

(обозначаемым

), удовлетворяющим неравенству

,                        (37)

где

- вероятность того, что количество отказов будет не менее

;

- функция распределения Пуассона с параметром

.

Предполагая, что значение

неизвестно,

определяют, используя следующую процедуру.

Шаг 1. Определяют нижнюю границу (

) одностороннего доверительного интервала для

с уровнем доверия

% в соответствии с формулами (5), (13) или (28).

Шаг 2. Подставляют это значение

в неравенство (37) и находят наибольшее целое число

, удовлетворяющее неравенству. Это целое число

является искомой нижней толерантной границей

.

Таким образом, изготовитель может утверждать с уровнем доверия

%, что "по крайней мере,

из

будущих периодов будет содержать не менее

отказов системы" или напротив "по меньшей мере

систем будет иметь не менее

отказов в период продолжительности

".

Приложение А

(справочное)

Примеры

Изготовитель хочет определить показатель надежности прибора, 3308 таких приборов находятся в эксплуатации в течение года. Для них зафиксировано 11 отказов. Предполагается, что:

- интенсивность отказов каждого прибора постоянна во времени;

- в кратчайшее время каждый отказавший прибор заменяют другим таким же работоспособным прибором;

- испытания завершаются при достижении заданной суммарной наработки.

А.1 Точечная оценка MTTF

В соответствии с уравнением (2)

А.2 Определение нижней границы одностороннего доверительного интервала с уровнем доверия 90% для средней наработки до отказа (MTTF)

В соответствии с уравнением (4)

.

Таким образом,

.

Кроме того, в соответствии с уравнением (23)

.

Таким образом, для заданной наработки, например для 10 лет,

.

А.3 Определение границ двустороннего доверительного интервала с уровнем доверия 90% для MTTF

В соответствии с уравнениями (8) и (9) нижняя граница двустороннего доверительного интервала для MTTF

.

Аналогично верхняя граница двустороннего доверительного интервала

.

Таким образом, нижняя и верхняя границы одностороннего доверительного интервала для истинного, но неизвестного MTTF - 182 года и 536 лет соответственно.

А.4 Определение двустороннего предикционного интервала с уровнем доверия 90%

Определим предикционный интервал для предстоящего года. Это означает, что

11,

1 и

1.

При использовании уравнений (34) и (35) можно показать, что

4 и

22.

Таким образом, можно утверждать с уровнем доверия 90%, что количество отказов в течение предстоящего года будет от 4 до 22.

А.5 Определение верхней границы 90%-ного толерантного интервала с уровнем доверия 95%

В соответствии с уравнением (3) верхняя граница одностороннего доверительного интервала с уровнем доверия 95% для интенсивности отказов

.

Таким образом,

(отказ/объект/лет). Поэтому верхняя граница одностороннего доверительного интервала с уровнем доверия 95% для математического ожидания количества отказов

за следующий год в 3308 раз больше этого значения, т.е. равна 18,2. Используя уравнение (36) и таблицу Е.1 (приложение Е), получаем, что наименьшим целым числом

, удовлетворяющим неравенству

, является

24.

Следовательно, верхняя граница одностороннего 90%-ного толерантного интервала равна 24 (

24).

Таким образом, можно утверждать с 95%-ным уровнем доверия, что, по крайней мере, в 90% будущих периодов продолжительностью в один год произойдет не более 24 отказов.

А.6 Определение нижней границы 90%-ного толерантного интервала с уровнем доверия 95%

В соответствии с уравнением (5) нижняя граница одностороннего доверительного интервала для интенсивности отказов с уровнем доверия 95%

.

Таким образом,

отказов на объект в год. Поэтому нижняя граница одностороннего доверительного интервала с уровнем доверия 95% для математического ожидания количества отказов

за следующий год в 3308 раз больше этого значения, т.е. равна 12,3/2=6,15. В соответствии с уравнением (37) и таблицей Е.1 (приложение Е) из этого следует, что наибольшее целое число

, удовлетворяющее неравенству

, определяется из уравнения

2.

Таким образом, нижняя граница 90%-ного толерантного интервала равна 3 (

3). Поэтому можно утверждать с уровнем доверия 95%, что, по крайней мере, в 90% будущих периодов продолжительностью в один год произойдет не более трех отказов.

Приложение В

(справочное)

Взаимосвязь доверительных, предикционных и толерантных интервалов

Введение

Основываясь на предварительном обсуждении значения истинной, но неизвестной средней наработки на отказ совокупности

, испытания необходимо спланировать так, чтобы суммарная наработка за время испытаний

* была достаточно большой относительно

(по крайней мере, в три раза больше

). Не следует излишне доверять испытаниям небольшого количества объектов, поскольку они могут быть не представительными для совокупности. Для дополнительной информации см. [1].

В.1 Доверительные интервалы

По сравнению с предикционными и толерантными интервалами доверительный интервал имеет дело со средним значением совокупности, хотя не всегда правильно описан в литературе. Описание принципа определения доверительного интервала в настоящем стандарте приведено на основе двустороннего доверительного интервала с уровнем доверия 90% для неизвестного среднего. Источником информации являются ресурсные испытания в соответствии с процедурой, описанной в настоящем стандарте. Исходные данные - выборка из

объектов.

Основанная на результатах испытаний статистическая процедура сводится к получению двух значений

и

. Первая величина

- это нижняя граница двустороннего доверительного интервала с уровнем доверия 90% для показателя надежности (например, MTTF, интенсивность отказов, коэффициент готовности), а последняя

- соответствующая верхняя граница этого интервала. Это означает, что если статистическая процедура определения границ доверительного интервала с уровнем доверия 90% будет повторена большое количество раз, то не менее 90% интервалов с полученными парами границ будут накрывать истинный, но неизвестный показатель надежности, а не накроет не более 10% интервалов. Таким образом, уровень доверия связан только с процедурой, используемой для построения интервала.

Из этого следует, что нельзя утверждать, что вероятность того, что среднее совокупности лежит между величинами

и

, равна 90%. Как только получены числовые значения для

и

, вероятность того, что эти границы содержат истинное среднее совокупности, равна нулю или единице. Это понятно, поскольку выборочное среднее является случайным, и таким образом, относительно него можно делать вероятностные утверждения. Однако среднее совокупности - постоянная величина и, следовательно, относительно него не может быть сделано никаких вероятностных утверждений.

В.2 Предикционные интервалы

Много типов предикционных интервалов описано в литературе. Например, существует предикционный интервал, содержащий единственное будущее событие, интервал, содержащий все события

или интервал, который содержит

из

событий. Существуют также интервалы для среднего (математического ожидания) или стандартного отклонения будущей выборки из

событий. Очевидно, что тип предикционного интервала зависит от вида исследуемой задачи.

В настоящем стандарте (см. раздел 6) рассматривается только предикционный интервал, для количества появления событий (например отказов) в заданном будущем промежутке времени, построенный на основе количества появлений этого события в предыдущем промежутке времени.

Особенностью всех типов предикционного интервала является то, что каждому интервалу соответствует уровень доверия, который также относится к процедуре построения предикционного интервала.

В.3 Толерантные интервалы

Предикционные интервалы, описанные в В.2, представляют интерес главным образом для изготовителей, которые хотят спрогнозировать выполнение требований к показателю надежности одного или нескольких будущих объектов. Изготовители, которые хотят сделать выводы о будущих показателях надежности большого количества будущих единиц продукции на основе данных случайным образом отобранной выборки из исследуемой совокупности, часто предикционные интервалы не используют. Такая ситуация возникает, например, если необходимо сделать выводы обо всем процессе производства. Для этой и других аналогичных задач более подходящим является использование прошлого опыта для выводов о будущей продукции. Это приводит к понятию толерантного интервала, который с заданным уровнем доверия будет содержать долю

будущей совокупности.

Например, изготовитель хочет построить интервал, который содержит 95% устройств, имеющих определенный параметр, с уровнем доверия 90%. Необходимо помнить, что для правильного применения статистических процедур выборочные данные (прошлый опыт) должны в действительности быть отобраны случайным образом из исследуемой совокупности. Это, конечно, создает некоторые трудности, когда большая часть совокупности состоит из "будущих" единиц продукции, а вся выборка состоит из "прошлых" единиц продукции. Одно событие, нарушающее принцип случайного выбора, может привести к возникновению дрейфа некоторого параметра, влияющего на будущую продукцию.

Необходимо заметить, что в толерантных интервалах используются две процентные величины, тогда как для доверительных и предикционных интервалов - только одна. Это не должно вызывать затруднений, так как одна из них (90%) относится к проценту исследуемой совокупности, а другая (95%) - к уровню доверия.

Приложение С

(обязательное)

Вычисление суммарной наработки за время испытаний

Рисунки С.1-С.3 поясняют три общих случая вычисления суммарной наработки за время испытаний для каждого вида плана испытаний.

С.1 Вариант 1. Один восстанавливаемый объект с постоянным параметром потока отказов

Суммарной наработкой за время испытаний одного восстанавливаемого объекта с постоянным параметром потока отказов является полное время его работы (исключая ремонт и другие периоды простоя); см. рисунок С.1.

- работоспособное состояние;

- простой, отказ, ремонт

Рисунок С.1 - Пример вычисления

для одного восстанавливаемого объекта

С.2 Вариант 2. Несколько восстанавливаемых объектов с одинаковыми постоянными параметрами потока отказов

Суммарная наработка за время испытаний нескольких восстанавливаемых объектов с одинаковыми постоянными параметрами потока отказов равна сумме всех полных наработок испытываемых объектов (исключая время ремонта и простоя); см. рисунок С.2.

Примечание -

- отказ;

5;

;

6

Рисунок С.2 - Пример вычисления

для пяти восстанавливаемых объектов (мгновенное восстановление)

Предполагается, что все объекты являются идентичными и работают в одинаковых условиях (окружающая среда и нагрузка).

С.3 Вариант 3. Невосстанавливаемые объекты

Для невосстанавливаемых объектов и каждого объекта существует единственный период испытаний (срок службы), равный наработке до первого события (отказа). Суммарная наработка за время испытаний в этом и только в этом случае, является суммой этих периодов (см. рисунок С.3).

Примечание -

- отказ;

- завершение наблюдений;

Рисунок С.3 - Пример вычисления

для пяти невосстанавливаемых объектов

Приложение D

(обязательное)

Таблица квантилей

-распределения

Таблица D.1 - Квантили распределения

	0,050	0,100	0,900	0,950
2	0,10	0,21	4,61	5,99
4	0,71	1,06	7,78	9,49
6	1,64	2,20	10,64	12,59
8	2,73	3,49	13,36	15,51
10	3,94	4,87	15,99	18,31
12	5,23	6,30	18,55	21,03
14	6,57	7,79	21,06	23,68
16	7,96	9,31	23,54	26,30
18	9,39	10,86	25,99	28,87
20	10,85	12,44	28,41	31,41
22	12,34	14,04	30,81	33,92
24	13,85	15,66	33,20	36,42
26	15,38	17,29	35,56	38,89
28	16,93	18,94	37,92	41,34
30	18,49	20,60	40,26	43,77
32	20,07	22,27	42,58	46,19
34	21,66	23,95	44,90	48,60
36	23,27	25,64	47,21	51,00
38	24,88	27,34	49,51	53,38
40	26,51	29,05	51,81	55,76
42	28,14	30,77	54,09	58,12
44	29,79	32,49	56,37	60,48
46	31,44	34,22	58,64	62,83
48	33,10	35,95	60,91	65,17
50	34,76	37,69	63,17	67,50
52	36,44	39,43	65,42	69,83
60	43,19	46,46	74,40	79,08
62	44,89	48,23	76,63	81,38
70	51,74	55,33	85,53	90,53
72	53,46	57,11	87,74	92,81
80	60,39	64,28	96,58	101,88
82	62,13	66,08	98,78	104,14
90	69,13	73,29	107,57	113,15
92	70,88	75,10	109,76	115,39
100	77,93	82,36	118,50	124,34
110	86,79	91,47	129,39	135,48
120	95,70	100,62	140,23	146,57
130	104,66	109,81	151,05	157,61
140	113,66	119,03	161,83	168,61
150	122,69	128,28	172,58	179,58
160	131,76	137,55	183,31	190,52
170	140,85	146,84	194,02	201,42
180	149,97	156,15	204,70	212,30
190	159,11	165,49	215,37	223,16
200	168,28	174,84	226,02	233,99
	-1,64	-1,28	+1,28	+1,64
Примечание - Линейная интерполяция промежуточных значений дает достаточно точные результаты.     Значения для различных целых включены. Для 200   ( - квантиль нормированного нормального распределения, приведенный в последней строке каждой графы).   Дополнительную информацию см. в [2].

Приложение Е

(справочное)

Интеграл вероятностей

-распределения и функция распределения Пуассона

Таблица Е.1

DF	J
		0,005	0,010	0,050	0,100	0,500	1,00	2,00	3,00	4,00	5,00	6,00	7,00
		0,0025	0,005	0,025	0,050	0,250	0,500	1,0	1,5	2,0	2,5	3,0	3,5
2	0	0,9975	0,9950	0,9753	0,9512	0,7788	0,6065	0,3679	0,2231	0,1353	0,0821	0,0498	0,0302
4	1	1,0000	1,0000	0,9997	0,9988	0,9735	0,9098	0,7358	0,5578	0,4060	0,2873	0,1991	0,1359
6	2	1,0000	1,0000	1,0000	1,0000	0,9978	0,9856	0,9197	0,8088	0,6767	0,5438	0,4232	0,3208
8	3	1,0000	1,0000	1,0000	1,0000	0,9999	0,9982	0,9810	0,9344	0,8571	0,7576	0,6472	0,5366
10	4	1,0000	1,0000	1,0000	1,0000	1,0000	0,9998	0,9963	0,9814	0,9473	0,8912	0,8153	0,7254
12	5	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	0,9994	0,9955	0,9834	0,9580	0,9161	0,8576
14	6	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	0,9999	0,9991	0,9955	0,9858	0,9665	0,9347
16	7	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	0,9998	0,9989	0,9958	0,9881	0,9733
18	8	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	0,9998	0,9989	0,9962	0,9901
20	9	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	0,9997	0,9989	0,9967
22	10	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	0,9999	0,9997	0,9990
24	11	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	0,9999	0,9997
26	12	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	0,9999
28	13	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
30	14	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
32	15	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
34	16	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
36	17	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
38	18	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
40	19	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
42	20	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
44	21	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
46	22	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
48	23	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
50	24	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
52	25	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
54	26	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
56	27	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
58	28	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
60	29	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
62	30	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
64	31	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
66	32	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
68	33	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
70	34	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
72	35	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
74	36	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
76	37	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
78	38	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
80	39	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
82	40	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
84	41	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
86	42	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
88	43	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
90	44	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
92	45	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
94	46	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
96	47	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
98	48	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
100	49	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
102	50	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
104	51	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000

 Продолжение таблицы Е.1

DF	J
		8,00	9,00	10,00	11,00	12,00	13,00	14,00	15,00	16,00	17,00	18,00	19,00
		4,0	4,5	5,0	5,5	6,0	6,5	7,0	7,5	8,0	8,5	9,0	9,5
2	0	0,0183	0,0111	0,0067	0,0041	0,0025	0,0015	0,0009	0,0006	0,0003	0,0002	0,0001	0,0001
4	1	0,0916	0,0611	0,0404	0,0266	0,0174	0,0113	0,0073	0,0047	0,0030	0,0019	0,0012	0,0008
6	2	0,2381	0,1736	0,1247	0,0884	0,0620	0,0430	0,0296	0,0203	0,0138	0,0093	0,0062	0,0042
8	3	0,4335	0,3423	0,2650	0,2017	0,1512	0,1118	0,0818	0,0591	0,0424	0,0301	0,0212	0,0149
10	4	0,6288	0,5321	0,4405	0,3575	0,2851	0,2237	0,1730	0,1321	0,0996	0,0744	0,0550	0,0403
12	5	0,7851	0,7029	0,6160	0,5289	0,4457	0,3690	0,3007	0,2414	0,1912	0,1496	0,1157	0,0885
14	6	0,8893	0,8311	0,7622	0,6860	0,6063	0,5265	0,4497	0,3782	0,3134	0,2562	0,2068	0,1649
16	7	0,9489	0,9134	0,8666	0,8095	0,7440	0,6728	0,5987	0,5246	0,4530	0,3856	0,3239	0,2687
18	8	0,9786	0,9597	0,9319	0,8944	0,8472	0,7916	0,7291	0,6620	0,5925	0,5231	0,4557	0,3918
20	9	0,9919	0,9829	0,9682	0,9462	0,9161	0,8774	0,8305	0,7764	0,7166	0,6530	0,5874	0,5218
22	10	0,9972	0,9933	0,9863	0,9747	0,9574	0,9332	0,9015	0,8622	0,8159	0,7634	0,7060	0,6453
24	11	0,9991	0,9976	0,9945	0,9890	0,9799	0,9661	0,9467	0,9208	0,8881	0,8487	0,8030	0,7520
26	12	0,9997	0,9992	0,9980	0,9955	0,9912	0,9840	0,9730	0,9573	0,9362	0,9091	0,8758	0,8364
28	13	0,9999	0,9997	0,9993	0,9983	0,9964	0,9929	0,9872	0,9784	0,9658	0,9486	0,9261	0,8981
30	14	1,0000	0,9999	0,9998	0,9994	0,9986	0,9970	0,9943	0,9897	0,9827	0,9726	0,9585	0,9400
32	15	1,0000	1,0000	0,9999	0,9998	0,9995	0,9988	0,9976	0,9954	0,9918	0,9862	0,9780	0,9665
34	16	1,0000	1,0000	1,0000	0,9999	0,9998	0,9996	0,9990	0,9980	0,9963	0,9934	0,9889	0,9823
36	17	1,0000	1,0000	1,0000	1,0000	0,9999	0,9998	0,9996	0,9992	0,9984	0,9970	0,9947	0,9911
38	18	1,0000	1,0000	1,0000	1,0000	1,0000	0,9999	0,9999	0,9997	0,9993	0,9987	0,9976	0,9957
40	19	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	0,9999	0,9997	0,9995	0,9989	0,9980
42	20	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	0,9999	0,9998	0,9996	0,9991
44	21	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	0,9999	0,9998	0,9996
46	22	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	0,9999	0,9999
48	23	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	0,9999
50	24	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
52	25	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
54	26	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
56	27	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
58	28	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
60	29	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
62	30	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
64	31	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
66	32	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
68	33	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
70	34	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
72	35	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
74	36	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
76	37	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
78	38	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
80	39	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
82	40	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
84	41	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
86	42	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
88	43	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
90	44	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
92	45	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
94	46	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
96	47	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
98	48	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
100	49	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
102	50	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000
104	51	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000	1,0000