明安圖和他的冪級數展開式

明安圖和他的冪級數展開式

羅見今

內容提要: 簡要介紹清代蒙古族科學家明安圖在天文、 地理測繪和數學方面的 貢獻。 重新表述明安圖在他的遺著 hh割圓密率捷法ii 中的三角函數冪級數公式, 特 別是係數中含卡塔蘭數 (Catalan numbers) 的無窮級數展開式, 指出在世界數學 史上明安圖是卡塔蘭數的首創者; 以一無窮級數 y10 自乘為例, 說明他在中國數學史 上奠定了無窮級數運算的基礎。 明氏的成就引起諸多國內外史學家和數學家的興趣, 例如當代電腦科學先驅者、 斯坦福大學的克努特 (D. E. Knuth) 教授的專著中、 英 國德爾比大學拉坎布 (P. J. Larcombe) 博士的多篇論文中均有論述。 國際天文學 聯合會 2002 年 5 月 26 日宣佈將新發現的 28242 號小行星命名為 “明安圖星”, 政府 也將他家鄉的城鎮命名為 “明安圖鎮”。

關鍵詞: 中國數學史, 無窮冪級數, 卡塔蘭數, 明安圖星。

1 、 明安圖的生活與工作

明安圖 (1692?−1763?), 蒙古族, “明安圖, 字靜庵, 蒙 古正白旗人。”¹ 他是中國著名的天文學家、 數學家、 測繪學 家。

根據康熙帝 1670 年的旨意: “天文關係重大, 必選擇 得人, 令其專心學習, 方能通曉精徹”, 禮部規定於官學生內

“每旗選十名, 交欽天監分科學習, 有精通者, 俟滿漢博士缺 補用²。” 1710年前明安圖被選入欽天監, 開始學習自然科學。

康熙帝喜歡學習自然科學和數學, 他親自教導學生, 是

明安圖的數學教師。 明安圖學習努力, 受到皇帝喜愛。 明安圖塑像

1見明安圖_,清史稿₄₆冊₅₀₆卷_,中華書局, 1977, 13961-13964。

2見東華錄·康熙十。

65

1.1. 明安圖是中國當時一流天文學家

1713 年明安圖完成學業, 就職於欽天監, 後來他成為當時國內一流天文學家。 他共參加了 三項重要天文書的編寫工作:

1、 hh律曆淵源ii 100卷 (1713−1723) 2、 hh曆象考成後編ii 10卷 (1737−1742) 3、 hh儀象考成ii 32卷 (1744−1752)

在 18 世紀中期他成為天文學的領導人, 任欽天監監正 (皇家天文臺臺長)。 他所接受的天 文學思想屬於第谷 (Tycho Brahe, 1546−1601) 體系。 現今保存在呼和浩特市五塔寺的石刻 蒙古文天文圖是那個時代的成果, 可能出自明安圖之手。

1.2. 明安圖完成了中國大地測量

1708−1716年清政府組織過一次中國國土測量, 未完成。 1755−1756年明安圖到新疆省 參加測量經緯度的工作。 1759年作為領導人率隊再赴新疆和塔什干等地測量經緯度。 根據幾次 測量的結果繪成中國全圖, 它成為後來中國大地圖的基礎。

hh割圓密率捷法ii 扉頁

hh割圓密率捷法ii 3卷 3頁

1.3. 明安圖對數學的貢獻

中國傳統數學到明代日漸式微, 大量古算書失傳。 與此 同時, 西方文藝復興, 出現了許多著名數學家, 中西數學間的 差距越來越大。 雍正之後形成了閉關自守的環境, 明安圖在 不能得到新的外部資訊的條件下付出巨大努力, 推進中國傳 統數學前進。

明安圖在 1730−1763年著書 hh割圓密率捷法ii³, 但生 前未能完稿, 臨終時囑咐他的兒子明新和學生陳繼新寫成此 書。 1774年完稿, 抄本曾經在數學家中流傳開, 幾經周折, 直 到 1839年才正式出版。 在清代他的數學工作產生了一百多年 的影響, 形成了 “明安圖學派”。 其主要數學貢獻是:

1、 他是中國數學史上的第一人, 計算並獲得了若干無窮級數。

2、 他是世界數學史上的第一人, 首先提出並應用了卡塔蘭數。

3、 他奠定了互反函數的理論和演算法, 並給出四組互反公式。

3明安圖_,割圜密率捷法_,道光己亥₍₁₈₃₉₎孟秋_,石樑岑氏校刊本。

4、 提出 “形數相生” 的理論, “堪與笛卡兒軔解析幾何媲美”(李儼) 5、 提出曲線和直線在無窮分割時可達到同一的極限理論。

6、 在中國首次將三角函數的無窮級數公式應用於天文計算。

本文以下討論上面第一 、 二項成就。

2 、 明安圖的冪級數展開式

18 世紀早期法國耶穌會傳教士杜德美 (P. Jartoux, 1668−1720) 到中國, 帶來三個級數 展開式⁴:

格列高裏 (J. Gregory) 1667年正弦、 正矢展開式:

r sina

r = a − a³

3!r² + a⁵

5!r⁴ − a⁷

7!r⁶ + · · · (1) rV ersa

r = a²

2!r − a⁴

4!r³ + a⁶

6!r⁵ − · · · (2) 牛頓 (I. Newton) 1676年 π 的展開式:

π = 3

1 + 1²

4 · 3! +1²· 3²

4²· 5! +1²· 3²· 5²

4³ · 7! + · · · ·

(3) 使用這三式計算的結果很精確, 中國傳統數學裏沒有這樣的公式, 但杜德美並未將三式的 證明帶到中國, 這引起了青年明安圖的懷疑, 認為西方人有 “金針不度之疑” 。 他決心做出這三 個公式的證明, 並獲得了若干三角函數其他的無窮級數展開式。

2.1. 明安圖獲得的 6個無窮級數展開式

令半徑為 r, 弧為 a, 圓心角為 α = a/r, 2a 的弦為 c, 2a 和 c 的中矢為 b, 明安圖獲得⁵ (現代形式):

c = 2r sin α = 2a − (2a)³

4 · 3!r² + (2a)⁵

4²· 5!r⁴ − (2a)⁷

4³· 7!r⁶ + · · · (4) 或 c = 2r sin α =

∞

X

n=0

(−1)ⁿ(2a)²ⁿ⁺¹ 4ⁿ(2n + 1)!r²ⁿ b = rV ersα = (2a)²

4 · 2!r − (2a)⁴

4²· 4!r³ + (2a)⁶

4³· 6!r⁵ − (2a)⁸

4⁴· 8!r⁷ + · · · (5)

4錢寶琮主編_,中國數學史_,科學出版社, 1964, 301。

5羅見今_,明安圖公式辨正_,內蒙古師大學報₍自然科學版_{), 1988,}第₁期_{, 42-48}。

或 b = rV ersα =

∞

X

n=1

(−1)ⁿ⁺¹(2a)²ⁿ 4ⁿ(2n)!r²ⁿ⁻¹ 2a = c + 1²· c³

4 · 3!r² + 1²· 3²· c⁵

4²· 5!r⁴ −1²· 3²· 5²· c⁷

4³· 7!r⁶ + · · · (6) 或 2a =

∞

X

n=0

[(2n − 1)!!]²c²ⁿ⁺¹ 4ⁿ(2n + 1)!r²ⁿ⁻¹ a = r sin α +1²(r sin α)³

3!r² +1²· 3²(r sin α)⁵

5!r⁴ + · · · (7) 或 a =

∞

X

n=0

[(2n − 1)!!]²(r sin α)²ⁿ⁺¹ (2n + 1)!r²ⁿ

a² = 2r²V ersα + 2 · 1²(2rV ersα)²

4! + 2 · 1²· 2²(2rV ersα)³

6!r + · · · (8) 或 a² =

∞

X

n=0

(n!)²(2rV ersα)ⁿ⁺¹ (2n + 2)!rⁿ⁻¹ =

∞

X

n=0

(n!)²(2b)ⁿ⁺¹ (2n + 2)!rⁿ⁻¹ (2a)² = (8b)r +2 · 1²(8b)²

4 · 4! + 2 · 1²· 2²(8b)³

4²· 6!r + · · · (9) 或 (2a)² = 2

∞

X

n=0

(n!)²(8rV ersα)ⁿ⁺¹ 4ⁿ(2n + 2)!rⁿ⁻¹ = 2

∞

X

n=0

(n!)²(8b)ⁿ⁺¹ 4ⁿ(2n + 2)!rⁿ⁻¹

上述 9個公式在中國數學史裏很有名, 在歷史上曾被不恰當地稱做 “杜氏九術”; 當代不少 數學史書中譯成現代形式的九術也出現了較多錯誤⁶, 這裏的寫法是正確的。

2.2. 無窮冪級數展開式中的卡塔蘭數

明安圖還獲得另外的重要的幾個無窮級數展開式⁶:

sin α 2

2

=

∞

X

n=1

Cⁿsin α 2

2n

(10)

sin 2α = 2 sin α −

∞

X

n=1

Cⁿ(sin α)²ⁿ⁺¹/4ⁿ⁻¹ (11)

sin 4α = 4 sin α − 10(sin α)³+

∞

X

n=1

(16Cⁿ− 2Cⁿ⁺¹)(sin α)²ⁿ⁺³/4ⁿ (12)

6明安圖原著_,羅見今譯注_,割圓密率捷法譯注_,內蒙古教育出版社_{, 1998}。

他還得到 sin 10α, sin 100α, sin 1000α, sin 10000α 的無窮級數。 英國數學家 P. J. Lar- combe 博士 2000 年將其推廣到 sin 2kα 的情況⁷。

在式中 Cⁿ 是卡塔蘭數 1, 1, 2, 5, 14, 42, 132, 429, . . . 。 卡塔蘭數在今天離散數學 (組 合數學、 圖論和數論) 中是一種著名的計數函數, 得到了廣泛的應用。 現在我們知道, 卡塔蘭數 有 6∼7個定義公式, 具有 50種組合背景的含義, 現代已發表了至少 600種論文或著作, 與著 名的 Fⁿ 斐波那契數 1, 1, 2, 3, 5, 8, 13, 21, 34, . . . 堪相媲美。

卡塔蘭 (E. Catalan, 1814−1894, 比利時) 1838年發表的一篇論文⁸ 中講到這種計數函 數, 後來該數就以他的名字命名。 後來發現大數學家歐拉 (L. Euler, 1707−1783) 在 1758年 曾經研究過它⁹。

卡塔蘭數的公式是

C₀ = 1, C₁ = 1, Cⁿ= 1 n + 1

2n n



(n ≥ 0) 卡塔蘭數的卷積遞迴公式是

C1 = 1, C2 = 1, Cⁿ=

n−1

X

k−1

Cn−kC^k (n ≥ 2)

牛頓二項式定理當指數為 1/2 時可得 (1 + z)¹² = 1 +

∞

X

n=1

(−1)ⁿ⁻¹Cⁿzⁿ

2²ⁿ⁻¹ (|z| < 1) 1839 年畢奈特 (J. Binet) 給出了卡塔蘭數的生成函數¹⁰

1 2 −1

2

√1 − 4z =

∞

X

n=1

Cⁿzⁿ (|z| < 1/4)

一些文獻¹¹介紹了它的發展史¹²。 但在這些文獻中, 西方學者不瞭解明安圖領先世界的工 作。 1730年代明安圖已遇到並應用了 Cn 數, 在他的著作 hh割圓密率捷法ii 卷三用三種方法都

7 Larcombe P. J., On Catalan numbers and Expanding the Sine Function, Mathematics Today, Bulletin of the Institute of Combinatorics and its Application(ICA), Vol. 28(2000), 39-47.

8Catalan E., Note sur une equation aux differences, J. Math. Pures Appl. 1838, 3(1), 508-516.

9Euler L., Novi Commentarii Acade Scientiarum Imperialis Petropolitanae, 7(1758-59), 13-14.

10Binet J., Reflexi ons sur le probleme de determiner le nomblre de manieres dontune figure rectiligne peut etre partagee en triangles au moyen de ses diagonals, J. Math. Pures Appl. 1839, 4, 79-91.

11Brown W. G., Historical note on a recurrent combinatorial problem, Ame. Math. Month. 72(1965), 973-977.

12Alter R., Some remarks and results on Catalan numbers, Proc. 2^nd Louisiana conf. on combinatorics, garph theory and computing(1971), 109-132.

獲得了卡塔蘭數列, 其中有兩種公式是 “新的” , 亦即現代數學界未知的, 1988年由筆者按照明 安圖原著的演算法表示為現代數學形式¹³:

公式A: 卡塔蘭數的遞推公式

C1 = 1, C2 = 1, Cⁿ+1 =X

k≥0

(−1)^kn − k k + 1



Cn−k (13)

公式B: 卡塔蘭數的有限多項式生成函數 M₁ = (1), M₂ = (0, 1),

M3 = (2M1+ M2)M2 = [2(1) + (0, 1)](0, 1)

= [(2) + (0, 1)](0, 1) = (2, 1)(0, 1) = (0, 0, 2, 1)

· · · · Mⁿ+1 =

 2

n−1

X

k=1

M^k+ Mⁿ

 Mⁿ

MCⁿ =

n

X

k=1

M^k = (1, 1, 2, 5, 14, 42, 132, 429, . . . , Cⁿ) (14)

這是一種表述方式複雜而特異的新函數, 1999年拉坎布給出它的一個證明¹⁴。

我們認為, 明安圖獲得的無窮級數中, 其係數包含有卡塔蘭數的那些結果是重要的, 以前 的研究沒有涉及, 因此這裏提出, 帶有卡塔蘭數的無窮級數應當引起更多的注意。

2.3. 明安圖使用的計算無窮級數的方法

中國傳統數學中沒有無窮級數, 明安圖必須建立一套計算無窮級數的新方法。 這種工作具 有原創的特點¹⁵:

1、 他建立了無窮級數的定義和符號系統。

2、 他創立了無窮級數的加法、 減法、 乘法的新演算法。

3、 他實際進行了運算並獲得若干無窮級數的正確結果。

本文以無窮級數 y₁₀ 為例, 介紹他的無窮級數自乘的方法, 並根據他的 hh割圓密率捷法ii 第 3 卷第 34−35頁的內容, 將明安圖的演算法翻譯成現代數學語言。

13羅見今,明安圖是卡塔蘭數的首創者_,內蒙古大學學報, 19卷2期, 1988, 239-245。

14Larcombe P. J., On a Finite Polynomial Generating Function for Catalan Subsequences: An 18th Century Observation Proved, Congresses Numerantium (Winnipeg Canada), 141(1999), 49-60.

15羅見今_,明安圖計算無窮級數的方法分析_,自然科學史研究_{, 9}卷₃期, 1990, 197-207。

令

x = 2 sin α (0 < α < π/4), yⁿm= (2 sin mα)ⁿ (m ≥ 1) 明安圖得到

y10= 5y2− 5y2³+ y₂⁵ = 10x − 165x³/4 − 3003x⁵/4³− 21450x⁷/4⁵

−60775x⁹/4⁷−

∞

X

n=1

16⁴Cⁿ− 8 · 16³Cⁿ+1+ 21 · 16²Cⁿ+2

−20 · 16Cⁿ⁺³+ 5Cⁿ+4

x²ⁿ⁺⁹/4²ⁿ⁺⁷

這裏 Cⁿ 是卡塔蘭數。 當 n = 1, 2, 3, 括弧內的係數是 41990, 22610, 29716。

現在我們將此無窮級數寫成如下的形式:

y10= a0x +

∞

X

n=1

aⁿx²ⁿ⁺¹/4²ⁿ⁻¹ (15) 明安圖必須計算 y₁₀ 的自乘積, 即 y₁₀· y¹⁰。 他的方法像這樣:

y₁₀² =

 a0x +

∞

X

n=1

aⁿx²ⁿ⁺¹/4²ⁿ⁻¹

2

= (a0x)²+ 2

∞

X

n=1

a0aⁿx²ⁿ⁺²/4²ⁿ⁻¹+

 ^∞ X

n=1

aⁿx²ⁿ⁺¹/4²ⁿ⁻¹

2

(16)

關鍵是計算式 (16) 中末項, 他發明的方法是

 ^∞ X

n=1

aⁿx²ⁿ⁺¹/4²ⁿ⁻¹

2

=

∞

X

n=2

n−1

X

k=1

4an−ka^k



x²ⁿ⁺²/4²ⁿ⁻¹ (17)

他最終得到:

y²₁₀= b₀x²+

∞

X

n=1

bⁿx²ⁿ⁺²/4²ⁿ⁻¹ (18) 這裏

b0 = a²₀, b1 = 2a0a1; when n ≥ 2 : bⁿ= 2a0aⁿ+

n−1

X

k=1

4an−ka^k 就這樣, 明安圖終於解決了無窮級數自乘的難題。

3、 學習和紀念

3.1. 史學家和數學家對明安圖的研究

日本著名學者三上義夫 (Mikami Yoshio, 1863−1950) 說 (1910): “圓理發達為最 緊要之事件, 可比西洋之定積分, 其演算法則始于所謂杜氏九術。 及蒙古族人欽天監監正明安 圖, 積三十餘年之辛勞, 始考出解析方法, 且別附以六術。”¹⁶ 李儼 (1892−1962) 、 錢寶琮 (1892−1974) 做了基本的研究工作。

李迪 (1927−2006) 發表了明安圖傳記¹⁷和 10多篇論文, 何紹庚、 羅見今、 特古斯做了深 入的工作。 法國學者詹嘉玲 (Catherine Jami) 說¹⁸: 明安圖數學工作的特點是綜合了西方和 中國兩種數學傳統。

程式設計的先驅者、 美國斯坦福大學克努特 (D. E. Knuth) 教授在他的名著 hh電腦程式 設計藝術ii 中提到卡塔蘭數的歷史, 他說: “蒙古族中國數學家明安圖在 1750年前研究無窮級 數時提出了卡塔蘭數· · ·· · · [見羅見今, 內蒙古大學學報, 19(1988), 239-245; hh組合數學和圖 論ii (世界科學出版社, 1993, 68-70)]”¹⁹。

德爾比大學的拉坎布 (P. J. Larcombe) 博士按照筆者對明安圖的解讀, 發表了 7篇論文, 研究明安圖和卡塔蘭數, 例如 “卡塔蘭數的歷史: 中國的最早紀錄”²⁰, “18世紀中國發現的卡塔 蘭數”²¹, “展開式中帶有卡塔蘭數的正弦函數: 超幾何函數研究的一個注記”²², 等。

3.2. 關於明安圖星

1992 年明安圖誕生 300 周年時, 在呼和浩特召開了紀念會, 北京天文臺台長、 中國科學院 王綬琯院士和 80多位學者 (其中包括 30多位外國學者) 參加了會議。

明安圖的天文學成就在國際天文學領域獲得承認。 1999 年中國天文學家發現了新的小行 星 28242 號, 國際天文學聯合會小天體提名委員會 2002 年 5 月 26 日宣佈將 28242 號命名為

“明安圖星”²³。

16三上義夫:中國數學之發展, 1910。 見中國算學之特色, 1925 (1929林科堂譯自日文₎。

17李迪,蒙古族科學家明安圖,內蒙古人民出版社, 1978。

18詹家玲_{, 18}世紀中國數學中的傳統工作以及所受的西方影響_,國際數學史雜誌_{, 1988}。

19 Donald E. Knuth, The Art of Computer Programming, Vol. 1, Fundamental Algorithms, Third Edition, Tsinghua University Press 2002 (English, for sale in mainland China only), 407.

20See Mathematics Today, Bulletin of the Institute of Mathematics and its Application (IMA), Vol. 35, No. 3, (1999), 89.

21See Mathematical Spectrum, Vol. 32, No. 1, (1999/2000), 5-7.

22 See The Journal of Combinatorial Mathematics and Combinatorial Computing (JCMCC), Vol. 37, (2001), 65-74.

23The document of “2002 May 26 M. P. C. 45750 New Name of Minor Planet” said: “Ming Antu (1692-1765?) was a Chinese astronomer and mathematician of Qing Dynasty. During the decades of his service in the Imperial Observatory, he participated in compiling and editing three very important astronomical works.”

2002 年是明安圖誕生 310 周年, 2002 年 8 月 “明安圖星” 命名大會和那達慕²⁴大會在明 安圖的故鄉舉行, 王綬琯院士和中國科學院地質地理科學研究所的滕吉文院士參加, 500多位代 表和 2 萬多本地居民參加了慶祝, 舉行了 “明安圖科學成果研討會”, 政府決定將明安圖故鄉城 鎮命名為 “明安圖鎮”, 宣佈建立 “明安圖科學技術博物館”。 王綬琯院士寫了一首紀念明安圖 的詩, 廣為流傳:

科學巨星明安圖 冉冉升起康乾初 精修曆象制皇輿 割圜密率冠前驅 碧野連空故上都 循公往跡纂公書 欲上青天攬公裾 今朝興國並興區

—本文作者任教內蒙古師大科學史研究所(呼和浩特), 浙江大學科技與文化研究所(杭州)—

The Symposium on Analysis and Probability

日 期 : 2010年8月10日 (星期二) ∼ 2010年8月12日 (星期四)

目 的 : The goal of this symposium is to expose young researchers and students to the latest developments in the field of anal- ysis and probability. The main topics consist of

• Stochastic Analysis and Its Applications

• Infinite dimensional Analysis

• Probability

• Nonlinear Analysis and Its Applications

• Harmonic Analysis

Please refer to http://www.math.sinica.edu.tw for further details.

24“那達慕_”是蒙古民族在草原上舉行的大型傳統聚會。