• J = 6 вимірювань у центрі пластини;

  • К = 6 днів;

  • L - 2 циклів.

Модель для трирівневого ієрархічного експерименту має такий вигляд:

^Ikj =Ц + 7/ + ^>Ік + slkj<

де / = 1,2, А-= 1,..., 6,/у = 16.

  1. Експеримент повторюють на кожній з М = 5 пластин, вибраних навмання з набору пластин. Ці контрольні еталони позначають в аналізі як пластини 138, 139, 140, 141 і 142. Експеримент також повторюють для Q = 5 зондів, що позначають в аналізі як зонд № 1, № 281, № 283, № 2062, № 2362. Дані включають:

— номер циклу;

  • ідентифікацію пластини;

  • місяць і день вимірювання;

  • оператора;

  • температуру;

  • середнє з шести вимірювань у центрі пластини;

  • короткотривале (повторюваність) стандартне відхилення з шести вимірювань.

  1. Аналізування повторюваності, щоденних і довготривалих складових

Точність зонда № 2362 звіряють з іншими зондами та перевіряють на відповідність позначен­ням на графіку стандартних відхилень повторюваності для кожного зразка та ідентифікаційного номера пластини. Визначають, що зонд № 2362 має найкращу прецизійність протягом циклу 1 і тому його вибирають для сертифікаційних вимірювань. Проте прецизійність не є ні єдиним, ні обов’язковим найважливішим внеском до невизначеності, у чому можна переконатися далі.

Вплив пластини розглядають як незмінний. Трирівневий ієрархічний експеримент із впливом перетину, обумовленим пластинами, має змішану модель і відрізняється від моделі у 5.2.2. У табли­ці 9 наведено оцінювання з ANOVA.

Таблиця 9 — Таблиця ANOVA для даних аналізу міри

Джерело

Ступені свободи, V

Середні квадрати, MS

Очікуваний середній квадрат

Цикл

1

MSR= 0,009 198

о2 +5сд +30од

День (Цикл)

10

MSD(R) = 0,003238

о2 + 5о q

Похибка

44

MSE = 0,0008046

о2

Примітка. У цьому прикладі перехресні впливи не розглядають. Такі складові додають, якщо очікується, що одне джере­ло може впливати на інше. Прикладом може слугувати аналіз складової двох медикаментів, в якому один медикамент може сповільнювати дію іншого.

Отримують також складові дисперсії, обумовлені складовими повторюваності, днів та циклів:



б2 = S2=MSE =0,0008046

MSRtRx—MSE
v2D=S2D= = 0,0004866

5

, „ MSp,-MSDlR}

dR=SR= — = 0,0001987.

R R ЗО

  1. Зміщення, пов’язане із зондом

    1. Графічний аналіз

Прецизійність зразка обговорено у 8.1. Підрозділ 8.4 присвячено зміщенню, пов’язане із зон­дом. Графічний аналіз дає відносні зміщення між п’ятьма зондами. Для кожної пластини відмінності від середнього пластини зображуються на графіку разом з ідентифікаційним номером пластини для кожного зразка. Графіки підтверджують, що зонд № 2362 (під індексом 5) є зміщеним відносно інших зондів, тобто завжди дає нижні значення. Зміщення стає більш очевидним, коли зонди по­чинають використовувати зонди (цикл 2).

Номер індексу зонда:

  • 1 = № 1;

  • 2 = № 281;

  • З = № 283;

  • 4 = № 2026;

  • 5 = № 2362.

  1. Кількісний аналіз

Поправки для п’яти зондів наведено у таблиці 10 як відмінності від середнього для кожної пластини. Для певного циклу bkj= Ykj-Y,, кожний елемент є, як визначено у 5.5.2 для k = 1,2,3,4,5 (зонд № 1, № 281, № 283, № 2026, № 2362) та і = 138, 139, 140, 141, 142.

Таблиця 10 — Зміщення для п’яти зондів

Пластина

Ідентифікаційний номер пластини

№ зонда

Цикл 1

Цикл 2

138

1

1

0,0248

-0,0119

281

0,0108

0,0323

283

0,0193

- 0,0258

2 062

- 0,0175

0,0561

2 362

-0,0372

-0,0507

139

2

1

-0,0036

- 0,0007

281

0,0394

0,0050

283

0,0057

0,0239

2 062

-0,0323

0,0373

2 362

-0,0094

- 0,0657

140

3

1

0,040 0

0,0109

281

0,0187

0,0106

283

- 0,0201

0,0003

2 062

- 0,0126

0,0182

2 362

- 0,0261

- 0,0398

141

4

1

0,0394

0,0324

281

-0,0107

- 0,0037

283

0,0246

- 0,0191

2 062

- 0,0280

0,0436

2 362

-0,0252

- 0,0534

142

5

1

0,0062

0,0093

281

0,0376

0,0174

283

- 0,0044

0,0192

2 062

- 0,0011

0,0008

2 362

- 0,0383

-0,0469



  1. Трактування зміщення, обумовленого зондом

Для сертифікації було вибрано зонд № 2362 через його кращу прецизійність, проте його зміщення відносно інших зондів створює проблему. Існує дві можливості для вирішення цієї про­блеми:

  1. виправити усі вимірювання, які проведено із зондом № 2362, відносно середнього зондів;

  2. увести в оцінювання невизначеності стандартне відхилення для різниці між зондами.

Кращим вибором є а), якщо можна припустити, що зонди в аналізі являють собою випадкову вибірку зондів цього типу. Це є особливо правильним, коли одиницю (питомий опір) визначають методом випробування. Отже, стратегія в цьому випадку полягає у введенні в усі вимірювання поправки на середнє зміщення зонда № 2362 та обчисленні стандартного відхилення поправки як оцінки невизначеності за типом А.

У таблиці 11 наведено поправки для зонда № 2362 з таблиці 10 відповідно до пластини. Усі поправки є від’ємними. На основі перевірки знака ми розглядаємо поправки з двох циклів як такі, що походять з однієї сукупності.

Таблиця 11 і— Зміщення для зонда N2 2362 відповідно до пластини

Пластина

цикл 1

цикл 2

Усього

138

-0,0372

- 0,0507


139

-0,0094

-0,0657


140

-0,0261

-0,0398


141

- 0,0252

- 0,0534


142

-0,0383

-0,0469


Середнє значення

- 0,0272

-0,0513

- 0,0393

Дисперсія (10 значень)

0,0002618

Стандартний відхил

0,01618



Зміщення зонда № 2362 оцінюють як мінус 0,0393. Невизначеність зміщення зонда № 2362 складає:

= 0101618=0,005117

°2362 J'] Q

  1. Зміщення, пов’язане із конфігурацією

    1. Вимірювання з двома конфігураціями зонда

Сертифікаційні вимірювання проводили із зондом № 2362, з’єднаними у певній конфігурації, яку називають конфігурація А, з багатьох можливих конфігурацій. Аналіз міри проводять із зондом № 2362, з’єднаним двома способами для ідентифікації великих розбіжностей. Необроблені дані про конфігурацію схеми з’єднання не наводять у цьому стандарті.

  1. Графічний аналіз

Різниці між вимірюваннями у конфігураціях А та В, проведених в один і той самий день, наве­дено на графіку для 6 днів для кожної пластини. На рисунку 5 представлено два цикли з проміжком приблизно два місяці. Лінія в центрі є нульовою лінією. Результати вимірювання достатньо випад­ково розкидані вище та нижче нульової лінії, не позначаючи жодної різниці між конфігураціями для зонда № 2362. Отже, внесок до невизначеності з цього джерела вважають незначним.

  1. Кількісне випробування для різниць між конфігураціями

Результати підтверджують (-випробуванням на основі таблиці 12, де середні різниці та стан­дартні відхилення обчислюють із вимірювань для 6 днів на п’яти пластинах. (-статистика для пере­вірки, чи відрізняються статистично засоби, приймає значення менше 2, що означає відсутність значної різниці. Можна зробити висновок, що зміщення, обумовлене конфігурацією схеми з’єднання та незначним внеском до невизначеності від цього джерела, немає.

Таблиця 12 — Різниці між двома конфігураціями

Статус

Середнє значення

Стандартне відхилення

Ступені свободи

t

Цикл 1

- 0,00858

0,0242

29

1,9

Цикл 2

- 0,0110

0,0354

29

1,7

Таким чином, невизначеність, обумовлена конфігураціями, Sconfjg, дорівнює нулю (Sconfig= 0).

Пояснення:

X — час, дні;

У — різниці питомих опорів між вимірюваннями, Ом • см;

1 5 — ідентифікаційні номери пластин (див, таблицю 10),


Рисунок 5 — Різниці між вимірюваннями у конфігураціях А та В відносно часу





































  1. Обчислювання невизначеності

    1. Загальні положення

Обчислювання невизначеності залежить від опублікованого значення. У цьому випадку серти- фікованим значенням питомого опору для кожної пластини є середнє з шести короткотривалих повторень, проведених у центрі пластини в один день, в один цикл.

З моделі в розділі 8, тобто:

Ylkj1 + 7/+5//< + £lkj-

Значення для певної пластини в день к та в цикл І дорівнює:

Y Ік. — ц + у'/ + 8lf<+ к/к*

  1. Бюджет невизначеності

Бюджет невизначеності наведено в таблиці 13.

Таблиця 13 — Бюджет невизначеності

Джерело

Тип

Дисперсія

Повторюваність

А

S2 = 0,0008046

Від дня до дня

А

Sq = 0,0004867

Від циклу до циклу

А

S^= 0,0001987

Зонд № 2362

А

Sprode= 0,00002618

Конфігурація схеми з'єднання

А

^config =$



  1. Сумарна стандартна невизначеність

Сумарну стандартну невизначеність, включаючи складові для повторюваності, днів, циклів та зміщення зонда, обчислюють із бюджету невизначеності так:

и = JS2+S?)+S^+S? =0,03894.

с V и н Ь2362

  1. Ефективні ступені свободи

Ефективні ступені свободи, пов’язані з ис, апроксимують за допомогою формули Велча- Саттерсвейта. Сумарну невизначеність переписують із 8.1 так:

24-MSe, MSD{R]MSrе2

c 5 6 30 ^2362

Відповідні ступені свободи для цих складових дорівнюють 44, 10, 1 і 9, відповідно. Ефективні ступені свободи відповідно до формули Велча-Саттерсвейта:

= 17


MS2tt о 4

ЗО2 І фзб2

1 9


цс
^^D(R)
б2
10


42-MSI

52

44



Критичне значення на рівні ймовірності 0,95 із 17 ступенями свободи з ї-таблиці дорівнює 2,11, і таким чином розширена невизначеність:

U = 2,11 с = 0,082 Ом ■ см

ДОДАТОК A
(обов’язковий)
ПОЗНАКИ

bk ■—оцінка зміщення для к-го ЗВТ;

S2 — оцінка дисперсії похибок;

S| — оцінка дисперсії зміщення для к-го ЗВТ;

S2onfg — оцінка дисперсії, обумовленої конфігурацією;

Sj -— оцінка складової дисперсії для складової дня;

S2un — оцінка складової дисперсії для складової циклу;

MS — середній квадрат;

MSD — середній квадрат для складової дня для дворівневого ієрархічного експерименту;