Документ взят из кэша поисковой машины. Адрес оригинального документа : http://lnfm1.sai.msu.ru/grav/russian/lecture/mon/5.doc Дата изменения: Mon Nov 4 17:50:24 2002 Дата индексирования: Mon Oct 1 23:33:57 2012 Кодировка: koi8-r Поисковые слова: п п п п п п п р п п р п п р п п р п п р п п р п п р п п р п п р п п р п

Глава 5

ОБРАБОТКА НАБЛЮДЕНИЙ

5.1. Прямые и косвенные наблюдения
Наблюдения называют прямыми в случае, когда результатом этих
наблюдений или измерений является сама искомая величина. В случае, когда
наблюдаемые величины являются исходными данными для вычисления других
параметров, зависящих от этих величин, то к таким наблюдениям применяют
термин косвенные (посредственные) наблюдения.
Однако, между прямыми и косвенными наблюдениями нет четкой границы.
Приведем простой пример. Допустим, что нам нужно знать площадь
прямоугольника. В нашем распоряжении есть линейка с нанесенной на ней
линейной шкалой. Мы можем с помощью этой линейки измерить стороны и
вычислить искомую площадь. Результатом измерений является численное
значение площади. Будут ли измерения прямыми, если для вычисления площади
приходится перемножать данные измерения? Хочется сказать, что нет, так как
в данном случае непосредственно площадь не измеряется, а вычисляется. Если
бы измерения были вычислены с помощью планиметра, то такие измерения были
бы прямыми.
С другой стороны, длина отрезка прямой определяется как разность двух
отсчетов по линейке. Можно ли такую математическую операцию как вычитание
не принимать во внимание и считать измерения сторон прямыми, а перемножение
дает результат косвенных наблюдений? Условность границы между понятиями
прямые и косвенные наблюдения очевидны. Автор известного учебника по
математической обработке астрономических наблюдений Б.М.Щиголев замечает:
"прямое измерение нужных величин заменяется измерением других величин...
Поэтому говорят иногда, что производятся косвенные или посредственные
измерения; эту терминологию, однако, нельзя считать удачной." [Б.М.Щиголев
"Математическая обработка наблюдений", М.,1969]
По-видимому, прямыми измерениями можно назвать такие измерения,
которые можно неоднократно повторить, и результаты которых могут быть
использованы при уравнительных вычислениях (в методе наименьших квадратов)
других величин, функционально связанных с результати прямых измерений.
Повторение измерений необходимо для того, чтобы образовать ряд наблюдений
для применения статистических методов обработки. Например, для вычисления
элементов орбит планет или астероидов требуются наблюдения их положений в
виде рядов ( и ( - координат на небесной сфере. По данным этих наблюдений,
которых должно быть заведомо больше определяемых элементов орбит,
организуют уравнительные вычисления поправок к элементам, в результате
получим посредственные наблюдения элементов. Однако, для определения ( и (
необходимо измерить линейные расстояния между изображениями выбранного нами
небесного тела и опорных звезд. Только после выполнения уравнительных
вычислений, снижающих влияние случайных ошибок измерений, коррекции за
неточность установки астронегатива в измерительном инструменте и влияние
оптических искажений получают искомые координаты. следовательно, в задаче
определения ( и ( (наблюдения астероидов, комет и т.п.) прямые измерения
дают лишь взаимые расстояния между изображениями небесного тела и опорных
звезд, а косвенные измерения дают ( и (. В другой задаче - задаче
определения орбит - прямыми наблюдениями можно считать координаты ( и (, а
косвенными - элементы орбиты.
Некоторая условность в разделении этих двух типов наблюдений
(измерений) никак не мешает применению статистических методов оценивания
искомых величин, т.к. наблюдения, которые называют косвенными или
посредственными, объединяют оба вида измерений, а прямые измерения можно
считать частным случаем косвенных.

5.2. Статистические методы оценивания
Оцениванием называют процедуру определения оценки - приближенной
величины - искомого параметра. Для оценивания, например, параметра ( по
ряду наблюдений [pic] необходимо выполнить вычисления, используя эти
наблюдательные данные. Запишем эту процедуру следующим образом:
[pic].
Функцию [pic] называют статистикой, крышка над (-обозначением параметра
будет обозначать оценку.
Оценку называют несмещенной, если она не содержит систематической
ошибки, т.е. [pic], где треугольные скобки, как и прежде, обозначают
осреднение по "бесконечному ансамблю реализаций".
Как правило, несмещенная оценка не единственная, их может быть
бесконечное множество. Возникает задача выбора наилучшей из этих оценок.
Кажется естественным отобрать ту из них, которая имеет минимальную
дисперсию погрешности
[pic].
Вычисление этого условия с учетом условия несмещенности часто приводит
к однозначности решения.
5.2.1. Линейная оценка параметра
Рассмотрим следующую задачу. Требуется определить наилучшую в средне-
квадратическом смысле оценку параметра Х, значения которого заданы рядом
прямых наблюдений [pic]. Линейной оценкой параметр Х будет равен линейной
статистике
[pic]
Эта оценка будет несмещенной, если [pic], другими слова, если
[pic].
Поскольку [pic] - точные константы, их можно вынести за скобки осреднения
[pic].
Примем гипотезу о случайности погрешностей измерений. Если [pic] содержат
лишь случайную ошибку, то средние их значения равны точному, т.е.
[pic].
Теперь условие несмещенности принимает вид
[pic], или при Х(0 [pic].
Это и есть условие несмещенности линейной оценки. Таким образом, имеем одно
уравнение для определения n неизвестных. Число решений - бесконечно, вернее
[pic], где (n-1) - число степеней свободы выбора коэффициентов [pic].
Действительно, если (n-1) коэффициентам [pic] задать любые произвольные
значения, то последний коэффициент определится однозначно
[pic].
Воспользуемся теперь вторым условием - минимума дисперсии погрешности
оценки. Рассмотрим два случая.
1) Наблюдения независимы и равноточны. В этом случае дисперсии ошибок
всех наблюдений равны
[pic].
Для независимых наблюдений дисперсия суммы равна сумме дисперсий. Поэтому
[pic].
Вследствие равноточности имеем [pic] .
Теперь нам нужно найти минимум функции [pic] при условии, что [pic].
Поступим следующим образом. Определим [pic] через остальные коэффициенты и
подставим в выражение для дисперсии [pic]:
[pic],
[pic] .
Теперь задачу можно решить как задачу поиска абсолютного минимума:[pic].
Отсюда [pic],
Итак, все коэффициенты равны одному и тому же числу. Легко понять, что это
число 1/n. Таким образом,
[pic].
Вывод: арифметическое среднее есть несмещенная и наилучшая в средне-
квадратическом смысле оценка величины, полученная по равноточным и
независимым наблюдениям.
2) Наблюдения независимы и неравноточны. Теперь будем считать, что
каждое наблюдение [pic] содержит погрешности, дисперсии которых не
обязательно равны. Пусть [pic] соответственно дисперсии погрешностей [pic].
Теперь дисперсия линейной оценки будет равна
[pic]
Определим минимум [pic] при условии, что [pic]. Поступим по аналогии с
предыдущими рассуждениями:
[pic],
[pic][pic],
[pic].
Поскольку [pic], то получим уравнение [pic]. Полученное равенство говорит о
том, что все произведения вида [pic] равны между собой и не зависят от
индекса к. Кроме того, сумма [pic] должна равняться единице.
Пусть [pic], где [pic] - некоторый нормировочный множитель, а [pic] -
константа, имеющая размерность дисперсии. Отсюда
[pic],
[pic].
Обозначим [pic], [pic]. Теперь определим [pic]:[pic]. Итак,
[pic].
Наилучшая линейная оценка в этом случае равна [pic], где [pic],
[pic].
Величины [pic] носят название веса наблюдения, ибо формула для определения
[pic] аналогична той, которую применяют для вычисления центра тяжести, если
"вес" точки с координатой [pic] будет равен [pic]. Постоянная [pic] может
быть произвольной, так как ее выбор не влияет на оценку [pic]. Ее можно
также интерпретировать как дисперсию наблюдения, вес которому мы
приписываем равным единице ("дисперсия единицы веса").

5.2.2. Дисперсия ошибок линейных оценок
Получим теперь формулы для дисперсии ошибок линейных оценок. Снова
возьмем два варианта: будем считать, что все наблюдения независимы, но они
либо равноточны (1), либо неравноточны (2).
1. Поскольку арифметическое среднее - оценка равноточных наблюдений -
является несмещенной, то погрешность арифметического среднего равна
арифметическому среднему погрешности измерений.
[pic],
а дисперсия ошибки среднего арифметического равна
[pic].
Таким образом, с увеличением числа наблюдений искомого параметра точность
оценки увеличивается, дисперсия ошибки арифметического среднего убывает
обратно пропорционально n.
Стандартное отклонение ошибки арифметического среднего [pic] убывает
как квадратный корень из n:
[pic].
Отсюда вывод: для того, чтобы увеличить точность оценки на один порядок (в
10 раз), нужно увеличить объем данных на два порядка (в 100 раз) и т.д.
Возникает естественно вопрос, можно ли таким способом увеличивать точность
беспредельно? Математика на этот вопрос дает утвердительный ответ. Да,
можно измерять отрезки прямой с точностью до одного микрометра линейкой с
миллиметровыми делениями. Для этого нужно выполнить 1 000 000 измерений,
при этом измерение должно содержать только строго случайные погрешности.
Последнее условие, конечно, не выполняется, так как полное отсутствие
систематических погрешностей на практике обеспечить нельзя.
2. Для оценки параметра Х при неравноточных измерениях мы имеем
формулу
[pic], где [pic].
Обозначая для краткости [pic], получим ошибку
[pic].
Все ошибки [pic] по условию независимы, поэтому
[pic].
Но [pic], отсюда [pic].
Итак, дисперсия ошибки среднего весового равна "дисперсии единицы веса",
поделенная на сумму весов наблюдений.
[pic] .
Учитывая, что дисперсия ошибки каждого наблюдения равна дисперсии единицы
веса, деленной на вес этого наблюдения (см. 5.2.1), [pic], делаем вывод,
что вес среднего весового равен сумме весов наблюдений.

5.2.3. Оценка дисперсии единицы веса по данным наблюдений. Средняя
квадратическая ошибка.
В разделе 5.2.1 мы показали, что произвольно взятая константа [pic]
есть не что иное, как дисперсия ошибки наблюдения, которому мы приписываем
вес единица. Произвольность этой константы определяется тем, что мы
произвольно приписываем вес единица какому-либо реальному или воображаемому
наблюдению. В случае равноточных наблюдений обычно вес, равный единице,
дают каждому наблюдению. В этом случае [pic] не может быть произвольной
величиной. Эта величина не может быть произвольной величиной и в том
случае, когда вес наблюдения известен и выбран обратно пропорционально
дисперсии ошибки.
Пусть [pic] - исходные данные измерений, [pic] - их веса. Дисперсии
ошибок отдельных наблюдений будем считать неизвестными. По определению
весов имеем
[pic].
Очевидно, справедливо равенство [pic].
Следовательно, [pic].
Таким образом, дисперсия единицы веса равна среднему значению
арифметического среднего произведений весов на квадраты погрешностей.
Операция, обозначенная треугольными скобками, требует бесконечного
повторения набора наблюдений [pic] и является лишь теоретической
процедурой. Поэтому в качестве оценки дисперсии единицы веса можно принять
величину [pic], полученную из единственного ряда наблюдений
[pic].
Как следует из приведенного выше равенства для [pic], эта оценка
несмещенная
[pic].
Для вычисления оценки дисперсии по полученной формуле требуется знать
точное значение параметра Х. На практике чаще всего это значение нам
неизвестно. Поэтому заменим Х на его оценку
[pic]
и проверим полученное равенство на несмещенность [pic].
Здесь и далее суммирование ведется по к от 1 до n. Для краткости пределы
суммирования мы будем опускать.
[pic].
Остается определить ковариацию: [pic].
Вследствие независимости измерений ковариации ошибок с разными индексами
равны нулю, а с одинаковыми - дисперсии:
[pic].
Учитывая, что [pic], получим [pic].
Таким образом, замена точного значения параметра Х на приближенное [pic]
при вычислении оценки дисперсии приводит к тому, что эта оценка оказалась
смещенной. Умножая ее на [pic], получим несмещенную оценку дисперсии
единицы веса
[pic].
Число в знаменателе (n-1) указывает на то, что наша статистика, принятая
для вычисления оценки дисперсии, [pic] имеет (n-1) степеней свободы. Одна
степень свободы из n потрачена на то, чтобы образовать оценку [pic].
При равноточных наблюдениях все дисперсии ошибок, а следовательно, и
веса наблюдений равны. Пусть [pic], тогда [pic]. Получим
[pic],
[pic].
Оценку стандартного отклонения для ошибок называют средней квадратической
ошибкой (СКО). Для этой величины, полученной a posteriory, т.е. на
основании опыта, мы будем применять обозначение (. Таким образом,
[pic] для равноточных наблюдений,
[pic] для неравноточных наблюдений.
В случае, когда величина Х известна, то в знаменателях приведенных формул
(n-1) нужно заменить на n.

5.2.4. О выборе весов
Правильный выбор весов наблюдений гарантирует корректность алгоритма
определения "наилучшего" значения параметра и его СКО. Однако, величины
дисперсий погрешностей каждого отдельного наблюдения, как правило,
неизвестны. Тем не менее, в отдельных случаях веса наблюдениям можно
назначить совершенно строго.
Например, пусть [pic] получено как арифметическое среднее из [pic]
равноточных измерений, [pic] - из [pic] таких же измерений и т.д. Тогда,
согласно изложенному ранее дисперсии ошибок [pic] будут соответственно в
[pic] раз меньше, чем дисперсия одиночного измерения. Если функция ошибок
одиночного измерения равна [pic], то
[pic].
Легко видеть, что в качестве весов можно взять число рядовых
измерений, из которых получили отдельное значение [pic]. поскольку веса
можно брать с точностью до постоянного множителя, то
p1 ( n1, p2 ( n2, .....
Второй способ - менее надежный. Он основан на замене дисперсии ошибки
n-го наблюдения [pic] на квадрат его средней квадратической ошибки [pic].
Недостаток этого метода заключается в том, что квадраты СКО могут
значительно отличаться от дисперсий, так как СКО вычисляют, как правило, на
основании внутренней сходимости данных при их сравнительно небольшом
количестве. Тем не менее, когда ничего другого не остается, полагают
p1 ( (1-2, p2 ( (2-2, .....
У астрономов существуют и субъективные методы назначения весов.
Например, если какое-то наблюдение всем не понравилось (этому должны быть
объективные причины!), вы можете назначить вес 1/2, когда как остальным
наблюдениям - вес 1. Возможна также субъективная оценка качества наблюдений
по 10-бальной шкале, что также можно использовать при назначении весов.
Приведем пример комплексного подхода при выборе весов. Задача состояла
в том, чтобы определить собственное движение звезд относительно
внегалактических туманностей. Собственное движение определяется как
изменение координат звезды за промежуток времени между двумя эпохами
наблюдений. В одной из работ, посвященных этой задаче, веса были выбраны
следующим образом:
pi ( [pic],
где [pic] - разность эпох, m - число внегалактических туманностей, [pic]-
качество изображения туманности, оцененное наблюдателем по 10-бальной
шкале. Здесь множитель [pic] указывает на то, что дисперсия ошибки
собственного движения обратно пропорциональна квадрату разности эпох. В
случае одинакового качества изображений, [pic], сумма [pic]пропорциональна
m - числу пар координат, из которых вычитается собственное движение как
среднее арифметическое. Однако, при [pic] больший вес придается таким парам
измерений, чье изображение туманности более звездообразное.

5.2.5. Обработка рядов равноточных и неравноточных наблюдений
Проиллюстрируем применение полученных формул на примерах обработки
рядов наблюдений.
Пусть наша задача состоит в получении наиболее точного значения
постоянной аберрации ( по 10 наблюдениям.
| ( | (| х| p | vI| vII|[pic], |
| | | | | | |[pic], |
| | | | | | |[pic], |
| | | | | | |[pic], |
| | | | | | |[pic], |
| | | | | | |[pic], |
| | | | | | |[pic]. |
| | ( | 45| 0.07| | -36| |
|20''.445|0.011 | | |-37 | | |
| | | 58| 0.14| | -23| |
|.458 |.008 | | |-24 | | |
| | | 98| 0.06| | +17| |
|.498 |.012 | | |+16 | | |
| | | 95| 0.05| | +14| |
|.495 |.013 | | |+13 | | |
| | | 83| 1.00| | +2 | |
|.483 |.003 | | |+01 | | |
| | | 112| 0.11| | +31| |
|.512 |.009 | | |+30 | | |
| | | 95| 0.05| | +14| |
|.495 |.013 | | |+13 | | |
| | | 91| 0.11| | +12| |
|.491 |.009 | | |+09 | | |
| | | 90| 0.07| | +9 | |
|.490 |.011 | | |+08 | | |
| | | 50| 0.11| | -31| |
|.450 |.009 | | |-32 | | |
| | |( 817 | | | | |
| | | | |-03 | | |
| | |[pic] | 1.77| | | |
| | |82 | |00 | | |


В первой колонке приведены значения ( в угловых секундах, во второй - их
ошибки по оценке наблюдателей. Поскольку все значения отличаются только в
двух последних знаках (кроме одной строки), именно эти два знака и возьмем
в качестве величины х, подлежащей обработке (третья колонка).
Обработаем двумя способами. Сначала пренебрежем сведениями об их
точности и будем считать эти данные равноточными. Тогда "наилучшей" в
среднеквадратическом смысле оценкой будет арифметическое среднее
[pic].
Под третьей колонкой приведено значение суммы и среднего с точностью до
целого. Для вычисления СКО вычислим разности [pic] (пятая колонка).
По формуле [pic] вычисляем ( - ошибку измерения х, а по формуле
[pic] - СКО среднего арифметического. В данном случае [pic].
Вспомним, что величина х - две последние цифры в тысячных долях секунды,
получим окончательный ответ: [pic].
Средняя квадратическая ошибка приписывается рядом с оценкой параметра
в виде [pic], указывая на интервал, внутри которого, по-видимому, находится
истинное значение.
Заметим, что ошибка единичного измерения в предположении, что все
наблюдения равноточны, у нас оказалась равной [pic], что куда больше, чем
все величины ( во втором столбце таблицы, полученные по внутренней
сходимости.
Теперь обработаем снова приведенный ряд наблюдений, принимая величины
( как реальные стандартные отклонения погрешностей. Иначе говоря, каждому
значению х мы припишем вес, равный [pic]. В качестве стандартного
отклонения [pic] возьмем наименьшую из приведенных СКО, т.е. припишем вес
единицу наиболее точному определению аберрации. Для весов достаточно
оставлять одну-две цифры (четвертая колонка). Сумма весов записана под
колонкой р=1.77. Вычислим среднее весовое
[pic].
Новые остаточные разности [pic] приведены в шестой колонке. Вычисляем СКО
единицы веса
[pic].
Наконец, вычисляем СКО весового среднего: [pic].
Таким образом, ответ будет выглядеть следующим образом: [pic]. Сравнивая
полученный результат с предыдущим, можно отметить, что они практически не
различаются.
В качестве априорного стандартного отклонения мы приняли [pic] (в
угловых секундах), а по экспериментальным данным мы получили оценку этой
величины ([pic]) равной [pic]. Столь большая разница говорит лишь о том,
что точность каждого определения аберрации сильно завышена. Любопытно
заметить, что наблюдение, оцененное средней квадратической ошибкой [pic]
мало отличается от среднего арифметического и среднего весового. Это
говорит о том, что действительно это наблюдение наиболее точное.
Приведем еще один пример, иллюстрирующий, как алгоритм обработки
неравноточных наблюдений автоматически отбраковывает "плохие" определения.
Определение гравитационной постоянной.
| Год | Автор |G(10-11 | p |
| | |((10-11 | |
| 1798 |Кавендиш |6.75 | 0.00|
| | |(0.05 | |
| 1852 |Райх |6.64 | 0.00|
| | |(0.06 | |
| 1891 |Пойнтинг |6.70 | 0.00|
| | |(0.04 | |
| 1896 |Этвеш |6.657 | 0.00|
| | |(0.013 | |
| 1942 |Хейл |6.673 | 0.03|
| | |(0.005 | |
| 1972 |Фаси |6.6714 | 0.69|
| | |(0.0006 | |
| 1978 |Сагитов |6.6745 | 0.39|
| | |(0.0008 | |
| 1982 |Лютер |6.6726 | 1.00|
| | |(0.0005 | |
| | | |p=2.11|


В приведенной таблице в первой колонке - год определения
гравитационной постоянной, во второй - автор этого определения. Третья и
четвертая колонки - результат определения постоянной G и ее средняя
квадратическая ошибка. В качестве стандартного отклонения ошибки единицы
веса примем СКО наиболее точного определения: [pic]. Тогда веса остальных
наблюдений определим по формуле [pic]. Соответствующие значения [pic]
приведены в пятой колонке. Веса, также как и средние квадратические ошибки,
выписывают с одной-двумя значащими цифрами. Однако, поскольку наблюдения
однородны (все наблюдения содержат только G), то необходимо сохранить и
число знаков после запятой. В результате получим, что веса наблюдений
прошлого века равны нулю, т.е. в обработку не включаются. Вычисление
среднего весового дает G=6.6726(10-11. СКО среднего весового можно получить
только по априорной оценке точности:
[pic].
Таким образом, G=(6.6726(0.0003)(10-11.
Разность значений по отдельным определениям относительно среднего весового
равны -0.0012 (Фаси), +0.0019 (Сагитов), 0.0000 (Лютер). Эти разности
превосходят по модулю их средние квадратические ошибки (кроме последнего),
что еще раз указывает на недостаточность оценки точности по внутренней
сходимости.