Вариант 5 (все задачи) (Задачи 1-6 вариант 5), страница 2
Описание файла
Файл "Вариант 5 (все задачи)" внутри архива находится в папке "Задачи 1-6 вариант 5". Документ из архива "Задачи 1-6 вариант 5", который расположен в категории "". Всё это находится в предмете "теория вероятностей и математическая статистика" из , которые можно найти в файловом архиве МГТУ им. Н.Э.Баумана. Не смотря на прямую связь этого архива с МГТУ им. Н.Э.Баумана, его также можно найти и в других разделах. Архив можно найти в разделе "курсовые/домашние работы", в предмете "теория вероятности и математическая статистика" в общих файлах.
Онлайн просмотр документа "Вариант 5 (все задачи)"
Текст 2 страницы из документа "Вариант 5 (все задачи)"
1) Метод моментов
Объем выборки .
Выборочная средняя
является точечной оценкой математического ожидания.
Математическое ожидание экспоненциального распределения выразим через его параметр.
Приравнивая мат. ожидание его выборочной оценке, получим оценку неизвестного параметра распределения.
2) Метод наибольшего правдоподобия
Функция правдоподобия ,
t i – варианты выборки, .
.
Логарифмическая функция правдоподобия:
.
Найдем первую производную по λ: .
Приравняв её нулю, найдем критическую точку: .
Найдем вторую производную по λ: .
Вторая производная отрицательна, следовательно, - точка максимума и она принимается в качестве оценки наибольшего правдоподобия неизвестного параметра λ.
Задача 4.
Пусть определяется величина H с помощью измерительного прибора. Среднеквадратичное отклонение погрешности измерения 0 равно 10 мм. По результатам измерений получена следующая выборка:
435.7 мм, 438.6 мм, 411.3 мм, 410.1 мм, 424.5 мм, 425.7 мм.
Считая, что H ~ N(, 0), определить доверительный интервал для параметра с уровнем доверия 0.9.
Решение:
Доверительный интервал для оценки параметра нормально распределенной случайной величины с известным среднеквадратическим отклонением имеет вид , - квантиль нормального распределения.
Объем выборки .
Выборочная средняя
По заданному уровню доверия определяем .
По таблице квантилей нормального распределения, находим квантиль порядка .
.
< 431,02 мм
C вероятностью 0,9 параметр находится в интервале 417,5 мм < <431,02 мм.
Задача 5.
Пусть измеряемая величина Y ~ N(, ) является пределом прочности при растяжении материала (Ст. 3). По результатам испытаний образцов получена выборка:
436.8 МПа, 410.2 МПа, 446 МПа, 422.3 МПа, 441.2 МПа, 436.7 МПа.
Установить доверительные интервалы для и с уровнем значимости 0.1.
Решение:
Объем выборки n=6.
Выборочная средняя
Выборочная дисперсия
Исправленное среднеквадратическое отклонение
1) Доверительный интервал для параметра случайной величины, распределённой по нормальному закону с неизвестным среднеквадратическим отклонением, имеет вид .
- квантиль распределения Стьюдента с степенями свободы порядка .
По таблице, .
432,2 – 11 = 421,2 ; 432,2+11=443,2.
C вероятностью 0,9 параметр находится в интервале 421,2 < <443,2.
2) Доверительный интервал для среднеквадратического отклонения случайной величины , распределённой по нормальному закону, имеет вид
.
По таблице находим квантили распределения с степенями свободы.
.
80,83< σ2 <781,19 откуда 8,99<σ<27,95.
C вероятностью 0,9 параметр σ находится в интервале:
8,99 МПа < σ< 27,95 МПа.
Задача 6.
Данные измерения величины y в зависимости от факторов m, b и s представлены в табл. 1 и 2.
Таблица 1
Значения y в зависимости от m, b и s
№ опыта | m, кг | b, Hc/м | s, H/м | y, H |
1 | 5 | 1 | 1 | 26.42 |
2 | 1 | 1 | 1 | 15.23 |
3 | 5 | 5 | 1 | 47.93 |
4 | 1 | 5 | 1 | 35.75 |
5 | 5 | 1 | 3 | 44.82 |
6 | 1 | 1 | 3 | 32.33 |
7 | 5 | 5 | 3 | 64.32 |
8 | 1 | 5 | 3 | 50.96 |
Таблица 2
Данные специальной серии опытов по измерению величины y
№ опыта | m, кг | b, Hc/м | s, H/м | y, H |
1 | 3 | 3 | 2 | 40.82 |
2 | 3 | 3 | 2 | 38.94 |
3 | 3 | 3 | 2 | 38.6 |
4 | 3 | 3 | 2 | 39.58 |
5 | 3 | 3 | 2 | 39.37 |
6 | 3 | 3 | 2 | 38.1 |
7 | 3 | 3 | 2 | 40.39 |
8 | 3 | 3 | 2 | 38.71 |
Используя данные из этих таблиц, найти экспериментальную зависимость y(m, b, s). В качестве первого приближения следует выбрать линейную регрессионную модель вида
y(m, b, s) = 01 + 1 m + 2 b + 3 s + e, (*)
где 0, 1, 2, 3 — коэффициенты регрессии; 1, m, b, s — базисные функции (F0 = 1, F1 = m, F2 = b, F3 = s); e — случайная величина.
Решение:
Чтобы получить точечные оценки параметров регрессионной модели (*) , используем данные табл. 1 и составим матрицы X и Y:
где хij – значение j-того фактора в i-том опыте, уi – значение отклика в i-том опыте. Численные значения всех базисных функций представим в матричной форме:
где fij – значение j-й базисной функции Fj в i-том опыте.
Чтобы определить точечные оценки b0 ,b1 ,b2 ,b3 коэффициентов регрессии, применим метод наименьших квадратов (МНК). Для этого найдём:
И по известной формуле получаем матрицу искомых точечных оценок (b0,.., b3):
Кроме точечных оценок коэффициентов βi необходима оценка дисперсии σe2 случайной величины е.
Так как вид адекватной модели заранее не известен, то по известной формуле, используя данные специальной серии опытов (при фиксированных значениях факторов), найдём (при n=N = 8):
где уi – значение отклика в i-том опыте, причём эти значения не используются для получения точечных оценок коэффициентов регрессии.
Затем проведём статистический анализ, состоящий из проверки значимости коэффициентов регрессии, проверки адекватности и работоспособности регрессионной модели.
При проверке значимости коэффициентов регрессионной модели выясним, обусловлено ли отличие bi от нуля чисто случайными обстоятельствами или же это отличие неслучайно и вызвано тем, что в теоретической регрессионной модели действительно присутствует соответствующий коэффициент регрессии. Проверка осуществляется путем вычисления статистик (t0,...,t3):
Если |ti| t* то соответствующий коэффициент регрессии полагаем незначимым и исключаем из регрессионной модели. Критическое значение t* равно значению t ,a/2 , которое является квантилем уровня 1-a/2 распределения Стьюдента, число степеней свободы равно 7 и уровень значимости a соответствует 0,05, т.е.:
t* = t7, 0,025 = 2,37.
После проверки значимости всех коэффициентов регрессии получим регрессионную модель, содержащую только значимые коэффициенты регрессии β(0) , β(1) , β(2) , т.е. :
y(m, b, s) = (0) m + (1) b + (2) s + e ,
где m, b, s – базисные функции (F(0) = m, F(1) = b, F(2) = s).
Затем численные значения всех базисных функций представим в матричной форме:
где fij – значение j-й базисной функции F(j) в i-том опыте.
Точечные оценки b(0) , b(1) , b(2) значимых коэффициентов регрессии определяем аналогичным образом , т.е. :
Тогда уравнение регрессии имеет вид:
y = b(0)F(0) + b(1)F(1) + b(2)F(2) . (**)
Проверка адекватности регрессионной модели возможна, так как получена точечная оценка S2e = 0,855 дисперсии σ2е и выполнено условие d+1 = 3 < N = 8. (Здесь d –больший индекс коэффициента регрессии). Для такой проверки используется статистика F_:
,
где
и
,
и где у_i – значения, предсказанные с помощью регрессионного уравнения (**).
Критическое значение F* равно значению F1,2,a , которое является квантилем уровня 0,95 распределения Фишера, 1 = N – d – 1 = 5 и 2 = 7 – число степеней свободы S2e , т.е. :
F* = F5,7, 0,05 =3,97.
Так как F_ не превосходит критического значения F* , то модель адекватна. Следовательно, для регрессионной модели имеет смысл рассматривать вопрос о её работоспособности.
Чтобы получить представление о точностных свойствах регрессионной модели, вычислим по формуле:
где
, и где
Тогда
Найденное значение R2 расположено близко к 1, что свидетельствует о хорошей точности регрессионной модели.
Таким образом, искомая экспериментальная зависимость у от массы механизма m, коэффициента затухания b демпфера и жёсткости s пружины имеет вид:
y = 2,973*m + 4,906*b + 8,111*s.