橢圓周長與橢球表面積的數(shù)值計算

四川省涼山州美姑縣中學 周鈺承

西南林業(yè)大學濕地學院2018 級環(huán)境科學 周園鈔

在科學與工程中，常常選定橢圓、橢球作為數(shù)學模型進行數(shù)值計算。辛普森數(shù)值積分、泰勒級數(shù)公式、龍貝格算法，都是對橢圓周長、橢球表面積進行數(shù)值計算的好工具。按不同精度要求，文章列出了橢圓周長和橢球表面積的數(shù)值計算公式，以供數(shù)學愛好者及普通工程師使用。

一、序言：橢圓與橢球是科學與工程常用的數(shù)學模型

星際的公轉(zhuǎn)和自轉(zhuǎn)，山水間的蜿蜒和盤旋，物種間的追逐和閃避……那拐彎處的美麗弧線，正是橢圓與橢球?qū)κ篱g萬物的慷慨饋贈。于是，在科學與工程計算中，選定橢圓和橢球作為數(shù)學模型，既是自然選擇，也是藝術惠顧。

然而，橢圓函數(shù)是一類高深莫測的超越函數(shù)。本文力圖利用辛普森數(shù)值積分、泰勒級數(shù)公式及龍貝格算法，對橢圓周長和橢球表面積提煉出相對簡潔而又精確的初等公式，以供數(shù)學愛好者及普通工程師使用。

二、用辛普森公式對橢圓周長和橢球表面積快速估值

1.橢圓周長的快速估值

與其他橢圓周長公式相比，這個公式最大的特點是沒有圓周率參與運算，減少了計算量，而當橢圓趨近于圓時，圓周率自動向近似值3.122 趨近，故公式（1）的誤差約為一千分之六，比辛普森公式直接結論，相對橢圓離心率誤差分布更均勻。如果橢圓兩個半軸長滿足勾股數(shù)，公式（1）完全可以口算：當a=4，b=3 時，周長L≈22（真值為22.103…）；當a=12，b=5 時，周長L≈56（真值為55.695…）；當a=15，b=8 時，周長L≈74（真值為73.939…）。

2.橢球表面積的快速估值

旋轉(zhuǎn)橢球體表面積是可以通過被積函數(shù)求原函數(shù)得到真值的，因此，可以用旋轉(zhuǎn)橢球體表面積的真值來驗證公式（2）的精度。驗證發(fā)現(xiàn)，盡管公式（2）沒有圓周率參與運算，但比用辛普森數(shù)值積分公式，不僅形式簡單得多，而且誤差更小，約為一千分之一。當a=1，b=0.5，c=0.5 時，S=5.375（真值為5.3696…）；當a=1，b=0.9，c=0.9 時，S≈10.936（ 真 值 為10.940…）；當a=1，b=1，c=0.7 時，S≈10.155（真值為10.144…）；當a=6378，b=6378，c=6357 時，S≈5.098 億（真值為5.100…億，地球表面積）。

三、用級數(shù)公式對橢圓周長與橢球表面積進行數(shù)值計算

在環(huán)境科學與工程的科學計算中，有時需要更高精度的數(shù)值計算。筆者曾在《數(shù)理化解題研究》發(fā)表《探索橢圓周長和橢球表面積的近似初等公式》（以下簡稱《探索》），被“知網(wǎng)”收錄。公式雖然當初發(fā)表在高中版，但因為利用了貝塞爾等人派生的橢圓周長級數(shù)公式（見《探索》），所以計算式簡單易算，而且精度很高。《探索》列出的橢圓周長公式和橢球表面積公式分別如下：

《探索》對橢球表面積的精度作了一個詳盡的統(tǒng)計表，當橢球三個半軸長度不超過2 倍關系時，誤差低于一萬分之一；當三個半軸長度趨近相等時，誤差趨近于0。《探索》因為編輯部的原因去掉了橢圓周長精度統(tǒng)計表，現(xiàn)列出，以饗讀者：（表中誤差率等于用近似值減去真值差，除以真值）

表1 橢圓周長公式（3）

四、用龍貝格算法對橢圓周長進行高精度數(shù)值計算

總之，科學與工程計算在選定數(shù)學模型后，應將超越函數(shù)初等化，將解析公式數(shù)值化，以供更多的數(shù)學愛好者及普通工程師使用。辛普森數(shù)值積分公式適用范圍廣，但若要提高精度，則表達式往往復雜。級數(shù)公式明確了函數(shù)變量與哪些自變量有關，故而更容易提高精度。龍貝格算法提高了橢圓周長公式的精度，讓誤差相對于離心率分布更均勻。

美國數(shù)學家、數(shù)學教育家波利亞曾指出：“類比是一個偉大的引路人，求解立體幾何問題往往有賴于平面幾何中的類比問題。”本文橢球表面積正是用了類比和驗證來代替復雜的數(shù)學推理和證明。