楊雙富

(云南一九八煤田地質(zhì)勘探隊(duì)，云南昆明 650208)

再議計(jì)算子午線弧長(zhǎng)的數(shù)值積分法

楊雙富?

(云南一九八煤田地質(zhì)勘探隊(duì)，云南昆明 650208)

給出了利用復(fù)化辛普森(Simpson)積分公式計(jì)算子午線弧長(zhǎng)的方法，分析了計(jì)算結(jié)果精度與積分區(qū)間的大小和區(qū)間等分?jǐn)?shù)的關(guān)系。

子午線弧長(zhǎng)；數(shù)值積分；辛普森(Simpson)積公式、復(fù)化辛普森(Simpson)積公式

1 前 言

《測(cè)繪通報(bào)》2006年第5期刊登的《計(jì)算子午線弧長(zhǎng)的數(shù)值積分法》一文中介紹了子午線弧長(zhǎng)計(jì)算的數(shù)值積分方法，遺憾是的沒(méi)有明確指出其計(jì)算結(jié)果源自于數(shù)值積分方法中的哪一種，但通過(guò)與經(jīng)典方法結(jié)果相對(duì)比，得出相互矛盾的結(jié)論:“數(shù)值積分法與經(jīng)典方法的計(jì)算結(jié)果之間存在著顯著的差異，數(shù)值積分法精度高?！盵6]，值得注意是經(jīng)典方法計(jì)算結(jié)果的精度優(yōu)于1 mm，而從文中表1～表3(僅摘錄表1)可以看出，緯度從10°～90°，子午線弧長(zhǎng)的數(shù)值積分法計(jì)算結(jié)果與經(jīng)典方法相差750多米到92 680多米，竟不如僅將積分區(qū)間二等分的辛普森(Simpson)積分計(jì)算結(jié)果，真讓人匪夷所思。

圖1 《計(jì)算子午線弧長(zhǎng)的數(shù)值積分法》一文表1摘錄

眾所周知，數(shù)值積分法精度高低與所采用積分方法、積分區(qū)間的大小、區(qū)間等分?jǐn)?shù)或節(jié)點(diǎn)數(shù)多少有密切的關(guān)系，下面以克拉索夫斯基橢球體及最常用的辛普森(Simpson)積分公式、復(fù)化辛普森(Simpson)積分公式為例，對(duì)采用數(shù)值積分法計(jì)算子午線弧長(zhǎng)作一簡(jiǎn)略的分析。

2 子午線弧長(zhǎng)計(jì)算公式及數(shù)值積分公式簡(jiǎn)介

2.1 子午線弧長(zhǎng)計(jì)算公式

設(shè)a為橢球的長(zhǎng)半軸，e為第一偏心率，則在子午圈上從赤道開(kāi)始到任一緯度B的弧長(zhǎng)計(jì)算公式為[2]:

2.2 經(jīng)典計(jì)算公式

式中系數(shù)a0、a2、a4、a6、a8詳見(jiàn)文獻(xiàn)[2]。

2.3 辛普森(Simpson)積分公式

2.4 復(fù)化辛普森(Simpson)積分公式

式中n為積分區(qū)間等分?jǐn)?shù)。

3 子午線弧長(zhǎng)計(jì)算結(jié)果對(duì)比

3.1 子午線弧長(zhǎng)計(jì)算函數(shù)的設(shè)計(jì)

為方便操作及使用，以Excel 2003為計(jì)算平臺(tái)，采用自定義函數(shù)的方式設(shè)計(jì)各計(jì)算函數(shù)。啟動(dòng)Excel 2003，選擇菜單工具(T)→宏(M)→Visual Basic編輯器(V)或按Alt＋F11進(jìn)入Visual Basic編輯器，再依次點(diǎn)擊插入(I )→模塊(M)，添加一個(gè)模塊，然后輸入如下程序:

Public Const pi As Double＝3.14159265358979'定義常數(shù)pi (即π)

Public Function jdhcjs(ByVal x As Double)

其中“jdhcjs”、“Simpson”、“Simpson0”三個(gè)自定義函數(shù)分別對(duì)應(yīng)于子午弧長(zhǎng)計(jì)算的經(jīng)典公式、辛普森(Simpson)積分公式和復(fù)化辛普森(Simpson)積分公式，而“Deg”是將形如“度.分分秒秒”格式的緯度轉(zhuǎn)換為弧度的函數(shù)。

3.2 子午線弧長(zhǎng)計(jì)算對(duì)比

(1)經(jīng)典計(jì)算公式與辛普森(Simpson)積分公式計(jì)算結(jié)果對(duì)比

建立如圖2所的表格，在A列3～12行填入緯度0°～90°，分別在B3、C3、D3單元格中填入公式“＝j(luò)dhcjs (A3)”、“＝Simpson(0，A3)”、“＝A3－C3”，選中B3～D3單元格，利用填充柄向下填至12行，結(jié)果見(jiàn)圖2。其中B列3～12行為子午線弧長(zhǎng)的經(jīng)典公式計(jì)算結(jié)果，C列3～12行為子午線弧長(zhǎng)的辛普森(Simpson)積分公式計(jì)算結(jié)果，D列3～12行為前面兩計(jì)算結(jié)果之差。

圖2 計(jì)算截圖

(2)經(jīng)典計(jì)算公式與復(fù)化辛普森(Simpson)積分公式計(jì)算結(jié)果對(duì)比

依照(1)操作步驟建立如圖3所示的表格，僅將C3單元格中的公式改為“＝Simpson0(66，0，A3)”，得到的結(jié)果如圖3。

圖3 計(jì)算截圖

(3)復(fù)化辛普森(Simpson)積分公式計(jì)算精度淺析

建立如圖4所示的表格，在第2行B列～J列填入緯度10°～90°，在A列3～14行中填入積分區(qū)間等分?jǐn)?shù)n＝20～100，在B3單元格中填入公式＝(jdhcjs(B＄2)－Simpson0(＄A3，0，B＄2))?1000”，選中B3單元格，利用填充柄向右填至J列；選中B3～J3單元格，利用填充柄向下填至14行，結(jié)果見(jiàn)圖4。表中B 3～J14區(qū)域內(nèi)單元格中的數(shù)值為子午線弧長(zhǎng)的經(jīng)典公式計(jì)算結(jié)果與復(fù)化辛普森(Simpson)積分公式計(jì)算結(jié)果之差的1000倍。從表中不難看出，當(dāng)積分區(qū)間相同時(shí)，隨著積分區(qū)間等分?jǐn)?shù)的增加，計(jì)算結(jié)果的精度提高；當(dāng)積分區(qū)間等分?jǐn)?shù)相同時(shí)，在緯度10°～60°范圍內(nèi)，隨著積分區(qū)間的增大，計(jì)算結(jié)果的精度降低，而在緯度60°～90°范圍內(nèi)，計(jì)算結(jié)果的精度卻呈現(xiàn)不規(guī)則變化。

圖4 計(jì)算截圖

4 結(jié) 語(yǔ)

通過(guò)上面簡(jiǎn)略分析，與經(jīng)典方法的計(jì)算結(jié)果相比，可以看出子午線弧長(zhǎng)數(shù)值積分法計(jì)算結(jié)果的精度不僅與所用積分方法有關(guān)，而且還與積分區(qū)間的大小、區(qū)間等分?jǐn)?shù)或節(jié)點(diǎn)數(shù)多少密切相關(guān)，這一點(diǎn)在實(shí)際運(yùn)用中應(yīng)加以區(qū)別對(duì)待，如復(fù)化辛普森(Simpson)積分公式要獲得任意積分區(qū)間的計(jì)算結(jié)果優(yōu)于1 mm，其區(qū)間等分?jǐn)?shù)n至少取66。

[1] 邢永昌，張鳳舉.礦區(qū)控制測(cè)量(上冊(cè))[M].北京:煤炭工業(yè)出版社，1987

[2] 孔祥元，郭際明，劉宗泉.大地測(cè)量學(xué)基礎(chǔ)[M].武漢:武漢大學(xué)出版社，2001

[3] 徐萃薇.計(jì)算方法引論[M].北京:高等教育出版社，1987

[4] 翟瑞彩，謝偉松.數(shù)值分析[M].天津:天津大學(xué)出版社，2000

[5] 數(shù)學(xué)手冊(cè)編寫(xiě)組.數(shù)學(xué)手冊(cè)[M].北京:人民教育出版社，1979

[6] 劉修善.計(jì)算子午線弧長(zhǎng)的數(shù)值積分法[J].測(cè)繪通報(bào)，2006(5):4～6

Reconsideration of the Numerical Calculation of Radial arc Length Integral Method

Yang ShuangFu
(198 coal geological exploration teams in Yunnan，Yunnan 650208，China)

This paper presents the use of complex－based Simpson(Simpson)integral formula meridian arc length method to analyze the accuracy of calculation results with the integral interval size and range of the relationship between such scores.

radial arc length；numerical integration；Simpson(Simpson)product formulas，complex of Simpson (Simpson)product formula

1672－8262(2010)06－140－03

P226.1

B

2010—02—05

楊雙富(1967—)，男，高級(jí)工程師，主要從事工程測(cè)量工作。