以方差分析法侦台湾一等一级水准网之残留系统误差.pptVIP

以方差分析法偵測台灣一等一級水準網之殘留系統誤差 Detect Remained Significant Errors In The First-Order ClassⅠLeveling Network of Taiwan Using Analysis of Variance 指導教授：許榮欣 學生：林曾進 Introduction 水準測量易受系統誤差的影響，如地球曲率、大氣折光、地面水蒸氣、重力變化、儀器沈陷、潮汐影響、水準尺系統性誤差等，因此在測量前後必須使用特定方式來移除系統誤差，但是系統性誤差很難完全移除而留下殘餘系統誤差，以往是將殘餘的系統誤差直接當成偶然誤差計算，但系統誤差仍影響到整個網形之精度。 此研究是利用統計方法中之方差分析法應用於偵測台灣一等一級(2001)水準網，找出網形是否仍有殘餘之系統誤差。 Analysis of Variance 假設有個n條測線所組成的水準網，其中每條測線中有Ki個測段數，N為水準網中所有測段數，Sij為各線段長度（km）。 其中 ij 表第 i 條測線中第 j 個線段 令Dij為為第i測線中第 j 測段的閉合差(mm)，則測段之每公里閉合差為dij（= Dij/Sij），各測線中之閉合差其平均值為di（= Σdij/Ki），則可將水準網中之各閉和差列如表1 假設所有的測段閉合差dij為互相獨立的隨機變數呈常態分佈，則每測線平均之閉合差d1 ，d2 ，d3 ，…，dn其方差期望值分別為 ， ，…， 。當所有測線檢驗出有明顯的週期，各相位上的理論值必不相同，也就是說各隨機變量具有不同的數學期望值，否則，他們的方差期望值應該相等[陶本澡，1984]。因此我們設定一零假說進行統計檢驗： H0： （1） against H1：not all are equal 全網方差（Grand mean squares）: （2） 每條測線（Li）之方差（Mean squares within lines) : （3） 測線與測線間（di）之方差（Mean squares between lines）: （4） 其中d為平均值 F-test 以F統計檢驗進行零假說測試，其自由度為n-1, N-n： （6） 在零假說 H0 下， 應該不會特別大過於 ，因其估值均來自於同一個母體方差σ2，在p%的信心水準下： 若 ，則H0將不被接受 （7） 當滿足此條件時則我們可判定此網形仍殘留系統性誤差。 Put the discrepancies to the normal distribution test 先假設所有的閉合差均為互相獨立不相關，設定間隔Ai(i=1,2,…,k) ，紀錄閉合差落入各區間之個數ni(i=1,2,…,k)，算出各區間理論期望值出現之機率p(Ai)= pi。當n (＝n1+ n2+…+ nk) 夠大時，則 （8） 應為一近似於自由度為k-1之χ2分佈[Hsu, 2002]。 若所有閉合差是從一常態母體隨機抽樣而來，會屬於一常態分佈N(0, σ2)，將所有閉合差標準化後得W，則W~N(0,I)。標準化之原因在於各觀測值可能由不同量測方式或儀器而得，資料將會有各種種類之觀測值。以下為統計檢定由一常態分佈母體之閉合差的步驟[Hsu, 2002]： 常態分佈檢定之步驟 將所有閉合差標準化 （9） i=1,2,…,n j=1,2,....,Ki 其中σ為所有dij之標準差 將所有標準化過之閉合差分成彼此獨立之k間隔，A1,A2,...,Ak。 統計落入各區間Ai (i=1,2,…,k)之個數oi。 算出出現個數期望值ei。 在信心區間α的條件下，若Dc，c由 查表求得，則否定零假說 H0: W ~ N ( 0 , I) 將上述成果繪製成統計分佈圖，以圖形來表示資料之分佈情形並檢驗數據是否成常態分佈[Wol