Поэтому в своем методе мы анализировали лишь ШИРОТЫ звезд «Альмагеста». Заранее было неясно — достаточно ли широт для датировки. Оказалось, что ответ положительный. Мы утверждаем, что «Альмагест» можно датировать, используя лишь сведения о ШИРОТНЫХ невязках.
Затем, работоспособность нашего метода была подтверждена анализом звездных каталогов Т. Браге, Улугбека, Гевелия и ряда, искусственно созданных нами каталогов, для чего использовался компьютер. Во всех случаях полученные нашим методом датировки каталогов совпали с заранее известными.
Предварительная работа по выявлению выбросов в «Альмагесте» во многом была уже проделана в более ранних исследованиях. См., например, [328]. Мы считали выбросами те звезды, у которых значение широтных невязок превосходило 1 градус. Кроме явных выбросов каталог содержит звезды, отождествление которых со звездами современного неба сомнительно. В упомянутой работе К. Петерса и Е. Кнобеля [328] такие случаи также отмечены. Один пример уже был приведен выше: это — звезда O2 Эридана. Поэтому для исключения всех таких сомнительных случаев необходимо было очистить каталог «Альмагеста» от неоднозначно отождествляемых звезд. Мы проверили список из более чем 80 быстрых звезд из современного каталога [281]. Из них в «Альмагесте», как выяснилось, отражено около 35 звезд. Затем мы выявили среди них звезды, имеющие неоднозначное, сомнительное отождествление. Таких звезд оказалось немного — всего три. Они были исключены из рассмотрения. Таким образом, наш анализ в основном подтвердил правильность отождествления подавляющего большинства звезд «Альмагеста», приведенного в труде [328].
Перейдем к анализу систематических ошибок. Если рассмотреть какую-нибудь совокупность звезд, то систематическая ошибка в положении этих звезд на небесной сфере может состоять только лишь в перемещении совокупности звезд как единого целого по небесной сфере. Такое перемещение имеет три степени свободы и, следовательно, может быть описано тремя параметрами. Однако поскольку мы интересуемся лишь широтными невязками, то достаточно рассмотреть только двухпараметрические вращения сферы. С вычислительной точки зрения удобно задать это вращение с помощью параметров φ и γ, где параметр φ задает ось, вокруг которой вращается сфера, а параметр γ задает угол поворота. См. рис. 1.13. А именно, мы выбираем в качестве φ угол между осью весеннего равноденствия, рассчитанной на какой-либо год t, и осью поворота, лежащей в плоскости эклиптики, также относящейся к году t.
Рис. 1.13. Вращение небесной сферы можно задать при помощи двух углов
Итак, если предположить, что звездный каталог составлялся в год t и истинные широта и долгота какой-либо звезды были равны b(t) и l(t) соответственно, то в результате ошибки в определении положения эклиптики, парамеризуемой γ=γ(t) и φ=φ(t), составитель каталога запишет в каталог координаты b'(t) и l'(t). С очень большой точностью можно считать, что
b'(t)=b(t)+γ sin(l(t)+φ)
Последняя формула справедлива при условии, что составитель каталога не делал никакой ошибки измерений. Если ошибка присутствовала, — а она присутствовала неизбежно, — и равнялась λ, то
b'(t)=b(t)+γ sin(l(t)+φ)+ λ
Последняя формула справедлива для всех звезд рассматриваемой совокупности, и, следовательно, можно поставить статистическую проблему оценки параметров γ и φ для данной совокупности звезд. Оценки параметров γ и φ можно найти, например, методом наименьших квадратов, когда γ и φ являются решением следующей задачи: