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

pdf.php@id=6180

.pdf
Скачиваний:
26
Добавлен:
14.05.2023
Размер:
7.56 Mб
Скачать

ПРАКТИЧЕСКОЕ ЗАНЯТИЕ №4 Практические кривые распределения статистических данных

4.1. Построение полигона распределения

Рассмотрим последовательность построения экспериментальной кривой распределения (рассеяния) размеров или погрешностей.

1.Изготовляется (обрабатывается) партия заготовок, деталей.

2.Измеряется каждая деталь (заготовка) обработанной партии по параметру, точность которого следует определить. Например, размер с заданным допуском 19,9+0,05.

Измерение деталей выполняют инструментом, цена деления

которого должна быть (1/6 1/10) , где допуск на измеряемый параметр. По результатам обмера составляется таблица (числовой ряд) распределения размеров (погрешностей), получают генеральную совокупность. Например:

19,93; 19,87; 19,97; 19,89; 19,95; 19,92; 19,89; 19,95; 19,93; …; 19,95; 19,88; 19,94; 19,93.

После получения выборки определяются возможные грубые ошибки (промахи). Их обнаруживают и исключают из расчетов (см. практическое занятие № 3).

При обнаружении грубой ошибкой этот результат исключают из расчетов и среднее значение х вычисляют заново для оставшихся достоверных результатов измерения.

Пользуясь данными таблицы или числового ряда, вычисляют практическое поле рассеяния (размах варьирования):

R = х max х min,

(4.1)

где х max, х min – максимальное и минимальное значения измеряемого параметра.

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

51

(интервалы). Число разрядов k должно быть увязано с количеством деталей. При N = 50 100 шт. k = 5 7, при N > 100 шт., k = 7 11. Для определения оптимального числа интервалов можно воспользоваться правилом Старджесса: k ≥ 1+3,3 lgN.

Когда число наблюдений N велико (например, более 200), то число интервалов приближенно можно найти по формуле

k = 4

0,75

(N 1)2

1/5 .

(4.2)

 

 

 

 

 

Число разрядов должно быть таким, чтобы цена разряда Cр = R/k была больше цены деления мерительного инструмента. Выполнение этого требования необходимо для того, чтобы уменьшить влияние погрешностей измерений. Например, для k = 7 расчетная цена разряда при R = 0,13 мм Ср = 0,13/7 = 0,0185 мм, принимаем Ср = 0,02 мм. Заметим, что округление цены разряда должно быть минимальным и допустимо в большую сторону.

Для удобства построения кривой результаты разбиений на интервалы и частоты, соответствующие этим интервалам, следует свести в таблицу (табл. 4.1).

 

 

 

 

Таблица 4.1

 

Определение частот и частостей

 

 

 

 

 

 

Номер

Интервалы, мм

Частота

Частость

интервала

 

 

mi

mi / N

от

до

 

 

 

1

19,85

19,87

3

0,03

2

19,87

19,89

16

0,16

3

19,89

19,91

22

0,22

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

7

19,97

19,99

2

0,02

 

 

 

∑ 100

 

Следующим этапом является построение кривой распределения. Она строится в координатах mi, R. Масштабы по осям выбираются произвольные, удобные для построения. По оси абсцисс

52

откладывается размах R (интервалы k), из середины интервала по вертикали откладываются соответствующие им значения чисел деталей, имеющих погрешности в пределах интервала. Полученные точки соединяются отрезками прямой (рис. 4.1).

Рис. 4.1. Полигон распределения

Практическая кривая (полигон) служит для первой приближенной оценки точности процесса и решения вопроса о выборе теоретического закона для характеристики данного распределения. Приближенной же потому, что форма практической кривой распределения зависит не только от объективных причин – характера распределения размеров, но и от случайных – числа интервалов, количества принятых для анализа деталей N. В связи с этим для объективной оценки точности обработки практическую кривую необходимо заменить теоретической кривой, изображающей вполне определенный закон распределения, описываемый математическим уравнением.

По накопленным частотам (частостям) можно построить ку-

муляту или огиву.

Для построения кумуляты из верхней границы каждого интервала на оси абсцисс восстанавливаются перпендикуляры, соответствующие по высоте накопленной частоте mi или частости mi/N, с начала ряда по данный интервал, а затем последовательно плавно со-

53

единяют вершины. Полученная кривая называется кумулятой. Она отражает характер нарастания частот (частостей) от группы к группе.

Если оси поменять местами, т.е. варианты откладывать на оси ординат, а накопленные частоты (частости) – на оси абсцисс, то построенная кривая называется огивой.

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

Закон нормального распределения (закон Гаусса) является одним из наиболее распространенных законов распределения погрешностей. Уравнение кривой нормального распределения имеет следующий вид:

y = f (x) =

 

1

 

e(xi x )2 /2 2 .

(4.3)

 

 

 

 

 

 

2

 

 

 

 

 

Функция распределения имеет вид

F (x) =

 

1

 

x

e( xi x )2 /2 2 dx.

(4.4)

 

 

 

 

 

 

 

2

 

 

График плотности нормального распределения называется нормальной кривой или кривой Гаусса (рис. 4.2).

Рис. 4.2. Распределение Гаусса

54

Отметим смысл характеристик этой кривой:

х – центр группирования, характеризует распределение размеров;

– характеризует кучность распределения размеров (погрешностей) около х ; чем меньше , тем кучнее распределяются размеры около х .

Кривая Гаусса имеет следующие особенности:

1.Кривая симметрична относительно х .

2.При xi = x кривая имеет максимум:

y =

 

1

 

0,4 ;

(4.5)

 

 

 

 

max

 

 

2

 

 

 

 

 

 

 

3. На расстоянии ± σ от вершины кривая имеет две точки перегиба А и Б, координаты которых соответственно

y

= y =

 

1

 

0,6y

 

0,24

;

(4.6)

 

 

 

 

 

А

Б

 

 

2 e

 

max

 

 

 

 

4.На расстоянии ± 3 σ от вершины кривой ее ветви так близки к оси абсцисс, что в пределах ± 3σ 99,73 % всей площади ограничивается кривой. Практически принято считать, что на расстоянии ± 3 σ от вершины кривой ее ветви пересекаются с осью абсцисс, и в этих пределах заключена вся площадь кривой, т.е. 100,0 %. Погрешность в этом случае составляет 0,27 %, что допустимо при решении многих задач производства.

5.σ – это мера рассеяния, мера точности. Кривые Гаусса при различных значениях средних квадратических отклонений представлены на рис. 4.3. На основании п.4 справедливо утверждение, что поле рассеяния

ω ≈ 6 σ.

(4.7)

При определении по данным непосредственных измерений заготовок и расчетов по формуле (2.1) погрешность определения среднего квадратического, обозначаемого в этом случае буквой S, зависит от общего количества N измеренных заготовок и в отдельных случаях весьма значительно. Учитывая это обстоятельство, для

55

предотвращения возможного появления брака целесообразно при использовании формулы (4.4) принять соотношение

σ = k S,

(4.8)

где k коэффициент, учитывающий погрешность определения среднего квадратического; S – среднее квадратическое, определяемое по формуле (2.1). Максимальная погрешность ( S) определения S выбирается по табл. 4.2.

Рис. 4.3. Нормальное распределение случайных погрешностей при различных значениях

 

 

 

 

 

Таблица 4.2

Значения максимальной погрешности S определения S

 

 

 

 

 

 

 

N, шт.

S, %

kσ

N, шт.

S, %

 

kσ

 

 

 

 

 

 

 

25

42,4

1,40

200

15,0

 

1,15

50

30,0

1,30

300

12,2

 

1,12

75

25,0

1,25

400

10,6

 

1,11

100

21,2

1,20

500

10,0

 

1,10

В тех случаях, когда поле рассеяния параметров (размеров) превосходит поле допуска (ω > δ), условие обработки без брака не выполняется и брак является возможным.

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

56

распределения Гаусса, оценка точности принимается с погрешностью не более 0,27 %. При этом считаем, что все детали партии имеют действительные размеры в пределах поля рассеяния:

σ = xmax xmin,

(4.9)

где xmax, xmin – максимальное и минимальное значения параметра (размера). При этом площадь, ограниченная кривой нормального распределения и осью абсцисс, равна единице и определяет 100 % заготовок партии. Площадь заштрихованных участков (рис. 4.4) представляет собой количество деталей, выходящих по своим размерам за пределы допуска.

Для определения количества годных деталей необходимо найти площадь, ограниченную кривой и осью абсцисс на длине, равной допуску δ. При симметричном расположении поля рассеяния относительно поля допуска следует найти значение интервала, определяющего половину площади, ограниченной кривой Гаусса и

абсциссой х1 (х2).

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

F (x) = x

y dx =

 

1

x

ex2 /2 2 dx.

(4.10)

 

 

 

 

 

 

2

 

 

Рис. 4.4. К определению количества годных деталей

57

Для случая, когда x = 0, =1, распределение называют

стандартным и функция распределения (4.4) имеет следующий вид

(рис. 4.5):

F (x) =

 

1

 

x

ex2 /2dx.

(4.11)

 

 

 

 

 

 

 

2

 

 

Таким образом, если случайная величина Х следует закону нормального распределения, то вероятность появления случайной погрешности определяется площадью, ограниченной кривой f (x) и ее частью и осью абсцисс:

 

 

1

 

x2

 

P x1 x x2 =

 

 

x ex2 /2 2 dx.

(4.12)

 

 

 

 

2

 

 

1

 

Подынтегральное значение есть элемент вероятности, равный площади прямоугольника с основанием dx и абсциссами x1 и x2, называемыми квантилями.

Произведем замену переменной: t = x / , dx = dt:

 

 

1

 

t2

et2 /2dx.

(4.13)

P x1 x x2 =

 

 

t

 

 

 

 

2

 

 

 

 

1

 

 

Рис. 4.5. Функция распределения F(x) и функция Лапласа Ф(x)

Представим правую часть в виде суммы двух интегралов:

58

 

 

1

0

 

 

 

 

 

1

 

t2

P x1 x x2 =

 

 

t

et2 /2dt +

 

 

0 et2 /2dt;

 

 

 

 

 

 

 

2

2

 

 

 

 

1

 

 

 

 

 

 

 

 

Интеграл вида

 

 

 

 

 

 

 

 

 

 

Ф(t) =

 

1

 

 

t

et2 /2dt;

(4.14)

 

 

 

 

 

 

 

 

 

 

2 0

 

 

 

 

 

носит название нормальной функции Лапласа. Значения этого интеграла сведены в табл. 4.3. Таким образом, указанная вероятность (4.11) сводится к разности нормальных функций Лапласа:

Р { x1 < x < x2 } = Ф (t2) – Ф (t1).

(4.15)

Расчет количества годных деталей сводится к установлению величины t и определению Ф(t) по табл. 4.3 с последующим пересчетом полученных величин в проценты или в число штук изделий.

В общем случае, когда x 0 , имеем следующую вероятность появления случайных погрешностей:

Р{x

x x }

 

x2

f (x)dx = Ф

 

x2

x

Ф

x1

x

.

(4.16)

=

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

 

 

Отметим свойства функции Лапласа: Ф(0) = 0; Ф(–х) = –Ф(х) (функция нечетная); Ф( ) = 1/2. Из рис. 4.5 видно, что кривые F(х) и Ф(x) эквидистантны.

Если в равенстве (4.13) положить х1 = – , то

 

1

x

x

 

 

P (x x2 ) =

2

+ Ф

2

 

 

,

(4.17)

 

 

 

 

 

 

 

 

так как Ф(– ) = –Ф( )= –1/2. Положив в соотношении (4.16) х2 = , находим

P (x x ) = 1

Ф

x1

x

.

(4.18)

 

 

1

2

 

 

 

 

 

 

 

 

59

60

Таблица 4.3

 

1

t

t2

Значения функции Лапласа при разных значениях t Ф(t)=

e

 

dt

2

 

 

 

0

 

 

 

 

 

 

t

Ф(t)

t

Ф(t)

t

Ф(t)

t

Ф(t)

t

Ф(t)

t

Ф(t)

0,00

0,0000

0,50

0,1915

1,00

0,3413

1,50

0,4332

2,00

0,4772

3,00

0,49865

0,01

0,0040

0,51

0,1950

1,01

0,3438

1,51

0,4345

2,02

0,4783

3,20

0,49931

0,02

0,0080

0,52

0,1985

1,02

0,3461

1,52

0,4357

2,04

0,4793

3,40

0,49966

0,03

0,0120

0,53

0,2019

1,03

0,3485

1,53

0,4370

2,06

0,4803

3,60

0,499841

0,04

0,0160

0,54

0,2054

1,04

0,3508

1,54

0,4382

2,08

0,4812

3,80

0,499928

0,05

0,0199

0,55

0,2088

1,05

0,3531

1,55

0,4394

2,10

0,4821

4,00

0,499968

0,06

0,0239

0,56

0,2123

1,06

0,3554

1,56

0,4406

2,12

0,4830

4,50

0,499997

0,08

0,0319

0,58

0,2190

1,08

0,3599

1,58

0,4429

2,16

0,4846

 

 

0,09

0,0359

0,59

0,2224

1,09

0,3621

1,59

0,4441

2,18

0,4854

 

 

0,10

0,0398

0,60

0,2257

1,10

0,3643

1,60

0,4452

2,20

0,4861

 

 

0,11

0,0438

0,61

0,2291

1,11

0,3665

1,61

0,4463

2,22

0,4868

 

 

0,12

0,0478

0,62

0,2324

1,12

0,3686

1,62

0,4474

2,24

0,4875

 

 

0,13

0,0517

0,63

0,2357

1,13

0,3708

1,63

0,4484

2,26

0,4881

 

 

0,14

0,0557

0,64

0,2389

1,14

0,3729

1,64

0,4495

2,28

0,4887

 

 

0,15

0,0596

0,65

0,2422

1,15

0,3749

1,65

0,4505

2,30

0,4893

 

 

0,16

0,0636

0,66

0,2454

1,16

0,3770

1,66

0,4515

2,32

0,4898

 

 

0,17

0,0675

0,67

0,2486

1,17

0,3790

1,67

0,4525

2,34

0,4904

 

 

0,18

0,0714

0,68

0,2517

1,18

0,3810

1,68

0,4535

2,36

0,4909

 

 

0,19

0,0753

0,69

0,2549

1,19

0,3830

1,69

0,4545

2,38

0,4913