Усе
.pdf
|
|
|
|
a j2a ji |
|
|
|
|
|
|
|
a jia ji 2 a jia ji 1 a j2a j2 |
a j2a ji ; |
i, j 3,4 , |
(14) |
||||||
|
|
|
|
a j2a ji |
|
|
|
|
||
|
|
a jil j 2 a jil j 1 a j2a j2 |
a jil j ; |
i, j 3,4 , |
(15) |
|||||
|
|
|
|
a j3a ji |
|
|
|
|||
|
a jia ji 3 a jia ji 2 a j3a j3 |
a j3a ji ; |
i, j 4 , |
(16) |
||||||
|
|
|
|
a j3a ji |
|
|
|
|||
|
|
a jil j 3 a jil j 2 a j3a j3 |
a jil j ; |
i, j 4 , |
(17) |
|||||
|
|
|
|
a j4l j |
|
3 |
|
|
||
|
|
|
ps a j4a ji |
3 ; |
|
(18) |
||||
|
|
1 |
|
|
|
|
|
|
|
|
|
|
pc |
a j3a ji 2 a j3l j 2 a jia j4 2 ps , |
(19) |
||||||
|
|
1 |
|
|
|
|
|
|
|
|
c2 |
|
a j2a ji 1 a j2l j 1 a jia j3 |
1 pc a jia j4 1 ps , |
(20) |
||||||
c1 a j11a ji a j1l j a jia j2 c2 a jia j3 pc |
a jia j4 ps . |
(21) |
За допомогою калькулятора обчислення за формулами (12) – (21) виконують у спеціальній таблиці (табл. 1), котра має назву: «Схема Гаусса-Дулітля».
Розв’язування системи нормальних рівнянь способом визначників.
Для зручності запису алгоритму цього способу перепишемо систему рівнянь похибок (6) у матричному вигляді [2]
A p L V , |
(22) |
де А – матриця коефіцієнтів рівнянь похибок, p – вектор параметрів трансформації, L – вектор вільних членів параметричного рівняння, V – вектор похибок. У випадку 3-ох пунктів матриця А має вигляд:
1 |
0 |
x1 |
y1 |
|
|
|
0 |
1 |
y |
x |
|
|
|
|
|
1 |
1 |
|
|
|
1 |
0 |
x2 |
y2 |
|
, |
(23) |
A |
|
|
|
|
||
0 |
1 |
y2 |
x2 |
|
|
|
1 |
0 |
x3 |
y3 |
|
|
|
|
|
|
|
|
|
|
0 |
1 |
y3 |
x3 |
|
|
|
7
Таблиця 1
Схема Гаусса-Дулітля
|
a] |
|
b] |
|
|
с] |
|
|
|
|
d] |
|
|
|
|
|
l] |
|
|
|
|
S] |
Контроль |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I |
[aa] |
|
[ab] |
|
|
[ac] |
|
|
[ad] |
|
|
|
[al] |
|
|
[aS] |
a]+b]+ +l] |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ab |
|
|
ac |
|
|
|
|
ad |
|
|
|
|
|
al |
|
|
|
|
ab |
|
|||||||||||||||
E1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
aa |
|
|
aa |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
aa |
|
|
|
|
|
|
aa |
|
|
aS |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
II |
|
|
[bb] |
|
|
[bc] |
|
|
[bd] |
|
|
|
[bl] |
|
|
[bS] |
b]+c]+ +l] |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I E1 |
|
|
ab ab |
|
ac ab |
|
ad ab |
|
|
al ab |
|
aS ab |
|
|||||||||||||||||||||||||
|
aa |
|
aa |
|
|
|
aa |
|
|
|
|
aa |
|
|
|
aa |
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
[bb.1] |
|
|
[bc.1] |
|
|
[bd.1] |
|
|
|
[bl.1] |
|
|
[bS.1] |
|
|||||||||||||||||||||
II |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
bc.1 |
|
|
|
bd.1 |
|
|
|
|
bl.1 |
|
|
|
bS.1 |
|
|||||||||||||||||
E2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
bb.1 |
|
|
bb.1 |
|
|
bb.1 |
|
|
bb.1 |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
III |
|
|
|
|
|
|
[cc] |
|
|
[cd] |
|
|
|
|
[cl] |
|
|
[cS] |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
ac ac |
|
|
ad ac |
|
|
|
al ac |
|
|
aS ac |
|
||||||||||||||||||||
I E1 |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
aa |
|
|
|
aa |
|
|
|
|
aa |
|
|
|
aa |
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
II E2 |
|
|
|
|
|
bc.1 bc.1 |
|
|
bd.1 bc.1 |
|
|
bl.1 bc.1 |
|
|
bS.1 bc.1 |
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
bb.1 |
|
|
|
bb.1 |
|
|
|
bb.1 |
|
|
|
bb.1 |
|
||||||||||||||||||
III |
|
|
|
|
|
|
[cc.2] |
|
|
[cd.2] |
|
|
|
[cl.2] |
|
|
[cS.2] |
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
8
продовження табл. 1
|
a] |
|
b] |
|
с] |
|
|
|
|
|
|
d] |
|
|
|
|
|
|
|
|
|
|
l] |
|
|
|
|
S] |
Контроль |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
E3 |
|
|
|
|
|
|
|
|
|
|
cd.2 |
|
|
|
|
cl.2 |
|
|
|
cS.2 |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
cc.2 |
|
|
|
|
|
cc.2 |
|
|
|
|
cc.2 |
|
|||||||||||||||
IV |
|
|
|
|
|
|
|
|
|
|
|
[dd] |
|
|
|
|
|
|
|
[dl] |
|
|
|
|
[dS] |
|
|||||||||||||
I E1 |
|
|
|
|
|
|
|
ad ad |
|
|
|
al ad |
|
|
|
|
aS ad |
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
aa |
|
|
aa |
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
aa |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
bd.1 bd.1 |
|
|
bl.1 bd.1 |
|
|
|
bS.1 bd.1 |
|
|
||||||||||||||||||||||
II E2 |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
bb.1 |
|
|
|
|
|
bb.1 |
|
|
|
|
bb.1 |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
cd.2 cd.2 |
|
|
cl.2 cd.2 |
|
|
cS.2 cd.2 |
|
|
|
||||||||||||||||||||||
III E3 |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
cc.2 |
|
|
|
|
cc.2 |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
cc.2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
IV |
|
|
|
|
|
|
|
|
|
|
|
[dd.3] |
|
|
|
[dl.3] |
|
|
|
[dS.3] |
|
||||||||||||||||||
E4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dl.3 |
|
|
|
|
dS.3 |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dd.3 |
|
|
dd.3 |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dl.3 |
|
|
|
|
|
|
|
|
|
|
|||||
pS |
|
|
|
|
|
|
ps |
|
dl.3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
dd.3 |
|
|
dd.3 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
cd.2 |
ps |
|
|
|
|
cl.2 |
|
|
|
|
|
|
|
|
|
||||||||||||||||
pC |
|
|
|
pc d ] l] |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
cc.2 |
|
|
|
cc.2 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
c2 |
|
c2 c] d ] l] |
|
bc.1 |
pc |
|
bd.1 |
ps |
|
|
|
|
bl.1 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
bb.1 |
bb.1 |
|
|
|
bb.1 |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c1 |
b] c] d ] l] |
|
ab c2 |
ac pc |
|
ad |
|
ps |
|
|
|
|
al |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
aa |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
aa |
|
aa |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
aa |
|
|
|
|
|
|
|
|
|
9
а cкладові вектора вільних членів L обчисляться як різниці відповідних координат певних пунктів, відомих в обох системах координат, тобто
xT1 x1yT1 y1 L xT 2 x2 .yT 2 y2xT 3 x3
yT 3 y3
Система нормальних рівнянь запишеться так:
AT A p AT L 0 . |
|
Введемо позначення |
|
N ATA ; |
W ATL . |
Тепер, з формулами (26) система рівнянь (25) прийме вид
N p W 0.
(24)
(25)
(26)
(27)
Обчислимо головний визначник і визначники для кожного невідомого ( с1,с2, pc, ps), маємо:
|
N11 |
N12 |
N13 |
N14 |
|
|
|
|
W1 |
N12 |
N13 |
N14 |
|
|
|
|
N21 |
N22 |
N23 |
N24 |
; |
c1 |
|
W2 |
N22 |
N23 |
N24 |
, |
(28) |
||
|
N31 |
N32 |
N33 |
N34 |
|
|
|
|
W3 |
N32 |
N33 |
N34 |
|
|
|
|
N41 |
N42 |
N43 |
N44 |
|
|
|
|
W4 |
N42 |
N43 |
N44 |
|
|
|
|
|
N11 |
W1 |
N13 |
N14 |
|
|
|
N11 |
N12 |
W1 |
N14 |
|
|
|
|
|
|
|
|
|
|
|||||||||
c2 |
N21 |
W2 |
N23 |
N24 |
; |
pc |
|
N21 |
N22 |
W2 |
N24 |
, |
(29) |
||
|
|
N31 |
W3 |
N33 |
N34 |
|
|
|
N31 |
N32 |
W3 |
N34 |
|
|
|
|
|
N41 |
W4 |
N43 |
N44 |
|
|
|
N41 |
N42 |
W4 |
N44 |
|
|
|
N11 |
N12 |
N13 |
W1 |
|
|
ps |
N21 |
N22 |
N23 |
W2 |
. |
(30) |
|
N31 |
N32 |
N33 |
W3 |
|
|
|
N41 |
N42 |
N43 |
W4 |
|
|
Для обчислення визначників 4-го порядку за допомогою калькулятора можна використати схему Гаусса-Дулітля і продовжити обчислення, але цей процес
10
трудомісткий. Якщо ж визначники обчислювати за допомогою електронних таблиць Microsoft Excel, тоді робота зовсім спрощується. Для цього сформуємо такі матриці, як у формулах (28) - (30), в середовищі електронних таблиць Excel. Встановлюємо курсор у вільній комірці (клітинці), вибираємо функцію «МОПРЕД», виділяємо матрицю для визначника, спочатку позначеного через, і натискаємо клавішу «Enter», або комп’ютерною мишкою на клавішу «OK», що на екрані монітора, і отримане число буде шуканим визначником. Виконуючи описані дії послідовно для кожного з визначників ( с1, с2, pc, ps), отримаємо їх значення.
Тепер поправки до початкових значень параметрів трансформації обчисляться за формулами:
c |
|
c1 |
; |
c |
2 |
|
c2 |
; |
p |
c |
|
pc |
; |
p |
s |
|
ps |
. |
(31) |
|
|
|
|
||||||||||||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Розв’язування системи нормальних рівнянь способом оберненої матриці. Обчислення цим способом також краще виконувати на комп’ютері, можна скористатися програмою Microsoft Excel. Алгоритм запишемо у матричній формі і продовжимо з системи нормальних рівнянь, формули (25). Звідси вектор поправок до початкових значень параметрів трансформації обчислимо, як
p Q AT L ,
де
Q AT A 1 N 1 ,
Тут Q – обернена матриця коефіцієнтів нормальних рівнянь, а вектор такий вид:
c1 |
|
|
|
p c2 |
. |
pc |
|
|
|
ps |
|
(32)
(33) p має
(34)
Таким чином, систему нормальних рівнянь виду (10) або (25) можна розв’язати одним із цих способів.
3.3. Контроль розв’язку системи нормальних рівнянь
Як правило, перевірку правильності розв’язання системи рівнянь виконують підстановкою знайдених коренів системи у кожне рівняння і отримують тотожності. Визначені невідомі с1, с2, pc і ps, тобто вектор p з
11
формули (34), підставимо в систему нормальних рівнянь і повинні отримати нулі з точністю обчислень.
3.4. Обчислення остаточних значень параметрів трансформації
Рівняння похибок (6) або (22) розв’язуємо методом ітерацій (послідовних наближень) Ньютона, тому що вони є лінійним відображенням нелінійного процесу. Це означає, що для трансформації координат нами прийнята певна логічна модель з певною кількістю параметрів, а всі можливі інші параметри, наприклад похибки координат вихідних пунктів, не враховуємо. Тому отримані параметри з розв’язку такої системи рівнянь залежать, у деякій мірі, від «близькості» початкових значень шуканих параметрів. Звідси, остаточний результат можна отримати послідовно наближаючись до істинного значення кожний раз (на кожній ітерації) розв’язуючи задачу з новими початковими значеннями параметрів.
Таким чином, шукані параметри обчислимо так:
c c(0) |
c , |
|
|
||
1 |
1 |
|
1 |
|
|
c2 c2(0) |
c2 |
, |
(35) |
||
pc pc(0) |
|
|
|||
pc , |
|
||||
ps ps(0) |
ps , |
|
де c1(0) , c2(0) , pc(0) , ps(0) – початкові значення параметрів трансформації.
3.5. Оцінка точності
Середня квадратична похибка одиниці ваги дорівнює:
|
6 |
|
|
|
|
|
|
|
|
v2j |
|
|
|
|
|
|
|
|
|
|
VT V |
|
|
|
||
|
j 1 |
|
, |
(35) |
||||
2n k |
2n k |
|||||||
|
|
|
|
|
|
де n - кількість пунктів, k - кількість шуканих параметрів трансформації. Оцінку точності кожного параметра з номером j з вектора p обчислимо так:
|
|
|
|
|
|
m p |
jj |
Q jj . |
(36) |
||
|
|
|
|
|
12
4. Числовий приклад Увага!!! Кожний пункт роботи обчислений за іншими вихідними даними. Зробити рисунок 1.
|
|
4.1. Вихідні дані |
|
|
|
||
1. Координати 3-ох пунктів X x |
y T |
на площині в проекції референц- |
|||||
еліпсоїда Красовського. |
|
|
|
|
|
||
2. Координати тих самих пунктів |
XT xT |
yT T на площині в проекції |
|||||
референц-еліпсоїда WGS-84 XT. |
|
|
|
|
|||
|
|
|
|
|
|
Таблиця 2 |
|
|
|
|
|
|
|
|
|
Пункти |
x, м |
y, м |
|
|
xT, м |
yT, м |
|
|
|
|
|
|
|
||
1 |
5527323,8576 |
215835,4775 |
5527095,8393 |
215894,7532 |
|
||
|
|
|
|
|
|
||
3 |
5528423,4113 |
229697,6169 |
5528195,2095 |
229756,8119 |
|
||
|
|
|
|
|
|
||
7 |
5512671,2450 |
220736,0700 |
5512443,2108 |
220795,3295 |
|
||
|
|
|
|
|
|
|
|
4.2. Обчислення коефіцієнтів і вільних членів рівнянь похибок
Таблиця 3
№ |
aj1 |
aj2 |
aj3 |
aj4 |
lj |
vj |
vj2 |
з/п |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
1 |
0 |
5527323,8576 |
215835,4775 |
228,0183 |
-0,0646 |
0,004175 |
2 |
0 |
1 |
215835,4775 |
-5527323,8576 |
-59,2757 |
-0,0047 |
0,000022 |
3 |
1 |
0 |
5529729,9610 |
223240,1882 |
228,2947 |
0,0190 |
0,000360 |
4 |
0 |
1 |
223240,1882 |
-5529729,9610 |
-59,1290 |
0,0765 |
0,005854 |
5 |
1 |
0 |
5523900,0377 |
222180,1830 |
228,2079 |
0,0457 |
0,002084 |
6 |
0 |
1 |
222180,1830 |
-5523900,0377 |
-59,1718 |
-0,0718 |
0,005162 |
|
|
|
|
|
|
vj2 = |
0,017656 |
4.3. Обчислення коефіцієнтів і вільних членів нормальних рівнянь
Таблиця 4
|
aj1] |
aj2] |
aj3] |
aj4] |
lj] |
|
|
|
|
|
|
[aj1 |
3 |
0 |
16580953,8563 |
661255,8487 |
684,5209 |
[aj2 |
|
3 |
661255,8487 |
-16580953,86 |
-177,5765 |
[aj3 |
|
|
9,17884793 1013 |
0 |
3744196090,046 |
[aj4 |
|
|
|
9,17884793 1013 |
1132344765,337 |
13
4.4. Розв’язування системи нормальних рівнянь способом Гаусса
Таблиця 5
|
aj1] |
|
aj2] |
aj3] |
aj4] |
|
lj] |
|
|
|
|
|
|
|
|
|
|
|
|
I |
|
3 |
0 |
16580953,86 |
661255,8487 |
|
|
684,5209 |
|
E1 |
|
|
0 |
-5526984,619 |
-220418,6162 |
|
-228,1736333 |
|
|
II |
|
|
3 |
661255,8487 |
-16580953,86 |
|
-177,5765 |
|
|
I E1 |
|
|
0 |
0 |
0 |
|
|
0 |
|
II |
|
|
3 |
661255,8487 |
-16580953,86 |
|
-177,5765 |
|
|
E2 |
|
|
|
-220418,6162 |
5526984,619 |
|
59,19216667 |
|
|
III |
|
|
|
9,17884793 1013 |
0 |
|
3744196090,046 |
|
|
I E1 |
|
|
|
-9,16426769 1013 |
-3,65475090 1012 |
|
-3783336486 |
|
|
II E2 |
|
|
|
-1,45753099 1011 |
3,65475090 1012 |
|
39141166,41 |
|
|
III |
|
|
|
49236173,82 |
0 |
|
770,9274191 |
|
|
E3 |
|
|
|
|
0 |
|
-1,56577443 10-5 |
|
|
IV |
|
|
|
|
9,17884793 1013 |
|
1132344765,337 |
|
|
I E1 |
|
|
|
|
-1,45753099 1011 |
|
-150881149,6 |
|
|
II E2 |
|
|
|
|
-9,16426769 1013 |
|
-981462584,2 |
|
|
III E3 |
|
|
|
|
0 |
|
|
0 |
|
IV |
|
|
|
|
49236173,81 |
|
1031,62182 |
|
|
E4 |
|
|
|
|
|
|
-2,0952518 10-5 |
|
|
ps |
|
|
|
|
-2,09525180 10-5 |
|
-2,0952518 10-5 |
|
|
pc |
|
|
|
-1,56577443 10-5 |
0 |
|
-1,5657744 10-5 |
|
|
с2 |
|
|
-53,1608 |
3,451258328 |
-115,8042449 |
|
59,19216667 |
|
|
с1 |
-137,0152 |
0 |
86,54011182 |
4,618325033 |
|
-228,1736333 |
|
||
|
4.5. Контроль розв’язку системи нормальних рівнянь |
|
|
||||||
|
|
|
|
|
|
|
Таблиця 6 |
||
-411,0456 |
|
+0,0000 |
-259,6203 |
-13,8550 |
+684,5209 |
=0,0000 |
|
||
0,0000 |
|
-159,4825 |
-10,3538 |
+347,4127 |
-177,5765 |
=0,0000 |
|
||
-2271842650,5 |
-35152903,1 |
-1437200536,4 |
+0,0000 +3744196090,0 |
=0,0000 |
|
-90602100,0 +881457102,2 +0,0000 -1923199767,5 +1132344765,3 =0,0000
4.6. Обчислення параметрів трансформації двовимірних координат пунктів
c1 c1(0) c1 = 0,0000 - 137,0152 = -137,0152 м , c2 c2(0) c2 = 0,0000 - 53,1608 = -53,1608 м ,
14
pc pc(0) pc = 0,0000 - 1,56577443 10-5 = -1,56577443 10-5 , ps ps(0) ps = 0,0000 - 2,09525180 10-5 = -2,09525180 10-5 ,
pc2 ps2 = 2,615670 10-5 ,
ps
arctg = 53,22936496 .
pc
4.7.Оцінка точності
Середня квадратична похибка одиниці ваги дорівнює:
|
|
0,017656 |
0,0940м. |
|
6 - 4 |
|
|||
|
|
|
|
Середні квадратичні похибки параметрів трансформації наступні: mc1 mc2 0,0940 621416,2756 74,067 м;
mpc mps 0,0940 2,03103 10 8 1,3 10 5 .
5.Література
1.Б. Гофманн-Велленгоф, Г. Ліхтенеггер, Д. Коллінз. Глобальна система визначення місцеположення (GPS). Теорія і практика.- Пер. з англ. під ред. Я.С.Яцківа.- Київ: Наук. думка, 1995.-380 с.
2.Маркузе Ю.И. Алгоритмы для уравнивания геодезических сетей на ЭВМ.- М.: Недра, 1989.- 248 с.
3.Цюпак І.М., Дульцев А.Т., Третяк К.Р., Савчук С.Г. Точність перетворення просторових координат пунктів //Сучасні досягнення геодезії, геодинаміки та геодезичного виробництва.- Львів, 2000.- С.
45-50.
4.Цюпак І.М., Я. Бойчук, В. Джаман. Визначення параметрів зв’язку між Державною та місцевою системами координат за допомогою GPS-вимірів//„Геоінформаційний моніторинг навколишнього середовища - GPS і GIS-технологій” (Збірник наукових праць ХІ міжнародного симпозіуму.- Алушта, вересень 2006 р.).- Львів, 2006.-
С. 36 – 39.
5.Цюпак І.М., Савчук С.Г., Бойчук Я.Д. Про точність двовимірного перетворення координат пунктів//“Геодезія, картографія і аерофотознімання”.- Міжвідомчий науково-технічний збірник.- Вип.68.- Львів:НУ ”Львівська політехніка”.-2007.- С.132-142.
15
НАВЧАЛЬНЕ ВИДАННЯ
МЕТОДИЧНІ ВКАЗІВКИ
до виконання лабораторної роботи з курсу “Космічна геодезія”
для студентів спеціальності 7.070901 “Геодезія”
Автори: Марченко Олександр Миколайович, докт.фіз.-матем.наук, проф. Цюпак Ігор Михайлович, канд.техн.наук, доц.
Редактор
Комп’ютерне складання
Підписано до друку Формат 70 1001/16 . Папір офсетний.
Друк на різографі. Умови друк. арк. 16. Обл.-вид. арк. Наклад 100 прим. Зам.
Поліграфічний центр Видавництва Національного університету «Львівська політехніка»
Вул. Ф. Колесси, 2, 79000, Львів
16