36
Практическая работа № 6
Контроль качества по методу однократной выборки
Контроль качества изделий предполагает проверку гипотезы о том, что качество изделий не ниже установленного уровня. При этом конечным результатом контроля является принятие одного из двух решений: принять партию изделий, считая качество изделий удовлетворительной, или забраковать контролируемую партию изделий как некачественную. При этом возможны два вида ошибок:
ошибка первого рода – когда хорошая партия бракуется, поставщик в этом случае рискует, и вероятность его риска обозначим буквой α;
ошибка второго рода – когда плохая партия принимается, рискует в этом случае заказчик, и вероятность его риска обозначим буквой β.
Одним из методов контроля качества является метод однократной выборки, основное достоинство которого заключается в том, что он легко планируется и осуществляется.
Метод однократной выборки заключается в том, что из контролируемой партии объема N изделий берется одна случайная выборка, объема n экземпляров. Определяются числа D0 и D1 – количество некачественных изделий во всей партии. При этом
D0 > D1 . |
(6.1) |
Если число дефектных изделий D<D0, в партии объемом N, то партия считается высоконадежной. Если число дефектных изделий D>D1, в партии объемом N, то партия считается дефектной. Если число дефектных изделий D0<D<D1, в партии объемом N, то партия считается неплохой и ее можно принять.
Исходя из следующих данных:
N – количество изделий в контролируемой партии; n – количество изделий в выборке;
d – количество бракованных изделий в выборке;
D0 – минимально допустимое число дефектных изделий в контролируемой партии;
D1 – максимально допустимое число дефектных изделий в контролируемой партии;
37
α – риск поставщика; β – риск заказчика,
определяются оценочные нормативные значения А0 и А1 для определения надежности контролируемой партии изделий.
Нормативные значения А0 и А1 могут быть определены из следующих соотношений
|
A |
d |
|
n−d |
|
|
|
||
α/ |
=1−∑0 |
CD0 CN −D0 |
, |
|
(6.2) |
||||
|
n |
||||||||
|
d =0 |
|
CN |
|
|
|
|||
|
A1 −1 |
d |
|
n−d |
|
|
|
||
β/ |
= ∑ |
CD1 |
CN −D1 |
|
, |
|
(6.3) |
||
|
|
n |
|||||||
|
d =0 |
|
CN |
|
|
|
|||
где α/ – риск поставщика близкий к заданному α; |
|
||||||||
β/ – риск поставщика, близкий к заданному β; |
|
||||||||
|
|
|
n |
|
N! |
|
|
|
|
|
|
CN = |
|
|
|
||||
|
|
n!(N −n)! |
|
||||||
В общем случае α/ ≠α и β/ ≠ β из-за дискретности значений получаемых по формулам (6.2) и (6.3), в которых определяется вероятность появления дискретной случайной величины, распределенной по гипергеометрическому закону. Поэтому должны выполняться следующие условия:
α′ ≤1,2 |
α |
(6.4) |
β′ ≤1,2 |
|
|
β |
|
Практическое использование формул (6.2) и (6.3) ограничено значениями выборки. При N >100 вычисление сочетаний в формулах (1.2) и (1.3) весьма затруднительно. Для приближенного вычисления n! в случае очень больших чисел n можно воспользоваться формулой Стирлинга
n n |
|
|
|
|
|
|
|
|
|
||
n!≈ |
|
2 π n . |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
||||
e |
|
|
|
|
|
|
|
|
|
||
При N ≤ 500 , q0 |
= |
D0 |
< 0.1 и q1 |
= |
D1 |
< 0.1 вместо формул (6.2) и (6.3) удобнее |
|||||
N |
|
||||||||||
|
|
|
|
|
|
|
|
N |
|
||
воспользоваться несколько упрощенными формулами |
|
||||||||||
|
|
|
|
|
A |
|
|
|
|
|
|
|
|
|
α/ =1−∑0 |
CDd |
0 |
f d (1− f )D0 −d , |
(6.5) |
||||
|
|
|
|
|
d =0 |
|
|
|
|
|
|
|
|
|
|
|
A1 −1 |
|
f d (1− f )D1 −d , |
|
|||
|
|
|
|
β/ = ∑CDd1 |
(6.6) |
||||||
|
|
|
|
|
d =0 |
|
|
|
|
|
|
38
где f = Nn .
Когда объем партии изделий N > 500 и n ≤ 0.1 N целесообразно использовать биномиальный закон распределения, в соответствии с которым
|
|
A |
Cnd q0d (1−q0 )n−d , |
|
|||||||
α/ =1−∑0 |
(6.7) |
||||||||||
|
|
d =0 |
|
|
|
|
|
|
|
|
|
|
|
A1 −1 |
|
|
|
|
|
|
|
|
|
β/ = ∑Cnd q1d (1−q1 )n−d , |
(6.8) |
||||||||||
|
|
d =0 |
|
|
|
|
|
|
|
|
|
Если выполняются условия: |
|
|
|
|
|
|
|||||
n ≤ 0.1 N; |
q0 |
≤ 0.1; |
q1 ≤ 0.1; |
(6.9) |
|||||||
то, пользуясь распределением Пуассона, получим |
|
||||||||||
|
|
∞ |
|
|
d |
|
|
|
|
||
α/ |
= ∑ |
a0 |
e−a0 |
, |
(6.10) |
||||||
|
|
||||||||||
|
|
d =A0 +1 d! |
|
|
|
|
|||||
|
/ |
|
|
∞ |
ad |
|
−a |
|
|
||
β |
|
=1 |
− |
∑ |
1 |
e |
1 |
, |
(6.11) |
||
|
d! |
||||||||||
|
|
|
|
d =A1 |
|
|
|
|
|||
где a0 =q0 n; a1 =q1 n .
Пример. Партия изготовленных изделий, качество которых необходимо проконтролировать, состоит из N=50 штук. Производитель и заказчик договорились, что если в изготовленной партии изделий содержится не более q0=0,1 дефектных изделий, то партия считается качественной. Если в изготовленной партии изделий содержится более q1=0,2 дефектных изделий, то партия считается бракованной. Если в изготовленной партии изделий содержится более q0=0,1 и менее q1=0,2 дефектных изделий, то партию можно считать удовлетворительного качества. Поставщик согласен на риск α=0,15, а заказчик согласен на риск β=0,15. Определить приемочное (А0) и браковочное (А1) числа дефектных изделий в выборке объемом n=20 изделий.
Решение. Партия изготовленных изделий не большая (N<100), а относительный объем выборки велик (n/N=0,4), то контроль необходимо проводить исходя из гипергеометрического распределения, т.е. расчеты проводить по формулам (6.2) и (6.3).
1. Определяются исходные данные, необходимые для решения задачи:
– объем изготовленной партии;
– объем выборки;
39
q0 = 0,1 – значение границы, определяющей изготовленную партию изделий, как качественную;
q1 = 0,2 – значение границы, определяющей изготовленную партию изделий, как дефектную;
D0 = n q0 =50 0,1 =5 – максимальное число дефектных изделий в качественной партии;
D1 = n q1 =50 0,2 =10 – минимальное число дефектных изделий в не качественной партии;
α= 0,15 – риск производителя;
β=0,15 – риск заказчика.
2. Для определения приемочного числа А0 дефектных изделий в выборке воспользуемся таблицей 6.1, из которой определим формулы, соответствующие диапазону значений исходных величин. Из таблицы 6.1 видно, что для представленных выше данных необходимо применить формулы (6.2) и (6.3). Для определения приемочного числа воспользуемся формулой (6.2) В этой формуле произведем суммирование вероятностей гипергеометрического распределения по тех пор, пока накопленная вероятность не приблизится к величине
P(d ≤ A0 ) =1−α =1−0,1 = 0,9 . (6.12)
Величины вероятностей для каждого d определится из следующих соотношений:
P(d = 0) = |
C0 |
C20−0 |
|
|
= |
1 3169870830126 |
= 0,067 |
(6.13) |
|||||||||||||
|
|
5 |
50−5 |
|
|
47129212243360 |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
C20 |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
50 |
|
|
|
|
|
|
|
|
|
|
|
|
|
P(d =1) = |
|
C1 |
C20−1 |
= |
5 |
2438362177020 |
= 0,258 |
(6.14) |
|||||||||||||
|
|
5 |
50−5 |
|
|
47129212243360 |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
C20 |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
50 |
|
|
|
|
|
|
|
|
|
|
|
|
|
P(d = 2) = |
|
C |
2 |
C20−2 |
|
|
= |
10 |
1715884494940 |
= 0,364 |
(6.15) |
||||||||||
|
|
5 |
50−5 |
|
|
|
47129212243360 |
||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
C20 |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
50 |
|
|
|
|
|
|
|
|
|
|
|
|
|
P(d = 3) = |
C3 |
C |
20−3 |
|
|
10 1103068603890 |
|
|
(6.16) |
||||||||||||
|
|
5 |
|
|
50−5 |
= |
|
|
47129212243360 |
= 0,234 |
|||||||||||
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
C20 |
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
50 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P(d = 4) = |
C4 |
C20−4 |
= |
5 646626422970 |
= 0,069 |
(6.17) |
|||||||||||||||
|
|
|
5 |
|
50−5 |
|
|
|
|
|
|||||||||||
|
|
|
|
C |
|
47129212243360 |
|||||||||||||||
|
|
|
|
|
|
|
20 |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
50 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
40
Таблица 6.1 Определение вида зависимостей для приемочного и браковочного чисел
Исходные данные |
|
|
|
|
|
|
Номера рекомендуемых формул |
|
||
|
|
|
|
|
|
|
|
|
|
|
N – объем изготовленной партии |
n – объем выборки |
|
– значение, определяющее |
изготовленную партию изделий, как качественную |
– значение, определяющее |
изготовленную партию изделий, как дефектную |
α – риск производителя |
β – риск заказчика |
|
|
|
0 |
1 |
|
|||||||
|
q |
q |
|
|||||||
100 |
< N |
|
|
0…1 |
|
0…1 |
0…1 |
0…1 |
(6.2); (6.3) |
|
|
|
|
|
|
|
|
|
|
|
|
500 |
< N |
|
|
< 0,1 |
|
< 0,1 |
0…1 |
0…1 |
(6.5); |
(6.6) |
|
|
|
|
|
|
|
|
|
|
|
500 |
< 0,1N |
|
|
0…1 |
|
0…1 |
0…1 |
0…1 |
(6.7); |
(6.8) |
|
|
|
|
|
|
|
|
|
|
|
500 |
> 50 |
|
|
< 0,1 |
|
< 0,1 |
< 0,1 |
< 0,1 |
(6.10); |
(6.11) |
|
|
|
|
|
|
|
|
|
|
|
Суммирование значений из (6.13)…(6.17) и сравнение со значением (6.12) проводим в следующем порядке:
P(d<=0) = 0,067 < 0,9
P(d<=1) = 0,067 + 0,258 = 0,325 < 0,9 P(d<=2) = 0,067 + 0,258 + 0,364 = 0,689 < 0,9
P(d<=3) = 0,067 + 0,258 + 0,364 + 0,234 = 0,923 > 0,9 P(d<=4) = 0,067 + 0,258 + 0,364 + 0,234 + 0,069 = 0,992 > 0,9
Принимая во внимание условие (6.4), определяем, что А0 = 3.
3. Для определения браковочного числа А1 дефектных изделий в выборке также воспользуемся таблицей 6.1, из которой видно, что для представленных выше данных необходимо применить формулы (6.2) и (6.3). Для определения браковочного числа воспользуемся формулой (6.3) В этой формуле произведем суммирование вероятностей гипергеометрического распределения до тех пор, пока накопленная вероятность не приблизится к величине
P(d A1) = β = 0,1 |
(6.18) |
Величины вероятностей для каждого d определится из следующих соотношений: