Контрольная подстановка полученных значений z в исходные уравнения показала правильность их вычислений (с учетом необходимых округлений промежуточных результатов).

§ 152. Коррелатный способ уравнивания

Приведенная выше система уравнений (16.8) имеет нелинейный вид. В математике не существует способов решения таких систем нелинейных уравнений. В связи с этим данную систему уравнений раскладывают в ряд Тейлора, ограничиваясь только первыми членами разложения с учетом того, что значения поправок vi достаточно малы (на основании выдержанных при измерениях допусков по точности), и вторые их степени будут весьма малыми, так что ими можно будет пренебречь. В результате уравнения (16.8) преобразуются к виду:

где i – номер измеренной величины (х); j – номер условного уравнения (или функции φ).

[arv] + Wr = 0

Равенства (16.82) и (16.83) называются условными уравнениями поправок.

Следует иметь в виду, что формулы (16.81) не используются, если известно, что система уравнений (16.8) имеет линейный вид, т.е. коэффициенты aij известны.

Для решения задачи уравнивания способом Лагранжа необходимо составить следующую функцию:

409

(16.85) Для определения поправок vi , при которых функция (16.85) достигает минимума, найдем частные производные по аргументам vi и приравняем их

нулю:

Уравнения (16.87) и (16.88) называют коррелатными уравнениями поправок.

Число коррелат всегда равно числу условий. В результате образуется система коррелатных уравнений поправок vi, содержащая n неизвестных поправок vi и r неизвестных коррелат kj, состоящая из n линейных уравнений.

Опуская промежуточные, хотя и важные преобразования (вывод можно посмотреть в [3, 8, 9, 40, и др.]), основанные на методе наименьших квадратов, приведем т.н. нормальные уравнения коррелат, в которых число неизвестных равно числу уравнений:

В уравнениях (16.89) неизвестными являются коррелаты ki , а свободными членами – свободные члены уравнений поправок (16.82) и (16.83).

Как видно, при каждом параметре (неизвестном) ki в уравнениях (16.89) стоит коэффициент в виде гауссовой суммы. Представим указанные коэффициенты в развернутом виде с помощью таблицы коэффициентов уравнений поправок (табл. 16.5).

410

Развернутый вид коэффициентов [qa1a1 ], [qa1a2 ],..., [qa r ar ] , в которых индекс при коэффициентах а – это второй индекс коэффициентов условных уравнений поправок:

………………………………………………

[qa r ar ] = brr = q1a1r a1r + q2 a2r a2r + q3 a3r a3r + ... + qn anr anr

Рассмотрим подробнее принцип вычисления коэффициентов bjj при коррелатах kj в нормальных уравнениях коррелат.

1-е уравнение коррелат.

Коэффициент при k1 равен сумме произведений обратных весов и квадратов коэффициентов первой строки матрицы.

Коэффициент при k2 равен сумме произведений обратных весов и коэффициентов 1-й и 2-й строк матрицы.

Коэффициент при k3 равен сумме произведений обратных весов и коэффициентов 1-й и 3-й строк матрицы и т.д.

………………………………………………………………………………….

Коэффициент при kj равен сумме произведений обратных весов и коэффициентов 1-й и j-й строк матрицы и т.д.

………………………………………………………………………………….

Коэффициент при kr равен сумме произведений обратных весов и коэффициентов 1-й и r-й строк матрицы.

2-е уравнение коррелат.

Коэффициент при k1 равен сумме произведений обратных весов и коэффициентов 2-й и 1-й строк матрицы.

411

Коэффициент при k2 равен сумме произведений обратных весов и квадратов коэффициентов 2-й строки матрицы.

Коэффициент при k3 равен сумме произведений обратных весов и коэффициентов 2-й и 3-й строк матрицы и т.д.

………………………………………………………………………………….

Коэффициент при kj равен сумме произведений обратных весов и коэффициентов 2-й и j-й строк матрицы и т.д.

………………………………………………………………………………….

Коэффициент при kr равен сумме произведений обратных весов и коэффициентов 2-й и r-й строк матрицы.

3-е уравнение коррелат.

Коэффициент при k1 равен сумме произведений обратных весов и коэффициентов 3-й и 1-й строк матрицы.

Коэффициент при k2 равен сумме произведений обратных весов и коэффициентов 3-й и 2-й строк матрицы.

Коэффициент при k3 равен сумме произведений обратных весов и квадратов коэффициентов 3-й строки матрицы и т.д.

………………………………………………………………………………….

Коэффициент при kj равен сумме произведений обратных весов и коэффициентов 3-й и j-й строк матрицы и т.д.

………………………………………………………………………………….

Коэффициент при kr равен сумме произведений обратных весов и коэффициентов 3-й и r-й строк матрицы.

…………………………………………………………………………………

j-е уравнение коррелат.

Коэффициент при k1 равен сумме произведений обратных весов и коэффициентов j-й и 1-й строк матрицы.

Коэффициент при k2 равен сумме произведений обратных весов и коэффициентов j -й и 2-й строк матрицы.

Коэффициент при k3 равен сумме произведений обратных весов и квадратов коэффициентов j-й и 3-й строк матрицы и т.д.

………………………………………………………………………………….

Коэффициент при kj равен сумме произведений обратных весов и квадратов коэффициентов j-й строки матрицы.

Коэффициент при kj+1 равен сумме произведений обратных весов и коэффициентов j –й и ( j+1)-й строк матрицы и т.д.

………………………………………………………………………………….

Коэффициент при kr равен сумме произведений обратных весов и коэффициентов j-й и r-й строк матрицы.

r-е уравнение коррелат.

Коэффициент при k1 равен сумме произведений обратных весов и коэффициентов r -й и 1-й строк матрицы.

Коэффициент при k2 равен сумме произведений обратных весов и коэффициентов r -й и 2-й строк матрицы.

412

Коэффициент при k3 равен сумме произведений обратных весов и коэффициентов r –й и 3-й строк матрицы и т.д.

………………………………………………………………………………….

Коэффициент при kj равен сумме произведений обратных весов и коэффициентов r-й и j-й строк матрицы и т.д.

………………………………………………………………………………….

Коэффициент при kr равен сумме произведений обратных весов и квадратов коэффициентов r-й строки матрицы.

Таким образом получают уравнения коррелат вида:

b11k1 + b12k2 + b13k3 + …+ b1jkj + …+ b1rkr + W1 = 0

b21k1 + b22k2 + b23k3 + …+ b2jkj + …+ b2rkr + W2 = 0

b31k1 + b32k2 + b33k3 + …+ b3jkj + …+ b3rkr + W3 = 0

………………………………………………………. (16.91)

bj1k1 + bj2k2 + bj3k3 + …+ bjjkj + …+ bjrkr + Wj = 0

……………………………………………………….

br1k1 + br2k2 + br3k3 + …+ brjkj + …+ brrkr + Wr = 0

Можно заметить, что коэффициенты b с обратными индексами равны между собой, т.е. b12 = b21, b35 = b53 и т.п. Т.н. диагональные коэффициенты bjj представляют собой сумму произведений обратных весов и квадратов коэффициентов а j–й строки, т.е. они всегда положительные. Коэффициенты b с обратными индексами располагаются с разных сторон от диагональной строки. В связи с этим достаточно вычислить диагональные коэффициенты и все коэффициенты, стоящие справа от диагонали. А далее дополнить уравнения недостающими коэффициентами b, записав их такими же, как и коэффициенты с обратными им индексами.

Решение систем линейных уравнений (16.91) выполняется различными способами, рассмотренными в § 151, но все они, как можно было убедиться из приведенных примеров, весьма громоздкие и требуют значительных затрат времени.

Полученные из решения уравнений (16.91) коррелаты kj используются для вычисления поправок vi по формулам (16.87) или (16.88). После введения поправок в измеренные величины получают уравненные значения измеренных величин (16.7).

При оперировании численными значениями коэффициентов условных уравнений, коррелат, весов (обратных весов) и т.п. необходимо иметь ввиду следующее:

-значения весов и обратных им величин вычислять до 0,01-0,001 единиц;

-значения коэффициентов a, b и коррелат k вычислять до 0,001-0,0001 единиц;

-чаще всего невязки W при обработке плановых построений выражают в дециметрах, в высотных сетях – в миллиметрах, угловые невязки и поправки выражают в секундах, десятых и сотых долях секунды.

Суммируя сказанное выше, приведем последовательность решения задачи уравнивания коррелатным способом.

413

Шаг 1. Для данного геодезического построения в системе n результатов xi , имеющих веса pi, определяют число k независимых и число r избыточных измерений.

Шаг 2. Составляют математические соотношения (условные уравнения) вида (16.5) с учетом следующих основных требований:

-все условные уравнения должны быть независимыми, т.е. ни одно из них не должно быть следствием другого (других);

-число уравнений должно быть равно числу избыточных измерений r;

-условные уравнения должны иметь возможно простой вид.

Шаг 3. Условные уравнения приводят к линейному виду, для чего выполняют их дифференцирование и находят коэффициенты aij (16.81) как частные производные функций φj по аргументам xi .

Находят свободные члены Wj уравнений, т.е. невязки в полученных уравнениях после подстановки в них измеренных значений xi.

Составляют таблицу (матрицу) коэффициентов aij и обратных весов qi (табл. 16.5).

Шаг 4. Находят коэффициенты bjj (16.90) нормальных уравнений коррелат (16.91) по алгоритму, изложенному выше, и решают полученную систему линейных уравнений.

После получения значений коррелат kj из решения уравнений обязательно необходимо выполнить контроль вычислений. Для этого значения коррелат подставляют в исходные уравнения (16.91) и проверяют выполнение условия. Незначительные отступления от условия уравнения допускаются, они возникают из-за округления результатов вычислений.

Шаг 5. Составляют, пользуясь табл. 16.5, условные уравнения поправок νi (16.82), (16.87) ( или (16.88). Для значений поправок, например, получим:

ν1 = q1(a11k1 + a12k2 + …+ a1jkj +…+ a1rkr) ν2 = q2(a21k1 + a22k2 + …+ a2jkj +…+ a2rkr) ν3 = q3(a31k1 + a32k2 + …+ a3jkj +…+ a3rkr)

Вычисляют поправки к измеренным величинам.

После вычисления поправок необходимо выполнить контроль вычислений. Для этого значения поправок следует подставить в условные уравнения поправок (16.82) и проверить выполнение указанного условия. Незначительные отклонения от указанного условия допускаются, они возникают изза округления результатов вычислений.

Шаг 6. Вычисляют уравненные значения xi' (16.7).

Контроль уравнивания осуществляют подстановкой xi' в условные уравнения (16.9). При правильном решении задачи все условные уравнения должны иметь указанное решение. Допускаются незначительные отклонения от указанных условий, они возникают из-за округления результатов вычислений.

После выполнения контроля значения xi' округляют с необходимой точностью и вычисляют искомые величины (координаты, высоты и т.п.).

414