- •Оглавление
- •Введение
- •Часть 1. Симуляция и сбор статистики на имитационной модели.
- •Часть 2. Анализ и обработка полученной информации.
- •Часть 2 лабораторных работ выполняется с целью ознакомления студентов с методами статистической обработки результатов испытаний изделий.
- •Лабораторная работа № 1. Обоснование статистических решений. Сравнение методов Байеса и минимакса.
- •Постановка задачи.
- •Описание методов : Метод Байеса :
- •Формирование расчетных алгоритмов Алгоритм на основе байесовского принципа.
- •Алгоритм на основе минимаксимального принципа
- •Рекомендации по выполнению работы.
- •Минимаксный принцип.
- •Варианты заданий
- •Пример оформления отчета по лабораторной работе №1
- •Обоснование статистических решений. Сравнение методов Байеса и минимакса.
- •Лабораторная работа № 2. Проверка статистических гипотез о законах распределения. Определение гарантированной скорости предела сквозного пробития Vпсп методом статистических испытаний.
- •Постановка задачи.
- •Описание метода.
- •Рекомендации по выполнению работы.
- •Теоретический материал
- •Исходные данные для лабораторной работы
- •Варианты заданий
- •Пример оформления отчета по лабораторной работе №2
- •Проверка статистических гипотез о законах распределения
- •Лабораторная работа №3. Определение характеристик случайных процессов.
- •Постановка задачи
- •Описание метода.
- •Рекомендации по выполнению работы.
- •2. Определение корреляционной функции Kx(t,t'):
- •Пример оформления отчета по лабораторной работе №3
- •Определение характеристик случайных процессов.
- •Исходные данные.
- •Нормированная корреляционная функция стационарного процесса
- •Лабораторная работа № 4 Планирование эксперимента. Регрессионный анализ
- •Постановка задачи.
- •Описание метода.
- •Рекомендации по выполнению работы.
- •Пример оформления отчета по лабораторной работе №4
- •Планирование эксперимента. Регрессионный анализ
- •Приложение 1
- •Приложение 2 Значения критерия Пирсона
- •Приложение 3
- •Литература:
Пример оформления отчета по лабораторной работе №2
ЛАБОРАТОРНАЯ РАБОТА N 2
Проверка статистических гипотез о законах распределения
гр. Е-341 Александров Александр Алексеевич
Цель работы: Ознакомление с методом статистических испытаний. Определение вида и параметров закона распределения случайной величины – скорости сквозного пробития преграды. Определение гарантированной скорости сквозного пробития.
Метод исследования: метод статистических испытаний заключается в воспроизведении исследуемого физического процесса при помощи вероятностной математической модели и вычислении характеристик этого процесса.
ИСХОДНЫЕ ДАННЫЕ
калибр [м] 0.0370 0.0002
масса [кг] 0.7200 0.0048
высота головной части [м] 0.0300 0.0015
толщина преграды [м] 0.0300 0.0002
предел текучести мат. преграды [Па] 910000000 9100
плотность мат. преграды [кг/м3] 7800.0000 100.0000
диаметр притупления головной части [м] 0.0200 0.0002
количество опытов 500
количество разрядов гистограммы 10
Результаты расчета
Vmin Vmax Vsr
511.21 543.00 524.45
Дисперсия Среднеквадратичное
отклонение
26.75 5.17
N gr V(i) V(i+1) N
--------------------------------------
1.00 511.21 514.39 12.00
2.00 514.39 517.57 35.00
3.00 517.57 520.74 72.00
4.00 520.74 523.92 111.00
5.00 523.92 527.10 108.00
6.00 527.10 530.28 89.00
7.00 530.28 533.46 52.00
8.00 533.46 536.64 18.00
9.00 536.64 539.82 1.00
10.00 539.82 543.00 1.00
Статистическая обработка результатов
Для аппроксимации полученного выборочного распределения с помощью нормального закона необходимо заменить теоретические параметры нормального распределения тх и х их выборочными значениями тх и x, полученными в результате расчета
тх = 524,47 м / с
х = 5,17 м / с
Для определения меры расхождения теоретического и эмпирического законов распределения необходимо вычислить вероятность:
Рi = Ф [(Vi+1 - mx ) / х] - Ф [( Vi - mx ) / х ] = Фi+1 - Фi
Ф - табулированная функция , (см. приложение 1) .
Результаты расчетов сведем в таблицу:
i |
V(i) , м / с |
V(i+1) , м / с |
(Vi - mx )/ х |
(Vi+1 - mx)/ х |
Фi |
Фi+1 |
Рi |
n Рi |
Ni |
2 |
1 |
511,21 |
514,39 |
-2,56 |
-1,95 |
0,0054 |
0,0256 |
0,0202 |
10 |
12 |
0,3574 |
2 |
514,39 |
517,57 |
-1,95 |
-1,33 |
0,0256 |
0,0885 |
0,0629 |
31 |
35 |
0,4007 |
3 |
517,57 |
520,74 |
-1,33 |
-0,72 |
0,0885 |
0,2266 |
0,1381 |
69 |
72 |
0,1260 |
4 |
520,74 |
523,92 |
-0,72 |
-0,11 |
0,2266 |
0,4404 |
0,2138 |
107 |
111 |
0,1572 |
5 |
523,92 |
527,10 |
-0,11 |
0,51 |
0,4404 |
0,6736 |
0,2332 |
117 |
108 |
0,6343 |
6 |
527,10 |
530,28 |
0,51 |
1,2 |
0,6736 |
0,8531 |
0,1795 |
90 |
89 |
0,0063 |
7 |
530,28 |
533,46 |
1,2 |
1,74 |
0,8531 |
0,9505 |
0,0974 |
49 |
52 |
0,2236 |
8 |
533,46 |
536,64 |
1,74 |
2,35 |
0,9505 |
0,9878 |
0,0373 |
19 |
18 |
0,0227 |
9 |
536,64 |
539,82 |
2,35 |
2,97 |
0,9878 |
0,9978 |
0,01 |
5 |
1 |
3,2000 |
10 |
539,82 |
543,00 |
2,97 |
3,58 |
0,9978 |
0,9997 |
0,0019 |
1 |
1 |
0,0026 |
=5,1309
На графике представлена сравнительная гистограмма, показывающая расхождение между теоретической частотой попадания значений в соответствующие разряды n Рi и полученной в численном эксперименте Ni
Рассчитаем критерий Пирсона 2 по формуле:
2 = ( Ni - n pi ) 2 / п р i = 5,1309
Оценочными значениями заменены 2 параметра нормального распределения , т (s=2), поэтому число степеней свободы:
k= n - 1 - s = 10 - 1 - 2 = 7
Из таблиц (приложение2) для уровня значимости =0,05 и k=7 находим критическое 2k .=14,1.
Сравнивая критическое значение с расчитаным выше, можем убедиться, в том, что гипотезу распределения скорости пробития по нормальному закону можно считать правдоподобной.
Среднее значение скорости пробития Vср=524,45 м/с, что соответствует вероятности 0,5. Поскольку скорость предела сквозного пробития соответствует вероятности 0,7 (т.е. два пробития из трех попыток), то для этой вероятности из таблиц находим значение =0,525.
Тогда скоростьVпсп будет вычисляться по формуле:
Vпсп=Vср + х*=524,45+ 0,525 5,17=527,16 м/с.
Определим доверительный интервал по формуле:
=t х* / N = 2,09 5,17/ 500=0,48
Таким образом, можем записать окончательный результат:
Vпсп=527,2 0,48 м/с.
Выводы:
Гипотезу распределения скорости пробития по нормальному закону можно считать правдоподобной так как значение параметра 2=5,13 2k . =14,1
В результате расчета была определена скорость и доверительный интервал скорости сквозного пробития, соответствующей вероятности 0,7.
