• 
    

    
    

      99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

      一種改進(jìn)的水準(zhǔn)網(wǎng)條件平差算法

      2021-02-05 11:43:46張廣興
      北京測(cè)繪 2021年1期
      關(guān)鍵詞:高差算例環(huán)路

      張廣興

      (河北省地礦局第八地質(zhì)大隊(duì)測(cè)繪處, 河北 秦皇島 066000)

      0 引言

      在高程控制測(cè)量中水準(zhǔn)網(wǎng)平差至關(guān)重要,通常使用條件平差和間接平差[1-3]。間接平差利用未知參數(shù)消除觀測(cè)值之間的不符值,通過(guò)自由極值求解未知參數(shù)的最或是值,達(dá)到平差的目的。間接平差參數(shù)設(shè)定簡(jiǎn)單,誤差方程在建立過(guò)程中唯一,但未知參數(shù)個(gè)數(shù)的選定要求高,一旦選定不合適就會(huì)產(chǎn)生較大誤差[4-10]。

      條件平差是一種經(jīng)典的平差方法,19世紀(jì)到20世紀(jì)之間發(fā)展的已經(jīng)相對(duì)成熟。條件平差是將全部觀測(cè)值的平差值作為未知數(shù),利用多余觀測(cè)數(shù)列方程,通過(guò)求條件極值求解改正數(shù),進(jìn)行平差。條件平差方程式直觀簡(jiǎn)單,計(jì)算過(guò)程規(guī)律性強(qiáng),但由于水準(zhǔn)網(wǎng)形多種多樣,導(dǎo)致條件平差方程不唯一[11]。自20世紀(jì)50年代開(kāi)始,隨著計(jì)算機(jī)的不斷發(fā)展,條件平差方法也得到了很大程度的發(fā)展,研究成果相當(dāng)普遍,但是由于水準(zhǔn)網(wǎng)形復(fù)雜多樣,條件平差方程列立不唯一,改正數(shù)的解存在多樣性,在使用計(jì)算機(jī)進(jìn)行處理時(shí)非常困難,由此本文提出了一種利用傳遞矩陣建立條件平差的算法。

      1 利用傳遞矩陣建立條件方程

      1.1 傳遞數(shù)矩陣的建立

      假設(shè)某一水準(zhǔn)網(wǎng)中水準(zhǔn)路線數(shù)為n,水準(zhǔn)點(diǎn)的個(gè)數(shù)為m,則該水準(zhǔn)網(wǎng)的傳遞數(shù)矩陣T為:

      (1)

      式(1)中水準(zhǔn)路線取值為1,2,…,n,水準(zhǔn)點(diǎn)從待定點(diǎn)開(kāi)始編號(hào),全部編號(hào)結(jié)束后緊接著對(duì)已知點(diǎn)編號(hào),直到m結(jié)束。

      1.2 傳遞數(shù)矩陣的運(yùn)算規(guī)則

      通過(guò)選取的水準(zhǔn)路線的起點(diǎn)和終點(diǎn)是否已知、未知水準(zhǔn)點(diǎn)的近似高程是否已經(jīng)計(jì)算,計(jì)算傳遞數(shù)矩陣[12],運(yùn)算規(guī)則為:

      (1)若水準(zhǔn)路線j的起始水準(zhǔn)點(diǎn)i已知,終點(diǎn)水準(zhǔn)點(diǎn)k未知,且該點(diǎn)近似高程未知,則有:

      T(k,j)=1,且T(k)其余分量值為0,式中,T(k)為第k號(hào)水準(zhǔn)點(diǎn)的傳遞向量。

      (2)若水準(zhǔn)路線j的起始水準(zhǔn)點(diǎn)i未知,且該點(diǎn)近似高程未知,終點(diǎn)水準(zhǔn)點(diǎn)k已知,則有:

      T(i,j)=-1,且T(i)其余分量值為0。

      (3)若水準(zhǔn)路線j的起始水準(zhǔn)點(diǎn)i未知,且該點(diǎn)近似高程已知,終點(diǎn)水準(zhǔn)點(diǎn)k未知,且該點(diǎn)近似高程未知,則T(k)由兩部分組成,一部分為該水準(zhǔn)路線對(duì)終點(diǎn)的影響:

      T(k,j)=1,且T(k)其余分量值為0。

      另一部分為起點(diǎn)傳遞數(shù)向量傳遞給終點(diǎn)的,為:

      T(k)=T(k)+T(i)

      (2)

      (4)若水準(zhǔn)路線j的起始水準(zhǔn)點(diǎn)i未知,且該點(diǎn)近似高程未知,終點(diǎn)水準(zhǔn)點(diǎn)k未知,且該點(diǎn)近似高程未知,則T(k)由兩部分組成,一部分為該水準(zhǔn)路線對(duì)起點(diǎn)的影響:

      T(i,j)=-1,且T(i)其余分量值為0。

      另一部分為終點(diǎn)傳遞數(shù)向量傳遞給起點(diǎn)的,為:

      T(i)=T(i)+T(k)

      2.3 利用傳遞數(shù)矩陣的條件平差系數(shù)矩陣的建立

      選取必要觀測(cè)數(shù)t條水準(zhǔn)路線計(jì)算待定點(diǎn)的近似高程,則多余觀測(cè)數(shù)r為水準(zhǔn)路線的總數(shù)n與必要觀測(cè)水準(zhǔn)路線數(shù)t的差值,利用r條水準(zhǔn)路線建立條件方程。

      假設(shè)系數(shù)矩陣A為:

      (3)

      假設(shè)A(s)為第s個(gè)條件方程的系數(shù)構(gòu)成的行向量,若第j條水準(zhǔn)路線未用于計(jì)算待定點(diǎn)的近似高程,起始水準(zhǔn)點(diǎn)點(diǎn)號(hào)為i,終點(diǎn)水準(zhǔn)點(diǎn)點(diǎn)號(hào)為k,則A(s)為:

      A(s)=T(i)-T(k)

      (4)

      常數(shù)項(xiàng)W為:

      W(s)=該點(diǎn)水準(zhǔn)路線起點(diǎn)高程+高差-終點(diǎn)高程

      (5)

      2 算例分析

      2.1 算例1

      水準(zhǔn)網(wǎng)1如圖1所示。已知高程為HX=237.483 m,水準(zhǔn)路線數(shù)n為5,水準(zhǔn)點(diǎn)總個(gè)數(shù)m為4,觀測(cè)高差與路線長(zhǎng)度如表1所示。

      圖1 水準(zhǔn)網(wǎng)1

      表1 水準(zhǔn)網(wǎng)1觀測(cè)高差和路線長(zhǎng)度

      水準(zhǔn)網(wǎng)1的水準(zhǔn)路線的編號(hào)如表1所示,水準(zhǔn)點(diǎn)的編號(hào)如表2所示。

      表2 水準(zhǔn)網(wǎng)1水準(zhǔn)點(diǎn)點(diǎn)名與對(duì)應(yīng)編號(hào)

      利用L1、L3、L5水準(zhǔn)路線計(jì)算待定點(diǎn)的近似高程,即第1、3、5號(hào)水準(zhǔn)路線,相應(yīng)確定傳遞數(shù)矩陣。第1、3、5號(hào)水準(zhǔn)路線為水準(zhǔn)路線起始水準(zhǔn)點(diǎn)已知,終點(diǎn)水準(zhǔn)點(diǎn)未知,且該點(diǎn)近似高程未知,則有T(1, 4)=1,T(2, 4)=1,T(3, 4)=1,其余分量值為0,所以傳遞數(shù)矩陣為:

      (6)

      第2、4號(hào)水準(zhǔn)路線未參與傳遞數(shù)矩陣的計(jì)算,可以用來(lái)計(jì)算條件平差的系數(shù)矩陣和常數(shù)項(xiàng)矩陣。以第2條水準(zhǔn)路線為例,起點(diǎn)點(diǎn)號(hào)為1,終點(diǎn)點(diǎn)號(hào)為2,所以可以求得:

      (7)

      同理可以計(jì)算出:

      (8)

      又因?yàn)锳(s)第j個(gè)分量的值為1,所以最終求得條件平差的系數(shù)矩陣為:

      (9)

      根據(jù)常數(shù)項(xiàng)W的計(jì)算公式,以W(1)為例,2號(hào)水準(zhǔn)路線起點(diǎn)Y近似高程為237.483+5.835=238.018 m,終點(diǎn)Z的近似高程為237.483+9.640=247.123 m,可以求得:

      W(1)=238.018+3.782-247.123=-0.023 m

      (10)

      同理可以求得W(2)=0.014 m。

      所以求得常數(shù)項(xiàng)矩陣為:

      (11)

      根據(jù)條件平差過(guò)程可以求得法方程聯(lián)系數(shù)K:

      (12)

      所以改正數(shù)V為:

      (13)

      這樣就可以求得觀測(cè)高差的改正數(shù),對(duì)水準(zhǔn)網(wǎng)1的原始測(cè)量高差進(jìn)行改正后如表3所示。

      表3 水準(zhǔn)網(wǎng)1觀測(cè)高差及其改正值 單位:m

      為了比較通過(guò)傳遞數(shù)矩陣改正后高程的精度,以環(huán)路閉合差進(jìn)行評(píng)價(jià),選取了水準(zhǔn)網(wǎng)1中的三個(gè)閉合環(huán)路,計(jì)算其改正前后環(huán)路閉合差,如表4所示。

      表4 改正后環(huán)路閉合差 單位:mm

      由表4可知,閉合環(huán)XYZ、XZQ以及XYZQ的環(huán)路閉合差在利用傳遞數(shù)矩陣改正后效果明顯變好,環(huán)路閉合差都為0 mm,說(shuō)明本文所用方法有效,可以在高程改正中應(yīng)用。為了避免數(shù)據(jù)的偶然性,利用算例2進(jìn)行驗(yàn)證。

      2.2 算例2

      水準(zhǔn)網(wǎng)2如圖2所示。已知高程為HA=23.0 m,HB=23.564 m,HC=23.663 m,水準(zhǔn)路線數(shù)n為6,水準(zhǔn)點(diǎn)總個(gè)數(shù)m為5,觀測(cè)高差與路線長(zhǎng)度如表5所示。

      圖2 水準(zhǔn)網(wǎng)2

      水準(zhǔn)網(wǎng)2的水準(zhǔn)路線的編號(hào)如表5所示,水準(zhǔn)點(diǎn)的編號(hào)如表6所示。

      表5 水準(zhǔn)網(wǎng)2觀測(cè)高差和路線長(zhǎng)度

      表6 水準(zhǔn)網(wǎng)2水準(zhǔn)點(diǎn)點(diǎn)名與對(duì)應(yīng)編號(hào)

      利用第3、4兩條水準(zhǔn)路線計(jì)算待定點(diǎn)的近似高程,相應(yīng)確定傳遞數(shù)矩陣。第3、4這兩條水準(zhǔn)路線為水準(zhǔn)路線起始水準(zhǔn)點(diǎn)已知,終點(diǎn)水準(zhǔn)點(diǎn)未知,且該點(diǎn)近似高程未知,則有T(1, 4)=1,T(2, 3)=1,其余分量值為0,所以傳遞數(shù)矩陣為:

      (14)

      第1、2、5、6條水準(zhǔn)路線未參與傳遞數(shù)矩陣的計(jì)算,可以用來(lái)計(jì)算條件平差的系數(shù)矩陣和常數(shù)項(xiàng)矩陣。以第1條水準(zhǔn)路線為例,起點(diǎn)點(diǎn)號(hào)為2,終點(diǎn)點(diǎn)號(hào)為1,所以可以求得系數(shù)矩陣A為:

      (15)

      根據(jù)常數(shù)項(xiàng)W的計(jì)算公式,求得常數(shù)項(xiàng)矩陣為:

      (16)

      根據(jù)條件平差過(guò)程可以求得法方程聯(lián)系數(shù)K:

      (17)

      所以改正數(shù)V為:

      V=QATK=
      [0.001 1-0.004 450.005 5-0.000 36
      0.011 60.012 4]T

      (18)

      對(duì)水準(zhǔn)網(wǎng)2的原始測(cè)量高差進(jìn)行改正后如表7所示。

      表7 水準(zhǔn)網(wǎng)2觀測(cè)高差及其改正值

      同樣,參考算例1,采用水準(zhǔn)網(wǎng)2中的三個(gè)閉合環(huán)進(jìn)行評(píng)價(jià),改正后效果,如表8所示。

      表8 環(huán)路閉合差 單位:mm

      由表8可知,閉合環(huán)ADE、BDE以及ABD的環(huán)路閉合差在改正前分別為10 mm、-7 mm、3 mm,利用傳遞數(shù)矩陣改正后效果明顯變好,環(huán)路閉合差分別為5 mm、-4 mm、1 mm,證明了利用傳遞數(shù)矩陣的條件平差的有效性。

      3 結(jié)束語(yǔ)

      基于條件平差和間接平差本身計(jì)算過(guò)程中存在的問(wèn)題,本文提出了一種利用傳遞數(shù)矩陣建立條件方程的算法,利用傳遞數(shù)的運(yùn)算規(guī)則,解決了條件方程建立過(guò)程中多樣性的問(wèn)題,采用兩個(gè)算例,對(duì)改正后結(jié)果的環(huán)路閉合差進(jìn)行比較,算例1通過(guò)改正之后三個(gè)環(huán)路閉合差全部為0,算例2在改正后三個(gè)環(huán)路閉合差為5 mm、4 mm、1 mm,效果明顯變好,兩個(gè)算例證明了本文所用方法的有效性,在實(shí)際應(yīng)用中具有參考價(jià)值。

      猜你喜歡
      高差算例環(huán)路
      高差影響下的城鎮(zhèn)燃?xì)夤艿浪τ?jì)算簡(jiǎn)化公式
      框架結(jié)構(gòu)梁板面鋼筋疊合產(chǎn)生的高差問(wèn)題探討
      同時(shí)對(duì)向間接高差精密測(cè)量技術(shù)應(yīng)用研究
      鐵道勘察(2018年5期)2018-10-22 06:12:00
      上海市中環(huán)路標(biāo)線調(diào)整研究
      上海公路(2018年4期)2018-03-21 05:57:46
      基于振蕩能量的低頻振蕩分析與振蕩源定位(二)振蕩源定位方法與算例
      互補(bǔ)問(wèn)題算例分析
      地形高差較大的別墅區(qū)排水設(shè)計(jì)要點(diǎn)分析
      基于CYMDIST的配電網(wǎng)運(yùn)行優(yōu)化技術(shù)及算例分析
      Buck-Boost變換器的環(huán)路補(bǔ)償及仿真
      燃煤PM10湍流聚并GDE方程算法及算例分析
      新乐市| 思茅市| 榆社县| 潞西市| 高青县| 汤原县| 盈江县| 汉川市| 清水河县| 宁乡县| 理塘县| 临漳县| 吉木萨尔县| 安国市| 潍坊市| 洱源县| 镇平县| 古田县| 防城港市| 崇阳县| 明溪县| 临朐县| 双峰县| 凤阳县| 广丰县| 邵阳县| 锦州市| 林西县| 五大连池市| 精河县| 手游| 博客| 亚东县| 保山市| 雷山县| 天峨县| 新巴尔虎右旗| 尚志市| 裕民县| 东源县| 富蕴县|