 |
Астронет: А. С. Расторгуев/ГАИШ Определение кривой вращения и шкалы расстояний в Галактике http://variable-stars.ru/db/msg/1172553/node13.html |
<< 11. Вычисление ошибок параметров | Оглавление | Литература >>
12. Метод решения
Максимально-правдоподобное решение для неизвестных параметров задачи (20)-(22) сравнительно слабо чувствительно к случайным ошибкам наблюдений, которые автоматически учитываются с помощью тензора ковариации (15). Тем не менее, оно также требует внимательного отношения к объектам, которые при малых случайных ошибках наблюдений сильно уклоняются от общего решения. Их не рекомендуется использовать для окончательного вычисления функции правдоподобия.
Для минимизации нелинейной функции (22) существует целый ряд эффективных алгоритмов, опубликованных, в частности, в монографии [19]. Это методы Гаусса - Ньютона и Левенберга - Маркардта (стр. 523), использующие интерполяцию и/или градиентный поиск, и метод деформируемых многогранников Нелдера - Мида (так наз. симплекс-алгоритм, стр. 289), более медленный, но не нуждающийся в вычислении производных от целевой функции. На их базе созданы готовые эффективные алгоритмы, написанные на языках FORTRAN, Pascal, C [19], а также включенные в пакеты MATLAB и IDL.
<< 11. Вычисление ошибок параметров | Оглавление | Литература >>