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

книги из ГПНТБ / Синдлер Ю.Б. Метод двухступенчатого статистического анализа и его приложения в технике

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

Г л a n a 6

ДВУХСТУПЕНЧАТЫЕ ПРОЦЕДУРЫ КОНТРОЛЯ КАЧЕСТВА И НАДЕЖНОСТИ ПРОДУКЦИИ

6.1. Предварительные замечания

Как известно, вопросы организации статистического приемоч­ ного контроля имеют важное значение для повышения качества и надежности продукции. Эти вопросы давно вызывают интерес исследователей. Имеется обширная литература в этой области. Библиография, приведенная в коице настоящей главы, содержит далеко не полный перечень наиболее существенных работ.

В настоящей работе мы рассмотрим некоторые типичные зада­ чи статистического контроля. Мы будем иметь дело с двумя типич­ ными объектами контроля в производстве:

партией сравнительно несложных изделий;

одиночным сложным изделием.

Контроль партий в пашем рассмотрении является выборочным. Рассматривая вопросы контроля партии продукции, мы примем одно весьма существенное допущение: контроль каждого изделия, выбранного для проверки, предполагается достоверным. Други­ ми словами, проверив изделие, мы можем достоверно судить о его качестве.

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

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

а) каждое изделие при проверке может быть квалифицировано либо как «годное», либо как «дефектное», такой вид называют конт­

ролем по альтернативному признаку.

V б) при проверке каждого изделия регистрируются и учиты­ ваются точные значения некоторых параметров; такой вид конт­

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

При выборочном контроле партий изделий могут быть при­ няты различные решения, в том числе:

решение о том, что партия изделий годная и ее следует принять;

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

5* 131

— решение о том, что следует провести сплошную, 100-про­ центную проверку партии и все дефектные изделия изъять.

Последнее

решение, естественно, может

быть принято только

в том случае,

если контроль не является

разрушительным.

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

Факт принятия любых решений, кроме решения о хорошем ка­ честве проверяемого изделия, можно рассматривать как сигнал о том, что произошла «разладка» производственного процесса и требуется принять необходимые меры для его наладки.

В настоящей книге, рассматривая процедуры контроля, мы будем опираться на теорию проверки двух статистических гипо­ тез, изложенную в предыдущих главах. Гипотезой На будем считать гипотезу о хорошем качестве контролируемого объекта, а гипотезой Нг — гипотезу о его плохом качестве. Более точные определения этих понятий мы приведем в каждой конкретной задаче.

Как мы увидим ниже, для весьма широкого круга задач ста­ тистического контроля данная теория является вполне адекват­ ной. Существует, однако, много задач контроля, точно не уклады­ вающихся в рамки этой теории, но близких к ней. Такой задачей является, например, задача контроля партии изделий, когда в результате контроля принимается одно из двух решений: ре­ шение о приемке партии и решение о сплошной проверке. Этот случай мы рассмотрим несколько ниже.

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

Впрактике контроля часто бывает, что объем выборок мал по сравнению с объемом партии, а возможными окончательными решениями при контроле являются решения о приемке и о бра­

ковке партии. Задачи для этого случая мы рассмотрим в разд. 6 . 3 - 6 . 5 .

В разд. 6.3 и 6.4 будет приведена задача контроля доли дефект­ ных изделий в партии. Задача решается в предположении, что число дефектных изделий, которое мы обнаруживаем в выборках, имеет пуассоновское или биномиальное распределение. При та­ ком допущении постановка задачи является приближенной. Эти распределения приближенно описывают закон распределения чис­ ла дефектных изделий, который, как известно [1—3], является ги­ пергеометрическим. Точная постановка задачи контроля доли де­ фектных изделий дана в разд. 6.6.

Вразд. 6.5 рассматриваются задача контроля доли дефектных изделий по количественному признаку и задача контроля одно­ родности изделий в партии.

Вразд. 6.7 решается задача контроля надежности сложных изделий по интенсивности потока отказов. В разд. 6.8 проведено сопоставление одного типа двухступенчатых процедур с опти­ мальной процедурой.

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

Рассмотрим, например, задачу контроля партии, содержащей N изделий, в случае, когда окончательным может быть решение о приемке партии или решение о сплошной ее проверке. В приведен­ ной выше теории мы определили показатели а и Е0 отдельно друг от друга. Напомним, что а — вероятность ошибки первого рода, Е0 — средний объем выборки при гипотезе Н0. В рассматривае­ мом случае связь между этими показателями имеет специфический характер. Ошибка первого рода здесь означает, что хорошую пар­ тию подвергают сплошной проверке. Ущерб, который мы несем, совершив ошибку первого рода, обусловлен тем, что мы должны произвести большие затраты на контроль. Средние затраты на контроль могут резко возрасти при увеличении а. Чтобы убедиться в этом, определим среднее число проверяемых изделий. Пусть на первой ступени производятся испытания выборки объема пх и пусть X — совокупность полученных данных. Конкретный смысл X уточнять не будем. Для определенности будем считать, что X может принимать дискретное множество значений, которое обозна­ чим Uх- Пусть щ — объем выборки на второй ступени и п2 (X) — детерминированное промежуточно^ правило. Допустим, что нам известны распределения Рв (X, Y) ж Рх (X, Y) при хорошем и пло­ хом качестве партии. Организуя контроль, мы стремимся к тому, чтобы хорошую партию принять, а плохую — подвергнуть сплош­ ной проверке. Пусть а у (X, 2) условная вероятность принятия решения о сплошной проверке при условии, что партия хорошая, а

на первой ступени получено данное значение X.

Непосредственно

при контроле мы проводим проверку пх + гс2 (X)

изделий. Кроме

того, с вероятностью а у (X, п„) (когда партия хорошая и X фикси­

ровано) мы проводим проверку оставшихся N — пх

— п« (X)

из­

делий.

 

 

 

 

 

 

 

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

фиксированном

X и при"~хорошем качестве партий

равно

+

+ пг

(X)

+ (N

— пх — пг

(X)) а у (X, п2).

Усреднив

это выраже- -

ние

по

всем

X е Ux,

получим среднее

число

проверяемых

изделий при хорошем качестве партии, которое мы обозначим через Е0:-

Е0 =

щ.+

2

n2(X)P0(X)

+

( N - n i ) a -

 

 

хеих

 

 

 

-

2

пг(Х)ау(Х,щ)Р0(Х),

 

(6.1)

где Ра (X) — распределение X при хорошем качестве партии. Вероятность ошибки первого рода можно выразить в виде

*=

2

* У ( Х . гк)Ро(Х).

(6 . 2)

 

хеих

 

 

Первые два слагаемых в правой части (6.1) определяют средний объем испытаний, проводимых непосредственно в процессе конт­ роля, Е0. Последние два слагаемых определяют средний объем до­ полнительных испытаний (при сплошной проверке). Обозначим

эту величину

символом

Ея.

Таким

образом, E Q

— Е0 -f-

Ея.

Определим пределы изменения величины Ея.

Пусть

есть

максимальное значение величины тг2. Тогда

 

 

 

г пх — па,) ос < £ д

^

(N — щ) ос.

 

 

(6.3)

Если

максимальное

значение

суммарного

объема

выборок

Щ + пш

невелико (пг +

гг2 м<^ N),

то величину

E R можно

приб­

лиженно оценить величиной Na. Тогда

 

 

 

E0^E0

+

Na.

 

 

 

 

 

 

Таким образом, показатель Е0 можно с некоторым приближе­ нием представить как линейную форму показателей а и Е0.

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

макс {я (s) | Е0

(s) <

Ъ0}.

(6.4)

С некоторым приближением ее можно заменить задачей

 

макс {я (s) | Е0

(s) +

Na (s) < 6 0 } .

(6.5)

Здесь b„ — заданная величина. Напомним, что я обозначает ве­ роятность правильного решения о плохом качестве партии. Функ­ ция Лагранжа для этой задачи представляет собой линейную фор­ му показателей я, а и Е0:

L = я - X (Е0 + Na) + %Ъ0.

(6.6)

Напомним, что функция Лагранжа для известной нам уже задачи

макс{іг(*)|а(в)<йв > / ? 0 ( . 9 )<6 0 }

(6.7)

может быть представлена в виде линейной формы тех же показате­ лей. Из того, что функции Лагранжа для задач (6.5) и (6.7) имеют одинаковую форму, следует, что двухступенчатые процедуры, при

которых обе функции Лагранжа достигают

экстремума, могут

быть описаны математически

одинаковыми

выражениями.

Таким образом, мы видим,

что в случае,

когда одним из воз­

можных решений является решение о сплошной проверке, зада­ ча (6.4) с некоторым приближением может быть заменена задачей (6.7).^Напомним, что такая замена возможна лишь при условии, что пг + пш N.

Вдальнейшем мы будем рассматривать процедуры с двумя возможными решениями: о приемке и о браковке контролируе­ мого объекта.

Втом случае, когда допущение пг + пт <^ N не выполняется, упомянутые выше задачи не соответствуют реальным условиям. Это связано с двумя обстоятельствами.

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

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

и f> не достаточ­

но полно характеризуют надежность контроля.

В случае, когда

проверке подвергается вся партия, эти показатели вообще те­ ряют смысл. В разд. 6. 6 мы рассмотрим задачу контроля доли дефектных изделий в партии; в этой задаче вместо а мы введем число бракуемых годных изделий при условии, что партия хоро­ шая, а вместо Р — среднее число принимаемых дефектных изделий при условии, что партия плохого качества.

Во-вторых, двухступенчатая процедура, имеющая в ансамбле окончательных решений решения о приемке и браковке партии, в рассматриваемом случае допускает фактически три исхода: прием­ ку партии, браковку партии и сплошную проверку. Последний ис­ ход может иметь место в случае, если пш = N — п-у. Подробно этот вопрос мы рассмотрим в разд. 6.6.

6.2.Выборочный контроль доли дефектных изделий

впартии. Постановка вопроса

Предположим, что производится

контроль

партии, состоящей

из ТУ изделий, по альтернативному

признаку.

Качество всей пар­

тии определяется долей 0 содержащихся в ней дефектных изде­ лий. Пусть D — число дефектных изделий в партии. Тогда их доля равна 0 = DIN.

Рассмотрим статистики, на которых может быть построена опти­ мальная двухступенчатая процедура. Очевидно, что статистикой первой выборки является число dx дефектных изделий, обнаружен­ ных в ней при контроле. Статистикой второй выборки является число d2 дефектных изделий, содержащихся в ней.

Найдем совместное распределение величин dx и d% при фикси­ рованных значениях объемов выборок % и п2 и при фиксирован­ ном значении D.

Как известно [13], если из партии объема N, содержащей D дефектных изделий, извлечь выборку объема ге, то с вероятностью

i?,tv мы обнаружим в выборке d дефектных изделий; величина

г>сШ

B n N равна

D\IN-D

 

 

 

 

d

)[

n—d

если

0 < | d < I мин [re, D ] ,

 

 

 

 

n

 

 

 

(o.«)

 

 

 

 

 

 

0

,

если

d ] > мин [re, £>].

 

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

(6.8)

называется гипергеометрическим. Величина

^!^ j обозначает количество сочетаний из г/ по а:. В выражении (6.8)

отражен тот очевидный факт, что величина d не может превышать ге и D . Такая запись нужна здесь потому, что ниже мы будем ис-

r>dD

пользовать величину nnN, ие оговаривая никаких ограничении сверху на d.

Найдем распределение величины d x . Очевидно, вероятность того, что в первой выборке объема пх окажется d x дефектных изде­ лий, равна

D \ I N — D

"1

Распределение числа дефектных изделий во второй выборке за­ висит от того, сколько дефектных изделий попало в первую выбор­ ку. Действительно, в тот момент, когда мы приступаем к извлече­ нию второй выборки, мы имеем партию, состоящую из N — гех изделий и содержащую D d^ дефектных изделий. Условная ве­ роятность того, что во вторую выборку попадет dn дефектных изде­ лий при условии, что в первую попало d x дефектных изделий, равна

/ D — d! \IN

— щ — D +

di \

 

R * . ( D - d . ) _ \ dj j \

m — d%

J

, f i . п .

Л2

Совместная вероятность того, что в первой и второй выборках окажется соответственно dx и d 2 дефектных изделий, равна произ­ ведению

P ( d 1 , d 2 ) = BdnfNB^fZdn'l

(6.11)

Алгоритм вычисления параметров двухступенчатой процедуры при гипергеометрическом распределении чисел d x и d 2 весьма сло­ жен (см. разд. 6.6).

Возможность упрощения алгоритма возникает в тех случаях, когда мы имеем дело с большими партиями изделий, в которых ве­ личина N исчисляется сотнями или тысячами, в то время как сум­ марный объем выборок 1 + п2 составляет небольшую часть от N. Тогда с достаточно хорошим приближением заменяют гипергео­ метрическое распределение более простыми, биномиальным пли пуассоновским.

Идея биномиального приближения основана на допущении о том, что при извлечении очередного изделия вероятность того, что оно окажется дефектным, постоянна для всей выборки и равна доле дефектных изделий в партии 0. Таким образом, мы имеем схему Бернулли [4]. Вероятность того, что в выборку из п изде­ лий попадает d дефектных изделий, выражается в виде

(6.12)

В тех случаях, когда 9 мало, применяют пуассоновское рас­ пределение; при этом выражение вероятности попадания d дефект­ ных изделий в выборку из їх изделий имеет вид

(6.13)

Математическое ожидание величины d для обоих распределе­ ний равно riQ.

При биномиальном и пуассоновском приближении величины di и d2 можно рассматривать как статистически независимые ве­ личины, если параметр 0 фиксирован. Совместное распределение величин dx и d2 можно записать в виде

P(d,, d2 ) = Р £ ( 0 ) Р £ ( 0 ) .

(6.14)

При контроле доли дефектных изделий обычно

рассматривают

две конкурирующие гипотезы, Н0 и Нх. Гипотеза

Н0 состоит в

том, что доля дефектных изделий меньше некоторой критической

величины 0 к р , а гипотеза

Их — в том, что доля дефектных изде­

лий не менее 8К р. Партия

полагается хорошей и подлежит прием­

ке, если 0 <С 0 к р , и плохой, подлежащей браковке, если 0 > 0 к р . Эти гипотезы являются сложными. Для того чтобы иметь дело с простыми гипотезами, пользуются различными искусственными приемами. Наиболее распространенный из них состоит в том, что

выбирают две характерные точки, 0О и 0^ при этом 0О <

Окр, 9i

>

> 0 к р . Тогда

вместо сложных

мы

имеем

две

простые

гипотезы:

гипотеза

На

о том, что 8 = 0О , и гипотеза Нх

о том, что 0 =

ву.

Совместное распределение

величин d x

и d2

при гипотезе

Hk

(ft = 0,1)

можно записать в виде

 

 

 

 

 

 

й2 ) = і>*(Є й )Р*(Є й ),

ft

= 0,1.

 

 

(6.15)

Рассмотрим статистику двух выборок, на основе которой дол­ жно быть построено оптимальное окончательное правило решения. Пусть в результате контроля в первой и второй выборках обна­ ружено соответственно dx и d 2 дефектных изделий. Естественно предположить, что окончательное решение о качестве партии должно приниматься в зависимости от суммарного числа обна­ руженных дефектных изделий dy -f- d 2 . Действительно, как было показано в главе 4, в оптимальном окончательном решающем пра­ виле в качестве статистики двух выборок должна быть использо­ вана величина отношения правдоподобия

1 ^ ^ = Ш % -

(6-16)

В случае биномиального распределения, используя формулы (6.12) и (6.15), имеем

Z(d1 ,d2 ) = (l r )

( т ^ )

.

(6.17)

В случае пуассоновского приближения, используя (6.13) и (6.15), получаем

I (dy, d2 ) = ( ~ ) d i + d ! exp [ - К + n8 ) (Є, - 80 )].

(6.18)

Таким образом, мы сформулировали две простые гипотезы о качестве партии изделий и определили для них совместные рас­ пределения величин dy и d 2 . Этих данных вполне достаточно для построения оптимальных двухступенчатых процедур контроля в случае биномиального или пуассоновского приближения.

6.3. Двухступенчатый контроль по альтернативному признаку

при постоянном объеме второй выборки

Допустим, что контроль доли дефектных изделий в партии про­ изводится с помощью процедуры, имеющей постоянный объем второй выборки. Напомним эту процедуру. Извлекают первую

выборку объема щ.

Если dy ^

Су, то партия признается

годной.

Если dy

с2 , то партия

бракуется.

Если Су <

dy ^

с3 , то извле­

кается вторая выборка

объема п2. Если dy + d 2

^ с3,

то

партия

принимается. Если dt +

d 2 ^> с3 , партия бракуется.

 

 

 

Таким образом, каждая процедура характеризуется пятью

параметрами: Пу, щ,

сх,

с2, с3. Сокращенно эти процедуры будем

записывать в виде s = (пх, п2, Су, с2,

с3). Например, запись (20, 25,

3,

5, 8) означает, что мы имеем процедуру s с параметрами Пу =

=

20, п2 =

25, сх =

3,

с„ = 5,

с 3 =

8.

 

 

 

 

Параметры Су, с2 и с3 обычно называют приемочными числами.

Приемочные числа должны удовлетворять соотношению

 

 

су < с2 <

сэ.

 

 

 

 

 

 

(6.19)

Действительно, в случае равенства сх = с, двухступенчатая про­ цедура вырождается в одноступенчатую. Это означает, что среди

значений dx

не найдется ни одного такого, при котором

должна

проводиться

вторая ступень контроля. Далее, если dx >

с3 , то

сумма dx -f- d2 обязательно будет превышать с3 . Поэтому при дан­ ном значении d± партия должна браковаться без второй ступени. Отсюда получаем условие сг < ; с 3 .

Распределение чисел dx и d2 будем считать биномиальпым или

пуассоновским. Показатели а,

В =

1 л, Е0, Ех

определяются

формулами (см. главу 3)

 

 

 

 

 

тії

 

 

Са

Пг

 

 

<*=

 

2

^ ' . ( в о ) +

2

2

рХл^Ркю,

(6.20)

 

d i = C a - i - l

 

 

d\=ci-\-x d a = c 8 — t i i + 1

 

 

 

ti[

 

 

Сз

Па

 

 

я =

 

2

РпЛ^)+

2

2

р'пЛ^г)pt(Qd,

(б.2і)

 

d i = c a + l

 

 

d . = r H - L d s = n a — d , + i

 

 

 

 

 

сз

 

 

 

 

#o

=

ih +

14

2

К (0O ),

 

 

(6-22)

 

 

 

d i = c , + l

 

 

 

 

Ді

=

» i +

na

2

pn\ ( 0 i) -

 

 

(6-23)

Предположим,

что требуется

построить таблицы

оптимальных

двухступенчатых

процедур (nx,

п2,

сг, с2 , с3 ) для различных зна­

чений показателей а и 6. При этом параметры ее и В должны лежать

в заданных пределах осн ^

a ^ а в , 6 Н <^ В

Вв . Под оптималь­

ной процедурой здесь

понимается процедура,

обеспечивающая

минимальное значение

Е0

при некоторых значениях а и р.

Метод построения таких таблиц состоит в следующем. Вводят два множителя Лагранжа, Ха и Х0. Величины \ а и %0 можно выби­ рать в достаточной степени произвольно. Порядки этих величин будут ясны из приведенных ниже примеров. После того как мно­ жители Лагранжа выбраны, ищут максимум функции Лагранжа

L = я — %ла — К0Е0

(6.24)

по всем пяти параметрам: щ, щ, ct, с2, с3. По параметрам пх и їй про­ изводят прямой перебор. Максимально допустимые значения nL и ?г2 предполагаются заданными. Иногда может быть задано лишь максимально допустимое значение суммы пх + п2. В процессе пе­ ребора при каждой паре значений пх и /г2 определяют параметры сх , с2 и с3 . Параметр с 3 определяют непосредственно из оптимального решающего правила (см. главу 4). В нашем случае это правило имеет вид:

если Z(rfb d2) ^ Ха, то партия принимается; если l(dx, d2) ^> Яа , то партия бракуется.

Для случая пуассоновского приближения (см. (6.13)) имеем

I {du d2 ) = ( i l ) d , + d ! e X p [ - (щ + ra2) (Єх - Є0 )].

(6.18)

Приравняем правую часть величине Ка и решим полученное урав­ нение относительно суммы dx + d 2 . Решение обозначим с3 . Имеем

' с', = l n ( ^ y 6 a ) [In Яа + (га, + на) (8Х -

60 )].

 

(6.25)

Очевидно, что параметр с 3 равен целой части числа с3 :

 

 

С з = Ы .

 

 

'

(6.26)

В случае биномиального приближения

 

 

Ийі, da) = (-Q^)

(ТТГТ^)

( 6 - 1 7 )

Отсюда аналогичным

образом находим

 

 

 

с3 = -g" I K + , г * ) + 1пЯа ],

 

(6.27)

где

 

 

 

 

Параметр с 3 по-прежнему равен с 3 = [с3 ].

 

 

Нетрудно видеть, что если 0О < ^ 1 и 0J

1, то

 

Т = In (1 - 8„) - In (1 - 0Х) з * 6t — 60 ,

 

(6.30)

и формула (6.27) переходит в (6.25).

 

 

Параметры сх и сг могут быть найдены следующим путем. Про­

изводят перебор значений величины а\ начиная с dx

0 в сторону

увеличения и для каждого dx вычисляют отношение правдоподобия

на первой ступени.

В случае

биномиального приближения

 

" • - ( £ ) * ( 4 = ! р .

 

 

<6 -3 1 >

а для пуассоновского

приближения

 

 

Zi = ( - ^ - ) d ' e x p [ r a 1 ( 0 o - 0 1 ) ] .

 

 

(6.32)

Далее вычисляют функцию A i^-f-,

(см. главы 4 и 5). В нашей

задаче эта функция

равна

 

 

 

А ( х " 'щ

)= а

( х 'n*)-bh=

 

2

^

( Є

і

) - х

2

^ n ; ( 0 o ) - ^ « 2 .

(6.33)

( і з = С з — ( I i + l

d j = c 3

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