В Харисове по Глонасс и некоторых других источниках приводится метод итерационного решения с использованием МНК
Все правильно. Записывая уравнения связывающее координаты потребителя, координаты спутника и измерения псевдодальности, для каждого спутника, вы получаете довольно сложную
нелинейную систему уравнений. Решать такую систему проще всего итерационным методом. Для этого нелинейную правую часть раскладывают в ряд Тейлора в окрестности некоторой точки и удерживают только линейные слагаемые. Получают таким образом систему
линейных уравнений, которую решают методом наименьших квадратов.
Полученное таким образом решение является приближенным, поскольку является решением линеаризованной системы, а не исходной. Его используют в качестве начальной точки для следующей итерации. Линеаризуют систему в окрестности этой точки и т.д. по циклу до тех пор, пока не будет достигнута нужная точность решения.
Метод наименьших квадратов,
основная идея. У вас есть избыточная система линейных уравнений вида
Fx=
H, где
F - матрица системы,
x - вектор неизвестных. Интересует приближенное решение системы, т.е.
Fx-
H=
g такое, что норма вектора ||
g|| = sqrt(sum(g
i2)) - минимальна.
П - вектор невязок. Пишете выражение для нормы вектора
G: ||
G||
2 = (
Fx-H)
T (
Fx-H). Получили квадратичную форму. Решением, минимизирующим невязки исходной системы будет решение, минимизирующее эту квадратичную форму. Ищете ее экстремум.
Решение, если я правильно помню, получается следующим:
x=(
FTF)
-1FTHЗабудте вы про баллистику! Читайте правильные книги - про математический метод следует читать в учебнике математики, про алгоритм работы приемника спутниковой радионавигации - книги по радиотехнике!