Файл: Кориков А.М. Математические методы планирования эксперимента учеб. пособие.pdf

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

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

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

Добавлен: 27.06.2024

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

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

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

- 1 6 2 -

Предположим сначала, что сделано К линейнс-независи- мюс наблюдений иХ = К . Тогда оценки коэффициентов можно получить из уравнения

^ -

F T

 

 

 

( 6 - 2 )

Решением'уравнения ( 6 . 2 ) является

 

 

1

т г ' -

 

'

( 6 , 3 )

 

 

 

 

 

Чувствительность решения ( 6 . 3 ) по отношению к неболь­

шим изменениям элементов вектора

обозначим через

 

 

 

 

 

 

( 6 . 4 )

s „ - j i e F " 1 4 5 =

F" 1 E=F"t

 

 

Для случая,

когда

J>f>k, умножив обе

части (6 . 2) на

Б Т , получим систему нормальных

уравнений

 

 

 

 

 

 

 

( 6 . 5 )

F T F f - F T ^ .

Система ( 6 . 5 ) аналогична системе нормальных уравнений ( 2 . 2 ) Системы уравнений ( 6 . 2 ) и ( 6 . 5 ) эквивалентны. Случайному векторуip системе ( 6 . 2 ) соответствует случайный вектор

в системе Х 6 .

5 ) . Учитывая, что

решением системы уравнений

( 6 . 5 )

является

 

'

• . 1 =

( F ' F j ' V u -

C I ,

 

 

 

( 6 . 6 )


• - 163 —

получим •

 

V £ - c 4 i - c .

( 6 . 7 )

Таким образом, чувствительность

решения ^о" по отноше­

нию к случайным возмущениям выражается ковариационной мат­

рицей С.

План эксперимента должен быть внбпан так, чтооы он обэспечнвал мпнтк:тьнуа чувствительность решения (6.6) к »лучайвш возьгугданаяи. Ясно, что чувствам-1 льность в таком

взде, как она определена в формуле (6.7) длох« подходит для этой цели, т-к как не дает возыогностл сравнивать в число­ вом виде чувствительность решения системы (6.6) при разных

планах F . Поэтому в качестве числовой характеристики чув­

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

Е - о п т и м а л ь н о с т ь . План называется Е -опти-

гальнш, если максимальное характеристическое число соответ­ ствуйте* ему ковариационной матрицы миншально.

Свойства «пенок параметров можно характеризовать свойст­ вами соответствующего эллипсоида рассеивания. Пусть величины

 

 

имеют распределение с ттоматгческими ожи­

даниями

j i , , {Ьг

jb^ « ковариациями

С ц ( i , j = 1 , 2 , . . .

Я. ) .

Рассмотрим Я -мерный эллипсоид, центр

которого сов­

падает с

точкой,

имеющей координата J * t ,

. .

. , ^

- 1 6 4 -

Допустам, что в этой области плотность вероятности посто­

янна. Ковариации

этого равномерного

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

можно

подобрать

гак, чтобы

они

совпадали о

Сц.

Величины

j » t ,

J a t , . . . ,

к и

C t - (

i . ,

j = 1 , 2 , . . . , К )

полностью

опре­

деляют параметры

эллипсоида, который

и называется

эллипсои­

дом рассеивания.

 

 

 

 

 

 

 

Отсюда и из

определения Е-оптимального плана

следует,

что Е -оптимальный план минимизируетмаксимальную

ось эллип­

соида рассеивания оценок

параметров.

 

 

 

 

А - о п т и м а л ь н о с т ь . План называется <А-опти- мальным, если его ковариационная матрица злеет наименьший след (сумму диагональных элементов).

© - о п т и м а л ь н о с т ь . План называется^ -опти- 7лальным, если ему соответствует ковариационная матрица с

наименьшим определителем (или, что то же, шбормационная мат­ рица с наибольшим определителем).

Определитель ковариационной матрицы с&1зан с эллипсоидом

рассеивания

следующим образом

Ссм.,иаприыер, книгу Г.Краме-

р а [ 1 б ] )

\

%

' .

'

где Г {\

+ \) - гамш-функция,

d a t

С- определитель ковариа­

ционной

матрица.

 

 

Таким образом, S -оптимальный

план минимизирует объем

зллшбоидА рассеивания оценок

параметров.


Критерии

Е - , Ai-оптимальности

оказываются экви-

валенткшли в

случае, если ковариационная матрица оценок

коэффициентов имеет наиболее просто! вид: С = 1 Е , где

^ - константа.

В случае, когда основная цель эксперимента - поиск

экстремума поверхности отклика, очень существеанш являет­

ся вопрос о величине дисперсии предсказанных значений рег­

рессионно?! функции в заданной области. Оказалось, что если

ввести обобщение понятия плава, то $ -отвмальввй план бу­

дет .совпадать с планом минимизирующим максимальную диспер­

сии

предсказанного

значения параметра оптимивации в задан­

ной

области планирования. Таков план называется G - о а -

т и м а л ь н ы м ,

ему соответствует наименьшее максималь­

ное

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

го значения параметра оптимизации.

Пусть регрессионная функция по-прежнему имеет вид:

а через

^ обозначим лучшую линейную оценку параметра •

( t

=

1 , 2 , . . . ,

Тогда для дисперсии предсказанных зна­

чений

получим

 

- 1 6 6 -

— E E b&Jii&Mh-

M V M -

При выводе этого соотношения ш воспользовались выражени­ ем ( 6 . 6 ) к ввели обозначение

Таким образом, & -оптимальным буде. план, мининизирупций'на множестве всех ШЕНОВ величину

 

m a x

Т . С ^ ( х ) ,

 

 

 

 

 

 

 

.

Х б Х

 

 

 

 

 

 

 

 

г д е . Х

- область

планирования.

 

 

 

 

 

 

§ 6 - 3 . О п о н я т и и

" н е п р е р ы в н ы й

 

 

 

 

 

 

 

п л а н "

 

 

 

 

 

Обычный, ЕЛИ "точнкч",

план представляет

собой мноаест-

во точек x . ^ ( i » =

1 , 2 , . . . ,

m,) в

области планирования с

за­

данным числом

п . - (1=

1 , 2 , . . . ,

т Л

наблюдений

в каядой

из

них. Общее количество

наблвдений

л =

п..^

. Обозначим

долю наблюдений из общего числа

Я

,

которая

 

приходился

 

i - у »

точку

I

. =J}±

. Д л я точного

плана

должны выпол-

няться

условия:

 

 

 

 

 

 

 

 

 


 

 

- 1 6 7 -

 

 

 

5 § ' ^ '

 

 

( 6 * 8 )

где л,^ -

целые числа.

 

 

 

 

Задачу построения точного плана можно, очевидно, рас­

сматривать

как задачу отыскания гн.* X

точек

х 4

( 1 = 1 , 2 , . . . ,

. . . . т,) и

соответствующих им значений

§ . ( . 1 =

1 , 2 , . ; . , т . ) ,

обладающих

свойствами ( 6 . 9 ) ,

( 6 . 1 0 ) . указывается, что зада­

ча построения точных""планов с

^я^тща

Jf

, удовлетворяю­

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

SD -оптимальных.точных планов, настолько сложна, что даже в

'простейшем случае - в случае полиномиальной регрессии на от­

резке - она пока еще очень датека от окончательного решения.

Задача построения об -оптимальных планов становится более дос-

тупной, если несколько обобщить понятие плана.

Американским статистиком Кифером введено понятие "непре­

рывного" плана. Оно не связано с заданным количеством наблю­

дений ;мало того, оно вообще не имеет непосредственной связи б обычным понятием наблюдения. Задание непрерывного плана

эквивалентно заданию некоторой неотрицательной функции / ь ^ ) ,

определенной на подмножествах множества X , такой, что

В частности, функция я- ( Л ) может быть такой, что для её ' полного определения достаточно задать её значения в конеч­

ном числе точек. В этом случае непрерывный план требует за-


 

 

 

 

 

- 1 6 8 -

 

 

 

 

 

 

 

 

дания только

множества

точек

 

х -

( i 1, , 2 . , ' . . . ,

m )

 

 

 

к доли набявдений в каждоЗ точке

Ц

.

,

причем

выполняется

свойство ( 6 , 1 9 ) . Вьтаолнения

свойства "(6.10)

не

требуется.

Таким образом, для некоторого заданного Л

 

непреравтй план

может требовать нецелого числа наблюдений в соответствующих

точках. Функция у.

(dt

) представляет

собой

вероятностную

ме­

ру, определенную на множестве

X

,

и задание непреравного

плана эквивалентно

заданию вероятностной

мера на X

. Оказа­

лось, что построение SD -оптимальных планов на множестве неп­

рерывных планов -

задала хотя

и слолшая,

но

в некоторых слу-

чаях разрешимая.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обозначим через Jlt("Jf§ )

информационную матрицу

непрерывно­

го плана с элементами

 

 

 

 

 

 

 

 

 

 

 

 

 

1

J

J H - t l

J

J

^

 

'

s

r

T

 

 

 

( 6

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

конечное

число

точек

т . , ЭТОТ интеграл превращается

в

сумму

 

 

 

 

 

 

i - ^ f v C ^ e ) f j C 2 c ) g t .

 

 

 

 

 

 

(

6 Л 2 )

В работах Кифера показано, что для области

X

,

обладаю­

щих некоторыми свойствами, которые практически всегда выпол­

няются, и для таких функций

 

^. ( х )

,

' j =

1,2, ... ,

К,

 

 

которые непрерывны и линейно-независимы, любую информационную

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

сосредоточенную

в конечном числе точек, т . е . как

матрицу

с

элементами

( 6 . 1 2 ) .