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

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

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

Если

ввести

дополнительным

условием db = 34,48 ч/ЮОО шт.,

то в примере 1 U х = 0,0290 X

34,48 =

1,0000 ч/ЮОО шт. Так как

К о =

1,0000,

экономия и потери в примере 1 уравновешиваются.

Итак,

db =

= 34,48 является

тем критическим значе­

нием db, которое отделяет операции, где выгодно полностью отка­ заться от приемочного контроля (при db < 34,48), от операций, где это невыгодно. В примере 2 таким критическим значением

является db — q = 135,14. Сопоставление S on для различ­

ных вариантов см. в табл. 13.

В связи с возможностью отказа от приемочного контроля эффективность увеличения точности настройки в примере 2 выгля­ дит в новом свете. До уровня db = 34,48 отказ от приемочного контроля выгоден в обоих примерах, но экономия выше при на­

стройках

с дополнительной проверкой

(пример 2). При 34,48 •<

< db <(

135,14 отказ от приемочного

контроля выгоден только

при настройке с дополнительной проверкой. Добавим, что в при­ мере 1 отказ от приемочного контроля является, в сущности, вве­ дением самоконтроля рабочего, а для этого необходимы определен­ ные условия. Иначе пришлось бы по меньшей мере прибегнуть к выборочному приемочному контролю. В примере 2 при некото­ рых мерах организационного обеспечения (см. гл. 11) вводится не самоконтроль, а совмещение функции ОТК с производствен­ ной функцией, благодаря чему отпадает необходимость в сплош­ ном и в выборочном приемочном контролях. Здесь на простом примере видно, что комплексная функция обеспечения качества действительно неделима при выборе эффективных вариантов СРК, что является основным исходным положением данной книги (см. предисловие и гл. 1).

Не следует думать, что оптимизация СРК выполняется путем частных сопоставлений, вроде сделанных выше. Методы выбора оптимального варианта описаны в гл. 8 и 9. В данной главе ста­ вится лишь ограниченная цель проиллюстрировать на близких к практике примерах смысл формул из гл. 3—6 и сделать попутные замечания имеющие общий интерес.

7.3. ПРИМЕР 3

Этот пример совпадает с примером 1 с той лишь разницей, что раз в час в сроки, не зависящие от начала технологического про­ межутка, в условиях примера 3 предусмотрены контрольные про­ верки отклонения у. н. v. План контрольной проверки предпо­ лагается таким же, как план дополнительной проверки в примере 2, т. е. с А, = 2,45, у+ = у~ = v5, что соответствует объему вы­ борки п = 6 и границам регулирования ±0,025 мм. Дополни­ тельная проверка настройки в примере 3 не предусмотрена, кон­ трольные данные при регулировках не используются.

141

Вычисление йероятности

 

 

 

 

 

 

 

oip-

+

 

 

 

 

 

 

 

 

 

CJ

 

 

 

II

 

 

 

о^-

сТ

о

 

 

 

 

 

 

О»

II

 

»а

 

 

й

 

 

 

CN

 

сч

 

 

 

 

 

 

1

 

 

 

 

а'

Ч*

° а

 

2

3

4

5

 

II-1

11-2

п - з

 

11-4

 

 

 

12

0,0000

0,0002

0.7G11

0,0002

0,0000

0,0000

0,0000

0,0000

— 11

0,0000

— 10

0,0008

0,6406

0,0005

0,0000

0,0002

0,0002

0,0001

0,0000

0,0022

0,5000

0,0011

0,0000

0,0004

0,0004

0,0002

— 9

— 8

0,0003

0,0048

0,3594

0,0017

0,0000

0,0009

0,0009

0,0003

0,0014

0,0104

0,2389

0,0025

0,0000

0,0021

0,0021

0,0005

— 7

0,0-101

0,0188

0,1423

0,0027

0,0008

0,0037

0,0045

0,0006

— 0

0.1908

0,0305

0,0764

0,0023

0,0059

0,0060

0 ,0 119

0,0009

— 5

0,5000

0,0458

0,0367

0,0017

0,0229

0,0090

0,0319

0,0012

— 1

0,8092

— 3

0,0620

0,0162

0,0010

0,0503

0,0122

0,0625

0,0010

0,9599

0,0788

0,0062

0,0005

0,0756

0,0155

0,0911

0,0006

2

0,9956

0,0920

—Т

0,0021

0,0002

0,0916

0,0181

0,1097

0,0002

0,9997

0 ,10 19

0,0007

0,0001

0 ,10 19

0

1,0000

0,104 8

0,0000

0,0201

0,1220

0,0001

0,0000

0,1048

0,0208

0,1256

0,0000

1

0,9997

0 ,10 19

 

0,0007

0,0001

0 ,10 19

0,0201

0 ,1220

0,0001

 

 

 

 

 

 

 

 

 

Таблица си

 

 

0 Г 0 ........................ I

0,0290 I

0,8028 |

0,1972 |

1,0000

|

0 ,0 114 I

 

 

 

 

<Э(2) =

10 ,19 72

 

Q(3) =

0,0988

 

 

 

V3 =

1 + 0,1972 +

0,0988 + 0,0688 + 0,0512 =

1,4160;

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

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

в соответствии с (5.5) и (5.6) и соотношениями, уже использован­ ными в предыдущих примерах.

В соответствии с общими условиями для всех примеров (см. п. 6.3) технологический промежуток длится 4 ч, причем в течение каждого числа операция повторяется 100 раз. Так как момент времени первой контрольной проверки не зависит от момента настройки, длительность первого межпроверочного промежутка T L является случайной величиной с равномерным распределением, принимающий значения от 0 до 100, с м. о. Т^ = 50. На основании аналогичных соображений для последнего, пятого межпровероч-

Таблица 14

брака q в примере 3

 

П овторяется предыдущий цикл

СО

°а

i y

°8

Ш - 1

III-2

Ш -З

1 1 1 -4

IV -1

IV-2

1V -3

1V-4

V -I

I V-2

V-3

V-4

0,0000

0,0000

0,0000

0,0000

0,0000

0,0000

0,0000

0,0000

0,0000

0,0000

0,0000

0,0000

0,0000

0,0001

0,0001

0,0001

0,0000

0,0001

0,0001

0,0001

0,0000

0,0000

0,0000

0,0000

0,0000

0,0002

0,0002

0,0001

0,0000

0,0001

0,0001

0,0001

0,0000

0,0001

0,0001

0,0001

0,0000

0,0005

0,0005

0,0002

0,0000

0,0003

0,0003

0,0001

0,0000

0,0002

0,0002

0,0001

0,0000

0,0010

0,0010

0,0002

0,0000

0,0007

0,0007

0,0002

0,0000

0,0005

0,0005

0,0001

0,0002

0,0002

0,0018

0,0020

0,0003

0,0001

0,0013

0,0014

0,0002

0,0001

0,0010

0,0011

0,0023

0,0030

0,0053

0,0004

0,0010

0,0020

0,0030

0,0002

0,0006

0,0016

0,0022

0,0002

0,0160

0,0045

0,0205

0,0008

0,0103

0,0031

0,0134

0,0005

0,0067

0,0024

0,0091

0,0003

0,0507

0,0061

0,0568

0,0009

0,0460

0,0041

0,0501

0,0008

0,0406

0,0032

0,0438

0,0007

0,0874

0,0078

0,0952

0,0006

0,0913

0,0053

0,0966

0,0006

0,0926

0,0040

0,0966

0,0006

0,1092

0,0091

0 ,118 3

0,0002

0 ,117 8

0,0061

0,1239

0,0003

0,1254

0,0047

0,1281

0,0003

0,1220

0,0101

0 ,13 2 1

0,0001

0 ,13 2 1

0,0068

0,1389

0,0001

0,1389

0,0052

0,1441

0,0001

0,1256

0,0104

0,1360

0,0000

0,1360

0,0070

0,1430

0,0000

0,1430

0,0054

0,1484

0,0000

0,12 20 0,0101

0 ,1321

0,0001

0 ,13 2 1

0,0068

0,1389

0,0001

0,1389

0,0052

0,1441

0,0001

мметричная

 

 

 

 

 

 

 

 

 

 

 

 

11,0000

10,0078

[0,9332

[0,0668

 

q (5 ) = 0,0512

 

 

 

 

 

 

 

 

 

 

 

 

0,0290

+ 0 ,0 114 + 0,0078 +

0,0064 +

0,0054 .) = 0,01095

 

 

 

4 - (

 

 

 

 

 

 

 

 

 

 

ного промежутка м. о. Ть — 50. Длительность остальных межпро­ верочных промежутков Т 2 Т 3 = Т4 = 100. В табл. 14 содер­ жатся вычисления, относящиеся к первому межпроверочному промежутку (гр. 3, 4 и 5) и основные циклы для всех последующих четырех межпроверочных промежутков (графы с двойной нуме­

рацией).

Втабл. 14 в гр. 2 оперативная характеристика L (у,-) совпадает

сзаписанной в гр. 2 табл. 12 в примере 2. В гр. 3 проставлены

значения функции вероятностей ф (иг) ошибки настройки, взятые из примера 1 (гр. 10 табл. 11). В гр. 4 записаны вероятности брака b (у,-) такие же, как в предыдущих примерах. В гр. 5 выполнено вычисление вероятности брака в течение первого_межпровероч-

ного промежутка с итогом 0,0290 совпадающим с q в примере 1. В гр. П-1, 11-2, П-З выполнен расчет перераспределения в соответствии с (5.5). Итогом гр. П-1 является вероятность ре­

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

межпроверочным

промежутком,

равная

Р (1) = 0,8028. Вероят­

ность настройки QO) в конце первого межпроверочного проме­

жутка равна

= 1 — Р Ю =

0,1972.

В гр. П-4 выполняется

142

расчет вероятности брака (/ \ Остальные циклы алгоритма вы­ полняются аналогично второму и в соответствии с (5 .6 ).

Ниже итоговой строки в .табл. 14 выполнены вычисления

v3 = Q(" + Q(2) -!-• • - + Q(5) = 1,4160

и 1 =

4 -

+ ? (2) + • • • +

+

- ^ - ) = 0,01095 .

Табл. 14

можно

упростить, если

при

вычислении ~q вынести

за скобку общий множитель b (и,-). Тогда вместо гр. П-4, Ш -4 и

1V-4

можно

ввести

в

конце таблицы графу

с суммами

s =

О

О

О

+

сВ4) (vt) в графу с

произведениями

= а (2) (у,-) +

а (3) (+•)

b (wf) s. Итогом последней графы можно заменить сумму гр.

П-4,

Ш-4, IV-4 при вычислении q, записанного внизу таблиц. Однако

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

ности брака <7 6) при последовательном повторении контрольных проверок, что иногда представляет интерес.

Вычисление показателей эффективности по данным примера 3 облегчается тем, что в этом отношении мало отличается от при­ мера 1. В частности, затраты на настройку R отличаются от при­ мера 1 только иным числом настроек в связи с контрольными

проверками

 

в

 

течение

технологического

промежутка,

именно

v3 =

1,416,

 

как это видно из табл.

14. В результате в соответствии

с (6 .2 ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

1,416-1,042(0,014 + 0,001)

 

 

 

 

 

 

К -------- 1

^

 

--- 1— 1— - =

0,0553 ч/ЮОО шт.

 

Затраты

 

на текущий

контроль K i

в течение технологического

промежутка

вычисляем

на

основании

(6.4)

с подстановкой J =

= 5;

 

0,029 +

6-0,007 = 0,071;

Т =

0,4:

 

 

 

 

 

 

 

 

 

■• Н И ’ 0 7 1 = 0 ,7 1 0 ч/ЮОО шт.

 

Затраты на приемочный контроль, как и в примере 1, равны

К о ~

1

ч/ЮОО

шт.

 

 

 

допусков

V в

соответствии

с (6.12)

Потери

 

на

нарушениях

с подстановкой

сь =

12,4

по условию

примера

и <7 =

0,01095

(см. табл.

14)

равны

 

 

 

 

 

 

 

 

 

 

 

 

 

F

= 0,0109

12,4 =

0,1360 ч/ЮОО шт.

 

В

итоге

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Son =

 

R + Кг +

К 2 + V

=

0,0553 + 0,7100 +

 

 

 

 

 

+

 

1,0000 +

 

0,1360 =

1,9013

ч/ЮОО шт.

 

Из

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

эффективности

S on

первых

трех примеров

(см. табл.

13)

ясно,

что

 

периодические

проверки

в течение тех­

нологического

 

промежутка

в данном

случае не

оправдываются.

144

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

Возвращаясь к примеру 3, найдем критическое значение dl потерь db на единице пропущенного в производство брака. Ана­ логично вычислениям в предыдущих примерах получим соотно­ шение

*# = -$?- = 91,7 ч/1000 шт., 9

т. е. соотношение более выгодное, чем даже в примере 1, который без отказа от приемочного контроля, как уже сказано, выгоднее варианта в примере 3. Способы оптимизации, позволяющие найти решение в подобных противоречивых ситуациях, рассмотрены в гл. 9.

Промежуточным решением между полным отказом от приемоч­ ного контроля и сплошным контролем является порядок, при котором сплошному контролю подвергается выработка за те межпроверочные промежутки, которые закончились нарушением границ регулирования. Вся остальная продукция считается при­ нятой, в том числе и за последний (в примере — пятый) техноло­ гический промежуток, в конце которого тоже требуется контроль­ ная проверка, хотя она не нужна для уточнения настройки и в при­ мере 3 не предусмотрена. Таким образом, при вычислении эффек­ тивности такого варианта надо в показатель затрат 5 0П включить затраты на выборочный контроль К.3 в конце последнего межпро­ верочного промежутка и потери из-за пропущенного брака и х. Кроме того, затраты на сплошной контроль /С2 уменьшаются, так как проверка производится только в случаях нарушения границ регулирования.

Показатель равен

Кг = сКг • 1000• Q = 0,001 • 1000*0,0674 = 0,0674 ч/1000 шт.,

где Q — ожидаемая доля продукции, подвергаемой сплошному контролю;

Q = 4

( " ^

+ Q‘3) +

+ Q(5) + Т - )

Значения Q6'),

j — 2,

3, 4,

5

взяты из табл. 14;

вычисление

Q(6) = J] а (5>(vt) L (V() = 0,0428

не показано.

 

10 В. В. ГоловннскнП

145

Затраты па выборочный контроль равны затратам на одну периодическую проверку в расчете на единицу продукции:

К3= °'029| ^ ° - — = 0,1775 ч/1000 шт

Способ вычисления Кз аналогичен К\-

Вычисление U 1 в данном случае мало отличается от вычисле­ ния V и выполняется с помощью следующих операций (см. алго­

ритм

в заголовках табл.

14; объяснение обозначений в

п. 5.2):

 

 

 

b\n =

а\'+1) b (щ);

 

 

 

 

 

 

 

 

 

 

г=12

V

 

 

 

 

 

 

 

- 7 . U )

_

 

Ш

-

 

 

 

 

Чпроп

Z J u i

>

 

 

 

 

 

 

 

 

 

 

I — —12

 

 

 

 

 

 

 

 

- ( 1 )

,

~ (5 )

|

4

 

 

 

 

 

 

' проп ~1

Ч проп

V

~ U )

 

 

 

 

vnpon

 

 

 

2

 

2 ПР0П’

 

 

где

— вероятность пропущенногобрака

за

/-й межпровероч­

ный

промежуток

при

v =

v(.

 

 

 

 

 

 

 

По данным примера 9проп =

0,0048

(вычисления не показаны);

сь — 91,7. Таким

образом

 

 

 

 

 

 

 

 

 

и г =

0,0048

91,7

= 0,4402

ч/1000

шт.

 

В

итоге

 

 

 

 

 

 

 

 

 

 

 

 

S on =

R +

7С1+7С2 +

з V

 

Ui =

 

 

= 0,0553 +

0,7100 + 0,0674 + 0,1775 -Ь

 

 

 

+ 0,1360 +

0,4402 =

1,5864.

 

 

Этот результат сопоставлен с остальными в табл. 13.

 

 

 

 

 

 

 

 

 

 

7.4.

ПРИМЕР 4

 

Как и до сих пор, в примере 4 предполагается, что ненормаль­

ности

технологического

процесса

наступают

настолько

редко,

что их опасностью можно пренебречь. В то же время в примере 4 имеет место неустранимый линейный износ настроенных элемен­ тов, что выражается в приросте w — 0,01 мм уровня настройки X в результате 100 повторений операции. Прирост отклонения у. и. v, выраженный в = 0,014 в примере равен w = 2/г, где /г = = 0,3571 — интервал округления. В технических единицах при­ рост равен 2/1(7* = 0,01 мм. Иначе говоря, статистическая законо­ мерность, в которой выражается износ настроенных элементов,

соответствует линейному

уравнению

 

v (т)

= v (0) + ах,

(7.5)

где v (т) — отклонение уровня настройки v при т-м повторении операции; т — число повторений операции после начала техно­ логического промежутка;

а = .2-°'о3о571 — 0,007142.

При неустранимом износе настроенных элементов нельзя обойтись без контрольных проверок отклонения у. и. v, кроме случаев, когда широкий допуск позволяет выделить часть его на покрытие смещения w за весь технологический промежуток. В рассматриваемом примере предполагается, что контрольные проверки выполняются раз в час. Технологически необходимые подналадки встречаются настолько редко, что длительность тех­ нологического промежутка можно без риска неправильных реше­ ний приравнять бесконечности. Настройки выполняются только в случае нарушения границы регулирования при очередных проверках. Схема настройки и все связанные с ней параметры при­ няты такими, как в примере 2.

Данные для вычисления S on в примере 4 частично заимство­ ваны из табл. 8, в которой описаны способы вычисления предель­ ных соотношений в случаях, аналогичных рассматриваемому здесь. Как видно из этой таблицы, составленной на основании исходных данных примера 4, предельная вероятность брака за

межпроверочиый промежуток равна qV]m=

0,0187. Предельная

вероятность нарушения

границ регулирования равна

QMm=

= 0,2663. Длительность

межпроверочного

промежутка

равна

Т ~ 100 операций.

 

 

 

При вычислении показателя R исходим из условий примера 2, но имея в виду, что настройка может потребоваться с вероятностью 0,2663 в конце любого межпроверочного промежутка длитель­ ностью Т = 100. Поэтому, пользуясь формулой (6.3), принимаем v3 = 0,2663 вместо v3 = 1 в примере 2 и Т = 100 вместо Т = 400.

Врезультате

R = °'2(^ 311'25 [1,042 (0,014 + 0,001) + 0,029 + 6 -0,007] =

= 0,2894 ч/ЮОО шт.

При вычисленни_потерь из-за нарушения допусков V поль­

зуемся значениями q = 0,0187 из табл. 6 и сь — 12,4 ч/ЮОО шт., как в примере 2. В результате

 

V = 12,4-0,0187 =

0,2319 ч/ЮОО шт.

 

Показатель затрат на контрольные проверки

вычисляем,

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

какие

приняты

в примере

2, т. е.

k =

0,071. Таким

образом,

K i =

= 0,071

: 0,1 = 0,710

ч/Ю00

шт.

 

 

 

10*

1 4 7

Затраты на сплошной контроль, как и в остальных примерах, приравняем 0,001 ч на штуку, что в расчете на единицу продук­ ции равняется: /С а — 1,000 ч/1000 шт.

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

5 0П -

R + Кг +

К * +

V = 0,2894 + 0,7100

+

 

-I- 1,000 +

0,2319

= 2,2313 ч/1000 шт.

 

Эту цифру

полезно сравнить со значением S on

в

примере 2,

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

примере, но

не было износа настроенных элементов технологической системы. В примере 2а (см. табл. 13) S on = 1,3627, иначе говоря, износ настроенных элементов привел к снижению производительности труда при выполнении функций, обеспечивающих качество про­ дукции в 1,7 раза (2,2313 : 1,3627). На операции без приемочного контроля 26 это соотношение достигло бы величины 2,4 [(2,2313 + + 1,0000) : 1,3627]. Конечно, эти данные соответствуют возмож­ ному частному случаю. Тем не менее, они дают некоторое пред­ ставление о том, во что может обойтись экономия па качестве инструмента, на форсировании режимов резания или на отказе от второго прохода, если при этом возникает заметное смещение уровня настройки вследствие затупления и износа резца.

Г л а в а 8

ПОИСК ОПТИМАЛЬНОГО ВАРИАНТА СРК ПРИ ЕДИНСТВЕННОМ УПРАВЛЯЕМОМ АРГУМЕНТЕ ЭФФЕКТИВНОСТИ

8.1.ПРЕДВАРИТЕЛЬНЫЕ ЗАМЕЧАНИЯ

Всвязи с приведенным в предыдущей главе примерами можно сделать вывод, что факторы эффективности системы статистиче­ ского регулирования и контроля разнообразны и многочисленны; соответственно показатель затрат 5 0П является функцией не­ скольких переменных. О многомерности показателя эффективности СРК говорилось и раньше, в первых главах, и этот вопрос будет рассматриваться в дальнейшем. Но сначала следует осветить ту особую, вполне возможную практически ситуацию, когда в ре­ зультате тех или иных ограничений число аргументов эффектив­ ности, которыми можно распорядиться, сведено к одному. Итак, стоит задача найти такое значение аргумента со — со*, при кото­ ром функция затрат S (со) достигает минимума. В предыдущих главах были показаны способы вычисления S on (со) для отдельных заданных значений со. Теперь следует решить, в какой последо­ вательности надо вычислять и как надо сравнивать значения S on(co) при различных значениях со для того, чтобы с наименьшими затратами и с достаточной точностью найти со = со*, при кото­ ром 5 0П(со) достигает минимума.

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

эффективности S on (со) в рассматриваемой задаче обладает свой­ ствами, при которых общеизвестные способы отыскания минимума неприменимы или невыгодны. Поэтому, прежде чем перейти к практическим примерам, необходимо хотя бы коротко остано­ виться на довольно разнообразных по процедуре и очень различ­ ных по эффективности методах отыскания экстремума одномер­ ной функции. Прежде всего назовем две особенности функции

эффективности 5 0П (со). Первая

особенность состоит в том, что

хотя показатель эффективности

S on (со) является дифференци­

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

ности 1 функции Son (со), когда, фигурально выражаясь, на ее

1 В работе [26] использован термин унимодальная функция.

149

«Скатах» нет «вмятин», в которых могла бы задержаться воДй (рис. 11). И если в случае одномерной S (со) существуют три точки coj, со о, соа, в которых

S 0n

(coi)

'S'on (^а)

5

(ыц);

 

ТО

со

со 2 <С ® з1

 

 

 

 

 

 

 

S 0n ((0) <

S

(соj) при всех

со •<

cot;

S 0n (ш) >

S

(со3) при

всех

со О

со3.

Рис. 11. Примеры одноэкстре­ мальной функции одного аргу­ мента

На этом свойстве одноэкстремальиой функции построены все эффективные способы отыскания экстремума.

8.2. МЕТОДЫ ПОЛНОГО И НАПРАВЛЕННОГО СПЛОШНОГО ПЕРЕБОРА

Самым простым и универсально применимым способом поиска максимума (минимума) функции / (х) одного аргумента л: является полный перебор, именуемый пассивным методом поиска или одновременным (параллельным) поиском («simultaneous search»). При небольшом числе /г возможных значений xit г = 1 , 2 , . . . , /г аргумента х полный перебор состоит в вычислении или в опреде­ лении путем эксперимента значений f (,v) при всех хс с последую­ щим сопоставлением результатов. Если х — непрерывная вели­ чина (или дискретная с малым интервалом между смежными значениями), то вычисления (эксперименты) выполняются в точ­ ках с интервалом 1га между ними (при дискретном аргументе /г0 кратно /г). Полный перебор уместен в случаях, когда число воз­ можных значений k мало.

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

эксперимент

в сельском

хозяйстве или

процесс изготовления

и испытания

механизма

с параметрами,

допускающими выбор

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

Так как несколько позже полный перебор будем сравнивать с направленными методами поиска, остановимся на оценке эффек-

150

тивности этого способа применительно к одноэкстремальным функциям непрерывного аргумента. В этой связи рассмотрим зависимость между точностью результатов полного перебора и интервалом /г0 между теми значениями х, для которых вычи­ сляется f (х). На рис. 12 сплошной линией показана одноэкстре­ мальная функция / (х) и ее значения в точках, попавших в полный перебор. На основании сопоставления полученных значений f (х) можно прийти к выводу, что min / (х) соответствует точке х = х°,

хотя, в действительности, он находится в точке х\. Следовательно, при определении точки минимума допущена ошибка б, равная б = х° х\. Но возмож­

но, что кривая f (х) ока­ жется такой, что точка

минимума совпадет с хг (нанесена пунктиром). Тогда ошибка при ее оп­

ределении

равнялась бы

с

о

*

О — X

Х2.

Таким образом, истин­ ной точкой х*, в которой f (х) достигает минимума, может оказаться любая точка в интервале (х°— /г0,

х° + /го), где 1г0— проме­ жутокмежду смежными

обсчитанными значениями Рис. 12. Остаточный интервал неопределен­ х. Этот интервал именует­ ности 3 ? ^ (т) при сплошном переборе

ся остаточным интервалом

неопределенности и при сплошном переборе обозначается (m), где т — число вычислений (экспериментов) при способе полного

перебора, после которых остался интервал (т), включающий то значение х* аргумента х, при котором f (х) достигает мини­ мума. Если приравнять единице промежуток значений х, суще­ ственный для решения задачи, то

(8 Л )

где т Г — число вычислений, включая граничные (в концах интер­ вала).

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

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

метод направленного сплошного перебора.

Поиск начинается

с вычисления / (х) в эвристической точке x lt

т. е. такой, которая

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

151

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