Опубликован: 11.08.2009 | Уровень: для всех | Доступ: платный
Лекция 14:

Принятие решений на основе моделей обеспечения качества

< Лекция 13 || Лекция 14: 1234567

Риск поставщика и риск потребителя, приемочный и браковочный уровни дефектности. С оперативной характеристикой связаны важные понятия приемочного и браковочного уровней дефектности, а также понятия "риск поставщика" и "риск потребителя". Чтобы ввести эти понятия, на оперативной характеристике выделяют две характерные точки, делящие входные уровни дефектности на три зоны (области) - А, Б и В. В зоне А почти всегда все хорошо, а именно - почти всегда экологическая обстановка признается благополучной, почти все партии принимаются. В зоне В, наоборот, почти всегда все плохо, а именно - почти всегда экологический контроль констатирует экологические нарушения, почти все партии бракуются. Зона Б - буферная, переходная, промежуточная, в ней как вероятность приемки, так и вероятность браковки заметно отличаются от 0 и 1. Для задания границ между зонами выбирают два малых числа-риск поставщика (производителя, предприятия) \alpha и риск потребителя (заказчика, системы экологического контроля) \beta, при этом границы между зонами задают два уровня дефектности - приемочный p_{пp} и браковочный p_{бр}, определяемые из уравнений

f(p_{пp}) = 1-\alpha,  f(p_{бр}) = \beta ( 6)

Таким образом, если входной уровень дефектности не превосходит p_{пp}, то вероятность забракования партии мала, т.е. не превосходит \alpha. Приемочный уровень дефектности выделяет зону А значений входного уровня дефектности, в которой нарушения экологической безопасности почти всегда не отмечаются, партии почти всегда принимаются, т.е. соблюдаются интересы проверяемого предприятия (в экологии), поставщика (при контроле качества). Это - зона комфортности для поставщика. Если он обеспечивает работу (входной уровень дефектности) в этой зоне, то его практически никогда никто не потревожит.

Если же входной уровень дефектности больше браковочного уровня дефектности p_{бр}, то нарушения почти наверняка фиксируются, партия почти всегда бракуется, т.е. экологи узнают о нарушениях, потребитель оказывается защищен от попадания к нему партий со столь высоким уровнем брака. Поэтому можно сказать, что в зоне В соблюдаются интересы потребителей - брак к ним не попадает.

При выборе плана контроля часто начинают с выбора приемочного и браковочного уровней дефектности. При этом выбор конкретного значения приемочного уровня дефектности отражает интересы поставщика, а выбор конкретного значения браковочного уровня дефектности - интересы потребителя. Можно доказать (см. ниже), что для любых положительных чисел \alpha и \beta и любых входных уровней дефектности p_{пp} и p_{бр}, причем p_{пp} меньше p_{бр}, найдется план контроля (n,c) такой, что его оперативная характеристика f(p) удовлетворяет неравенствам

f(p_{пp})\ge 1 -\alpha,  f(p_{бр})  \le\beta

При практических расчетах обычно принимают \alpha = 0,05 (т.е. 5%) и \beta = 0,1 (т.е. 10%).

В качестве примера вычислим приемочный и браковочный уровни дефектности для плана (n,0) . Из формул (4) и (6) вытекает, что

(1 - p_{пp})^n = 1 - \alpha, p_{пp} = 1 - (1 -\alpha)^{1/n}.

Поскольку риск поставщика \alpha мал, то из известного соотношения математического анализа

(1-\alpha)^{1/n}=1-\frac{\alpha}{n}+O\left( \frac{\alpha^2}{n^2}\right)

вытекает приближенная формула

p_{пp}\approx \frac{\alpha}{n}

Для браковочного уровня дефектности имеем

p_{бр} = 1 - \beta^{1/n}.

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

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

При среднем входном уровне дефектности р и применении контроля с разбраковкой с вероятностью f(p) партия принимается (и уровень дефектности в ней по-прежнему равен р ) и с вероятностью (1 - f(p)) бракуется и подвергается сплошному контролю, в результате чего к потребителю поступают только годные изделия. Следовательно, по формуле полной вероятности средний выходной уровень дефектности равен

f_1(p)= pf(p) +0(1 - f(p)) = pf(p).

Средний выходной уровень дефектности f_1(p) равен 0 при р=0 и р=1, положителен на интервале (0;1), а потому достигает на нем максимума, который в теории статистического контроля называется пределом среднего выходного уровня дефектности (сокращенно ПСВУД):

ПСВУД = max_{0 \le p \le 1}f_1(p)

Пример. Рассмотрим план (n,0) . Для него f(p) = (1 - p)^n и f_1(p) = p(1-p)^n. Чтобы найти ПСВУД, надо приравнять 0 производную среднего выходного уровня дефектности по среднему входному уровню дефектности:

\frac{df_1(p)}{dp}=(p(1-p)^n)\prime=(1-p)^n+pn(1-p)^{n-1}(-1)=(1-p)^{n-1}(1-p-pn)=(1-p)^{n-1}(1-(n+1)p)=0

В полученном уравнении корень р = 1 соответствует минимуму, а не максимуму. Поскольку непрерывная функция на замкнутом отрезке достигает максимума, то максимум достигается при

p_n=\frac{1}{n+1}

Следовательно,

ПСВУД = p_n(1-p_n)^n=\frac{1}{n+1}\left(1-\frac{1}{n+1}\right)^n ( 7)

По выражению (7) могут быть проведены конкретные расчеты. Однако оно довольно громоздко. Его можно упростить, используя один замечательный предел из курса математического анализа, а именно:

\varlimsup_{n \to \infty} \left(1-\frac{1}{n+1}\right)^{n+1}=e^{-1}=\frac{1}{2,718281828}\approx0,368 ( 8)

Сравнивая соотношения (7) и (8), видим, что

ПСВУД = \left(\frac{\frac{1}{n+1}}{1-\frac{1}{n+1}}\right)\left(1-\frac{1}{n+1}^{n+1}\right)

Первая скобка равна 1/n, а вторая согласно соотношению (8) приближается к 0,368 при росте объема выборки. Поэтому получаем простую асимптотическую формулу

ПСВУД \approx\frac{0,368}{n}

Для более сложных планов ПСВУД рассчитывают с помощью более или менее сложных компьютерных программ.

При проведенном выше рассмотрении основ статистического контроля расчетные формулы удалось получить лишь для простейших планов, в основном для планов вида (n,0) . Если ослабить требования и рассчитывать не на точные формулы, а на асимптотические, при n\to \infty, то можно справиться и с одноступенчатыми планами вида (n, c) .

Асимптотическая теория одноступенчатых планов. Пусть Х - число дефектных единиц продукции в выборке объема n. Как уже отмечалось, распределение Х является биномиальным и имеет вид

Р(Х = k) = C_n^k p^k (1-p)^{n - k},

где C_n^k - число сочетаний из n элементов по k, а p - входной уровень дефектности.

Пусть используется одноступенчатый план контроля (n, c) . Тогда оперативная характеристика этого плана имеет вид

f_p)=\sum_{1\le k \le c}P(X=k)=\sum_{1 \le k \le c}C_n^k p^k(1-p)^{n-k}

Пусть n\to \infty Тогда по Закону Больших Чисел теории вероятностей (по теореме Бернулли)

\frac{X}{n}\to p

(сходимость по вероятности). Значит, если с/n окажется заметно больше входного уровня дефектности р, то партии будут почти всегда приниматься, а если с/n окажется заметно меньше входного уровня дефектности р, то партии будут почти всегда отклоняться. Ситуация будет нетривиальной только там, где величины с/n и р близки друг к другу.

Хотя оперативная характеристика рассчитывается с помощью сумм биномиальных вероятностей, анализировать эти суммы затруднительно. Поэтому целесообразно найти для нее приближение с помощью теоремы Муавра-Лапласа. Имеем цепочку тождественных преобразований:

f(p)=P(X \lec)=P \left(\frac{X-np}{\sqrt{np(1-p)}}\le\frac{c-np}{\sqrt{np(1-p)}}\right)

Справа стоит именно то выражение, которое участвует в теореме Муавра-Лапласа. Воспользовавшись равномерной сходимостью в этой теореме, можно записать, что

\varlimsup_{n \to \infty} f(p)=Ф\left(\frac{c-np}{\sqrt{np(1-p)}}\right)

где ?(х) - функция стандартного нормального распределения с математическим ожиданием 0 и дисперсией 1.

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

Ф\left(\frac{c-np_{np}}{\sqrt{np_{np}(1-p_{np})}}\right)=1-\alpha,\\
Ф\left(\frac{c-np_{bp}}{\sqrt{np-{bp}(1-p-{bp})}}\right)=\beta ( 9)

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

p_{np}=\frac c n-\frac{1}{\sqrt n}\sqrt{\frac c n(1-\frac c n)}Ф^{-1}(1-\alpha),\\
P_{bp}=\frac c n-\frac{1}{\sqrt n}\sqrt{\frac c n(1-\frac c n)}Ф^{-1}(\beta) ( 10)

Так как величины с/n и р близки друг к другу, то при переходе от формулы (9) к формуле (10) в подкоренных выражениях приемочный и браковочный уровни дефектности заменены на c/n (с точностью до бесконечно малых более высокого порядка).

Поскольку при практическом применении статистического приемочного контроля, как уже отмечалось, принимают \alpha= 0,05, \beta=0,10, то в предыдущие формулы следует подставить Ф^{-1}(0.95)=1,64 и Ф^{-1}(0,10)=-1,28 Итак, итоговые формулы для приемочного и браковочного уровней дефектности имеют вид

p_{np}=\frac c n - \frac{1.64}{\sqrt n} \sqrt{\frac c n \left(1-\frac c n \right)},\\
P_{bn}=\frac c n - \frac{1.28}{\sqrt n} \sqrt{\frac c n \left(1-\frac c n \right)}

Перейдем к задаче синтеза. Пусть заданы приемочный и браковочный уровни дефектности. Требуется построить одноступенчатый план, имеющий эти характеристики. Из формул (9) следует, в частности, что

c-np_{np}=Ф^{-1}(c-\alpha)\sqrt{np_{np}(1-p_{np})},\\
c-np_{bp}=Ф^{-1}(c-\beta)\sqrt{np_{bp}(1-p_{bp})} ( 11)

Вычитая из первого уравнения второе, получаем, что

np_{bp}-np_{np}=\sqrt n\{Ф^{-1}(1-\alpha)\sqrt{p_{np}(1-p_{np}})-Ф^{-1}(\beta}\sqrt{p_{bp}(1-p_{bp})}\}

Следовательно, оценка n* необходимого объема выборки имеет вид

n*= \left(\frac{ Ф^{-1}(1-\alpha)\sqrt{p_{np}(1-p_{np}}-Ф^{-1}(\beta}\sqrt{p_{bp}(1-p_{bp})}}{p_{bp}-p_{np}} \right)^2

Для стандартных значений рисков \alpha = 0,05, \beta = 0,10 имеем:

n*= \left(\frac{1,64\sqrt{p_{np}(1-p_{np})}+1,28\sqrt{p_{bp}(1-p_{bp})}}{p_{bp}-p_{np}} \right)^2 ( 12)

С помощью уравнений (11) нетрудно найти оценку с* приемочного числа, заменив неизвестный объем выборки на его оценку n*. Будем использовать оценку

c*=n*p_{bp}+Ф^{-1}(\beta)\sqrt{n*p_{bp}(1-p_{bp})}

Для стандартного значения \beta = 0,10 имеем

c*=n*p_{bp}+1,28\sqrt{n*p_{bp}(1-p_{bp})} ( 13)

Итак, по формуле (12) можно рассчитать оценку объема выборки, затем по формуле (13) найти оценку приемочного числа. Необходимо отметить, что результаты расчетов по рассматриваемым асимптотическим формулам отнюдь не всегда дают целые числа, поэтому необходима корректировка полученных результатов.

Полученные формулы позволяют решить сформулированную выше задачу - по заданным приемочному и браковочному уровням дефектности подобрать такой одноступенчатый план контроля, что его оперативная характеристика f(p) удовлетворяет неравенствам

f(p_{пp}) \ge 1 - \alpha,  f(p_{бр})  \le \beta.

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

Пример. Пусть p_{пp} = 0,02, p_{бр} = 0,09. Тогда по формуле (12) оценка объема выборки равна

n*= \left(\frac{1,64\sqrt{0,02(1-0,02)}+1,28\sqrt{0,09(1-0,09)}}{0,09-0,02}\right)^2= \left(\frac{1,64*0,14+1,28*0,286}{0,07} \right)^2=8,51^2=72,42

Полученное число не является натуральным, поэтому вполне естественно откорректировать объем выборки до ближайшего целого, т.е. до n* = 72.

Оценку приемочного числа находим по формуле (13):

c*=72*0,09-1,28\sqrt{72*0,09*0,91}=6,48-1,28*2,428=3,37

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

Если объем выборки округлить до 73, то аналогично получим

c**=73*0,09-1,28\sqrt{73*0,09*0,91}=6,57-1,28*2,445=3,44

При округлении снова получаем 3.

С помощью первого из уравнений (11) можно построить оценку с* на основе приемочного уровня дефектности:

c*=n*p_{np}+Ф^{-1}(1-\alpha)\sqrt{n*p_{np}(10p_{np})}=n*p_{np}+1,64\sqrt{n*p_{np}(1-p_{np})}

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

c*= n*p_{np}+1,64\sqrt{n*p_{np}(1-p_{np})}=72*0,02+1,64\sqrt{72*0,02*0,98}=3,39

Итак, в результате асимптотических расчетов найден одноступенчатый план (72, 3).

< Лекция 13 || Лекция 14: 1234567
Михаил Агапитов
Михаил Агапитов

Не могу найти  требования по оформлению выпускной контрольной работы по курсу профессиональной переподготовки "Менеджмент предприятия"

Подобед Александр
Подобед Александр

Я нажал кнопку "начать курс" и почти его уже закончил, но для получения диплома на бумаге, нужно его же оплатить? Как оплатить? 

Александр Кулаков
Александр Кулаков
Россия, Челябинск
Максим Овцин
Максим Овцин
Россия, Санкт-Петербург, НИУ ИТМО