如何用EXCEL進行附和導(dǎo)線近似平差及計算中若干問題的解決方法
地測隊:楊沖
2011年11月
如何用EXCEL進行附和導(dǎo)線近似平差及
計算中若干問題的解決方法
在煤礦測量工作中,我們經(jīng)常需要進行單一附和導(dǎo)線近似平差。傳統(tǒng)的手工計算,不但效率低,而且算錯之后不易改正(因為后一步的起算數(shù)據(jù)由前一步計算得到,一旦算錯,后續(xù)步驟都需重新計算),運用EXCEL強大的函數(shù)計算功能,就可以解決這個難題。
用EXCEL進行導(dǎo)線近似平差的流程大致為:a.輸入各站水平角、豎直角、起終端坐標方位角、起終端坐標及高程、各站斜距、儀器高、棱鏡高。b.計算坐標方位角及其閉合差,并進行閉合差分配。c.計算各導(dǎo)線點平距、坐標和坐標閉合差。如閉合差不超限,分配坐標閉合差。d.計算導(dǎo)線點高程和高差閉合差,如不超限,分配高差閉合差。(具體的計算方法及格式見附表算例)
但是,在計算過程中,我們會遇到諸如角度計算、數(shù)值取位等問題。
首先是角度計算問題。我們知道,在EXCEL中,不能對角度值直接進行三角函數(shù)運算,需要將角度值換算成弧度。這里有兩個問題,第一是角度輸入問題。在EXCEL表格中,不能直接輸入度、分、秒。這里,筆者提供一種輸入角度的方法:首先自定義單元格格式,將其定義為[h]"°"mm"′"ss"〃"的形式。角度值輸入時,度、分、秒間以冒號隔開,例如“23:30:01”,輸入完成后將自動變成“23°30′01〃”,這樣輸入的角度可直接進行加減運算;第二是角度的三角函數(shù)運算問題,需運用角度轉(zhuǎn)弧度函數(shù)“RADIANS()”。這里需要注意的是,先前輸入的角度值需要乘以24,例如:COS(RADIANS(A1*24))(假設(shè)輸入的角度值位置在A1上)。經(jīng)過以上兩個步驟,便可將斜距、水平角輕松轉(zhuǎn)換為平距、坐標方位角。
再談數(shù)值取位問題,我們知道,在角度計算中,遵循“四舍六入,五前奇進偶不進”原則。在坐標增量(Δx=S×cosα、Δy=S×sinα)計算時,我們也會遇到取位問題。如果一律按四舍五入方法計算,勢必造成坐標增量偏大。故坐標計算時,也應(yīng)遵循角度計算的取值方法,以避免因取值方法不當造成的誤差累積。我們可以采用下面的方法實現(xiàn)這一過程。
需要用到EXCEL中的幾個函數(shù):取絕對值函數(shù)ABS、嵌套條件函數(shù)IF、取小數(shù)位數(shù)函數(shù)ROUND、取余數(shù)函數(shù)MOD。以中間計算過程保留四位小數(shù)為例。需先將要計算的數(shù)值取絕對值,最后帶上正負號。例如:ABS(A1)/A1*K。計算中需要檢測小數(shù)點后第四位和第五位的數(shù)值(A1為要保留四位有效數(shù)字的數(shù)據(jù))。過程為:a、如果MOD(A1*10000,1) <0.5,則直接保留四位有效數(shù)字輸出結(jié)果,即:ROUND(A1,4);b、否則,如果MOD(A1*10000,1)<0.6,需要檢測小數(shù)點后第四位(保留四位小數(shù)后)的數(shù)值,如果是偶數(shù),則輸出保留四位小數(shù)后的數(shù)值ROUND(A1,4),否則,輸出ROUND(A1,4)-0.0001。c、如果MOD(A1*10000,1)>=0.6則輸出ROUND(A1,4)。具體實現(xiàn)算法為:
“=A1/ABS(A1)*IF(MOD(ABS(A1)*10000,1)<0.5,ROUND(ABS(A1),4),(IF(MOD(ABS(A1)*10000,1)<0.6,(IF(MOD(ROUND(ABS(A1),4)*10000,2)=0,ROUND(ABS(A1),4),ROUND(ABS(A1),4)-0.0001)),ROUND(ABS(A1),4))))”
解決了上述兩個難題,再運用簡單的三角函數(shù)及加減運算,便可輕松在EXCEL上進行導(dǎo)線近似平差。當然,上述方法仍有一些瑕疵,如:負值角度在EXCEL中不能正確顯示,雖然并不影響計算結(jié)果,但還是希望大家在以后的研究中能夠進一步改進;其次,由于不能在一個單元格中進行復(fù)雜的角度運算,表格中存在部分冗余數(shù)據(jù)(如90°00′00〃、180°00′00〃),使得表格整體感覺龐大;再次,由于存在取位(近似值)問題,計算結(jié)果與理論值仍存在毫米級偏差,當然以目前礦山控制測量要求的精度來說,完全可以忽略不計。
地測隊:楊沖
2011年11月