Н.Н. Калиткин - Численные методы (1133437), страница 106
Текст из файла (страница 106)
Сравним дисперсию совокупности измерений с' зз с дисперсиями отдельных серий с) =з'. Если окажется, что при выбра1»- ном уровне достоверности р, можно считать о) о, для всех (, то зависимость г от х имеется. Если достоверного превышения ") С такой сиязыо мы уже ястречались я стохастических задачах нахождения корня уравнения (гл. Ч,$2, п.4) и минимумз функции (тл.
ЧП, Е!, п. 4). 49б СТАТИСТИЧЕСКАЯ ОБРАБОТКА ЭКСПЕРИМЕНТА !ГЛ. ХР нет, то зависимость не поддается обнаружению (при данной точности эксперимента и принятом способе обработки). Дисперсии сравнивают по критерию Фишера (30). Поскольку стандарт з определен по полному числу измерений А!, которое обычно достаточно велико, то почти всегда можно пользоваться коэффициентами Фишера г" (со, т; р,), приведенными в таблице 25. Второй способ анализа заключается в сравнении средних г! при разных значениях х; между собой. Величины г! являются случайными и независимыми, причем их собственные стандарты выборки равны ББ! =з!6 а!.
Поэтому их сравнивают по схеме независимых измерений, описанной в п. 3. Если различия г! значимы, т. е. превышают доверительный интервал, то факт зависимости г от х установлен; если различия всех г! незначимы, то зависимость не поддается обнаружению. Многофакторный анализ имеет некоторые особенности. Величину г (х, у) целесообразно измерять в узлах прямоугольной сетки (хь у,) „чтобы удобнее было исследовать зависимость от одного аргумента, фиксируя другой аргумент.
Проводить серию измерений в каждом узле многомерной сетки слишком трудоемко. Достаточно провести серии измерений в нескольких узлах сетки, чтобы оценить дисперсию единичного измерения; в остальных узлах можно ограничиться однократными измерениями. Дисперсионный анализ при этом проводят по первому способу. 3 а м еч а н и е 1. Если измерений много, то в обоих способах отдельные измерения или серии могут с заметной вероятностью довольно сильно отклониться от своего математического ожидания.
Это надо учитывать, выбирая доверительную вероятность ро достаточно близкой к 1 (как это делалось в п.2 при установлении пределов, отделяющих допустимые случайные ошибки от грубых). А н а л и з р е г р е с с и и. Пусть дисперсионный анализ указал, что зависимость г от х есть. Как ее количественно описать? Для этого аппраксимируем искомую зависимость некоторой функцией г(х) 1(х, а), а=[а„аэ ..., а ',.
Оптимальные значения параметров аА найдем методом наименьших квадратов, решая задачу ~ч, '!Е(х!) [г! — 1(х!, а)1'=ппп, (34) где и!(х,) — веса измерений, выбираемые обратно пропорционально квадрату погрешности измерения в данной точке (т, е.
и!! (Ог!)-!). Эта задача была разобрана в главе 11, 9 2. Остановимся здесь лишь на тех особенностях, которые вызваны присутствием больших случайных ошибок. гл, хч) СТАТИСТИЧЕСКАЯ ОБРАБОТКА ЭКСПЕРИМЕНТА 497 Вид 1(х, а) подбирают либо из теоретических соображений о природе зависимости г(х), либо формально, сравнивая график г(х) с графиками известных функций, Если формула подобрана из теоретических соображений и правильно (с точки зрения теории) передает асимптотику г (х), то обычно она позволяет не только неплохо аппроксимировать совокупность экспериментальных данных, но и экстраполировать найденную зависимость на другие диапазоны значений х.
Формально подобранная функция ):(х, а) может удовлетворительно описывать эксперимент, но редко пригодна для экстраполяции. Проще всего решить задачу (34), если ) (х, а) является алгебраическим многочленом ~ч~~ аех". Однако такой формальный выбор функции редко оказывается удовлетворительным. Обычно хорошие формулы зависят от параметров нелинейно (транс!(едентная регрессия). Трансцедеитную регрессию наиболее удобно строить, подбирая такую выравнивающую замену переменных У (г), Х (х), чтобы зависимость 3(Х) была почти линейной (см. гл.
11, 3 1, п. 8). Тогда ее нетрудно аппроксимировать алгебраическим многочленом; Х=Р(Х, а). Выравнивающую замену переменных ищут, используя теоретические соображения и учитывая асимптотику г (х). Дальше будем считать, что такая замена уже сделана. Замечание 2, При переходе к новым переменным задача метода наименьших квадратов (34) принимает вид ~х~ ~))У!(Х! — Р (Хь а))'=ш(п, (35) где новые веса связаны с исходными соотношениями )р'! = ( — „) и!!. (Зб) Нетрудно показать, что всегда выполняется соотношение )р)( 1. Поэтому, даже если в исходной постановке (34) все измерения имели одинаковую точность, так что !в!= — 1, то для выравнивающих переменных веса не будут одинаковыми.
Корреляционный анализ. Надо проверить, действительно ли замена переменных была выравнивающей, т. е. близка ли зависимость г, (Х) к линейной. Это можно сделать, вычислив коэффициент парной корреляции ;э, ег! (х! — х!) (г! — г!) М Я вЂ” М1) (~ — Мй)1 ! = ! (37) ',~„эг; (х! — х!) ~ч; ег! (х! — 2)' !=! != ! стАтистичнскАя ОБРАБОткА экспвпимвнтд [Гл. хч Если зависимость Я(Х) строго линейная (и не содержит слу- чайных ошибок), то р= + 1 или р= — 1 в зависимости от знака наклона прямой. Чем меньше [р[о тем менее зависимость л (Х) похожа на линейную.
Поэтому, если )р) 1, а число измерений [[[ достаточно велико, то выравнивающие переменные выбраны удовлетворительно. Подобные заключения о характере зависимости по коэффици- ентам корреляции называют корреляционным анализом, При корреляционном анализе не требуется, чтобы в каждой точке х! проводилась серия измерений. Достаточно в каждой точке сделать одно измерение, но зато взять побольше точек на исследуемой кривой, что часто делают в физических экспериментах.
Замечание 3. Существуют критерии близости [р ~ к 1, поз- воляющие указать, является ли зависимость А (Х) практически линейной. Мы на них не останавливаемся, поскольку далее будет рассмотрен выбор степени аппроксимирующего многочлена. За,мечание 4. Соотношение [р) — О указывает на отсутствие линейной зависимости Е (Х), но не означает отсутствия какой- либо зависимости. Так,. если к = Х' на отрезке — 1 Х = 1, то р= О. Оптимальная степ ен ь ми огоч лен а.
Подставим в за- дачу (35) аппроксимирующий многочлен, степени тл! т Р(Х) = да а»Х». (38) »-о Тогда оптимальные значения параметров а» удовлетворяют системе линейных уравнений (2.43): »1 ,У', А»[а! — В„О == й = пт, т=е Аы= ~ )(У[Х';+', Д = У )(У[~[Х'!, !=1 !=! и найти их нетрудно. Но как выбрать степень многочлена? Для ответа на этот вопрос вернемся к исходным переменным и вычислим дисперсию аппроксимационцой формулы с найденными коэффициентами.
Несмещенная оценка этой дисперсии такова*): ь) =).)[)(х, а„а„..., а ))в (39) ~~ се[[э! — у(х[, а))з/ ~~ [а!. (40) а=! а= 1 ') В формулах типа (32) имелся делитель вида Лl — [; он связан с тем, что при вычислении стандарта выборки используется одна величина к, определенная по той же выборке. В (40) используется аз+1 козффипиентов а», определенных по выборке, позтому появляется делитель й[-(пт+!). гл.
хтч стлтистичвскля овелаоткл зкспввиментл 499 Очевидно, при увеличении степени многочлена т дисперсия (40) будет убывать: чем больше взято коэффициентов, тем точней можно аппроксимировать экспериментальные точки. Сравним г?ж с дисперсиями з) единичных измерений (32), определенными по небольшим сериям экспериментов хотя бы в нескольких точках хо Если О ->з1 для всех (, то погрешность аппроксимации больше погрешности, с которой измерены значения гь Надо увеличивать т до тех пор, пока отличие В от з) хотя бы для одного 1 не перестанет быть значимым по критерию Фишера (30).
Наоборот, если В (з), то надо уменьшать и. Если полученная таким образом оптимальная степень т, удовлетворяет условию ть((6(, то выравнивающие переменные выбраны удачно; если т, )У, то следует подобрать другую замену переменных. 3 а м е ч а н н е 5, Описанный способ нахождения оптимального числа параметров т, можно применять при произвольном виде функции 1(х, а); но сами коэффициенты в этом случае вычисляются не по формулам (39).
Точность коэффициентов. Коэффициенты аь определяются по случайным величинам г; и поэтому сами являются случайными величинами. Какие их значащие цифры достоверны, а какие можно отбросить? На первую половину вопроса ответить нетрудно. Проведем математический эксперимент. Зная дисперсию единичных измерений з), искусственно внесем в величины г, случайные ошибки бг„ распределенные по нормальному закону с дисперсиями з) (это делается методами Монте-Карло), и вычислим соответствующйе аь.
Повторим эту процедуру многократно. Для каждого аь получим набор случайных значений, по которому вычнслил1 среднее аь и стандарт з(а,). Отсюда по критерию Стьюдента (9) найдем для аь доверительный интервал и, тем самым, выясним, какие значащие цифры коэффициента достоверны, Однако, вообще говоря, недостоверные цифры коэффициента аь нельзя отбрасывать, Коэффициенты аь можно округлять только все одновременно, меняя их на взаимно согласованные величины Лаь. Для такого округления многочлен Р (Х) представляют в виде линейной комбинации: Р (Х) = — Ьи аьХь = г, уД, (Х), (41) ь=ь г=о где (~, (Х) — алгебраические многочлены, ортогональные на системе точек Х, (1 «1«Ж) с весами яг; *). Коэффициенты у, можно округлять независимо друг от друга в пределах их доверительных интервалов.
") См. Приложение. баО СТАТИСТИЧЕСКАЯ ОБРАБОТКА ЭКСПЕРИМЕНТА [ГЛ. ХУ ЗАДАЧИ !. Для примера,<b>Текст обрезан, так как является слишком большим</b>.