誤差理論與測量平差基礎(chǔ)間接平差PPT課件

1、 7-1 間接平差原理1、函數(shù)模型間接平差的函數(shù)模型就是誤差方程，其一般形式為式中：且nntnnnndLdLdLlxxxxtbatbatbaBvvvV22112122211121,111nttnnlxBVtBrk)(第1頁/共48頁2、隨機模型間接平差的隨機模型與條件平差的隨機模型相同，即3、基礎(chǔ)方程及其解 誤差方程的個數(shù)為觀測值的個數(shù)n，而未知數(shù)的個數(shù)為n+t n。所以誤差方程有無窮組解。而滿足 解只有一組。由于向量V是向量 的函數(shù)，按數(shù)學上求自由極值的方法有：nnnnLLnnLLPQD12020minPVVTx 022PBVxVPVxVVPVVxPVVTTTT第2頁/共48頁轉(zhuǎn)置后得：將此

2、式與誤差方程聯(lián)立，得間接平差的基礎(chǔ)方程為：基礎(chǔ)方程的個數(shù)與未知數(shù)的個數(shù)相等，故有唯一解。為解此基礎(chǔ)方程，將第二式代入第一式，消去V，得因為 ，所以上式有唯一解。令則0PVBT0PlBxPBBTTlxBVPVBT0tBrkPBBrkT)()(PlBWPBBNTTBB,WNxBB1第3頁/共48頁由上式解出參數(shù) 后，代入誤差方程可得到改正數(shù)V。進而可求得觀測值的平差值：間接平差的計算步驟1、根據(jù)平差問題的性質(zhì)，選擇 t 個獨立量作為參數(shù)；2、列出誤差方程；3、組成法方程；4、解算法方程；5、計算改正數(shù)V；6、計算觀測值的平差值x VLLVLL教材：71，72習題：7.1.04，7.1.05第4頁/

3、共48頁7-2 誤差方程 間接平差的關(guān)鍵是列誤差方程，而列誤差方程的關(guān)鍵是選擇待估參數(shù)（未知數(shù)）。1、未知數(shù)的個數(shù) 在間接平差中，未知數(shù)的個數(shù)等于必要觀測數(shù)t。2、未知數(shù)的選擇 選擇原則：a、所選取t個待估參數(shù)必須相互獨立； b、所選取t個待估參數(shù)與觀測值的函數(shù) 關(guān)系容易寫出來。第5頁/共48頁3、不同情況下未知數(shù)的選擇及誤差方程的列立（1）、水準網(wǎng) 在水準網(wǎng)平差中，通常選t個待定點的高程平差值作為待估參數(shù)。這樣選 既足數(shù)，又獨立， 而且容易寫出參數(shù) 與觀測值之間的函 數(shù)關(guān)系。如圖，選 FEHXHX,21第6頁/共48頁于是有：DCBAHhXvHhXvhXXvHhXvHhXv525424321

4、3212111，則有令iiixXX05254243213212111lxvlxvlxxvlxvlxv0255024402013301220111,XHhlXHhlXXhlXHhlXHhlDCBA其中：第7頁/共48頁例：水準網(wǎng)如右圖所示，已知 =5.000m， =3.953m， =7.650m。各點的近似高程為：觀測值見下表，試列出誤差方程。(m)ih(m) 12345670.0501.1002.3980.2001.0003.4043.452AHBHCHm450. 7m452. 8m053. 5407020321hHHhHHhHHCABppp第8頁/共48頁解：設(shè)于是誤差方程為：0520103

5、27216325343131211xvxxvxxvxvxxvxvxv321,321pppHXHXHX00ipiHX習題：7.2.15第9頁/共48頁（2）、GPS網(wǎng)三維無約束平差 在GPS網(wǎng)三維無約束平差中，常常選某點 i作為參考點，則該點在WGS84系下的三維坐標 、 、 可看作已知數(shù)據(jù)，其余各點作為待定點。在WGS84系下，要確定一個點的空間位置，需要X、Y、Z三個坐標分量，設(shè)GPS網(wǎng)中的總點數(shù)為m個，則必要觀測數(shù)為 ，因此，可選 個點的坐標平差值作為參數(shù)。 如圖，以A點為參考點，即 已知，則t個參數(shù)為：iXiYiZ) 1(3mt1mAAAZYX,DDDCCCBBBZXYXXXZXYXXX

6、ZXYXXX,987654321第10頁/共48頁于是，誤差方程為：CDZADAZCDYADAYCDXADAXBDZACAZBDYACAYBDXACAXBCZABAZBCYABAYBCXABAXZXXvZZXvYXXvYYXvXXXvXXXvZXXvZZXvYXXvYYXvXXXvXXXvZXXvZZXvYXXvYYXvXXXvXXXvCDADCDADCDADBDACBDACBDACBCABBCABBCAB969858747936825714633522411第11頁/共48頁（3）、三角網(wǎng) 在三角網(wǎng)平差中，通常選m個待定點的坐標平差值作為待估參數(shù)，即t=2m 。 這樣選，既足數(shù)，又獨立， 而

7、且容易寫出參數(shù)與觀測 值之間的函數(shù)關(guān)系。一般 地，角度觀測值可由右圖 表示，于是有： LikjhijhjhjkjkjhjkiLXXYYXXYYvarctanarctan教材：74習題：7.2.14第12頁/共48頁 例如右圖所示的大地 四邊形，其必要觀測 數(shù)為4，圖中待定點坐 標也是4，故選：DDCCYXXXYXXX,4321第13頁/共48頁于是，誤差方程為：834871234761212651213245431423443343432341221121arctanarctanarctanarctanarctanarctanarctanarctanarctanarctanarctanarct

8、anarctanarctanarctanarctanLXXYXXXYYvLXXYXXXYXvLXXXYXXXYvLXXXYXXXXvLXXXXXXXYvLXXXYXXXYvLXXYXXXYXvLXXYXXXYYvBBBABABDBABBBBBDBCBBAACBCAAACACDBBDCDBBBAADBDAAAAAADACAAABABACAB第14頁/共48頁（4）、三邊網(wǎng) 有足夠起算數(shù)據(jù)的三邊網(wǎng)與三角網(wǎng)一樣，也是選m個待定點的坐標平差值作為待估參數(shù)，即t=2m 。一般地，邊長觀測值可由下圖表示，于是有：SijkijkjkiSYYXXv22)()(第15頁/共48頁例如在下圖，我們選DDCCYXX

9、XYXXX,4321第16頁/共48頁 于是，誤差方程為：52221542221432422313224232124231)()()()()()()()()()(LXYXXvLXYXXvLXXXXvLXYXXvLXYXXvBBAABBAA教材：75習題：7.2.16第17頁/共48頁（5）、導(dǎo)線網(wǎng) 導(dǎo)線網(wǎng)為特殊的邊角網(wǎng)，其必要觀測數(shù)t=2m（m為待定點個數(shù)），其觀測值為角度觀測值和邊長觀測值兩類。所以誤差方程也是角度誤差方程和邊長誤差方程兩類?？梢韵攘薪嵌日`差方程： 再列邊長誤差方程。iijijikikijikiLXXYYXXYYvarctanarctanijijiiLYYXXv22)()(第

10、18頁/共48頁（6）、擬合模型 a、曲線擬合 如圖，觀測了很多散點，要求將其擬合成一條曲線。設(shè)此曲線為：332210iiiixaxaxaay第19頁/共48頁 由于觀測值 y 有誤差，故由上式可得曲線擬合的誤差方程為： b、曲面擬合 曲面擬合在DEM、GPS水準等工作中常常用到。將地面視為一個連續(xù)的曲面，則高程可表達為平面坐標的函數(shù)，且可用多項式表達為： 由于觀測值H有誤差，故由上式可得曲面擬合的誤差方程為：iiiiiyxaxaxaav332210iiiiiiiyxayaxayaxaaH52423210iiiiiiiiHyxayaxayaxaav52423210第20頁/共48頁 c、標準曲

11、線擬合 對于標準曲線，由于其方程已知，其擬合方法有所不同。如圖所示，測得m個點的坐標，要求擬合圓曲線。由于圓曲線的參數(shù)方程為：式中： 為圓心坐標，R為半徑，這三個參數(shù)是圓的基本參數(shù)， 為第i點矢徑的方位角。所以確定一條圓曲線 的必要觀測數(shù)為t=3+m。在圓周上觀測了n=2m個點的坐標，則r=m-3（ ）。于是誤差方程為： ),(00yxiiiiRYYRXXsincos00iiiyiixyRYvxRXviisincos003m第21頁/共48頁（7）、坐標變換 不論是GPS， 還是GIS，還是 RS，都會經(jīng)常用 到坐標變換。測 量中的坐標變換， 一般采用如圖所 示的相似變換。efigyyqhik

12、xxii00第22頁/共48頁 由于兩坐標系不是用同一個長度基準定義的，所以長度基準不一定嚴格相等，即兩坐標系的單位長度之比可能為： 于是坐標系 中的長度變換到坐標系 中時應(yīng)乘以尺度比m。于是： 式中， 為待定參數(shù)。由于坐標觀測值有誤差，于是坐標變換的誤差方程可寫為：1mSSyoxxoysincossincos00iiiiiixmymyyymxmxx,00myxiiiyiiixyxmymyvxymxmxviisincossincos00第23頁/共48頁上圖表示某相機在空中對地面進行攝影，uvw為相機本身的坐標系，XYZ為地面坐標系，S為相機鏡頭中心，SM為相機成像面中心與鏡頭中心連線。成像數(shù)

13、學關(guān)系可以用下式表示：第24頁/共48頁)()()()()()()()()()()()(333222333111ssssssssssssZZcYYbXXaZZcYYbXXafyZZcYYbXXaZZcYYbXXafx攝影測量的核心共線方程coscoscossincossinsinsinsincoscossinsincoscossincoscossincossinsinsincossinsinsincoscos321321321cccbbbaaaR標在影像上的對應(yīng)像點坐），為地面某點（為相機的焦距，的方向，向量下的位置和在地面坐標系為攝影中心其中ZYXyxfSMXYZZYXSSS,S,第25頁/

14、共48頁相對定向相對定向第26頁/共48頁XYZa1(x1,y1)x1y1z1S1A(X,Y,Z)a2(x2,y2)z2y2x2S2前方交會前方交會第27頁/共48頁影像點坐標地面點坐標x(mm)y(mm)X(m)Y(m)Z(m)1-86.150-68.99036589.4125273.322195.172-53.40082.21037631.0831324.51728.693-14.780-76.63039100.9724934.982386.50410.46064.43040426.5430319.81757.31已知4對點的影像坐標和地面坐標：f=153.24mm，試編程求參數(shù) 的最小二

15、乘解,SSSZYXXs，Ys近似值取地面點重心，Zs近似值取7500m,近似值均取零第28頁/共48頁7-3 非線性誤差方程的線性化 由以上所列誤差方程知，角度觀測值的誤差方程： 邊長觀測值的誤差方程： 圓曲線的誤差方程： 以及坐標變換的誤差方程都是非線性誤差方程。都必須線性化。下面介紹線性化的方法。iiiyiiixyxmymyvxymxmxviisincossincos00iiyiixyRYvxRXviisincos00ijhjhjkjkjhjkiLXXYYXXYYvarctanarctanijkjkiSYYXXv22)()(第29頁/共48頁1、變量代換法對于坐標變換的誤差方程：令則有：上

16、式即為坐標變換的線性誤差方程。iiiyiiixyxmymyvxymxmxviisincossincos00sin, cos,00mdmcybxaiiiyiiixydxcybvxdycxavii第30頁/共48頁2、線性近似 對于角度觀測值的誤差方程、邊長觀測值的誤差方程和圓曲線的誤差方程一般都是采用線性近似的方法線性化。角度觀測值的誤差方程：令：將hhhkkkjjjhhhkkkjjjyYYyYYyYYxXXxXXxXX,000000ijhjhjkjkjhjkiLXXYYXXYYvarctanarctan第31頁/共48頁在按臺勞級數(shù)展開，取至一次項，得式中：ihjhjhhjhjhkjkjkkj

17、kjkjjhjhjkjkjjhjhjkjkilySXxSYySXxSYySXSXxSYSYv )()()()()()()()(20020020020020020020020000000000arctanarctanjhjhjkjkiiXXYYXXYYLl第32頁/共48頁注意：上式是相對與右圖中三點均為代定點導(dǎo)出的。1、當圖中j點為已知點時，由于已知點的改正數(shù)為零，即于是，誤差方程變?yōu)椋?、當h、k兩點為已知點時，由于Likjh0,0jjyxihjhjhhjhjhkjkjkkjkjkilySXxSYySXxSYv )()()()(2002002002000hhkkyxyx第33頁/共48頁則誤

18、差方程變?yōu)椋?、當h或k點為已知點時，誤差方程變?yōu)椋篿jjhjhjkjkjjhjhjkjkilySXSXxSYSYv )()()()(200200200200ikjkjkkjkjkjjhjhjkjkjjhjhjkjkilySXxSYySXSXxSYSYv )()()()()()(200200200200200200第34頁/共48頁或邊長觀測值的誤差方程：令：將ihjhjhhjhjhjjhjhjkjkjjhjhjkjkilySXxSYySXSXxSYSYv )()()()()()(200200200200200200kkkjjjkkkjjjyYYyYYxXXxXX,0000ijkjkiSYYX

19、Xv22)()(第35頁/共48頁按臺勞級數(shù)展開，取至一次項，得式中：注意：1、若j點為已知點，則上式變?yōu)椋?、若k點為已知點，則：20020000)()(,jkjkjkjkiiYYXXSSSlikjkjkkjkjkjjkjkjjkjkilySYxSXySYxSXv00000000ikjkjkkjkjkilySYxSXv0000ikjjkjkjjkjkilySYxSXv0000第36頁/共48頁圓曲線的誤差方程令將按臺勞級數(shù)展開，取至一次項，得式中iiyiixyRYvxRXviisincos000000000000000,iirRRyYYxXXiiiiyiiyxiixlryvlRrxv 000

20、0000cossinsincos000000000000sin,cos,iiiiyiiiixRYYYYlRXXXXlii第37頁/共48頁例：水準網(wǎng)如下圖，觀測高差和路線長度為 ：已知點高程分別為：用間接平差求 、點高程平差值。標準答案: kmSSkmSSmhmhmhmh2,1505. 0,503. 0,501. 0,003. 132414321mHmHmHCBA008.12500.11000.111P2PmHmH5083.12,0047.1221習題：7.2.24， 7.2.25第38頁/共48頁7-4 精度評定1、單位權(quán)方差的估值2、 的計算直接計算：用常數(shù)項計算： PVVT2222211

21、nnTvpvpvpPVVxWPlllxBPlPVlPVBxPVlxBPVVTTTTTTTT)()(rPVVT20第39頁/共48頁3、基本向量的協(xié)因數(shù)矩陣間接平差中，基本向量為觀測向量L，參數(shù)向量 ，改正數(shù)向量V和觀測值的平差值向量 。令 =XLLVXLZ0111111LPBBNEPBBNPBNELPBBNEPBBNPBNETbbTbbTbbTbbTbbTbb第40頁/共48頁由協(xié)因數(shù)傳播律得：TbbTbbbbLLTbbTbbTbbLLLVLXLLLVVVVXLVLXVXXXXLLLLVXLLLZZBPBNEBPBNPBNEQPBBNEPBBNPBNEQQQQQQQQQQQQQQQQQ111111第41頁/共48頁展開得：于是：TbbbbTbbTbbLLLLTbbTbbbbTbbTbbLLTbbbbLLZZBBNBNBBNBBNQQBBNBNNBNBBNQBBNBNQQ111111111110000ttttttXXXXXXXXXXXXXXXXXXbbXXQQQQQQQQQNQ1212222112111第42頁/共48頁4、