Файл: Кориков А.М. Математические методы планирования эксперимента учеб. пособие.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 ) . |