Файл: Бошняк, Л. Л. Измерения при теплотехнических исследованиях.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 130

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Последнее выражение аналогично (XIV.4) и, следовательно, сви­ детельствует (с учетом сделанных выше допущений) о правомерности использования для оценки cLметода наименьших квадратов, изло­ женного в предыдущем параграфе. Многочлен

Yo + Yi* + • ’ + Yтх т,

сглаживающий экспериментальные данные, проведенный методом наименьших квадратов, называется линией параболической регрес­ сии порядка т. Интуитивно ясно: чем меньше т, при котором откло­ нения наблюдений от многочлена можно считать чисто случайными, тем лучше. Надо начинать с небольшой степени (т — 1; 2, вряд ли более 3) и увеличивать т лишь при наличии достаточных экспери­ ментальных оснований. Если т = 1, то прямая линия у 0 + у гх называется линией регрессии или линейной регрессией. Если пере­ менная х, характеризующая условия опыта, меняется в небольших пределах, то многие зависимости можно приближенно считать ли­ нейными. Поэтому проведение линии регрессии является одним из наиболее распространенных приемов первичного осмысливания ре­ зультатов эксперимента. Широкое применение линия регрессии находит, в частности, при градуировке шкал измерительных при­ боров.

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

ний у ъ .

. ., уп взаимно

независимые

случайные

величины,

под­

чиняющиеся нормальному распределению с параметрами

М [У;] =

= ту (xt),

D (Y l ] = a 2,

i = 1, . . ., п;

значения

x lt

. . ., хп

не­

случайные известные величины, функционально не связанные между

собой; линия

регрессии определяется

на

отрезке c ^ x ^ . d , где

с и

d — заданные числа;

коэффициенты

регрессии

неизвестны.

Обычно полагают у 0 = а,

у г =

b ах

и ищут оценку для линей­

ной

функции

у = а (х х) +

Ь. В

этом случае

М (УД =

= a (xt х) + Ь.

Согласно методу наименьших квадратов, эффективные оценки

для линейной функции у имеют вид

[16]

 

 

у — <х(х х)

13,

(XIV.9)

где

 

 

 

П

 

' _

п

^ y i i x i —x)

* = Ч

. Р = у = — Ъ У г

s

to — х)2

 

 

1=1

 

 

 

Здесь а и р — соответственно оценки а и Ь. Доверительные интер­ валы для а и b определяются уравнениями:

аа =

а t2sa\ 1

Ьн== р — t2Sfi\ |

ав =

а + t2sa- J

(XIV.10)

Ьв = р + t2sр, |

431


где

12 — табулированное

значение аргумента

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

Стью-

дента, определяемое

по

таблицам

(см.

 

[16]),

имеющим

два

входа

v =

п — 1 и Q =

0,5 (1 — р 2) 100%,

р 2 — доверительная

вероят­

ность.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Доверительная

область функции

у

= а (х — х) + 6

с

коэффи­

циентом доверия

р 3 ограничивается графиками функций

 

 

 

 

 

а (х — х) - f р ±

ип_2 (ра. A,) s

У

 

 

(х — х)2

+ т

(XIV.11)

 

 

X, (*i *)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

/ Ц '

 

 

 

1 +

пСР

 

__

 

 

 

 

 

 

 

 

 

К О +

пС2) (1 +

nD2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С ==

 

 

С— JC

 

 

D = -

 

d х

 

 

 

 

 

 

 

 

 

 

(* ;- * )* ’

 

 

£

 

(Xi-X)2

 

 

 

 

 

 

 

 

 

 

 

 

 

у

 

 

 

 

 

un_2 {pz,X) — табулированная функция (см. [16]);

 

 

 

 

 

 

 

 

 

 

 

:= ]'

/

 

 

(=1

 

 

 

 

 

 

 

 

 

 

•— оценка

квадратичного

 

отклонения.

 

2 (р2,

A,) s

заменяется

на

«со

Если

параметр

а

известен,

 

то

 

 

(р г, А) о. Если нужно построить не доверительную зону,

а лишь

два

совместных доверительных интервала на отрезке

с ==s;

х ^

d

для у (с)

и у (d),

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

соединяющие точки с координа-

 

Л

 

 

 

 

 

А

 

 

R

(d) ]. Тогда у

 

 

_

Ь

тами [с, у (с) ± R

(с) ] и [d, у (d) ±

 

= а (х х) +

будет заключена в полосе,

ограниченной этими прямыми.

Величина

 

 

R (х) =

и„_2 (pa,

X)s

 

 

 

(х — х)2

 

 

 

(XIV. 12)

 

 

 

 

 

 

(Xi — X)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

( р А- ) — табулированная

функция

 

(см. [16]).

 

нужно

Если

параметр

а

известен,

 

то

в

 

формуле

 

(XIV. 12)

Vn- 2

2.

A,) s заменить на

Va, (р8,

А,)

о.

 

 

 

 

 

 

 

 

 

 

Квадратичное отклонение а и р

оценивается следующими выра­

жениями:

 

 

 

_____ S______ .

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

s0

 

 

 

(XIV. 13)

 

 

 

 

 

у ^ { Х1- х ) 2

 

Vri

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Использование изложенного метода рассмотрим на следующем

примере

[16].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

432


 

Пусть n = 20,

a (хп

t/f) заданы:

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

1

 

2

 

 

3

 

4

5

 

6

 

7

 

8

 

9

 

10

 

 

Xi

 

 

 

0

 

0

 

 

0

 

0

0,1

 

0,1

 

0,1

0 ,5

 

0 ,5

1,0

 

 

Ui

 

 

 

3 ,3 2

 

4 ,4 9

 

2,01

5 ,0 7

0 ,2 7

3 ,5 9

 

3 ,4 9

5 ,0 0

 

0 ,4 4

5 ,4 2

 

 

i

 

 

 

 

11

 

12

 

13

 

14

15

 

16

 

17

18

 

19

 

20

 

 

XI

 

 

 

 

1,0

 

1,5

 

1,5

1,9

1,9

1,9

 

2 ,0

2 ,0

 

2 ,0

2 ,0

 

 

Ус

 

 

 

7 ,8 7

 

5 ,6 4

1,27

7 ,4 5

7 ,5 8

0 ,0 0

8 ,3 2

10,25

6 ,5 6

4,21

С о гл асн о

ф о р м у л а м

(X IV .9 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ов =

т 1 ж == 1-599^

1’60;

Р=:У = - ^ Г

 

= 4.609 ^

4,61;

 

 

 

 

 

 

 

 

 

 

 

 

у = 1,60 (л: — 1) -)- 4,61

 

 

 

 

 

 

 

п о стр о и м п о в е р и т ел ь н у ю

зо н у д л я

л и н е й н о й ф у н к ц и и

у = а (х 1) +

b на

о т р езк е

— 3

^

л

: ^

5

.

 

С о гл а с н о

ф о р м у л а м

( X I V .11)

о ц е н к а

к в а д р а т и ч н о го

о т к л о н е н и я а

 

 

S =

i

у

f

 

123,23

 

 

0 й , л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

----j-|— == 2,616;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

— с = D =

 

 

4

.

 

 

 

 

 

 

 

 

 

 

 

 

^

=

~ У

13,86 +

32 0

~

°> 9 7 9 0 -

 

 

 

 

 

 

 

 

 

 

 

П о л а г а я р 2 = 0 ,9 по т а б л . 4 .6 а [1 6 ],

 

 

 

 

 

 

 

 

 

 

 

н ах о д и м

ы 18(0,9;

0 ,9 7 9 ) =

2 ,2 8 4 -0 ,4 2 +

 

 

 

 

 

 

 

 

 

 

 

+ 2,291

-0 ,5 8

 

=

2 ,2 8 8 ;

п о это м у с

к о эф ­

 

 

 

 

 

 

 

 

 

 

 

ф и ц и ен то м д о в е р и я р 2 — 0 ,9 м о ж н о

 

 

 

 

 

 

 

 

 

 

 

у т в е р ж д а т ь , что н а о т р е зк е

— 3 ^

 

5

 

 

 

 

 

 

 

 

 

 

 

н е и зв е ст н а я п р я м а я л и н и я р е гр е сс и и

 

 

 

 

 

 

 

 

 

 

 

за к л ю ч е н а м еж д у г р аф и к а м и ф у н к ц и й

 

 

 

 

 

 

 

 

 

 

 

 

 

1,60 (х— 1) 4- 4,61 +

 

 

Р и с .

173. Д о в е р и т е л ь н ы е

зо н ы

д л я

л и н и и

±

5,99 V 0,05 +

0,072 — I)2.

 

 

 

 

р е гр е сс и и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С о о тветств у ю щ ая д о в е р и т е л ь н а я

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

на р и с .

173.

Е сл и

бы в

этом п р и ­

м ере

б ы л о и зв е с тн о ,

что

о

=

2 ,

то

п р о и зв ед ен и е ип — 2

( р 2,

X) s — 5 ,9 9

с л е д о в а л о бы

за м е н и т ь н а «ш

(р2,

X) а = 4 ,2 9 .

 

 

 

 

 

 

 

 

 

 

 

 

 

П о с тр о и м т е п е р ь д о в е р и т е л ь н у ю о б л а с т ь с п р я м о л и н ей н о й гр а н и ц е й . Т а к к а к

п о т а б л . 4 .6 . [1 6 ]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vv (p2,

А) =

 

i/18 (0,9;

0,979) =

vv (р2;

V

1 — А2 ) =

 

 

 

 

 

 

 

=

 

% (0,9;

0,204)= 1,888-0,92+ 1,920-0,08= 1,891,

 

т о с о гл а с н о

 

ф о р м у л е

(X IV . 12)

с

к оэф ф и ц и ен том

д о в е р и я

0 ,9

м о ж н о

у т в е р ж д а т ь ,

ч то

н а о т р е з к е

— 3 ^

х

 

5

н е и зв е с т н а я п р я м а я за к л ю ч е н а

в полосе

 

 

 

 

 

 

 

1,60

1) +

4,61

±

1,891

-2 ,6 1 6

]А о ,05 +

0 ,0 7 2 2 -1 6 .

 

 

 

 

28

 

Л. Л.

Бошняк

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

433


Доверительная зона с прямолинейными границами изображена на рис. 173 штриховыми линиями. Если бы было известно, что а = 2, то произведе­ ние уи—2 (р2, Я) s = 4,947 следовало бы заменить на Усо (р2, Я) а — 3,57. Квад­ ратичное отклонение а и (3 вычисляется по формулам (XIV.13):

 

%

2,616

 

«3 =

2,616

 

 

 

 

0,71;

 

0,58.

 

 

 

 

1/13786

 

 

1/20

 

 

По формулам (XIV.10) определим доверительные интервалы для а

и b при р 2 =

= 0,90 (Q =

5%):

 

 

 

 

 

 

 

 

ан =

1,60— 1,729-0,71

= 0 ,3 7 ;

1

Ьи =

4,61

— 1,729-0,58 =

3,61;

|

ав =

1,60 + 1,729-0,71

= 2,83;

|

&в =

4,61

+ 1,729-0,58 =

5,61.

J

4. Оценка точности прибора при линейной градуировочной зависимости

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

Рис. 174. Схема включения и характеристики приборов при их проверке по методу трех приборов

Для этого случая в настоящее время известны две методики обработки данных, отвечающие условиям рассматриваемой задачи; метод трех приборов [18] и метод Вальда [79], [180]. Оба метода позволяют оценить дисперсии нескольких приборов, измеряющих одновременно один и тот же параметр х, по сигналам этих приборов, зарегистрированным независимым образом. Математическая модель, используемая в обоих методах, одинакова. Предполагается, что систематические погрешности полностью исключены (или учтены), а полученные отклонения — случайны.

Метод трех приборов заключается в испытании трех приборов, одновременно измеряющих случайную величину X (рис. 174). Для

434


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

У2 = А12Ух “Ь В12', у 3 = АхзУх В 19.

(XIV.14)

Так как коэффициенты уравнений (XIV. 14)

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

тодом наименьших

квадратов и знак равенства справедлив в сред­

нем, то для конкретных реализаций у и , у 21,

у ы уравнения (XIV.14)

принимают вид:

 

 

У и = АтУи

В 12 -f- А 2it» Уы = АхзУх1

“Ь 5 i3 ~Ь Asi/i

где А2 х,-, Азх; — соответственно разница между показаниями при­ боров. Выразим величину показания прибора через ее математи­ ческое ожидание и погрешность: у { = М (г/() + е,, где М (г/,) — математическое ожидание г/(., et — погрешность. Тогда последние уравнения запишутся так:

 

М (у2)

+

8 2

А 12[М (уг)+ 8 Х]

+ В 12 +

Д21; |

 

М (Уз) +

е 3 =

Л х з

(г/х)+

е х]

+

В 13 +

А 31;

J

 

А4 (у2)

А 12М

(г/х)

В 12=

 

е2

+ A 12&i + А21;1

 

А4 (Уз)

 

A 1SM (г/х)

^хз

=

 

 

ез ~Ь ^хзех “Ь А31.J

С

учетом

(XIV. 14)

левая

часть

последней

системы

уравнений

равна

нулю

и,

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

 

 

 

 

 

 

 

 

 

 

 

 

 

А 2 х

=

ег

А 12&1,

1

 

 

 

 

 

 

 

 

^31

=

83 ^хз81- J

 

 

 

Вычислим начальные моменты второго порядка

М (Д 2 1 ),

М 31) и

М 2 1 А31),

учитывая,

что

погрешности

 

 

е2

и

е3 независимы и

поэтому М [ех8 2] = М [8 x8 3 ] = М

2е3]

= 0.

 

 

 

Получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

М (All)

=

М [(е -

А 12ех)2]

= D (е2) +

A2nD (ej;

 

М (Ali)

=

М [(е3 -

 

А гзе,)2]

=

D (е8) +

А \ф (ej;

А4 (А 21 Азх)

=

А4 [(е 2

А 12^х) (ез —

А 1381) ]

А 12А 13D (sjJ

• С другой стороны статистические оценки тех же моментов равны

S Д2и

S Лз1г

М ( All) = ^ —

; М ( & ) = ! = L - — ,

Tj ^21^31

1=1

М(А21Ьз1) п

28*

435