Файл: Шаракшанэ, А. С. Испытания сложных систем учеб. пособие.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 21.10.2024
Просмотров: 143
Скачиваний: 0
При натурных испытаниях сложных систем обычно привлекают большое количество измерительной техники, обслуживающий пер сонал которой имеет различный уровень подготовки. Вследствие этого, а также из-за случайной неоднотипности условий проведения натурных экспериментов порой возникают значительные ошибки в выходных результатах. При ограниченном объеме реальной инфор мации эти ошибки могут привести к неправильным статистическим выводам. Чтобы избежать подобные явления, при анализе сложных систем большое внимание уделяют определению условий, при кото рых можно исключить из рассмотрения резко выделяющиеся на блюдения.
При выявлении подобных аномальных результатов наблюдений прежде всего необходимо провести тщательный анализ измерений и проверить, не являются ли эти результаты -следствием грубого промаха или нарушения условий измерений. Если подобный анализ не устранит появившихся сомнений, то прибегают к помощи статис тических методов выявления грубых ошибок, позволяющих произ водить целесообразную отбраковку аномальных данных. Для этого можно использовать ряд различных критериев: Колмогорова; Пир сона; Шовенэ и др. В частности, если допустить предположение о нормальном распределении погрешностей измерений, то можно при бегнуть к-правилу «трех-сигма».
Учитывая, что появление значительной по абсолютной величине ошибки маловероятно, так как
Р (— Зо.г < х — т*х < За*} 0,997,
где Зет* — предельная ошибка измерений, то все |х*—т х*|, превы шающие величину За*, можно отнести к.категории грубых и из даль нейшего рассмотрения их исключить. Однако подобный подход к отработке результатов измерений требует особо тщательного ана лиза, так как при малой выборке оценить правильно предельную погрешность достаточно трудно.
Для отбраковки грубых ошибок молено также воспользоваться критерием Смирнова. При этом необходимо подсчитать среднюю
арифметическую величину
п
и эмпирическое среднее квадратическое отклонение
Затем необходимо найти отношение абсолютной величины раз ности между сомнительным результатом измерений х* и средним арифметическим значением к величине S„, т. е. рассчитать статис
124
тику:
£==| |
х п|/S ii. |
Бели для данного числа измерений п и выбранной надежности Р= 1—р величина | превосходит критическое значение £П(Р), то с вероятностью, большей Р, можно считать сомнительный результат измерения грубой ошибкой.
|
|
|
|
|
|
|
|
Т а б л и ц а |
7.2.3 |
|
|
|
Критические значения ^ |
|
|
|
|||
|
|
|
|
|
9 |
|
|
|
|
п |
0,0005 |
0,001 |
0,002 |
0,005 |
0,01 |
0,02 |
0,05 |
0,10 |
0,20 |
|
|||||||||
3 |
1,414 |
1,414 |
1,414 |
1,414 |
1,414 |
1,414 |
1,414 |
1,412 |
1,406 |
4 |
1,732 |
1,732 |
1,731 |
1,730 |
1,728 |
1,723 |
1,710 |
1,689 |
1,645 |
5 |
1,996 |
1,994 |
1,990 |
1,982 |
1,972 |
1,955 |
1,917 |
1,869 |
1,791 |
6 |
2,219 |
2 , 2 1 2 |
2,203 |
2,183 |
2,161 |
2,130 |
2,067 |
1,996 |
1,894 |
7 |
2,408 |
2,395 |
2,377 |
2,344 |
2,310 |
2,265 |
2,182 |
2,093 |
1,974 |
8 |
2,568 |
2,547 |
2,521 |
2,476 |
2,431 |
2,374 |
2,273 |
2,172 |
2,041 |
9 |
2,704 |
2,677 |
2,643 |
2,586 |
2,532 |
2,464 |
2,349 |
2,238 |
2,097 |
10 |
2,822 |
2,788 |
2,747 |
2,680 |
2,616 |
2,540 |
2,414 |
2,294 |
2,146 |
11 |
2,925 |
2,884 |
2,837 |
2,760 |
2,689 |
2,606 |
2,470 |
2,343 |
2,190 |
12 |
3,015 |
2,969 |
2,915 |
2,830 |
2,753 |
2,663 |
2,519 |
2,387 |
2,229 |
13 |
3,098 |
3,044 |
2,984 |
2,892 |
2,809 |
2,713 |
2,563 |
2,426 |
2,264 |
14 |
3,167 |
3,111 |
3,046 |
2,947 |
2,859 |
2,759 |
2,602 |
2,461 |
2,297 |
15 |
3,232 |
3,171 |
3,102 |
2,997 |
2,905 |
2,800 |
2,638 |
2,494 |
2,327 |
16 |
3,290 |
3,225 |
3,152 |
3,042 |
2,946 |
2,837 |
2,670 |
2,523 |
2,354 |
17 |
3,343 |
3,274 |
3,198 |
3,083 |
2,983 |
2,871 |
2,701 |
2,551 |
2,380 |
18 |
3,392 |
3,320 |
3,240 |
3,120 |
3,017 |
2,903 |
2,728 |
2,577 |
2,404 |
19 |
3,437 |
3,361 |
3,278 |
3,155 |
3,049 |
2,932 |
2,754 |
2,601 |
2,426 |
20 |
3,478 |
3,400 |
3,314 |
3,187 |
3,079 |
2,959 |
2,779 |
2,623 |
2,447 |
2 2 |
3,552 |
3,469 |
3,378 |
3,245 |
3,132 |
3,008 |
2,823 |
2,664 |
2,486 |
24 |
3,616 |
3,529 |
3,434 |
3,295 |
3,179 |
3,051 |
2,862 |
2,701 |
2,521 |
26 |
3,673 |
3,582 |
3,483 |
3,340 |
3,220 |
3,089 |
2,897 |
2,734 |
2,553 |
28 |
3,724 |
3,629 |
3,528 |
3,380 |
3,258 |
3,124 |
2,929 |
2,764 |
2,582 |
30 |
3,769 |
3,672 |
3,567 |
3,416 |
3,291 |
3,156 |
2,958 |
2,792 |
2,609 |
35 |
3,866 |
3,762 |
3,652 |
3,494 |
3,364 |
3,224 |
3,022 |
2,853 |
2 ,6 6 8 |
40 |
3,943 |
3,835 |
3,720 |
3,557 |
3,424 |
3,281 |
3,075 |
2,904 |
2,718 |
45 |
4,007 |
3,896 |
3,778 |
3,610 |
3,474 |
3,329 |
3,120 |
2,948 |
2,762 |
50 |
4,062 |
3,948 |
3,827 |
3,556 |
3,518 |
3,370 |
3,160 |
2,987 |
2,800 |
В табл. 7.2.3 приведены критические границы для различных |3. Если л>50, то верхние р-процентные точки статистики | (тх*, ох*) можно рассчитывать при j3=SC0,2 по формулам [39]:
1р = к | / |
-------------------— — |
+ --------------------2и'|)(1/ 6 |
(2л — 5) |
; |
У |
2 л — 5 + м 2 + (3 + и2 |
|
u=W (l-P/2n),
где lF(p, п ) — функция, обратная функции нормального распреде ления.
125
Пример 1. Пусть в результате натурных испытаний получена выборка .v,, х2, .. . . хп, численные значения элементов которой совпадают с первой группой
нормально распределенных чисел, приведенных |
в работе [39]. Указанные числа |
||||||||
выписаны и представлены в табл. 7.2.4. |
|
|
|
|
|
||||
Докажем, |
что в результатах 19-го эксперимента не допущено грубых ошибок. |
||||||||
Для решения поставленной |
задачи |
рассчитаем |
выборочные |
значения: |
|||||
а) |
математического ожидания |
|
|
|
|
|
|||
|
|
|
|
/п* =0,211; |
|
Т а б л и ц а 7.2.4 |
|||
|
|
|
|
|
|
|
|||
1 |
Х1 |
1 |
xi |
1 |
x i |
1 |
xi |
1 |
x i |
1 |
0,464 |
11 |
2,455 |
21 |
—0,068 |
31 |
—0,288 |
41 |
0,241 |
2 |
0,060 |
12 |
—0,531 |
22 |
0,543 |
32 |
0,187 |
42 |
0,022 |
3 |
1,486 |
13 |
—0,634 |
23 |
0,926 |
33 |
0,785 |
43 |
—0,853 |
4 |
2,022 |
14 |
1,279 |
24 |
0,571 |
34 |
0,0194 |
44 |
—0,501 |
5 |
1,394 |
15 |
0,046 |
25 |
2,945 |
35 |
—0,258 |
45 |
0,439 |
6 |
0,137 |
16 |
—0,323 |
26 |
0,296 |
36 |
1,298 |
46 |
—0,957 |
7 |
—2,526 |
17 |
—0,194 |
27 |
—1,558 |
37 |
—1,190 |
47 |
0,525 |
8 |
—0,354 |
18 |
0,697 |
28 |
1,375 |
38 |
—0,963 |
48 |
—1,865 |
9 |
—0,472 |
19 |
3,521 |
29 |
—1,851 |
39 |
1,192 |
49 |
—0,273 |
10 |
—0,555 |
20 |
0,321 |
30 |
1,974 |
40 |
0,412 |
50 |
—0,035 |
б) |
дисперсии |
|
о*2 = |
1,3482. |
|
|
|
|
|
|
|
|
|
|
|
|
|
На основании этих значений величина статистики:
£so ( m * . ° х ) = 2 , 8 5 1 .
Если принять вероятность ложной отбраковки анализируемых результатов, равной 0,05, то по табл. 7.2.3 при заданном а=50 нетрудно найти ей соответст вующее значение £ р =о,о5=3,16.
Сопоставляя выборочное значение | 50 (тх*, |
ох*) |
с критическим | р, приходим |
|
к следующему выводу: так как £5о<£ р = 0 ,0 5 , |
то |
нет |
никаких оснований считать |
результаты 19-го эксперимента аномальными |
среди |
всех остальных результатов. |
Нетрудно убедиться в том, что этот вывод справедлив и для (3<0,10.
Пример 2. Пусть для п = 15 независимых равноточных измерений некоторой величины получено среднее арифметическое значение, равное тх*= 10,17 и эмпи рическое значение среднего квадратического отклонения сг**= 0,73.
Пусть значение jci5= 11,5 будет сомнительным. Необходимо, используя ста~и- стические методы, решить, можно ли отбраковать этот результат?
Находим значение
S i 5 ( m * , а * ) = 1 , 8 2 .
Для и=15 и 1—Р—0,95 критическим значением £ р является величина 2,638,
что свидетельствует о том, что данный результат не следует относить к числу грубых ошибок.
§7.3. ОЦЕНКА ДИСПЕРСИЙ И МАТЕМАТИЧЕСКИХ ОЖИДАНИЙ
Врезультате измерений получают некоторую выборку, состоя щую из совокупности значений хи х2, х3, ..., хп. В общем случае каж дый замер отличается от других замеров и их отклонение от истин
126
ного значения искомого параметра носит случайный характер (при этом предполагается, что результаты измерений уже свободны от грубых и систематических ошибок).
Чтобы по данной совокупности измерений получить приближен ное значение величины х, необходимо найти центр группирования последней и оценить рассеивание всех рассматриваемых результа тов относительно этого центра. В качестве числовых характеристик центра группирования обычно принимают математическое ожида ние, медиану и моду. Правда, последние две характеристики в этих целях используют сравнительно редко. На практике для характе ристики центра группирования определяют эмпирическую число вую характеристику, среднюю арифметическую величину, получен ную по совокупности значений хь х2, ..., хп. Среднее арифметическое значение х* в случае равноточных измерений находят по формуле:
П
Как следует из обоснований, приводимых в математической ста тистике, такая оценка является несмещенной и состоятельной. Если ошибки измерения подчинены нормальному закону распределения, то эта оценка будет к тому же и эффективной. Последнее предполо жение о нормальном характере распределения ошибок измерения на практике в большинстве случаев оправдано.
Если же распределение ошибок измерений подчиняется другому закону, то эффективной оценкой может быть другая статистика. Так, для равномерного распределения ошибок измерений медиана является более эффективной оценкой, чем арифметическая средняя.
В случае неравноточных измерений в качестве оценки искомой величины х принимают взвешенное среднее арифметическое значе ние:
|
P i x |
i |
х = P l x l + Р 2 Х 2 + . . • |
+ Р п х п ___ /=1 |
(7.3.2) |
PY + Р 2 + . . . |
+ . Р п |
|
|
2 Pi |
|
|
i = l |
|
где pi —веса измерений, принимаемые обычно обратно пропорцио нальными значениям дисперсий ошибок, т. е.
А = 1/°/, ^=1, 2,-.., л. |
(7.3.3) |
Вследствие влияния различных случайных факторов данные на блюдений можно охарактеризовать определенным рассеиванием результатов измерений. В качестве характеристики рассеивания случайных величин при обработке результатов испытаний обычно
127
принимают эмпирическую дисперсию:
п
(7.3.4)
Приведенное выражение для эмпирической дисперсии позволяет получить несмещенную, 'состоятельную и асимптотически эффектив ную оценку. Последнее означает, что при п, стремящемся к беско нечности, отношение дисперсии D{SX2} к минимально возможной ее величине, неограниченно приближается к 1.
Для равноточных и независимых измерений дисперсию оценки можно охарактеризовать выражением:
D |
= |
'■'7^1 |
I |
'fSi |
= |
(7‘3'5) |
|
|
|
|
Таким образом, средняя квадратическая ошибка средней ариф метической величины может быть представлена в виде:
/■}/«, |
(7.3.6) |
т. е. средняя арифметическая величина х* в У п раз точнее отдель ного измерения Xi (i= 1, 2, ..., /г). В случае неравноточиых измере ний необходимо учитывать веса последних. Взвешенное эмпириче ское значение среднего квадратического отклонения величин хи Л'о, ..., хп от их среднего значения х* при этом можно определить по формуле:
2 Pi'x > -л;*)2 |
(7.3.7) |
L i = i |
|
где Pi — веса измерений. |
|
Если же имеют к серий наблюдений, в которых было получено соответственно пи п2, ..., па количество измерений, то оценку диспер сии можно получить в виде средней взвешенной величины из эмпи рических дисперсий по сериям:
(щ — 1) + («2 — |
1) S'j.a + . . . + (« й— 1) |
(щ — 1) + («2 |
— 1) + ••• + ( п к — 1) |
Для эмпирической оценки среднего квадратического отклонения величины х* в случае неравноточных измерений можно воспользо ваться формулой:
р Лх . - х У |
(7.3.8) |
|
/=1 |
Рассуждая аналогично, находят погрешность, допускаемую при оценке среднего квадратического отклонения ах. В математической
128
статистике доказывается, что эмпирическая дисперсия распределе на по нормальному закону и характеризуется средним квадратиче ским отклонением:
а {S2x)— <sxJ y 2 (и— 1).
Рассмотренные эмпирические характеристики представляют со бой оценки искомых параметров в виде чисел, характеризующих распределение случайных величин в данной совокупности наблюде ний. Такую оценку параметров называют точечной. Кроме этого способа оценки параметров на практике прибегают к другому спо собу, основанному на определении интервала, накрывающего ис тинное значение искомого параметра с заданной вероятностью. Эту оценку называют интервальной, или доверительной.
Например, доверительная оценка для математического ожида ния в симметричном'случае может быть представлена в виде:
е<[.х<С-х:*-|-е или | х — х* | s. |
(7.3.9) |
Доверительный интервал 1$ = (х*—е, х* + г) определяет область возможных значений полученной при обработке средней арифмети ческой величины х* для данного параметра х. Причем попадание х* в эту область гарантируется с заданной доверительной вероятно стью, т. е.
(|JC — JC* | < ep)= P,
где х* — среднее арифметическое значение, полученное в результа те обработки данных наблюдений; х — истинное значение искомого параметра; ер — наперед заданная положительная величина, опре деляемая по доверительной вероятности (надежности оценки), при нимаемой равной величине р.
В случае равноточных измерений, когда заранее известна точ ность измерений ах, доверительную оценку математического ожида ния можно представить через функцию Лапласа, предположив, что ошибки измерения подчинены нормальному закону. Тогда, прини мая во внимание, что
из уравнения
можно наити ер:
ер .= У1>аг* Ф -1ф ) = У 2 -Ц з - Ф -1 (р),
Уп
где Ф-1(р) — функция, обратная функции Лапласа.
129