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

 

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

(14.78)

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

Введем обозначения:

(i = 1, 2, …, n; j = 1, 2, …, r), (14.79)

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

С учетом введенных обозначений получим:

а₁₁v₁ + a₂₁v₂ + … + a_n1v_n + W₁ = 0

а₁₂v₁ + a₂₂v₂ + … + a_n2v_n + W₂ = 0 (14.80)

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

а_1rv₁ + a_2rv₂ + … + a_nrv_n + W_r = 0

В обозначениях гауссовых сумм:

[ a₁v ] + W₁ = 0

[ a₂v ] + W₂ = 0 (14.81)

………………

[ a_rv ] + W_r = 0

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

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

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

, (14.82)

где - неопределенные множители Лагранжа.

Обозначим: , где - коррелаты. Тогда функцию (14.82) можно записать со значениями коррелат:

. (14.83)

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

(14.84)

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

Из полученной системы уравнений следует, что

(i= 1, 2, …, n) (14.85)

или

, (14.86)

где - обратный вес измерения с индексом i.

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

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

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

(14.87)

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

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

В уравнениях (14.87) неизвестными являются коррелаты k_i, а свободными членами – свободные члены уравнений поправок (14.80) и (14.81).

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

Развернутый вид коэффициентов , в которых индекс при коэффициентах а – это второй индекс коэффициентов условных уравнений поправок:

Таблица 14.3

Матрица коэффициентов линейных уравнений

i   j				…	i	…	n
	a11	a21	a31	…	ai1	…	an1
	a12	a22	a32	…	ai2	…	an2
	a13	a23	a33	…	ai3	…	an3
…	…	…	…	…	…	…	…
j	a1j	a2j	a3j	…	ajj	…	anj
…	…	…	…	…	…	…	…
r	a1r	a2r	a3r	…	ajr	…	anr
qi	q1	q2	q3	…	qi	…	qn

 

……………………………………………… (14.88)

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

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

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

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

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

Далее производятся последовательные действия для всех остальных строк до последней строки. Так, коэффициент при k_r равен сумме произведений обратных весов и коэффициентов 1-й и r -й строк матрицы.

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

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

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

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

Аналогично выполняются вычисления для остальных строк.

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

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

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

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

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

b₁₁k₁ + b₁₂k₂ + b₁₃k₃ + …+ b_1jk_j + …+ b_1rk_r + W₁ = 0

b₂₁k₁ + b₂₂k₂ + b₂₃k₃ + …+ b_2jk_j + …+ b_2rk_r + W₂ = 0

b₃₁k₁ + b₃₂k₂ + b₃₃k₃ + …+ b_3jk_j + …+ b_3rk_r + W₃ = 0

………………………………………………………. (14.89)

b_j1k₁ + b_j2k₂ + b_j3k₃ + …+ b_jjk_j + …+ b_jrk_r + W_j = 0

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

b_r1k₁ + b_r2k₂ + b_r3k₃ + …+ b_rjk_j + …+ b_rrk_r + W_r = 0

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

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

Полученные из решения уравнений (14.89) коррелаты k_j используются для вычисления поправок v_i по формулам (14.85) или (14.86). После введения поправок в измеренные величины получают уравненные значения измеренных величин (14.7).

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Шаг 5. Составляют, пользуясь табл. 14.3, условные уравнения поправок ν_i (14.80), (14.85) (или (14.86). Для значений поправок, например, получим:

ν₁ = q₁(a₁₁k₁ + a₁₂k₂ + …+ a_1jk_j +…+ a_1rk_r)

ν₂ = q₂(a₂₁k₁ + a₂₂k₂ + …+ a_2jk_j +…+ a_2rk_r)

ν₃ = q₃(a₃₁k₁ + a₃₂k₂ + …+ a_3jk_j +…+ a_3rk_r)

………………………………………….. (14.90)

ν_n = q_n(a_n1k₁ + a_n2k₂ + …+ a_njk_j +…+ a_nrk_r)

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

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

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

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

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

Если условные уравнения изначально существенно нелинейны и при разложении в ряд Тейлора, вообще говоря, недостаточно ограничиваться первыми членами разложения, то условия (14.9) могут не выполниться. В этом случае производят второе приближение уравнивания, считая уравненные из первого приближения значения x_i^' измеренными, а свободными членами W_j – остаточные невязки в уравнениях (14.9).

В § 137 рассмотрены примеры уравнивания различных геодезических построений коррелатным способом.