- •Федеральное агентство по образованию московский государственный технический университет «мами»
- •Б.А.Дидусёв
- •Часть 2. Математический аппарат исследования надёжности технических
- •Содержание
- •1.Математический аппарат исследования
- •1.1.Основные предпосылки.
- •1.2.Качество статистической информации.
- •2.Первичная обработка экспериментального
- •2.1.Упорядочение выборочных наблюдений.
- •2.2.Вычисление частостей (относительных частот).
- •2.4.Графическое представление статистических выборок.
- •2.5.Пример первичной обработки экспериментальных данных.
- •3.Вероятностные распределения, используемые
- •3.1.Общие замечания.
- •4. Точечные и интервальные оценки
- •5.Анализ однородности исходного
- •5.1.Цель анализа и общие подходы.
- •5.2. Дисперсионный анализ при оценке однородности выборок.
- •5.3.Критерий Вилкоксона.
- •5.4.Критерий Смирнова – Колмогорова.
- •6.Оценка параметров распределения
- •6.1.Метод моментов.
- •6.2.Метод максимального правдоподобия (метод наибольшего правдоподобия, метод максимума правдоподобия).
- •6.3.Метод квантилей.
- •6.4.Графический метод.
- •7.Проверка согласия эмпирического
- •7.1.Общие принципы оценки расхождения.
- •7.2.Критерий согласия λ.
- •7.3.Критерий согласия χ2.
- •7.4.Критерий ω2.
- •8. Представление эмпирического распределения
- •9.Корреляционный (регрессионный) анализ
- •10.Дисперсионный и последовательный анализы
- •10.1Дисперсионный анализ как факторный анализ.
- •10.2.Последовательный анализ.
8. Представление эмпирического распределения
В ВИДЕ СУММЫ РАСПРЕДЕЛЕНИЙ
При обработке экспериментальных данных для получения характеристик надёжности нередко бывает удобным представить исследуемую выборку как смесь (сумму, суперпозицию) нескольких распределений.
К таким распределениям могут привести различные причины изготовления и эксплуатации – изготовление одних и тех же деталей на различном оборудовании или по разной технологии, изменение конструкции детали или узла, различия в условиях эксплуатации.
По внешнему виду распределения, например, наличию двух и более максимумов, можно судить о целесообразности использования суперпозиции распределений. В этом случае плотность эмпирического распределения представляют в виде суммы k
f(x) = ∑ cifi(x) (52)
i=1
где ci – коэффициенты весомости i – го распределения (доля реализаций i –го распределения в смешанной выборке), fi(x) –плотность i – го распределения определённого вида, k
∑ сi = 1.
i=1
Оценку математического ожидания x и дисперсии s суперпозиции производят по формулам: k
x = ∑ ci xi; (53)
k k i=1
s = ∑ cisi2 + ∑ ci (xi – x)2, (54)
i=1 i=1
где xi и si - оценки математического ожидания и дисперсии i- го распределения.
Например, если смесь распределений можно представить в виде k
нормальных распределений, то функция распределения суперпозиции
k x k
Fc(x1) = ∑ ci (1/√2π si) ∫ exp [ – (x – xi)2 / 2s2i ]dx = ∑ ci Ф[(x1 – xi) / si], (55)
i=1 – ∞ i=1
а плотность k
fc(x) = ∑( ci/si) φ[(x–xi)/si]. (56)
i=1
При этом функция распределения суперпозиции в общем случае уже не будет иметь нормального распределения.
9.Корреляционный (регрессионный) анализ
ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ
При анализе результатов экспериментальных данных по надёжности часто приходится рассматривать распределение не только одной случайной величины (например, наработки на отказ,), но и влияние на эту случайную величину другой случайной величины, то есть приходится устанавливать – есть ли и какая взаимная связь двух случайных величин.
Обычно говорят, что величины Х и Y связаны функционально, если каждому значению х1 соответствует определённое (может быть не одно) значение у1. Между случайными величинами, кроме функциональной, может существовать ещё и стохастическая связь. При стохастической связи одна из случайных величин реагирует на изменение другой или других случайных величин изменениями параметров закона распределения. Исследованием круга вопросов, связанных с использованием стохастической связи, занимается теория корреляции (теория стохастической регрессии).
При оценке стохастической связи между двумя или несколькими случайными величинами определяют форму связи (криволинейная или прямолинейная) и силу (тесноту) связи. Метод определения тесноты связи (степени взаимосвязи) между случайными величинами называют корреляционным анализом.
При рассмотрении системы двух случайных величин Х и Y необходимо ввести новые статистические характеристики, устанавливающие их взаимную связь.
Стохастическую зависимость Y от Х описывают условным математическим ожиданием ∞
Y(x) = M[Y / (X=x)] = ∫ y f(Y/x) dy. (57)
–∞
В механической аналогии распределения, если единичная масса распределена на плоскости хоу с плотностью f(x,y), то Y(x) есть ордината центра тяжести массы, распределённой на прямой Х = х.
Дисперсия M[(Y – α)2] минимальна при α = M[Y / (X=x)].
Поэтому линия (57) даёт наилучшее предсказание значения величины Y по значению Х = х и называется линией регрессии.
При исследовании линии регрессии (57) и определяют форму связи случайных величин. Таким образом, регрессия Y по Х определяет изменение математического ожидания Y при изменении Х.
Форма связи исследуемых величин зависит от физической сущности явления (физической сущности отказов) и может иметь вид прямой линии, параболы и т.д. Определение формы линии регрессии, соответствующей реальной форме зависимости, составляет принципиальный момент в изучении корреляции. Глубокое знание исследуемых явлений, влияющих на отказы, зависимостей между ними, является существенным элементом для правильного выбора линии регрессии.
При линейной зависимости
Y(x) = a + bx. (58)
Задача определения линии регрессии сводится в этом случае к определению коэффициентов a и b. Эта задача решается методом наименьших квадратов, согласно которому необходимо, чтобы математическое ожидание квадратов отклонений условных средних от расчётных значений было минимально, т.е. надо найти a и b так , чтобы обеспечить
min M[y(x) – y(x)] = min M[y(x) – a – bx].
Тесноту связи между случайными величинами при линейной корреляции оценивают корреляционным моментом (ковариацией) Kxy и коэффициентом корреляции rxy.
Корреляционный момент представляет собой математическое ожидание произведения отклонений X и Y от их математических ожиданий
Kx = M[(X – MX) (Y – MY)] (59)
Так как Kxy зависит от единиц измерения и сам по себе не может служить показателем связи, то рассматривают связь нормированных отклонений
(X – MX)/σx и (Y – MY)/σy ,
в результате чего получают коэффициент корреляции
rxy = Kxy/σxσy . (60)
Коэффициент корреляции может находиться в пределах – 1≤ rxy ≤ 1.
Чем ближе к нулю │rxy│, тем слабее линейная связь между величинами, чем │rxy│ , ближе к единице, тем связь сильнее. При │rxy│=0 (или близком к нулю), линейная корреляционная связь отсутствует. При │rxy│=1 (или близком к единице) статистическая линейная связь становится функциональной.
Не следует думать, что если линейная корреляционная связь отсутствует, то нет и других форм связи. Криволинейная связь может существовать и даже быть функциональной.
При криволинейной корреляционной связи теснота связи характеризуется корреляционным соотношением ρxy
√ 1 – ρxy = ( √ M[{Y – y(x)}2]) / σy = { σ [Y – y(x)]} / σy . (61)
Из (61) также следует, что – 1 ≤ ρxy ≤ 1.
