Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Анализ точности

.pdf
Скачиваний:
40
Добавлен:
23.02.2015
Размер:
2.24 Mб
Скачать

 

 

 

 

 

 

 

 

 

 

Таблица 2.4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Параметры

Метод получения

Алгоритм вычисления

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

распределе-

случайного числа x

случайного числа

 

ния

 

 

 

 

 

 

 

 

 

 

 

 

Равномерное распреде-

a

Функциональное

a b a

ление на интервале

преобразование слу-

b

(a…b)

чайной величины

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Экспоненциальное

 

Функциональное

 

1

 

ln

 

преобразование слу-

 

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

 

 

 

чайной величины

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функциональное

 

 

 

 

 

 

 

 

 

 

 

k

1

ln

Распределение Вейбулла

преобразование слу-

 

k

чайной величины

 

 

 

 

 

 

 

 

 

 

 

 

 

Нормальный закон

M

Метод

суммирова-

 

 

 

 

12

 

 

M ( i

6)

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

 

ния

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

n – число

 

 

n

 

 

 

 

 

12

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

степеней

Метод

суммирова-

2j , где i

i 6

 

свободы

ния

 

j 1

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

программой, показан на рис. 2.13.

Рис. 2.13. Окно программы с результатами статистического моделирования

31

Вторая группа моделей

Ниже рассматриваются математические модели, построенные исходя из предположения наличия вероятностной (стохастической) зависимости ме-

жду входными и выходными параметрами. Из всего многообразия таких мо-

делей ограничимся двумя частными случаями – математическими моделями с регрессионной и корреляционной зависимостями между входными и вы-

ходными параметрами. В этих случаях модели представляют собой аналити-

ческие выражения, устанавливающие связи между числовыми характеристи-

ками (математическими ожиданиями и дисперсиями) с помощью коэффици-

ентов, определяемых эмпирически методами регрессионного и корреляцион-

ного анализа.

Основная задача регрессионного анализа заключается в установлении

зависимости ~ , где W – неслучайная входная переменная, которую мы

Y f {W}

можем изменять по своему усмотрению. Дальнейший порядок следующий:

1. Проводится эксперимент, в котором для каждого значения wi

( i 1, m) получаем n независимых значений выходного параметра Y. Результа-

ты эксперимента заносим в табл. 2.5.

 

 

 

 

 

Таблица 2.5

 

 

 

 

 

 

Вход

Выход

Повторы

yi

i2{y}

 

 

 

 

 

 

 

 

w1

y11,…, y1n

n

y1

2

{y}

 

 

 

 

 

1

 

 

w2

y21,…, y2n

n

y2

 

2{y}

 

 

 

 

 

 

2

 

 

...

 

 

 

 

 

 

wm

ym1,…, ymn

n

ym

 

2 {y}

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

2.Построчно вычисляем арифметические средние yi и дисперсии

i2{y} .

3.Проверяются гипотезы.

Согласно правилам теории вероятностей и математической статистики

для регрессионной зависимости выходные данные каждой строки (yi1,…,yin)

32

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

наковую дисперсию, но разные математические ожидания. Чтобы считать

исследуемую зависимость регрессионной, необходимо доказать следующее:

а) выходные данные каждой строки принадлежат генеральной сово-

купности с нормальным распределением;

б) дисперсии i2{y} однородны (т.е. некоторые различия в их значениях объясняются только малым размером выборок).

Доказательства (проверки гипотез) осуществляются путем вычисления соответствующих статистических критериев (Стьюдента, Колмогорова, Кох-

рена, Пирсона и др.) с последующим сравнением с их табулированными зна-

чениями. В зависимости от задаваемого уровня вероятности ошибки гипоте-

за принимается или отклоняется. Поскольку математическая модель пред-

ставляет приближенное описание того или иного процесса, то для решения

практических задач проверка подобных гипотез мало что значит.

Например, если экспериментальные данные yi1,…, yin графически будут

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

при не-

котором уровне

 

вероятности ошибки

можно считать, что они взяты из гене-

ральной совокупности с нормальным распределением.

 

 

 

 

 

 

4. Наносим на график (рис. 2.14) совокупность точек yi ( i 1, m ), для

которых

исходя из теоретико-профессиональных соображений подбираем

функцию

~

 

 

 

 

 

 

 

 

 

 

 

 

Y f {W}. В качестве последней часто принимают полином вида

~

 

2

... (где b0, b1, b2, … – постоянные коэффициенты).

 

 

Y b0 b1W b2W

 

 

 

5. Оцениваем степень соответствия выбранной функции и опытных

данных. Рассчитываем величину

 

 

 

 

 

 

 

 

 

 

 

 

 

m

~

 

 

 

 

 

 

 

 

 

 

 

N m n yi

}

2

 

 

 

 

 

 

v

 

 

y{wi

,

(2.7)

 

 

 

2

i

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m n

 

2

 

 

 

 

 

 

 

 

 

 

m 2

yij yi

 

 

 

 

i 1 j 1

33

где yi – средние арифметические выхода в точке wi;

~ – значение функции в точке w ; y{wi } i

yij – экспериментальные значения i 1, m , j 1, n ;

N – общее число точек наблюдений (N=mn).

Рис. 2.14. Построение линии регрессии

Числитель выражения (2.7) характеризует разброс экспериментальных групповых средних yi относительно выбранной линии регрессии (подобно методу наименьших квадратов).

Для проверки гипотезы статистическими методами необходимо задать-

ся малым уровнем значимости α (где α – вероятность принятия неправильно-

го решения). Далее по таблицам математической статистики для F-критерия Фишера со значениями (m-2, N-m) находят 100α-процентную точку v 2 . Если v2 v2 , то гипотеза о линии регрессии не противоречит экспериментальным данных.

Если влиянием факторов Z можно пренебречь (например, если считать

их только погрешностями измерений), то уравнение регрессии ~ пе-

Y f {W}

реходит в функциональную зависимость Y= f {W}.

Последовательность построения уравнения регрессии для случая не-

скольких входных параметров приведена в главе 4.

Основная задача корреляционного анализа заключается в исследовании наличия, формы и силы связи между двумя или несколькими случайными ве-

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

Корреляционная связь может существовать и между входами «черного ящи-

34

ка». Последовательность корреляционного анализа (на примере случайных величин X и Y).

1. Данные пассивного эксперимента

x1, x2,…, xN – для входной переменной; y1, y2,…, yN – для выходной переменной

заносим в корреляционную таблицу (табл. 2.6).

Таблица 2.6

x

y

y

min

y

min+1

...

y … y

max

y

i

 

2{y}

 

 

 

 

 

max-1

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

xminxmin+1

 

 

n11

 

...

n1l

 

l

y1

12{y}

 

 

 

 

 

 

 

 

 

 

n1 j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xmax…xmax-1

 

 

nk1

 

nkl

 

l

yk

 

k2 {y]

 

 

 

 

 

 

 

 

 

 

nkj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

k

 

 

N

 

 

 

 

ni1

 

nil

 

 

 

 

 

 

 

 

 

i 1

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x j

 

 

x1

 

xl

 

 

2

{x}

 

2{x}

2

{x}

 

 

j

 

 

 

1

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П р и м е ч а н и я: 1. Число интервалов определяют по формуле Стерджесса с округлением до ближайшего целого числа

 

 

 

 

 

 

 

k

xmax xmin

 

;

 

l

ymax

ymin

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 3,32 lg N

 

 

 

 

1 3,32 lg N

 

 

 

 

 

 

 

 

 

 

 

2. nij i

 

 

j

 

– число значений, попав-

 

 

 

 

 

 

 

 

1, k,

1, l

ших в соответствующий интервал.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее для строк и столбцов рассчитываем средние арифметические ( x j

и yi ) и дисперсии ( 2j {x} и i2{y} ) по формулам

 

 

 

 

 

 

 

 

 

 

k

 

 

 

l

 

 

k

xi

 

 

2

 

 

l

y j yi 2

 

 

 

 

nij xi

 

 

 

nij y j

 

nij

x j

 

nij

 

x

 

 

i 1

, y

 

 

j 1

, 2 {x}

 

i 1

 

 

 

 

 

 

, 2 {y}

j 1

 

,

 

j

 

k

 

i

 

l

j

 

 

 

k

 

 

 

 

i

 

l

 

 

 

 

nij

 

 

 

nij

 

 

 

 

nij

 

 

 

 

 

nij

 

 

 

 

i 1

 

 

 

j 1

 

 

 

 

i 1

 

 

 

 

 

 

 

j 1

 

где xi и yj – середины интервалов случайных величин X и Y соответственно.

35

2.На основании данных табл. 2.6 вычисляем:

a)коэффициент корреляции

 

 

 

 

k

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nij xi

 

 

 

y j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

Y

 

 

 

 

 

 

 

 

 

rxy

 

i 1

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

2

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ni xi

X

 

 

 

 

n j

y j

 

Y

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

k

 

 

 

 

 

l n j x j

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

n y

 

 

 

ni nij ,

n j nij ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

X

 

 

 

 

 

 

 

,

 

 

 

Y

i

 

i

;

 

 

 

N

 

 

 

 

N

 

 

j 1

i 1

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

б) корреляционное отношение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

ni

 

Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

y j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n j

 

Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1

3.По величине коэффициента корреляции и корреляционного отноше-

ния судят о форме и о силе связи (табл. 2.7).

Таблица 2.7

rxy

ηy

Характеристика связи

 

 

 

rxy = ±1

Линейная функциональная

rxy = 0

ηy =1

Криволинейная функциональная

rxy = 0

ηy =1

Связь отсутствует

rxy = 0

ηy <1

Криволинейная корреляционная

|rxy| = ηy

Точная линейная корреляционная

|rxy| > 0

ηy <1

Линейная корреляционная

 

 

 

Последний случай на практике встречается чаще всего.

4. Рассчитываем коэффициенты для проверки гипотез:

a) о нормальном законе распределения данных по строкам и по столб-

цам;

б) об однородности дисперсий по строкам и по столбцам;

в) о статистической значимости коэффициента корреляции.

Первые две гипотезы проверяются так же, как и в регрессионном ана-

лизе. Что касается коэффициента корреляции, то обычно проверяется обрат-

36

ная гипотеза – гипотеза об отсутствии корреляционной связи. Если неравен-

ство

rxy

 

N 2

 

t0,05 N 2 ,

 

 

 

 

 

 

1 rxy

 

 

 

где t0,05 n 2 – пятипроцентная точка распределения Стьюдента с N-2 степенями свободы, выполняется, то гипотеза об отсутствии корреляционной связи принимается. В противном случае с вероятностью 0,95 можно утверждать,

что корреляционная связь действительно имеет место.

5. Как и в случае регрессионного анализа наносим на график точки yi

 

 

~

f {X}).

(рис. 2.15) и соединяем их линией (обычно прямой Y

При линейной корреляционной зависимости уравнение регрессии име-

ет вид

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

Y a bX .

 

 

Коэффициенты a и b определяются из выражений

b rxy

{X }

 

 

 

 

 

 

a Y b X .

{Y },

Рис. 2.15. Построение линии регрессии

2.2.4. Проверка адекватности математической модели

Для оценки адекватности математических моделей и реальных техно-

логических процессов используют:

а) визуальное сопоставление соответствующих теоретических и экспе-

риментальных распределений, точностных диаграмм и т.п.;

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

экспериментальных распределений; чаще всего в качестве таких показателей

37

используют сумму квадратов различий в средних значениях (в регрессион-

ном и корреляционном анализе) или в количествах значений, соответствую-

щих каждому интервалу (например, при проверке гипотез по 2 -критерию Пирсона).

В последнем случае рассчитывается величина

2

эi

 

2

 

тi ,

k

 

 

 

 

 

i 1

 

тi

 

 

где эi – число значений экспериментального распределения, попавших в i

интервал;

т i – число значений теоретического распределения, попавших в i-й ин-

тервал, причем

где k – количество интервалов (при объединяют).

k

k

эi тi ,

i 1

i 1

эi

или т i меньше 5 соседние интервалы

3. МАТЕМАТИЧЕСКИЕ МОДЕЛИ Н.А. БОРОДАЧЕВА

Рассмотренные выше математические модели строились на основе классических методов теории вероятностей и математической статистики.

Они пригодны для описания в основном стабильных технологических про-

цессов (см. рис. 1.2 ). У нестабильного технологического процесса со време-

нем изменяется распределение выходного параметра (вид распределения, чи-

словые характеристики). Исследователей в подобных случаях интересует суммарное распределение выходного параметра за определенный период времени. Это необходимо для принятия своевременных мер по корректиров-

ке (подналадке) технологического процесса, для того чтобы значения выход-

ного параметра не вышли за пределы поля допуска. Рассмотренные ранее математические модели для решения подобных задач не годятся.

В теории точности существует раздел, в котором рассматривается класс распределений со следующими свойствами [13]:

38

а) мгновенное распределение – всегда нормальное;

б) математическое ожидание и дисперсия этих распределений является функцией времени, т.е.

 

 

 

1

 

 

t

{y}

 

 

 

 

exp

 

 

 

 

t

 

2

 

 

 

 

 

( y Ct )2 .

2 2

t

Впервые этот класс распределений был исследован профессором Боро-

дачевым Н.А. Для описания процессов возникновения случайных погрешно-

стей им был предложен ряд схем.

Схема 1

YX i .

i 1n

Погрешность Y рассматривается как функция воздействия n (где n

достаточно большое число) случайных факторов – причин, также имеющих какие-то распределения. При n распределение Y будет стремиться к га-

уссовому распределению. Поскольку число факторов и степень их влияния считаются относительно постоянными, то распределение Y может быть опи-

сано законом Гаусса с числовыми характеристиками m0 и ζ0. Таким образом,

для стабильного технологического процесса мы получаем закон распределе-

ния, не зависящий от времени (t). Он же будет и суммарным распределением за период T:

 

 

 

1

 

 

{y} 0

{y}

 

 

 

 

exp

 

 

 

 

0

 

2

 

 

 

 

 

( y m0 )2 .

2 2

0

Теоретическая точностная диаграмма такого процесса приведена на рис. 3.1.

Рис. 3.1. Теоретическая точностная диаграмма с m0 и ζ0

39

Если же технологический процесс таков, что на величину Y , помимо независимых или слабо зависимых факторов Xi, действует еще группа факто-

ров Cs(t ) , которые не являются случайными, причем число или значения этих факторов систематически изменяются в зависимости от времени, то получим следующую схему 2.

Схема 2

n

 

Yt X i Ct ,

 

i 1

 

где C – значение суммы неслучайных слагаемых

C(t ) , отвечающих моменту

t

s

времени t. В этом случае мгновенное распределение будет гауссовым, но ма-

тематическим ожиданием, зависящим от t. Следовательно,

 

 

 

1

 

 

( y Ct

)2

 

t

{y}

 

 

 

 

exp

 

 

 

 

,

 

 

 

 

2

2

 

0

 

2

 

 

 

 

 

0

 

 

 

где Ct m0 f {t}.

В отличие от регрессионной и корреляционной зависимостей, рассмот-

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

ления за весь период T tk t0 . Для этого применяется следующий искусст-

венный прием (правда, не совсем справедливый с позиций классической тео-

рии вероятностей): параметр t рассматривается как случайная величина, рав-

номерно распределенная в интервале от tk до to:

 

 

1

tk

 

1

 

 

( y Ct

)2

 

{y}

 

 

 

 

 

 

 

exp

 

 

 

 

dt .

 

 

0

 

 

 

 

2

 

 

 

 

 

 

 

tk

t0 t0

 

2

 

 

2

0

 

 

 

Теоретическая точностная диаграмма такого процесса приведена на рис. 3.2.

40