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

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

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

Величина А щ) возрастает по мере увеличения dx. Параметр сх определяют как наибольшее значение dx, при кото-

dx = 0 должна быть отнесена к области продолжения испытаний. Это означает, что решение о приемке партии никогда не прини­ мается после первой ступени. На практике такие процедуры не встречаются; обычно сх 2> 0.

После того как значение параметра сх найдено, продолжают перебор dx в сторону увеличения и для каждого dy вычисляют, как

и ранее, функцию А

Щ . Но при этом функцию сравнивают

уже не с нулем, а с величиной, равной 1

^ - . Вначале имеет

место соотношение

 

 

 

 

(6.34)

затем это неравенство меняет знак.

Параметр с2 определяют как наибольшее значение dx, при ко­ тором условие (6.34) выполняется.

После того как найдены параметры сх, с2 и с3, определяют ве­ личины а, я и Ей и, подставив результаты вычисления в (6.24), находят функцию Лагранжа L .

Перебор по riy и « 2 продолжают до тех пор, пока не будет най­ ден максимум функции Лагранжа. Таким образом определяют оп­ тимальную процедуру типа (пх, п2, сх, с2 , с3 ), соответствующую данным Ха и Х0.

Напомним, что вначале мы задались значениями Ха и Х0 про­ извольно. От выбора этих величин зависят значения параметров а и В оптимальной процедуры. Если выбор Ха и Х0 оказался удачным, мы получим значения а и 6 внутри заданных интервалов а н

— а в и 6 Н — р„. Если выбор оказался неудачным, то значения пока­ зателей а и В выходят за эти пределы. Если эти пределы достаточ­ но широки, то практически нетрудно подобрать нужные значения Ха и Х0. Выбрав их, мы получим одну из оптимальных процедур для нашей таблицы. Варьируя далее параметры Ха и Ха, можно

получить много оптимальных процедур, у

которых параметры

а

и В «заполняют» заданные интервалы а н

— Ов и 6Н — 6В .

 

В табл. 1 приведены двухступенчатые процедуры типа (пх,

п2,

с 1> c 2i с3 ), построенные этим методом. Предполагалось, что распре­ деление чисел dx и di является пуассоновским. Доля дефектных

изделий в хорошей партии предполагалась равной

60 = 0,1,

а

в плохой партии — Вх

= 0,2.

В табл. 1 против значений I n Ха

и

1пХа указаны восемь

чисел,

характеризующих

оптимальную

Таблица 1

 

 

 

 

 

 

 

 

 

 

 

— 6 , 5

 

— 6,0

 

— 5,5

 

75

100

33,3

50

70

26

35

35

13,9

— 0 , 5

8

15

24

5

10

16

3

6

9

 

0,047

 

0,059

0,088

 

0,1007

0,164

 

0,143

 

80

120

33,8

65

80

25,6

40

50

17,5

—0,25

ї : 9

17

28

7

13

20

4

8

12

 

0,0291

0,0622

0,056

 

0,0858

0,106

 

0,149

 

і

 

 

 

 

 

40

60

21,8

 

85

140

33

70

100

26,8

0

10

20

32

8

15

24

4

9

14

 

0,018

 

0,0656

0,0339

 

0,0897

0,0714

 

0,162

процедуру. Эти числа расположены согласно следующей схеме:

12 ч

С1 С2 сз

а1 3

Напомним, что Е\ = Еа — п^.

6.4. Двухступенчатый контроль

по альтернативному признаку при переменном объеме второй выборки

Напомним, что процедуры контроля с переменным объемом

второй выборки отличаются от процедур типа (пх,

п2,

сх, с2, с3)

тем, что когда dx изменяется в пределах сх <С dx

с2 ,

величины

п2 и с 3 имеют не постоянное значение, а являются функциями ве­

личины dx. Эти функции будем обозначать пг

(dj) и с 3 (dj). Функция

с 3 (dj) определяется на интервале сх <

dx ^

сг , а функция п2

(dx)

определяется при всех значениях dy =

0,1,..., пх. При этом п2

(dx)

положительна при сх < с\ ^ сг и равна нулю вне этого интервала.

Каждая

процедура

характеризуется одним

параметром Пу

и двумя функциями, пг

(dx)

и с3

(dy). Для краткого обозначения про­

цедуры введем запись s =

(пх,

п2

(dx), с3

(dy)). В этой записи мы не

указываем

параметров

сх

и

с2.

Эти

параметры

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

функцией да (dy) или, точнее, нижней и верхней границами области, где эта функция положительна.

А

Рис. 24

Рис. 25

Методика

вычисления таблиц процедур типа s = (пх, пг (dy),

с 3 (dy)) во многом аналогична методике вычисления таблиц проце­

дур типа (riy, п2, Су, с2 , с3). Отличие состоит

лишь в

следующем.

При вычислении таблиц процедур типа (пх, п2,

сх, сг, с3)

цикл пере­

бора по параметрам пх и п2 является внешним по отношению к цик­ лу перебора но dy. Здесь же цикл перебора по Пу является внешним

по отношению к циклу перебора

по ^ , а цикл перебора по dy —

внешним к циклу перебора по п2.

Другими словами, сначала мы

производим перебор по пх, при каждом пх

— перебор по dx и при

каждых

Пу ж dy — перебор по п2.

 

 

 

 

Перебор п2 (при фиксированных пх и dx)

производят до тех пор,

пока не будет найден максимум функции A (kjlx, п2)

(см. (6.33)).

и

Вычисление оптимальной функции п2

(dy) можно

производить

графоаналитическим методом,

пользуясь графиками функции

a

(%Jly,

п2) (см. главу 5). Напомним, что графики функции а (ка/

Цх, п2) строят с помощью семейства рабочих характеристик одно­ ступенчатой процедуры. В случае, когда наблюдаемая величина (в данном случае число дефектных изделий d) дискретна, рабочие характеристики представляют собой последовательности изоли­ рованных точек на плоскости а 0 я. Каждая точка рабочей харак­ теристики определяется двумя параметрами: щ (объемом выборки) и с (пороговым числом); абсцисса а точки (п2, с) равна вероятности того, что d > с при условии, что верна гипотеза Н0: 8 = 0О ;

п2),

ордината я точки (?і2, с) равна вероятности того, что d > с, если верна гипотеза Нх: 9 = Qx.

На рис. 24 показано семейство рабочих характеристик для слу­

чая, когда

d имеет биномиальное распределение, причем 0О =

= 0,1; Qx =

0,25. Возле каждой точки указано соответствующее

число с. Точки, соответствующие одинаковым щ, соединены пунк­ тиром.

Дискретный характер рабочих характеристик приводит к то­ му, что функция а (ка/1х, 2) имеет «волнистый» характер. На рис. 25 в качестве примера приведено несколько кривых a (kjlx, построенных по рабочим характеристикам рис. 24. На рис. 25

функции

a (kjlx, п2) для X J l x =

0,2; 1; 5 изображены волнистыми

рядами

точек — каждая точка

соответствует целочисленному

значению 7г2. При графических приближенных построениях удобно производить аппроксимацию рядов точек непрерывными кривы­ ми линиями. На рис. 25 такие аппроксимирующие кривые даны

для

приведенных

выше

значений X J l x , а

также для значений

%Jh

= 0,05; 0,5;

2; 8.

 

 

 

Волнистый характер

функций

а (kjlx,

п2 ) можно объяснить

следующим образом. Построение

графиков

функции a (XJlx , п2)

производится с помощью наклонных прямых, проходящих под углом arctg X J l x к оси абсцисс и совпадающих с рабочей характе­ ристикой либо в одной, либо в двух точках. При этом все осталь­ ные точки рабочей характеристики должны быть ниже наклонной прямой. Так, например, прямая АВ на рис. 24 и рабочая характе­ ристика п2 = 20 имеют одну общую точку с = 5. Значение функ­ ции a (kjlx, п2) мы определяем как ординату точки пересечения наклонной прямой с осью ординат. Если наклонная прямая не параллельна ни одному пунктирному отрезку рабочей характери­ стики, то их совпадение происходит в одной точке. Если же хотя бы один отрезок рабочей характеристики параллелен наклонной прямой, то совпадение происходит в двух точках (на концах этого отрезка).

Пусть заданы произвольные значения X J l x и ?г2. Допустим, что при этом совпадение происходит в одной точке («г, £)• Найдем с. Пусть / (d2, п2) — функция, связывающая отношение правдоподо­ бия на второй ступени 12 и d2. l2 — f (d2, п2 ), и пусть f"1 (Za, "г) —

— d2 — обратная функция. Если принять во внимание, что угол, между осью. 0 а и отрезком, соединяющим точки (п2, с) и (п2, с -f- 1),

равен

arctg

/

(с, п2) (см. рис. 24), то

становится очевидным, что

величина

с

равна

ближайшему

большему целому (ББЦ)

от

/ - 1

(Klk>

"г):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.35)

В случае, если

функция / 1 а/1Х,

пг) точно равна некоторому

целому числу

с, то

совпадение

происходит в двух точках,

(п2,

с)

и (п2,

с +

1).

 

 

 

 

 

 

Пусть отношение KJlx фиксировано (фиксирован угол между наклонной прямой АВ и осью Оа (см. рис. 24)). Будем изменять п2, и для каждого п2 вычислим с. На рис. 26 изображен качествен­ но вид зависимости с от п2. Как видно из рис. 26, на оси щ суще­ ствуют некоторые отрезки, внутри которых значение с является неизменным. Найдем размеры этих отрезков для случая, когда ве­ личина dz имеет биномиальное и пуассоновское распределение. Функция / (d2, пг) в случае биномиального распределения равна

(см. (6.17))

U

 

Г M l - В о )

(6.36)

f(d

г.)

i * / 1 - 6 , у

а обратная функция (при фиксированном 1г = % J l x )

 

 

К

/ 1 - Є о

(6.37)

Єї (1 - Bp) In бо (1 - Qi)

Следует подчеркнуть, что переменная щ здесь рассматривает­ ся как величина непрерывная. Мы видим, что в случае, когда di изменяется на единицу, величина тг2 изменяется на

Ди =

1

+

In Єі/Є0

 

 

 

 

 

 

(6.38)

In

( 1 - Є о )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 1 - Є і )

 

 

 

 

 

 

 

Аналогично

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

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

полу­

0

 

 

 

 

 

чаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

=

("gj") 2 e x p [ - r a 2 ( 0 1

•бо)],

(6-39).

4 -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- З

 

 

 

 

 

An

In (ЄХ/ЄО)

 

 

 

(6.40)

г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Волнистый

характер

функций

S i Ю

15

20

\2S

л2

a (XJl t ,

щ) объясняется тем, что при уве­

 

личении га2 точки совпадения наклонной

A'n

иг

7n

A'n

 

прямой и рабочих

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

Рис.

26

 

 

 

 

ходят с одного значения

с

на другое.

 

 

 

 

Для кривых а /^і,

изображенных

 

 

 

 

 

 

на рис.

2 5 , 0 О =

0 , 1 ; 0J. =

0,25;

вели-

чина An («период волнистости») равна

 

 

 

 

Ап = 1 + • In0,25/0,1

=

6,02 s ;

6.

 

 

 

 

 

 

1п-

•0,1

 

 

 

 

 

 

 

 

 

•0,25

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функция

с3

(dx) при

вычислении процедур ( « l t

n 2

(dj), с 3 (d^

определяется по тем же формулам (6.25) — (6.29), что и при вы­ числении процедур (%, « 2 , сх , са , с3 ). Отличие состоит лишь в Т0М (

что вместо параметра ?г2 в эти формулы подставляется функция

Приведем в качестве примера результаты вычислений опти­

мального

двухступенчатого плана при

следующих условиях:

V = 1; Я0

= Ю - 2 ; 60 = 0,1; 0Х = 0,25.

Полагаем, что числа dy и

d2 имеют биномиальное распределение. При этих условиях опти­ мальный план имеет следующие характеристики:

щ =

20;

Су =

2;

 

с2

= 6;

 

 

 

правила n2(dy)

 

и

cs(dy):

 

 

 

если

dy — 3,

то

щ =

14,

с3

—• 5;

если

dy =

4,

то

2

=

31,

с3

=

8;

еслп

dy =

5,

то га2

=

46,

с3

=

10;

если

сД =

6,

то

« 2

=

60,

с3

=

12.

Показатели плана: а = 0,064, [3 = 0,147, 2?0 = 27,4.

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

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

Контроль доли дефектных изделий

Рассмотрим случай, когда контролируемая партия состоит из большого числа (сотен или тысяч) изделий и качество каждого из­ делия определяется значением некоторого параметра и. Величина и различна для разных изделий. Пусть Ф (и) — функция, которую будем называть распределением величины и в партии. Это означает, что доля изделий, у которых а < и Ь, равна Ф (Ь) — Ф (с). Здесь а и Ъ — произвольные числа, причем Ъ^> а. В литературе часто распределение Ф (и) полагают равным

du,

(6.41)

—оо

причем параметры и0 и сг неизвестны.

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

Г' f Изделие считается дефектным, если

и >

u , t p , где и к

р

неко­

торый заданный порог.

 

 

 

 

 

 

 

 

 

Без ущерба для общности положим,

что г/.Іір = 0. Если и к р =f=

=j= 0, то мы можем вместо и рассматривать разность и

и к р .

 

но,

Нас интересует доля дефектных изделий в партии 8. Очевид­

что

 

 

 

 

 

 

 

 

 

Є =

Вер > 0} = 1 -

Ф

) = Ф

 

 

.

 

(6.42)

 

Мы видим, что величина 0 зависит от

отношения параметров

иа

и о, но не от их значений. Обозначим это отношение у

=

и0/а.

 

Параметр у однозначно связан с интересующим нас парамет­

ром©:

0 = Ф (у), -или у

= Ф'1

(0). Здесь

Ф - 1

(а:) — функция,

об­

ратная

Ф (а:).

 

 

 

 

 

 

 

 

 

Сформулируем две простые

гипотезы

относительно

парамет­

ра у. Допустим, что качество партии считается удовлетворитель­ ным, если 0 <^ 81 ( р , где 0цр — заданная норма. Выберем два чис­

ла 0О

и 0! (Оо <

 

0«р, 0Х > е к р ) .

 

 

Пусть у0

=

Ф 1 (0О ), 7а =

Ф 1 (0а). Тогда наши гипотезы можно

записать

в

следующем виде: гипотеза Н0:

у = у0; и

гипотеза

Нх: у

=

ух.

 

 

 

 

 

Рассмотрим двухступенчатую процедуру для проверки этих

гипотез. Пусть

X = (хг,...,

хП1) — совокупность измеренных зна­

чений параметра

и в первой

выборке, a Y

= (г/х ,..., упг)

— такая

же совокупность

во второй

выборке.

 

 

Процедуру построим на следующих статистиках. Для первой

выборки вычисляется статистика

 

 

-

-

^

S1

 

 

 

 

(6.43)

 

 

 

 

 

 

 

 

где

Sпі і=і.

выборочное среднее,

выборочная дисперсия. Аналогично для второй выборки вычис­ ляется статистика

т2 = J ^ L ,

(6.46)

S2

 

где

П2 -12 ( г / 1 - у ) 2 -

( 6 - 4 8 )

Статистики тх и т2 имеют нецентральное і-распределение [5]. Плотность вероятности статистики Xj (/ = 1,2), когда верна ги­ потеза Як . (/с = 0,1), может быть записана в виде [5, стр. 302]

со

 

 

i

 

 

/ , (т,-) = С^х1 *ч-*>

ехр

{ -

4" [х + (*} Y-фт

~

V^r,f]}dx,

 

 

 

 

 

(6.49)

где С — постоянная; /

=

1,2.

При этом функция Д (т,-) не зависит

от параметров и0 и а, а зависит от их отношения, которое в рас­

сматриваемом случае может принимать два

 

значения, у0

и yv

 

Поскольку

данные измерений

в

первой

и

второй

выборках,

X

и Y, в каждой партии изделий

статистически независимы,

ве­

личины ТІ

E T S также взаимно независимы. Поэтому их совместная

плотность

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

 

 

 

 

 

 

 

 

 

 

 

 

Д ( т 1 , т 2 ) = Д.(т1 )Д.(г2 ).

 

 

 

 

 

 

 

 

 

(6.50)

 

Пусть

заданы число пг и функции п2

х)

и с4

х ).

При

 

этом

пг

(ti) = 0, если t j <

Cj и

>

с2 , где ^

и с2

— некоторые числа

2

> сх).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Процедура контроля построена следующим образом. Извлека­

ют первую выборку X и вычисляют статистику хх. Если хх

^

сх,

то партия принимается. Если хх

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

Если

c i

< т і ^

с 2> то извлекается вторая выборка

Y

объема щ

 

х).

Далее вычисляют статистику т2 . Если т2

^> с4

х), то партия бра­

куется. Если т2

^

с4

х), то партия принимается.

 

 

 

 

 

Отношения

правдоподобия 1Х х ) = ^ ^ |

и Z2

2 ) =

^

 

яв-

ляются монотонно возрастающими

функциями статистик хх

и т2

соответственно. Оптимальное

окончательное

решающее

правило

имеет вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

если

<С_ Яа ,

то

принимается

гипотеза

Н0;

 

(6.51)

 

если ^ 2 ^ > Я а ,

 

то

принимается

гипотеза

 

Нг.

 

 

 

 

 

 

 

 

 

 

Оптимальная функция с4 х ) может быть

найдена с помощью

уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

К.

 

 

 

 

 

 

 

 

 

 

 

(6.52)

Для того чтобы найти с4 х ), необходимо решить это уравнение относительно т2 . Величина с4 равна корню этого уравнения.

Заметим, что уравнение (6.52) имеет единственный корень. Это следует из того, что функции 1г х ) и 12 2 ) являются монотонно возрастающими. Функция с4 (Ху) — монотонно убывающая.

Для вычисления оптимальной функции п2 (ху) можно исполь­ зовать графоаналитический метод, изложенный в главе 5.

Показатели а, я и Е0 можно записать в виде:

оо Сг со

а

=

J /о (Ті) dxx

+ J

J

/о ( t j /о ( т а ) d t i d t 2 ,

 

 

с.

с, с4 (т,)

 

 

 

оо

Сз

оо

 

я = S Л ( T i ) d T i +

$

/ і ( t i ) / і ( т г ) d t j d T 2 ,

 

 

СI

C j

C i ( - C i )

 

 

 

 

 

 

# 0

=

« 1 + S " 2

(Ті) / О ( Т х )

d t i .

 

 

C l

 

 

 

(6.53)

(6.54)

(6.55)

Вычисление этих величин следует проводить путем численного ин­ тегрирования.

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

Мы рассмотрели выше случай контроля партии по количест­ венному признаку, когда контролируемым параметром считалось отношение у = щіа. Теперь рассмотрим случай, когда контроли­ руемым параметром партии является а. Параметр о* характеризует степень «неоднородности» параметра и в партии. Распределение

параметра и будем по-прежнему считать нормальным N {и0,

ст2).

Допустим, что качество партии считается удовлетворительным,

если а <^ а к р , и неудовлетворительным, если

а ^> сткр, где

сткр

заданная постоянная. Выберем два числа, ст0

и ау, такие, что сг0

<

< о"Кр, cTj. J> crKp. Гипотеза Н0 состоит в том, что о* = а 0 , а гипо­ теза Ну — в том, что 0 = а у.

Процедура контроля в данной задаче имеет такую же структу­ ру, как процедура контроля доли дефектных изделий по количест­ венному признаку. Отличие состоит лишь в том, что в рассматри­

ваемой процедуре

мы будем использовать статистики s\ и s\

вместо статистик тг

и

т2 .

 

Распределения статистик s\ и s\ выражаются через

х 2 _ Р а с п Р е ~

деление. Известно

[5], что для выборки обьема п из нормальной

совокупности N (и0,

о- 2 ) величина s2 (п — 1)/ст2 имеет

%2 -распре-

деление сп 1 степенями свободы (здесь s2 — выборочная диспер­

сия).

Таким образом, величина si

(пу 1)/CTQ при

гипотезе Я 0

(о =

ст0) и величина si (tiy — 1)/стї при гипотезе Ну(а

= ау) имеют

одинаковое распределение %2 с п

1 степенями свободы. Анало-

гичпо величины s2 2 — 1)/ст2 и si (щ — 1)/ст2 имеют ^-распределе­ ние с по, — 1 степенями свободы (при гипотезах Н 0 и Н І соответ­ ственно).

Для краткости запишем выражение для плотностей вероятно­ сти статистик s\ и s\ в виде одной формулы. Плотность вероятности

статистики

s) (/

=

1,2) при гипотезе

Hh (к — 0,1) равна

U (А) =

-

Ф

1

є' т

> 4 1 4 ,

/ = 1.2,

к = 0, 1 (6.56)

 

 

( 4 ) ' Г

(иг,-)

 

 

 

где ???7- = (7ij

— 1)/2, Г (jn-j)

гамма-функция.

 

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

мы. Из этого

следует, что

величины si и si при

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

значении ст (ст = ст0 или ст = CTl) также статистически

независимы

Приведем

выражение

для отношений правдоподобия

lx (s2)

и l2 {si). Величина Ц (s)), j

= 1, 2, равна [5]

 

 

 

 

 

/

=

1,2.

(6.57)

Отсюда видно, что функции lx (s\) и Z2 (s\) являются монотонно воз­ растающими. Двухступенчатая процедура в рассматриваемом слу­ чае строится таким же способом, как и процедура контроля доли дефектных изделий по количественному признаку (см. первую часть настоящего раздела).

6.6. Двухступенчатый контроль

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

Рассмотрим снова случай контроля доли дефектных изделий по альтернативному признаку; при этом будем считать, что объем выборки составляет значительную часть объема партии. В этом случае уже нельзя использовать биномиальное и пуассоновское приближения. Действительно, если контролю подвергается зна­ чительная часть партии, то в процессе извлечения изделий из партии' существенно изменяется вероятность того, что очередное изделие окажется дефектным. Эту вероятность нельзя считать по­ стоянной (какэто мы делали, применяя биномиальное и пуассо­ новское распределения).

Как было отмечено в разд. 6.1, числа дефектных изделий в вы­ борках имеют гипергеометрическое распределение (см. (6.8)). В данном разделе мы приведем соотношения, необходимые для построения оптимальной двухступенчатой процедуры контроля при гипергеометрическом распределении чисел dy и а^.

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