Z Wyrównanie obserwacji bezpośrednich, Geodezja i Kartografia, Rachunek Wyrównawczy


Wyrównanie obserwacji bezpośrednich. Podstawowe definicje i oznaczenia

Zdefiniujemy pewne pojęcia, którymi często będziemy się posługiwali.

Błędem prawdziwym 0x01 graphic
, i-tej obserwacji 0x01 graphic
nazywamy różnicę między wartością prawdziwą L mierzonej wielkości, a wartością zaobserwowaną:

0x01 graphic

Błędem pozornym 0x01 graphic
albo poprawką i-tej obserwacji nazwiemy różnicę między wartością najbardziej prawdopodobną 0x01 graphic
mierzonej wielkości, a wartością zaobserwowaną 0x01 graphic

0x01 graphic

Z powyższych określeń wynika, że błędami prawdziwymi można operować w tych przypadkach, gdy znana jest wartość prawdziwa mierzonej wielkości. Błędami pozornymi posługujemy się wtedy, gdy wartość prawdziwa nie jest znana, lecz potrafimy obliczyć wartość bliską prawdziwej, którą nazywa się powszechnie wartością najprawdopodobniejszą. W przypadku, gdy dokonujemy wielokrotnych pomiarów tego samego elementu geometrycznego, np. kąta, długości za wartość najprawdopodobniejszą uważamy średnią arytmetyczną uzyskanych wyników.

Przykład. Wykonano trzykrotny pomiar trzech kątów trójkąta otrzymując następujące wyniki:

1

70°02'

50°18'

59°36'

2

70°01'

50°20'

59°37'

3

70°04'

50°19

59°58'

Obliczyć błąd prawdziwy sumy kątów trójkąta dla każdej z trzech serii.

Ponieważ wartość prawdziwa sumy kątów w trójkącie wynosi L = 180°, natomiast z obserwacji otrzymaliśmy odpowiednio 0x01 graphic
; 0x01 graphic
, 0x01 graphic
, więc otrzymamy następujące wartości dla błędów prawdziwych:

0x01 graphic
; 0x01 graphic
, 0x01 graphic

Przykład. Pomiar kąta, którego wartości prawdziwej nie znamy, dał następujące wyniki: 37°05', 37o05'; 37°01'. Obliczyć błędy pozorne obserwacji.

Jako wartość ostateczną mierzonego kąta przyjmiemy średnią arytmetyczną z uzyskanych wyników. Wynosi ona 0x01 graphic
skąd, mamy:

0x01 graphic
, 0x01 graphic
, 0x01 graphic

Jak widać każdy z obliczonych błędów prawdziwych, czy pozornych, dotyczy odpowiedniej i tylko jednej obserwacji. W pierwszym przykładzie obserwacją jest każda z trzech zmierzonych sum kątów trójkąta; w przykładzie drugim za obserwację uważamy każdy z trzech uzyskanych wyników.

Przykład wprowadzający. Załóżmy, że oddając do tarczy serię 100 strzałów, każdy z dwóch strzelców uzyskuje przeciętnie następujące wyniki:

I strzelec

II strzelec

Liczba uzyskanych punktów

5

2

1

3

2

1

Liczba trafień

70

15

15

60

40

0

Postawmy sobie zadanie wyznaczenia kryterium liczbowego, które pozwoliłoby dać odpowiedź na pytanie, który z nich jest lepszy? Za liczbę charakteryzującą jakość każdego ze strzelców przyjmiemy średnią ilość punktów przypadającą na jeden strzał każdego z nich. Łatwo, bowiem można policzyć liczbę punktów uzyskiwanych średnio w 100 strzałach. Wynosi ona:


dla I strzelca PI = 5 • 70 + 2 • 15 + 1 • 15 = 255;

dla II strzelca PII =3 • 60 + 2 • 40 + 0 • l = 260.

Stąd już widać, że pierwszy strzelec jednym strzałem uzyskuje przeciętnie 2,55 punktu, zaś drugi 2,60 punktu. Uznamy, więc, że drugi jest lepszy.

Rozpatrzmy teraz zadanie formalnie podobne do powyższego, lecz dotyczące teorii błędów. Załóżmy, że pewien kąt pomierzono w dwóch seriach wykonu­jąc po 100 obserwacji w każdej z nich, a następnie obliczono błędy prawdziwe uzyskanych wyników i zestawiono je w dwóch tabelach, oddzielnie dla każdej serii, otrzymując:

I seria

-1"

0"

+1"

II seria

-10"

0

+10"

Ilość wystąpień

25

50

25

Ilość wystąpień

25

50

25

Jak na tej podstawie ocenić, który z tych dwóch zespołów obserwacji jest lepszy? Gdybyśmy chcieli rozumować analogicznie jak w zadaniu poprzednim to należałoby obliczyć średnią wartość błędu prawdziwego dla każdej serii i dwie otrzymane liczby porównać. Tą serię, której średnia wartość błędu prawdziwego byłaby większa uznalibyśmy za gorszą. Zauważmy, że kryterium to nie nadaje się jednak do charakterystyki jakości tych dwóch zespołów obserwacji. Średnia wartość błędu prawdziwego v; obu przypadkach jest równa zeru, chociaż widać wyraźnie, że drugi pomiar dał wyniki gorsze. Przypadek ten jest charakterystyczny dla błędów pomiarów geodezyjnych. Zestawiając szeregi błędów prawdziwych wielkości zmierzonej n-krotnie, widzimy, że przy dużych wartościach n = 500; 400; 500 obserwacji wystąpi, w przybliżeniu, ta sama liczba błędów dodatnich i ujemnych mających tę samą wartość bezwzględną. Zmusza nas to do wybrania innego kryterium liczbowego charakteryzującego jakość zespołu obserwacji. Musimy je zdefiniować w ten sposób by zlikwidować wpływ różnych znaków na ostateczny wynik. Obliczmy np. wartość średnią bezwzględnej wartości błędu prawdziwego dla każdej serii. Otrzymamy:

0x01 graphic
0x01 graphic

Liczby dI i dII nazwiemy błędami przeciętnymi. Jak widać można je przyjąć za liczbowe kryterium jakości obserwacji, choć w geodezji jest ono rzadko stosowane, z następującego powodu.

Często zachodzi konieczność scharakteryzowania jakości pewnej funkcji obserwacji

0x01 graphic

gdy błędy obserwacji są znane. Zwykle mówi się wtedy, że chcemy obliczyć dokładność funkcji y, znając dokładność argumentów 0x01 graphic
. Gdybyśmy za miernik dokładności przyjęli błąd przeciętny, okazałoby się, że obliczenie błędu przeciętnego funkcji y na podstawie błędów przeciętnych obserwacji byłoby uciążliwe. Rezygnujemy, więc z błędu przeciętnego i wprowadzimy inaczej zdefiniowaną liczbę jako miernik dokładności zespołu obserwacji. Jest nią wartość średnia kwadratu błędu prawdziwego, którą nazywamy kwadratem błędu średniego i oznaczamy symbolem 0x01 graphic
, Dla naszego przykładu możemy, więc obliczyć:

0x01 graphic

skąd otrzymujemy następujące błędy średnie poszczególnych serii:

0x01 graphic

Możemy teraz powiedzieć, że druga seria jest mniej dokładna, gdyż jej błąd średni jest większy. Okazuje się, że znając błędy średnie obserwacji można obliczyć błąd średni funkcji obserwacji, bowiem, jak to dalej zobaczymy, błąd średni łatwo przenosi się na sumę.

Możemy stwierdzić, że w geodezji, jako kryterium dokładności obserwacji przyjmuje się błąd średni. W przypadku, gdy mamy duży szereg błędów prawdziwych tej samej wielkości, jej błąd średni wyznaczamy z wzoru

0x01 graphic
.

gdzie n jest liczbą wszystkich błędów prawdziwych rozważanego szeregu.

Nie zawsze istnieje możliwość obliczenia błędów prawdziwych; na ogół nie są nam znane wartości prawdziwe mierzonej wielkości. W tych przypadkach wielokrotny pomiar tego samego elementu geometrycznego pozwala obliczyć błędy pozorne v oraz wyznaczyć błąd średni ze wzoru

0x01 graphic
.

W niektórych przypadkach, na obliczenie błędu średniego istnieją wzory empiryczne, np.

0x01 graphic
na wyznaczenie błędu średniego długości.

Prawo przenoszenia się błędów średnich

Wyznaczymy błąd średni funkcji obserwacji

0x01 graphic

zakładając, że błędy średnie obserwacji są dane i wynoszą odpowiednio 0x01 graphic
.

Dla funkcji

0x01 graphic

gdzie a jest stałą, zaś m dane np. obliczone na podstawie szeregu błędów prawdziwych dla konkretnych wartości 0x01 graphic
funkcja y przybierze wartości 0x01 graphic
z błędami prawdziwymi równymi odpowiednio 0x01 graphic
. Błąd średni funkcji y można obliczyć ze wzoru 0x01 graphic
. Mamy, więc

0x01 graphic
.

Natomiast błąd średni 0x01 graphic
funkcji

0x01 graphic

przy założeniu znajomości błędów średnich 0x01 graphic
i 0x01 graphic
obserwacji 0x01 graphic
i 0x01 graphic
określony jest wzorem: 0x01 graphic
.

Wykorzystując powyższe zależności błąd średni funkcji liniowej obserwacji 0x01 graphic
wyznacza się z zależności:

0x01 graphic
.

Natomiast błąd średni dowolnej funkcji 0x01 graphic
wyznaczamy przy założeniu, że jest ona rozwijalna w szereg Taylora. Przyjmując bowiem przybliżone wartości 0x01 graphic
obserwacji x, y ,..., t , można przedstawić każdą wielkość zaobserwowaną 0x01 graphic
, jako sumę odpowiadającej jej wartości przybliżonej 0x01 graphic
i niewielkiego przyrostu df uzależnionego od przyrostów obserwacji 0x01 graphic
tj.: 0x01 graphic
. Z wzorów tych wynika, że błąd średni 0x01 graphic
obserwacji x jest równy błędowi średniemu przyrostu dx, bowiem błąd średni stałej 0x01 graphic
jest, oczywiście, równy zeru. Stąd, gdy rozwiniemy funkcję f w szereg Taylora w sąsiedztwie punktu 0x01 graphic
zachowując tylko wyrazy pierwszego rzędu, otrzymamy:

0x01 graphic
.

(Zauważmy, że wartości pochodnych cząstkowych wyznaczono w punkcie 0x01 graphic
.)

Jest to przybliżenie liniowe funkcji f o zmiennych 0x01 graphic
. Ponieważ błędy średnie 0x01 graphic
tych zmiennych są znane, a błąd średni stałej 0x01 graphic
jest równy zeru, to można zastosować wzór na błąd średni funkcji liniowej. Otrzymujemy, więc związek

0x01 graphic
.

Przykład. Obliczyć błąd średni różnicy wysokości w ciągu niwelacyjnym złożonym z 0x01 graphic
stanowisk, jeśli błąd średni odczytu łaty 0x01 graphic

Ponieważ przewyższenie na i-tym stanowisku obliczamy jako różnicę: 0x01 graphic
odczytów wstecz i w przód, dlatego dla n stanowisk będzie: 0x01 graphic
.

Stosując tu wzór na błąd funkcji liniowej mamy: 0x01 graphic
.

Przykład. Obliczyć błąd średni boku L = 180 m zmierzonego taśmą 20 m, gdy błąd średni jednego przyłożenia taśmy wynosi m = 1 cm.

Ponieważ uzyskany wynik jest rezultatem dziewięciokrotnego odłożenia taśmy, więc możemy napisać 0x01 graphic
, gdzie 0x01 graphic
- i-te przyłożenie taśmy, a 0x01 graphic
.

Stąd mamy 0x01 graphic
.

Przykład. Obliczyć błąd średni powierzchni prostokąta o bokach 0x01 graphic
, 0x01 graphic
, wiedząc, że 0x01 graphic
, 0x01 graphic

Ponieważ 0x01 graphic
mamy: 0x01 graphic

Wyrównanie obserwacji bezpośrednich jednakowo dokładnych

Załóżmy, że z tą samą dokładnością pomierzyliśmy n razy pewien element geometryczny otrzymując wyniki: 0x01 graphic
i chcemy wyznaczyć jeden ostateczny rezultat L tej wielokrotnie zaobserwowanej wielkości. Żądamy przy tym, aby był on w pewnym sensie najlepszy, najbardziej zbliżony do prawdziwego. Należy zwrócić uwagę, że wyrażenie "wynik najlepszy" jest pozbawione ścisłego znaczenia dopóki nie ustali się kryterium tej "najlepszości". Oznaczmy przez v różnicę między nieznaną na razie, wartością najwłaściwszą liczby L, a wartością i - tej obserwacji, czyli

0x01 graphic
.

Przyjmiemy, że znane jest przybliżenie 0x01 graphic
wielkości wyznaczanej L, różne od niej o niewielką wartość x. Będzie wtedy

0x01 graphic
.

Jeżeli oznaczyć 0x01 graphic
to otrzymamy n równań

0x01 graphic
,

z których nie da się jednoznacznie wyznaczyć 0x01 graphic
niewiadomych 0x01 graphic
bez poczynienia dodatkowych założeń. Wprowadźmy warunek, który pozwoli rozwiązać zadanie jednoznacznie, precyzując tym samym pojęcie najlepszego wyniku L. Wyznaczymy taką wartość 0x01 graphic
, aby błąd średni pojedynczego spostrzeżenia obliczony ze wzoru

0x01 graphic
(m1)

przyjął wartość najmniejszą. Wtedy, oczywiście, musi być 0x01 graphic
minimum, ponieważ n - 1 jest stałe przy ustalonej liczbie obserwacji n. Sprowadziliśmy, więc zadanie do znalezienia minimum funkcji:

0x01 graphic
,

gdzie: 0x01 graphic
tj.

0x01 graphic
.

Skąd, wyznaczając pochodną i przyrównując ją do zera otrzymamy

0x01 graphic

czyli 0x01 graphic
, a więc ostatecznie

0x01 graphic

i

0x01 graphic
.

Za szukaną wartość L wyznaczaną w oparciu o szereg obserwacji jednakowo dokładnych przyjmiemy, więc średnią arytmetyczną uzyskanych wyników. Uważamy ją za najlepszą w tym sensie, że suma kwadratów różnic 0x01 graphic
jest najmniejsza.

Różnice te nazwaliśmy błędami pozornymi. Łatwo jest uzasadnić, że suma błędów pozornych jest równa zeru, czyli

0x01 graphic
.

Ponieważ

0x01 graphic

skąd po zsumowaniu otrzymamy

0x01 graphic
.

Błąd średni M średniej arytmetycznej L oblicza się ze następującego wzoru (Dlaczego ? Dowiesz się w dalszej części wykładu.)

0x01 graphic
.

Należy nadmienić, że przy niewielkiej liczbie obserwacji n, wartość m obliczona z wzoru (m1) ma cechę dużej przypadkowości. Lepiej jest wyznaczyć m ze znanego wzoru empirycznego, co w geodezji jest bardzo często możliwe.

Przykład. Zmierzono bok sześciokrotnie otrzymując wyniki zestawione w tabeli. Obliczyć średnią arytmetyczną, jej średni błąd oraz błąd średni pojedynczego spostrzeżenia.

Nr obs.

0x01 graphic
[m]

0x01 graphic
[cm]

0x01 graphic
[cm]

0x01 graphic

1

214,27

+ 5

- 2,5

6,25

2

,22

0

+ 2,5

6,25

5

,22

0

+ 2,5

6,25

4

,28

+ 6

- 5,5

12,25

5

,25

+ 5

- 0,5

0,25

6

,25

+ 1

+ 1,5

2,25

0x01 graphic
[t] = 15 [v]= 0,0 [vv] = 35,50

(Aby operować mniejszymi liczbami od wszystkich wyników można odjąć, ten którego wartość jest najmniejsza. W naszym przypadku będzie to 0x01 graphic
)

0x01 graphic
;

0x01 graphic
;

0x01 graphic
;

0x01 graphic
.

Wyrównanie obserwacji bezpośrednich o różnej dokładności

Niech będzie dany szereg obserwacji tej samej wielkości 0x01 graphic
o błędach średnich równych 0x01 graphic
. Obliczymy najwłaściwszą wartość mierzonej wielkości, jej średni błąd oraz błąd średni tzw. "typowego spostrzeżenia" (tj. spostrzeżenia o wadze p = 1).

Powyższe zadanie można sprowadzić do znanego nam przypadku obserwacji jednakowo dokładnych. Zauważmy, że każdą wartość 0x01 graphic
można traktować jako średnią arytmetyczną z liczby 0x01 graphic
jednakowo dokładnych obserwacji o błędzie średnim spostrzeżenia równym 0x01 graphic
. Będzie, więc

0x01 graphic
, a stąd 0x01 graphic

Obserwację 0x01 graphic
traktujemy, więc jako średnią arytmetyczną z 0x01 graphic
fikcyjnych, jednakowo dokładnych obserwacji, 0x01 graphic
będzie średnią arytmetyczną z 0x01 graphic
fikcyjnych obserwacji o tym samym błędzie średnim 0x01 graphic
itd. Wynika stąd, że wartość 0x01 graphic
można wyznaczyć ze wzoru

0x01 graphic

gdzie 0x01 graphic
, jest nieznaną nam sumą wartości fikcyjnych, jednakowo dokładnych obserwacji, potrzebnych do wyznaczenia średniej. Poszukiwaną wartość średnią obliczymy, więc jako średnią arytmetyczną ze wzoru:

0x01 graphic
,

gdzie w liczniku występuje suma wartości wszystkich fikcyjnych obserwacji, natomiast w mianownika liczba tych obserwacji. Otrzymany wynik można też napisać w postaci:

0x01 graphic

Liczbę L nazywamy ogólną średnią arytmetyczną a 0x01 graphic
wagą i-tej obserwacji.

W zastosowaniach praktycznych stosuje się inną postać powyższego wzoru. Ponieważ wartości 0x01 graphic
niewiele odbiegają od, siebie, przeto każdą z nich można przedstawić jako sumę pewnej stałej L oraz niewielkiego przyrostu 0x01 graphic
, czyli: 0x01 graphic
. Podstawiając prawą stronę powyższego związku do wzoru otrzymamy:

0x01 graphic
.

Zwykle za 0x01 graphic
obiera się najmniejszy z zaobserwowanych wyników.

Z uzyskanych wzorowych wynika, że użycie w rachunku k - krotnych wartości wag nie zmienia wartości ogólnej średniej arytmetycznej. Ponieważ błąd średni typowego spostrzeżenia liczymy ze wzoru: 0x01 graphic
, gdzie 0x01 graphic
więc błąd średni M ogólnej (ważonej) średniej arytmetycznej L przedstawia następujący związek:

0x01 graphic
.

Wynika to z tego, że ogólną średnią arytmetyczną traktujemy jako zwykłą średnią z 0x01 graphic
jednakowo dokładnych obserwacji. Ze wzoru wynika, że użycie w rachunku k - krotnych wartości wag nie zmienia błędu średniego ogólnej średniej arytmetycznej. Analogicznie jak w przypadku obserwacji jednakowo dokładnych można udowodnić, że [pv] = 0, oraz że błędy pozorne obliczone na podstawie ogólnej średniej arytmetycznej spełniają warunek 0x01 graphic
.

Przykład. Pomierzono bok poligonowy otrzymując wyniki, które wraz z błędami średnimi zamieszczono w poniższej tabeli. Obliczyć L, m, M.

Nr obs.

0x01 graphic

m [cm]

0x01 graphic

0x01 graphic
[cm]

0x01 graphic

0x01 graphic

pv

0x01 graphic
[cm2]

1

128,36

4

1

3

3

- 2,1

- 2,1

4,41

2

,34

1

16

1

16

- 0.1

- 1,6

0,16

3

,33

2

4.

0

0

+ 0,9

+ 3,6

3,24

L0 = 128

[p] = 21

[pt] = 19

[pv] = -0,1

[pvv]=7.81

0x01 graphic
, 0x01 graphic
,

0x01 graphic
,

0x01 graphic
.

Przykład. Pomierzono kąt otrzymując następujące wyniki;

68°32' jako średnia z 2 obserwacji,

68° 34' jako średnia z 4 obserwacji;

68°35' jako średnia z 6 obserwacji.

Obliczyć wartość kąta, błąd średni pojedynczego pomiaru i błąd średni średniej arytmetycznej. Przyjmujemy, że wagi poszczególnych wyników są równe liczbie obserwacji.

Nr obs.

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

pv

0x01 graphic
[cm2]

1

68,32

2

0

0

+2,2

+4,4

9.68

2

,34

4

2

8

+0,2

+0,8

0,16

5

,35

6

3

18

-0,8

-4,8

3.84

0x01 graphic
0x01 graphic
0x01 graphic
0x01 graphic
0x01 graphic

0x01 graphic
, 0x01 graphic
, 0x01 graphic
, 0x01 graphic
.

Metoda spostrzeżeń pośredniczących. Uwagi ogólne.

Załóżmy, że chcemy wyznaczyć wysokości punktów A i B (rys. m2) w odniesieniu do pewnego poziomu, przy czym znane są wysokości punktów R1 i R2 w odniesieniu do tego poziomu (tj. w tym samym układzie odniesienia).

0x01 graphic

Rys. m2.

Ponieważ bezpośredni pomiar wysokości punktów A i B na ogół nie jest możliwy, to wyznacza się je na drodze pomiaru pośredniego na przykład mierząc różnice wzniesień między którymkolwiek z punktów R1, R2 a punktami A i B. Pomiar dwóch ciągów (np. R1A i AB) pozwala pośrednio obliczyć wysokości punktów A i B, a tym samym rozwiązać postawione zadanie.

W geodezji rzadko stosuje się konstrukcje dające matematycznie jednoznaczne rozwiązanie. Jest niemal regułą, że materiał polowy zawiera obserwacje nadliczbowe, a to z następujących powodów;

1) można wtedy stwierdzić, czy wyniki pomiarów nie są obarczone błędami grubymi,

2) możliwe jest obliczenie dokładności wyznaczanych wielkości,

5) dokładność takich konstrukcji jest większa.

Te względy powodują, że oprócz obserwacji koniecznych do jednoznacznego wyznaczenia wysokości punktów, dokonuje się zwykle szeregiem dodatkowych obserwacji pozwalających obliczyć te wysokości na innej drodze. Gdybyśmy w rozważanym przykładzie oprócz ciągów R1A i AB rozważyli jeszcze ciąg BR2 mielibyśmy trzy wartości pomierzonych przewyższeń. Mając dane wysokości punktów R1, R2 i wybierając z tego zespołu po dwie dowolne obserwacje, zawsze jesteśmy w stanie obliczyć poszukiwane wysokości. Wartość wysokości każdego punktu będzie na ogół zależna od tego, jakie obserwacje uwzględniono w rachunku. Różnice wyników tych obserwacji pozwalają zwykle upewnić się, że w pomiarze nie ma omyłek, a ponadto obliczyć dokładność wyznaczanej wysokości.

W rozważanej przez nas sieci niwelacyjnej, mierząc wszystkie trzy ciągi R1A, AB, BR2, mamy o jedną obserwację więcej niż to jest potrzebne do jednoznacznego określenia wysokości punktów A i B. Będziemy, więc mieli jedną obserwację nadliczbową. Obserwacjami nadliczbowymi będziemy dalej nazywali te obserwacje, których wykonanie nie jest konieczne do jednoznacznego rozwiązania danej konstrukcji geodezyjnej. Aby ustalić ich ilość należy przede wszystkim ustalić ilość obserwacji koniecznych. Różnica między wszystkimi wykonanymi obserwacjami, a ilością obserwacji koniecznych, daje ilość obserwacji nadliczbowych. Można to ująć wzorem: 0x01 graphic
, gdzie nn - ilość obserwacji nadliczbowych, nw - ilość wszystkich obserwacji, nk - ilość obserwacji koniecznych.

Rachunek wyrównawczy zajmuje się matematyczno rachunkowym opracowaniem wyłącznie takich konstrukcji geodezyjnych, w których występują obserwacje nadliczbowe. Jako cel stawia sobie obliczenie najbardziej prawdopodobnych wartości wielkości poszukiwanych oraz wyznaczenie ich średnich błędów. To matematyczne opracowanie nosi nazwę wyrównania, które można przeprowadzić na różnych formalnie drogach. Sposoby, jakimi się tego dokonuje, noszą nazwę metod. Jedną z nich jest omawiana metoda pośrednicząca. Celem wniknięcia w jej istotę posłużymy się prostym przykładem rachunkowym. Trzeba zaznaczyć, że ten sposób rozwiązania nie jest stosowany w praktyce, gdzie podobne zadania opracowuje się w na innej drodze.

Weźmy pod uwagę sieć niwelacyjną przedstawioną na rysunku m2. Opierając się na dwóch punktach stałych (reperach) oznaczonych na szkicu literami R1 i R2, chcemy wyznaczyć wysokości punktów A i B mając zaobserwowane przewyższenia h1, h2, h3 (różnice wysokości między punktami R1 i A itd. Strzałki oznaczają kierunki wzniesień ( mogą też oznaczać spadki), a długości ciągów są w przybliżeniu równe.

Niech wysokość punktu R1, H1 = 0,000 m, natomiast punktu R2, H2 = 6,000 m, a zaobserwowane różnice wzniesień między punktami sieci wynoszą:

dla ciągu R1A: h1= 0,996 ,

dla ciągu AB: h2 = 1,999 ,

dla ciągu BR2: h3 = 2,999 .

Widać, że posiadamy tu obserwacje nadliczbowe (jedną!), trzeba więc będzie przeprowadzić wyrównanie tej sieci. Aby pokazać, że jest to konieczne, obliczymy wysokość jednego tych punktów, np. punktu A, opierając się na wielkościach pomierzonych:

0x01 graphic
, 0x01 graphic
.

Zauważmy, że w zależności od tego, jaką kombinację zaobserwowanych przewyższeń uwzględnimy w rachunku, otrzymamy różne wartości dla wysokości punktu A. Dla punktu B rzecz przedstawia się podobnie. Aby każdej wysokości zapewnić jednoznaczną wartość, niezależną od wyboru ciągów, jakimi się posłużymy do jej wyznaczenia można zniekształcić nasze (trzy) obserwacje przyporządkowując im poprawki. Poprawione wyniki pomiaru przyjmą, więc wartości:

0x01 graphic
, 0x01 graphic
, 0x01 graphic
.

Oznaczając poszukiwane wysokości punktów A i B odpowiednio przez X1, X2, możemy ułożyć następujące trzy równania wiążące nieznane poprawki 0x01 graphic
i wartości X1, X2

0x01 graphic

0x01 graphic
(w1)

0x01 graphic
.

Aby uprościć dalszy rachunek dokonujemy pewnych przekształceń. Wynikają one stąd, że zwykle znamy przybliżone wartości X10, X20 poszukiwanych rzędnych X1, X2 różniące się od nich odpowiednio o 0x01 graphic
. Będzie, więc:

0x01 graphic
, 0x01 graphic
.

(Takie postępowanie pozwala operować, co jest korzystne z obliczeniowego punktu widzenia, liczbami o małych wartościach)

Równania (w1) można, więc zapisać następująco:

0x01 graphic
;

0x01 graphic
;

0x01 graphic
.

czyli

0x01 graphic
;

0x01 graphic
;

0x01 graphic
.

Składniki ujęte w nawiasy są to przewyższenia odpowiednich ciągów obliczone z wysokości przybliżonych. Nazywamy je przewyższeniami przybliżonymi i oznaczamy przez 0x01 graphic
. Związki można zapisać, więc w postaci:

0x01 graphic
;

0x01 graphic
;

0x01 graphic
.

Przyjmując, że w naszym zadaniu: 0x01 graphic
, 0x01 graphic
i podstawiając te oraz pozostałe

( H1 = 0,000 m, H2 = 6,000 m) wartości liczbowe do powyższej zależności otrzymamy równania:

0x01 graphic
;

0x01 graphic
; (w2)

0x01 graphic
.

w których wyrazy wolne wyrażone są w milimetrach.

Po dokonaniu powyższych przekształceń zadanie nasze polega na wyznaczeniu przyrostów0x01 graphic
, które dodane do rzędnych przybliżonych X10, X20, dadzą rzędne poszukiwane X1, X2. Łatwo zauważyć, że otrzymany układ nie może dać jednoznacznego rozwiązania bez postawienia dodatkowych założeń. Zawiera on, bowiem trzy równania i pięć niewiadomych 0x01 graphic
; niewiadomych jest więcej niż równań. Dodatkowego warunku, pozwalającego rozwiązać to zagadnienie, dostarcza nam teoria błędów, według, której należy przyporządkować obserwacjom takie poprawki, aby suma ich kwadratów przyjęła wartość minimalną. Stwierdziliśmy już, że układ nie ma rozwiązania jednoznacznego. Posiada on, jednakże, rozwiązania wieloznaczne, które łatwo otrzymać z przytoczonych związków. Wystarczy w tym celu do równań podstawić dowolne liczby dla x1 oraz x2 i wyznaczyć odpowiednie wartości 0x01 graphic
. I tak, jeśli np. 0x01 graphic
, 0x01 graphic
to mamy: 0x01 graphic
. Natomiast przy 0x01 graphic
, 0x01 graphic
otrzymujemy:0x01 graphic
itd.

Stawiając warunek, aby suma kwadratów poprawek była najmniejsza, co zapisujemy w postaci

0x01 graphic

żądamy, aby ze wszystkich możliwych trójek wartości 0x01 graphic
wybrać tą trójkę, dla której suma 0x01 graphic
jest najmniejsza.

Zadanie to można, więc sformułować następująco: znaleźć minimum funkcji

0x01 graphic
(f1)

gdy zmienne v, spełniają zależności:

0x01 graphic
;

0x01 graphic
; (w3)

0x01 graphic
.

Jak wiemy funkcja f przyjmuje wartość najmniejszą w punkcie, dla którego jej pochodne cząstkowe względem zmiennych niezależnych są równe zeru. Obliczając te pochodne należy pamiętając, że poprawki v są funkcjami zmiennych x1 i x2 i wobec tego sama funkcja f jest też funkcją tych zmiennych. Wykonując różniczkowanie względem x1 i x2 otrzymamy:

0x01 graphic
, 0x01 graphic

Podstawiając tu pochodne cząstkowe obliczone ze wzorów (f1) dla funkcji f i (w3) dla poprawek otrzymamy dodatkowe związki:

0x01 graphic
.

W rezultacie otrzymujemy do rozwiązania układ pięciu równań z pięcioma niewiadomymi,

0x01 graphic
;

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic
,

który posiada już rozwiązanie jednoznaczne.

Obliczenie niewiadomych znacznie upraszcza, gdy wielkości v wzięte z pierwszych trzech równań podstawimy do dwóch równań ostatnich. W ten sposób mamy:

0x01 graphic

Rozwiązując ten układ równań z dwiema niewiadomymi uzyskujemy wartości niewiadomych: 0x01 graphic
, 0x01 graphic
, które podstawione do pierwszych trzech równań układu pozwalają obliczyć wartości poprawek: 0x01 graphic
.

Suma kwadratów tych poprawek jest najmniejszą ze wszystkich możliwych. Wynosi ona w naszym przypadku [vv] =12. Dla przyjętych powyżej wartości zmiennych niezależnych 0x01 graphic
oraz 0x01 graphic
stwierdzimy, że 0x01 graphic
, a więc

0x01 graphic
.

Mając wartości poprawek możemy łatwo obliczyć poprawione wartości obserwacji (obserwacje wyrównane) pamiętając, że rozwiązania otrzymaliśmy w milimetrach mamy:

0x01 graphic

0x01 graphic

0x01 graphic
.

2

15



Wyszukiwarka

Podobne podstrony:
S 6 Spostrzeżenia bezpośrednie, Geodezja i Kartografia, Rachunek Wyrównawczy
S 7 Równania obserwacji 3, Geodezja i Kartografia, Rachunek Wyrównawczy
S 7 Równania obserwacji 4, Geodezja i Kartografia, Rachunek Wyrównawczy
S 7 Równania obserwacji 1, Geodezja i Kartografia, Rachunek Wyrównawczy
S 7 Równania obserwacji 2, Geodezja i Kartografia, Rachunek Wyrównawczy
S 6 Wyrównanie spostrzeżeń bezpośrednich jednakowo dokładnych, Geodezja i Kartografia, Rachunek Wyró
3.Wyrównanie obserwacji bezpośrednich, Geodezja, rachunek wyrówmawczy
Opracowanie bezpośrednich wyników pomiarów, Geodezja i Kartografia, Rachunek Wyrównawczy
Z Obliczenia dla sieci kątowej, Geodezja i Kartografia, Rachunek Wyrównawczy
Sieci płaskie, Geodezja i Kartografia, Rachunek Wyrównawczy
Obliczenia na liczbach przybliżonych, Geodezja i Kartografia, Rachunek Wyrównawczy
Ćw. 1 Zastosowanie form rachunkowych Hausbrandta, Geodezja i Kartografia, Rachunek Wyrównawczy
Wyrównanie parametryczne - metoda macierzowa, Geodezja i Kartografia, Rachunek Wyrównawczy
Wagi i błędności, Geodezja i Kartografia, Rachunek Wyrównawczy
Wyrównania korelat, Geodezja i Kartografia, Rachunek Wyrównawczy
Ściaga RW, Geodezja i Kartografia, Rachunek Wyrównawczy
Równania ogólne poprawek, Geodezja i Kartografia, Rachunek Wyrównawczy

więcej podobnych podstron