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

книги из ГПНТБ / Головинский В.В. Статистические методы регулирования и контроля качества. Расчет оптимальных вариантов

.pdf
Скачиваний:
15
Добавлен:
24.10.2023
Размер:
10.75 Mб
Скачать

На основании теоремы умножения вероятностей и пользуясь обычными в таких случаях способами комбинаторики, можно за­ писать

& (Щ = 0) = (1 — р')п;

 

& ( m i = l ) = C ' np ' ( l -

р')»"1-

р' (1 - РУ1- 1; (3• 16)

^ (mi = b) = Cba (p')b(1 -

р ) п~ ь = 6|(; 1

6)| ( рУ (1 - РУ " .

Пользуясь введенными обозначениями и соотношениями, можно

предложить алгоритм вычисления L?. 2 (v), который в сокращен­ ной записи выглядит следующим образом:

- р т ~ а ;

(злу»

т — О

 

р' = Ф (ик----п); п = 2Ь-\-\.

 

Ниже изложен тот же алгоритм в развернутом виде.

 

Алгоритм вычисления оперативной характеристики планов Г. 2 при п — 5

I.

Входные данные:

значений о,- отклонений у. н. о; i = /т щ

1.

Последовательность

— 3, — 2, — 1, 0, 1, 2, 3 ............

(щах-

2.uK_. ик + .

3.п.

4.Ф (О-

5 . С}, = 5: С2„ = 10.

6. Л.

И. Операции (при всех i):

1. t . =

u ^ - v r

2-

Pi =

Ф Hi).

3-

p2i=PiPi-

4-

Pi =

1— Pi-

5.

q2i =

qi ql -

6 . 9? =

9?9(.

7. q\ = qzqL-

8 . 9-=<?4--

9. a. =

Cnpi .

70

10. ff>^=a.qj.

1 2.

?><.2>= £м?)>.

13.

с£ = 9? +

<?f>.

14.

L ? . 2 (o) =

c, + s f > -

15- 4 . 2 («K_ - i/ i) = L ; 2 (ик_ + £Л).

16.

d. — LT 2

(»,.) -f- I r 2 (“ ,•)■

17.

L (o,-) = d[ — 1.

Частная оперативная

характеристика слева Lj\ 2 (и) является

функцией распределения

А

вероятностей выборочной

ме­

дианы с центром в ик~.

 

 

 

Так как выборочная

медиана х является при симметричном

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

качества

х несмещенной оценкой

для

м. о. х, точки ик- и ик+ соответствуют, как и в планах класса А, точкам равновесия:

Lp. 2 (ик~) =

Ly. 2

(ык+) =

0,5.

 

(3.18)

Если при одинаковом объеме выборки план Г.2 меняется только

в отношении границ ик-

и

« к+,

новые

Lf, 2. 1 (у)

и Lr.

2, 1 (и)

можно получить, сместив

старые L r . 2 ( у)

и Lr.2 ( у)

таким

обра­

зом, чтобы значения их, равные 0,5, совпадали с новыми ик- и и'к + (так же, как при планах класса А). На этом основании и с уче­

том (3.13) в алгоритме (см. операцию 15) Lr) 2 (у) вычисляется на основании Lr. 2 (у).

3 .5 . О С Т А Л Ь Н Ы Е РА ЗН О ВИ Д Н О СТИ П Л А Н О В К Л А С С А Г

Планы класса Г.З, первоначально именовавшиеся методом ка­ либров распределения и позже — методом группировки, были раз­ работаны и впервые опубликованы в 1944 г. В. Н. Гостевым (СССР) и несколько позже, независимо от него Л. Типпеттом (Англия). Так же, как уже рассмотренные разновидности Г.1 и Г.2, в планах Г.З решающая функция определена в трехмерном

пространстве на множестве целочисленных векторов т = т ъ т 2, т 3, где nij — число наблюденных выборочных значений признака качества х, попавших в интервал / = 1 , 2 , 3.

Единственным отличием вероятностных схем планов Г.З от схем планов Г.1 и Г.2 является то, что в соответствии с решающим

71

правилом решение не вмешиваться в технологический процесс, принимается, если

т 3 —

т 1

с; т 1 m3s^ c ,

(3.19)

где с — фиксированное

для

каждого объема выборки

п число.

Это отличие не существенно ни в отношении способов измере­ ния (при планах Г.1 и Г.2 применим и шкальный, и предельный инструмент, как при планах Г.З), ни в отношении оперативной характеристики в том особом случае, на который рассчитаны пла­ ны Г.З (в этом мы убедимся позже). Существенное отличие пла­ нов Г.З от планов Г.1 и Г.2 заключается в условии, не имеющем отношения к математической статистике *. Условие это заклю­ чается в следующем. В то время, как границы интервалов груп­ пировки при планах Г.1 и Г.2 можно выбрать в зависимости от экономических или технических соображений, аналогичные гра­

ницы планов Г.З являются

функциями

единственной

величины,

а именно — технического допуска.

 

 

 

Обозначим границу между первым и вторым интервалами

3)

и границу между вторым и третьим интервалами

При

пла­

нах Г.З **, если допуск б ^

6 стЛ.,

 

 

 

ЩГ. 3)

б

< г . 3)

б

 

(3.20)

4а*

4а* *

 

К

ик+

 

Если б - < 6 0*, то считается,

что операция вышла из состояния

контроля (доказательством этого служит соотношение т 1 -f- /п2 >

>

с2). Тогда, прежде чем проверять и уточнять настройку, надо

устранить ненормальность.

Переходя к вычислению оперативной характеристики Lr. 3 (у)

планов Г.З, следует учесть, что вследствие несовместимости нару­ шения соотношений (3.9) в данном случае, как и в остальных рас­

смотренных, выполняется

равенство

 

L r.3(v) =

L7.3 (v) + L t.з (о) — 1,

(3.21)

где LF. з (у) — частная оперативная характеристика слева, рав­ ная вероятности, что при данном v не будет нарушено соотноше­

ние (3.19): т х tn3

с; Lr. з (у) — частная

оперативная харак­

теристика справа, равная вероятности, что т

3 т х ^ с.

Прежде чем изложить алгоритм вычисления LF з (v) и Lr. з (у), введем следующие обозначения для группировок выборочных дан­ ных в общем случае. Общее число групп обозначим через J, каж­ дую группу — через G,-, / = 1, 2, . . ., /, границы между интер-

* Л. Типпетт предложил аналогичный метод для классической контрольной карты У. Шыохарта, на которой границы не зависят от допуска. Поэтому ска­ занное здесь не имеет отношения к предложению Л. Типпетта.

* * Напомним, что началом координат всюду принят заданный уровень на­ стройки SB. В планах Г.З. предполагается, что SB соответствует середине поля допуска.

72

валами,

соответствующими

группам, через

и\1

j — 1 ,

2 , . .

., /,

причем

Gx, G2, . . .,

Gj

состоят из

« выб

соответственно

в интер­

валах.

Цвыб

«к I

 

^

Нвыб

j

• •

Uk

=

1,

Нвы<5*,

ЧИ­

СЛО значения

мвыб,

попавших в группу Gy,

my, j

2, . .

., /;

вектор

группировки т

— т 1, т 2,

. . ., т у

вероятность того,

что

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

Gy,

Ру, / =

1, 2 , . . . , / .

В

условиях

рассматриваемой

задачи

 

 

 

р/ =

ф (и,(/> — и) — Ф (а */"1» — и),

 

 

(3.22)

где Ф (/). — функция гауссова распределения; вероятность век­ тора группировки с компонентами ту, на основании теоремы умножения и пользуясь известными способами комбинаторики, можно записать в виде

 

 

&

(mi, m 2, ..., my,

..., m j )

=

 

 

 

 

_________л!________

pT 'P i\

. . . . P?J .

(3.23)

 

 

mL! . .

п и ! . . m.j\, . . m j\

 

 

 

 

 

При

вычислении

Lt. з (у) число групп

J

= 3. На основании

теоремы

сложения вероятностей

можно вычислить

L r)3

(y) как

сумму

вероятностей

векторов,

удовлетворяющих

требованию

т 3 mi < с.

В сокращенной записи алгоритм Lp(3 (y)

таков:

 

 

т , = 0

 

 

 

 

 

 

 

•*

»*— 3

л!

 

рпирп-пи-т,рт,

(3.24)

 

 

 

 

 

 

 

! (« — /лх — т 3)! т 3!

Приравняв

середине поля допуска, для вычисления Lp. з (v)

можно

воспользоваться соотношением

 

 

 

 

 

 

 

L7.3 (v) =

L t.3

(— v).

 

 

 

(3.25)

Разновидностью Г.4 планов класса группировки является комбинация метода медиан (Г.2) и метода крайних значений (Г.1). Когда на заводе электроосветительной арматуры (Москва) в 1947 г. Московским инженерно-экономическим институтом был впервые применен так называемый метод медиан, и позже, когда он был усовершенствован на заводе «Калибр», возник вопрос — как быть, если границы выборочной медианы не нарушены, а в выборке от­ мечено нарушение допуска? Вопрос этот в свое время не был до конца исследован. На практике в инструкциях на некоторых заво­ дах предусматривалось, что наличие в выборке брака даже без нарушения границ регулирования выборочной медианы влечет вмешательство в технологический процесс, что приводит к ино­ му плану, чем Г,2, Обобщая это решение в том смысле, что про­

73

межуток между границами можно по усмотрению уменьшить сравнительно с допуском, приходим к планам Г.4.

В процедурно-вычислительном отношении планы Г.4 можно описать следующим образом. На основании выборки нечетного

объема п = 2b + 1 вычисляется выборочная медиана х. Решение о невмешательстве в технологический процесс принимается, если

хк-

4

«с

Хц+,

где л'к-, хк+ — границы

для

выборочной

медианы; хк- х,<+ —

границы для крайних выборочных значений.

В статистических понятиях планы Г.4 формулируются сле­

дующим образом.

2b +

1 независимых вероятностных ис­

1. Выполняется п =

пытаний, каждое из которых состоит в изготовлении и измерении

одного изделия с результатом

в виде отклонения

ишб =

-----

выборочного значения хвЬ1б от

заданного уровня

а х

настройки

S6.

2.Группировкой значений ивый определяются компоненты т и

т2, rn3, mit пц вектора группировки т . Границами интервалов группировки являются

и

~(2) _ а'к-

&

Ии—— "

ГГ

К-

* '

(1

„ ( 2 )

(3.26)

к+ =

ик+ =

3.Решающее правило: решение о невмешательстве в техноло

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

вектор

т =

т 1,

т 2, т а,

т 4, т ъ обладает следующими

свойствами:

 

 

 

 

mL+ m2 <

т 4 +

т 5

b\

т х =

т ъ =

0,

(3.27)

При вычислении

оперативной

характеристики

Lr. 4

(v) пла­

нов Г.4 предполагается, что интервал («£4),

и ^ ) достаточно

ве­

лик, чтобы вероятностью совмещения неравенств т 1

)> 0

и т 5>

О

в одной выборке можно было пренебречь. Поэтому при вычисле­

нии

L r .4 (v) пользуемся соотношением

 

 

Lr. 4 (v) — L r.4 (у) +

Lt. 4 (у) — 1,

(3.28)

где

Lp. 4 (v) — частная

оперативная

характеристика слева,

рав­

ная

суммарной вероятности

векторов

т , для которых

 

 

т х +

т 2 с

Ь; /пх

=

0 ;

(3.29)

74

Lf. 4 (v) — частная оперативная характеристика справа, равная суммарной вероятности векторов т , для которых

 

 

 

 

ш4

-|- т 5 ^

Ь\

т 5

=

0.

(3.30)

Алгоритм

вычисления

Z.F. 4 (0 )

и

 

4

(v) имеет

в основе

соотношения

(3.29)

и (3.30);

 

 

 

 

 

^

<(“) “

t

0

 

to +

р<+

(3.31)

 

 

 

ш2=

 

 

 

 

 

 

где т и р

— см.

(3.22) и

(3.23).

 

выбор границ

интервалов

Планы

Г.4 допускают

свободный

группировки в соответствии с экономическими и технологическими

соображениями с

единственным ограничением

 

 

 

 

SP {т х= \ ,

/п2 + т 3 +

т 4 =

п — 2,

т ъ — 1)^ е ,

(3.32)

где

е — пренебрежимо

малая

вероятность.

 

 

При

п — 5

для

соотношения

(3.32)

достаточно,

чтобы

4 5) -

4 1) > 2,5.

 

 

 

 

 

 

 

Он

Метод индивидуальных значений обозначается здесь через Г.5.

был

разработан

в

1947 г.

и внедрен

на одном из

заводов

Л.

М.

Брагинским.

Основой этих

планов

является оригиналь­

ная идея использования для проверки уровня настройки по возможности большего числа членов вариационного ряда. Способ

применяется при условии, что допуск 6

6 о.(.

В процедурно-вычислительном отношении планы Г.5 можно описать следующим образом. На контрольной карте наносятся две пары границ для индивидуальных (наблюденных) выборочных значений хпЫб признака качества: внешние нижняя и верхняя

границы 4 -> 4 +( и внутренние нижняя и верхняя границы

х<2> х<2>

ЛК-> лк+- При выборке фиксированного объема п — 5 решение не вме­

шиваться в технологический процесс принимается, если за внеш­ ними границами 4 -\ 4 + не окажется ни одного наблюдения хвыб,

а за каждой из внутренних границ 4 -> 4 + — не больше одного наблюдения. Положение границ фиксировано на поле допуска, причем внешние границы совпадают с границами поля допуска, внутренние смещены к середине поля допуска на 1 / 1 0 последнего.

В статистических понятиях планы Г.5 формулируются сле­ дующим образом:

1. Выполняется п — 5 независимых вероятностных испыта­ ний, каждое из которых состоит в изготовлении и измерении одного

изделия с результатом в виде отклонения ивыб — выб-~~— .

Ох

2 . Группировкой значений пвыб определяются компоненты т х, т 2 > tn3, /п4, т ь вектора группировки т .

75

Считая

равным Середине поля допуска 6 ,

границами интер­

валов группировки

являются

 

 

 

 

 

 

ul" =

— 0 ,5 6 .

~ 0 ''

, ( 3

)

0,46 .

 

 

и ‘ 2 ) -

° х

 

 

Щ

 

° х

 

 

 

 

 

(4)

0,5

6

 

 

 

 

 

 

к ~

о *

 

 

 

 

3.Решающее правило: решение о невмешательстве в техно

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

группировки т —

= т х, т 2, т 3,

т 4,

т 5 обладает следующими свойствами:

т

1 =

т ъ = 0; т 2 < 2 ; т 4 < 2 .

(3.33)

Так как по условию применимости планов Г.5 нарушение верхних и нижних границ в одной выборке практически несов­ местимо, для вычисления оперативной характеристики можно

пользоваться соотношением

 

 

 

^ г.б (у) = C f. 5 (v) +

L t .s (t>) - 1,

(3.34)

где

L^,5 (v) — частная оперативная

характеристика

слева, рав­

ная

суммарной вероятности векторов группировки т

с т г = 0;

/«.,

<(2; Lt. 5 (н) — частная оперативная характеристика справа,

равная суммарной вероятности векторов группировки т

с rni <(

<

2, т ъ = 0.

 

 

 

 

 

На основании соотношений (3.33) для планов Г.5 алгоритм

вычисления в

сокращенной записи

сводится к следующему:

 

-г.5м

- 2

п\

(pi +рг-\-рэ)я т ‘ рТ'.

(3.35)

 

-(п — т 4)\ т А!

 

 

ШЯ4 =0

 

 

 

При планах Г.5 в том виде, в каком они получили практиче­ ское применение, оперативная характеристика зависит только от допуска, и выбор варианта по технологическим или экономи­ ческим соображениям невозможен.

3.6. АППРОКСИМАЦИЯ ПЛАНОВ КЛАССА Г НОРМАЛЬНЫМ РАСПРЕДЕЛЕНИЕМ

В гл. 2 в связи с вопросом экономической оптимизации СРК отмечено, что многочисленность методов статистического регули­ рования, отличающихся в процедурно-вычислительном отношении, и сложность формул их оперативных характеристик крайне за­ трудняет выбор лучшего варианта для того или иного класса технологических операций. Там же упоминалось, что это затруд­ нение можно снять, так как оперативные характеристики всех рассмотренных выше планов можно приближенно выразить с помощью одной и той же формулы оперативной характеристики планов А с различными значениями трех параметров %, у~ , у+ для каждого из планов.

76

Следует добавить также на основании сказанного в этой главе, что достаточно найти общие формулы для частных оперативных характеристик планов класса группировки, так как переход от них к полным оперативным характеристикам не представляет никаких затруднений. Частные оперативные характеристики планов А имеют только два параметра (X, у~ или у+) и, как из­ вестно (см. и. 3.2), совпадают с функцией нормального распре­

деления

вероятностей [Lr

(о)] или

выражаются через него

[L+ (w)l.

Речь

идет, таким

образом,

об аппроксимации

частных

оперативных

характеристик

планов типа

группировки

функции

нормального

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

 

 

 

 

Отыскание для данной частной оперативной характеристики

плана Г

L f

(у)

наиболее близкой к

ней

оперативной

характе­

ристики плана A L J (у) можно выполнить различными способами, причем обоснование последних [включая вычисление линейной

регрессии для выпрямленной характеристики LF (у)] в большей или меньшей степени включает произвольные постулаты и носит интуитивный характер. Как увидим, это ие существенно для выво­ дов, которые будут представлены ниже. Не вдаваясь в элемен­ тарные мотивировки и детальные пояснения, перейдем к изложе­ нию использованного способа, который можно назвать аппрокси­ мацией через функцию нормального распределения по двум точкам. Этот способ по идее совпадает с выпрямлением кривых накопленных частостей на «вероятностной бумаге», отличаясь от него большей объективностью и удобством (по крайней мере, при отсутствии хорошей вероятностной сетки).

Рекомендуемая аппроксимация функции L~ (у) функцией нормального распределения Ф (t) выполняется с помощью таб­

лицы

функции

ф (/), обратной относительно Ф (t)

(см. табл. II

приложения 1).

Если левая частная оперативная характеристика

L~ (у) вычислена для интервалированных значений

у ,-, то пара­

метр крутизны X находим следующим образом:

 

 

а)

взять два

значения LF ( у)

при у (1) и у (2) таких,

чтобы

L w =

L f ( у (1))

было ближайшим к 0,158 значение LF ( у);

L (_) =

= LF ( у (2)) — ближайшим к 0,841

значение LF ( у);

 

 

б)

определить по таблице обратной функции ф (t) нормального

распределения (см. приложение

II) два значения:

 

 

 

 

фФ = ф (L*1));

ф<2>=,ф(А<2>);

 

 

в)

вычислить

 

 

 

 

 

 

 

,

уЫ -уЫ

 

(3.36)

 

 

 

ф(2! — т])Ф ’

 

 

 

 

 

 

г)

определить параметр

положения слева

 

 

 

 

у- =

у( 0

— ф'А,.

 

(3.37)

77

Вычисления для правой частотной оперативной характери­ стики аналогичны. Точки аппроксимации у(1) и г/2) при /_Г(и(1)) =

=

0,158

и LJr

(u(2)) =

0,841 соответствуют значениям ф (/) =

=

—1; ф

(t) =

1. Они

выбраны такими потому, что находятся

в интервале значении v, над которыми погрешности аппроксима­ ции могут привести к ощутимым отклонениям оценок эффектив­ ности статистического регулирования в отличие от периферийных участков.

Втабл. 3 и 4 частные оперативные характеристики L? (v) планов Г.1, Г.2, Г.З, Г.4 сопоставлены с ближайшими к ним пла­ нами А.

Вприведенных ниже примерах приняты объемы выборок и границы, указанные в табл. 3. В той же табл. 3 в гр. 8 показано

отношение

допуска к среднему квадратическому

отклонению,

а в гр. 5,

б и 7 — параметры аппроксимирующей

оперативной

характеристики плана А. В гр. 9 проставлены отношения объема выборки па плана А к объему выборки пг планов Г при равной

Условныечения ­обозна

I

Г.1

Г.2

Г.З

Таблица 3

 

 

Параметры оперативных характеристик

 

 

 

 

и статистическая

полезность одного наблюдения

 

 

 

 

 

при способах группировки

 

 

 

Статистическаяпо­ одноголезность наблюденияН

 

Обычные

выборкиОбъем п

Удалениеграниц интерваловгруппи­ серединыотровки

допускаполяв до­ допускалях

А,

V-

Y+

вДопусках

 

названия

 

 

 

 

 

 

 

2

3

А

 

5

6

7

8

9

Способ

крайних

10

—0,5;

0,5 1,75

-1 ,4 7

1,47

6

0,31

значений

 

 

 

 

 

 

 

 

 

Способ медианы

5

1 .

1

1,87

— 2,00

2,00

6

0,70

3

3

 

 

 

 

 

 

 

 

 

Способ

калиб­

6

1 .

1

2,05

1,31

1,31

6

0,70

4

4

ров

распределе­

 

 

 

 

 

 

ния;

с =

2

 

 

 

 

 

 

 

 

Г.З

То же

 

6

1 .

1

2,0 0

— 2,10

2,10

9

0,67

 

4

4

 

 

 

 

 

 

 

 

 

Г. 4

Способ

медиан

5

о , 5

; ф ; 1,71

— 1,67

1,67

6

0,58

и

крайних

значе­

 

 

 

 

 

 

 

 

ний

 

 

l b

0,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

Способ

средней

4 '

 

2,00

1,30

1,30

6

1,00

арифметической

 

 

 

 

 

 

 

 

78

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 4

 

 

Сопоставление оперативных характеристик L~ (и), вычисленных по своим формулам

 

 

 

 

 

 

с вычисленными по параметрам у - и Я.

 

 

 

 

 

 

 

 

Оперативные характеристики L ~ (v)

при способах статистического регулирования

 

 

 

 

г л

 

 

Г . З при допуске,

Г.З при допуске,

Г .4 при допуске,

А

i

 

 

 

Г . 2

равном бегv

равном 9ал

равном ба

vi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

по фор­

по пара­

по фор­

по п ара­

по фор­

по п ара­

по фор­

по п ара­

по фор­

по п ар а ­

по п ара­

 

 

мулам

метрам

мулам

метрам

мулам

метрам

мулам

метрам

мулам

метрам

метрам

1

2

3

4

5

6

7

8

9

10

11

12

13

— 3

- 1 , 5

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

— 2

— 1,0

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

— 1

— 0 ,5

0,9980

0,9998

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

1,0000

0 ,9 9 98

0

0

0,9871

0,9952

1,0000

1,0000

0,9959

0,9972

1,0000

1,0000

0,9934

0 ,9978

0,9953

1

0 ,5

0,9382

0 ,9564

0,9974

0 ,997 5

0,9486

0 ,9505

1,0000

1,0000

0,9 674

0,9770

0,9 45 2

2

1,0

0,7948

0,7948

0,9688

0 ,9683

0,7334

0,7340

0,9851

0,9821

0,8722

0,8729

0,7257

3

1,5

0,5000

0,4880

0,8256

0 ,8256

0,3463

0,3372

0,8643

0,8643

0,6276

0,6123

0,3 4 4 6

4

2,0

0,1778

0,1778

0,5000

0,5000

0,0775

0,0778

0^5421

0,5398

0,2835

0,2843

0,0 80 8

5

2 ,5

0,0250

0,0351

0 ,1746

0,1746

0,0074

0,0071

0,1809

0,1841

0,0632

0,0738

0 ,0082

6

3 ,0

0 ,0036

0,0000

0,0310

0,0307

0,0003

0,0002

0,0271

0,0288

0,0058

0,0104

0 ,0003

7

3 ,5

0,0000

0,0000

0,0027

0,0025

0,0000

0,0000

0,0015

0,0013

0,0002

0,0009

0,0 00 0

8

4 ,0

0,0000

0,0000

0,0000

0,0000

0,0000

0,0000

0,0000

0,0000

0,0000

0,0000

0,0000

-j

га

Соседние файлы в папке книги из ГПНТБ