約略算術及其應用

約略算術及其應用

劉源俊

我們自小學算術, 不管是加、 減、 乘、 除、 開方, 通常講究的是精確。 中學階段又學到不等 式, 當時腦裡仍然存的是精確的概念。 其實, 有許多情況, 近似算法不祇有用, 甚至我們也習慣 於近似; 特別在音樂的算學裡, 還必須用到約略的觀念才行。 我稱這類算術為 「約略算術」 (翻譯 為英文, 可作 approximate arithmetic), 它講究的是 「扼重要」, 有別於講究 「求精確」 的 「通 常算術」 (ordinary arithmetic)。

一般講 「科學精神」 時總不忘提到 「求精確」, 其實 「扼重要」 更是 「求精確」 之本。 《禮 記·大學》 裡 「物有本末, 事有終始, 知所先後, 則近道矣。」 一句可為註解。

以下介紹約略算術, 大部分算式用 「約等於」 符號 (≈) 而不用 「等於」 符號 (=)。

一、 從圓周率 π 與閏年說起

大家熟知的圓周率 π = 3.14159265 · · · , 近年用計算機已可算至約十兆 (10¹³) 位小數。

有些印度孩童作興即時背出幾百位小數。 其實算出小數位十幾位以後的數字, 在數學上唯好玩 而已, 在實驗科學裡著實沒用 ! 古書有言 「圓徑一而周三。」 (《周髀算經》 上卷 : 「勾股圓方 圖。」 漢 · 趙爽注), 意思就是圓周率 π 約等於 3, 當然太粗略。 西元前約 250 年, 阿基米德已知

223

71 < π < ²²₇ , 其後知道 ²²₇(≈ 3.142857) 其實是圓周率相當好的近似數, 誤差約萬分之四, 平常很夠用了。 到了約西元 480 年, 南北朝時期南朝宋人祖沖之除了給出 「約率」π ≈ ²²7 之外, 更得出 「密率」 π ≈ ³⁵⁵113(3.14159292 · · · ) — 此一約略分數更準確到千萬分之一。

利用連分數的表示法 π = 3 + ₇₊ ¹1 15+ 1

1+···

, 可以列出 π 的近似分數, 依次是 3, ²²₇ , ³³³₁₀₆,

355 113, · · · 。

大家也都知道每四年要 「閏年」 一次, 但到了第一百年要 「不閏年」 一次。 為什麼? 因為一

「回歸年」 為 365.2422 · · · 日, 換言之每一年多出 0.2422 · · · 日, 而 0.24 · · · ≈ ¹4−100¹ 。 然而 這樣只照顧到前面兩位小數 0.24。 若要再照顧到後面的小數 0.0022 · · · , 國際的約定是: 第四 百年要閏年, 第八個四百年則又不要閏年; 為什麼呢? 原來 0.2422 · · · ≈ ¹4−100¹ +₄₀₀¹ −3200¹ 。

56

然而約略等式當然不只一種, 所以也不妨規定第五百年要閏年, 第五千年也要閏年; 其所依據的 近似式則是 0.0022 · · · ≈ 500¹ + ₅₀₀₀¹ 。

中國曆法是陰陽合曆, 所以一定要照顧到陰曆月。 陰曆閏月的道理類似: 因為一 「朔望月」

為 29.5306 日, 而 365.2422 ÷ 29.5306 ≈ 12.3683 ≈ 12 + 19⁷, 所以每 19 年需閏月 7 次。

因為 ₁₉⁷ = 0.3684 · · · , 與 0.3683 相差無幾, 所以此一約略分數式相當準確, 可以沿用歷久。

二、 自然對數底 e 與 2 及 10 的關係

數學裡有兩個最重要的常數, 一個是上節講的圓周率 π, 另一個是所謂 「自然對數的底」 e。

根據指數函數的定義 e^x = 1+_1!^x+^x_2!²+^x_3!³+· · · , 於是 e = 1+1!¹+_2!¹+_3!¹+· · · = 2.71828 · · · 。 現在介紹兩個關於 e 的約略等式。 首先, e⁷ = 1, 096.3 · · · ≈ 2¹⁰, 誤差約 7%, 大了些!

但若兩邊取對數, 則得 ln 2 ≈ 0.7, ln 2 = 0.693 · · · , 此一約略等式的兩邊誤差減至 1%, 可以 接受了! 換言之, 2 ≈ e^0.7。

其次, ln 10 ≈ 2.30 ≈ ⁷3, 或 10 ≈ e^2.30。

以下是兩個應用實例: 其一, 在原子物理、 固態物理、 核物理或天體物理裡常遇到 e 的 冪次, 需換算為 10 的冪次, 以利認知。 例如問 : e⁻⁴² 以 10 的冪次記, 為何? 答 : e⁻⁴² ≈ 10^−242.30 ≈ 10⁻¹⁸。

其二, 在複利 (按月) 計算中, 常需估算本金經多少年能翻倍。 有關的公式是 n ≈ ^{ln 2r} , 式 中 r 是年利率, n 是翻倍的年數。 根據前面所述 ln 2 ≈ 0.7 的近似值, 便可輕易算出: 如果年 利率是 2%, 要 35 年; 如果利率達 18%, 則只需要約 4 年。

三、 黃金比 ϕ

有個出名的數列, 名費氏數列 (Fibonacci sequence) : 1, 1, 2, 3, 5, 8, 13, 21, 34, · · · , 除第一、 二兩數為 1 外, 其後每一數都是前兩數的的和。 不難證明, 後一數與前一數的比 ²₁, ³₂,

5

3, ⁸₅, ¹³₈ , · · · 。 這一數列趨近 ^√5+12 ≈ 1.61803 · · · ; 此極限值就是所謂的 「黃金比」 ϕ (the golden ratio)。

取上列第五個數 ¹³₈(= 1.625) 為例, 與黃金比的誤差約為 0.4%。 取第六個數 ²¹₁₃, 則誤差 減小為約 0.16%。 其後每一分數都可當作是 ϕ 的約略數, 越來越接近。

費氏數列在天地間常常可看到, 例如在螺旋生長的機制裡 : 向右 (左) 的螺旋數與向左 (右) 的螺旋數的比一定是費氏數列 — 觀察鳳梨、 桑椹、 松果、 向日葵花, 都可發現。 生物界裡 應用的實例還很多, 在此不贅。

四、 約略等式 5

≈ 2

現在一說簡單的近似算式 5³ ≈ 2⁷。 大家都會算 : 5³ = 125, 2⁷ = 128, 兩數差 2.4%。

利用這一近似算式, 我們來推演一下 : 兩邊各乘上 2³ 再開 10 次方, 於是得到 10^0.3 ≈ 2 ⇒ 10^0.7≈ 5, 10^0.1 ≈ 2

1 3。

但 10^0.1 或 2¹³ 約等於多少呢? 將 2⁷ ≈ 5³ 兩邊各三次方, 再稍微移項, 不難得到 2 · 5 = 10 ≈ ⁵2¹⁰²⁰ = (⁵₄)¹⁰, 於是知道 : 10^0.1 ≈ ⁵4 = 1.25。 10^0.1 的精確值是 1.2589 · · · , 2¹³ 的精確 值是 1.2599 · · · , 與 ⁵4 這一約略式的誤差值各約為 0.7% 與 0.8%。

10³ ≈ 2¹⁰ 這個式子常用到。 計算機科學裡常講的的 1K 約略等於 10³, 相當於二進位 10 個位元的消息量, 精確值為 1,024。 以此類推, 則 1M 相當於 20 個位元, 1G 相當於 30 個位 元, 1T 相當於 40 個位元。

以下舉兩個簡單的應用實例。

其一 : 某種細菌每 20 分鐘增生一倍, 問一天共增生若干倍? 因為一天 24 小時, 一小時 是 20 分鐘的 3 倍, 所以答案為 2^3×24 = 2⁷² ≈ 4 × 10²¹。 一分鐘呢? 答案是 2²⁰¹ ≈ e

0.7 20 = e^0.035 ≈ 1.035。

其二 : 一個 32 位元電腦所處理的數字, 最多有效位數以二進位言是 31 位, 問相當於十 進位數的幾位數? 答案是 31 × 0.3 ≈ 9。

常聽說音量高若干分貝, 卻很少人清楚 : 1 分貝 (decibel) 是什麼意思? 原來, 人感覺響 度的大小與實際聲音功率強度呈對數關係, 於是定義 B= log10

P

P₀, 式中 P 是測到的功率, P0

是某一標準功率, B 是響度的 「貝數 (bel)」。 所以, 響度高 1 貝指功率大 10 倍。 1 分貝則是 1 貝的十分之一, 定義為 dB=10 log10

P

P₀, 所以響度高 1 分貝意指輸出功率大 10¹⁰¹ 倍, 亦即約 為 1.26 倍。 準此, 則響度高 3 分貝指輸出功率大約大 10^0.3 ≈ 2 倍。

天文學裡, 「星等」 的定義也類似。 星等定義為 m = −2.5 log^10II0, 式中 I 為亮度, I0 為 某標準亮度, m 為星等。 因此, 一等星比六等星亮 10^25.5 ≈ 100 倍, 而 1.7 等星比 2.2 等星的 亮度則大 10^02.5.5 = 10^0.2 ≈ (1.25)² ≈ 1.56 倍。

在地震學裡有所謂芮氏規模 (Richter scale), 大家都一定耳熟, 但未必能詳。 例如 1999 年 臺灣集集 「九二一」 地震的規模是 7.3, 2008 年四川汶川地震的規模是 8.0, 試問汶川地震釋放 出的能量比集集地震大多少倍? 原來芮氏地震規模與釋放能量的關係可寫為 M = ²₃log_{10 E}^E

0, 式中 E 為釋放能量, E0 為某標準能量, M 為規模大小。 所以規模若差 1, 釋放能量差 10³² = 31.62 · · · ≈ 32 倍。 而汶川地震釋放的能量比集集地震大 10^0.7×32 ≈ 5

3

2 ≈ 11 倍。

五、 約略等式 3

≈ 2

再來看 3¹² ≈ 2¹⁹ 這一約略等式。 3¹²= 531, 441, 2¹⁹ = 524, 228; 兩者約略相等, 誤差

-1.4% 。 大家一定很難想到這一約略不等式是音樂裡 「十二平均律」 的基本原理。

從物理觀點看, 一個 「高八度」 的音與原音 (例如高音 C 與 C) 之差在哪裡呢? 差在高 音 C 的頻率是 C 頻率的 2 倍。 現在如果我們把兩者之間的這一頻率比 2 依等比級數平均分 為 12 等分, 則每兩相鄰音 (兩者差 「半音」, 例如 C# 與 C, 或 F 與 E) 之間的頻率比便是 2¹²¹ = 1.0594630 · · · ≈ 1.06, 這就是所謂 「十二平均律」; 於是, G (Sol) 與 C (Do) 之間的 頻率比便是 2¹²⁷ 。 鋼琴就是依據十二平均律調音的典型樂器。

音樂學界把每個 「半音」(semitone) 又分為 100 分, 稱為 「音分」(cents), 於是高一音分 相當於頻率高 2¹²⁰⁰¹ = 1.000577 · · · 。 音感佳者的耳朵可以分辨出四分之一個半音, 亦即 25 音分的差別 — 相當於比差 2¹²⁰⁰²⁵ = 2⁴⁸¹ ≈ 1.015; 據說稟賦特優的耳朵可以分辨出八分之 一個半音 — 相當於比差 1.007。 十二平均律音階的 G 與 C 能和諧嗎? 讓我們來算算看 : 2¹²⁷ = 1.49830 · · · ≈ ³2, 兩者相比約為 1.001。 這樣小的頻率差, 再好的耳朵也很難辨別! 所以 結論是: 十二平均律的 G 與 C 相當和諧!

其間的奧秘說穿了就是 2¹²⁷ ≈ ³2 這一約略等式; 將此式兩邊各乘方 12 次, 再略加整理, 就得到 3¹² ≈ 2¹⁹!

以下列表簡單說明天然律 (純律) 與十二平均律裡諸音的頻率差:

音名 頻率比 (天然律) 頻率比 (十二平均律) 兩者比值

C(Do) 1 1 1

D(Re) ³₂²3 = 1.125 2¹²² ≈ 1.1225 0.9978 E(Mi) ³₂⁴6 = 1.2656 · · · 2¹²⁴ ≈ 1.2599 0.9955 F(Fa) ²₃² = 1.3 2¹²⁵ ≈ 1.3348 1.0011 G(Sol) ³₂ = 1.5 2¹²⁷ ≈ 1.4983 0.9987 A(La) ³₂³4 = 1.6875 2¹²⁹ ≈ 1.6818 0.9966 B(Si) ³₂⁵7 = 1.8984 · · · 2¹¹¹² ≈ 1.8877 0.9944

高音 C 2 2 1

從上表列十二平均律諸音頻率與天然律諸相關音頻率的比值可以看出, 最多差不過 10 音 分 (比差 ≈ 1.006)。 測量儀器固然可以分辨, 經過特殊訓練的耳朵或許能夠分辨, 但普通人的 耳朵是絕對聽不出差異的。

其實, 每種樂器發出的各個基本頻率本身, 總附帶一定的 「頻寬」; 換言之, 發出的頻率並 不十分準確, 出入比差 1% 是正常的。 換一個角度說, 還正是因為我們的聽覺並不十分敏銳, 我

們才不會太挑剔, 而生美感, 也才會接受十二平均律的 「和諧」。 而人的美感也是可以培養的 — 本來, 習慣成自然!

六、 約略速算簡介

通俗科學書 《別鬧了, 費曼先生!》 中有一節介紹物理學家費曼自述其約略速算的本領。 例 如他一下子可以心算出 e^3.3 ≈ 27.11 與 e^3.0 ≈ 20.08。 為什麼呢? 首先要熟知: 當 x 小時, e^x ≈ 1 + x。 於是, 從本文前幾節所述, 不難推得:

e^3.3 = e¹e^2.3 = e¹e^2.3026e^−0.0026 ≈ 2.718 × 10 × (1 − 0.0026) ≈ 27.11; e^3.0 = e^2.3e^0.7 = e^2.303e^−0.003e^0.693e^0.007= e^2.303e^0.693e^0.004≈ 10 × 2 × 1.004 = 20.08。

費曼所舉其他速算的例子, 好些是與二項式定理有關; 在此要介紹一下。 二項式定理最早是 巴斯卡 (Pascal, 1623∼1662) 提出的, 當時只用在整數冪次情形。 牛頓在大學四年級時 (∼1664) 將它推廣到任何冪次情形:

(1 + x)^a= 1 + ax + a(a − 1)

2! x² + · · ·^當x小時≈ 1 + ax.

這一定理在約略算術裡極為有用, 以下舉幾個例子, 讀者當可舉一反三。

例一: 28¹³ = (27+1)¹³ = 3(1+₂₇¹)¹³ ≈ 3+27¹ = 3+₂₅¹(1+₂₅² )⁻¹ ≈ 3+0.040(1−0.080) ≈ 3.037。 如果熟練, 的確可以在兩秒鐘內給出答案, 準確到四位數。

例二: 40¹² = (36 + 4)¹² = 6(1 +¹₉)¹² ≈ 6 + ¹3 ≈ 6.33。

例三: 2¹² = (1.96 + 0.04)¹² = 1.4(1 + ^0.04_1.96)¹² ≈ 1.4(1 + 0.01) = 1.414。

例四 (費曼舉的): _1.73¹ = 7 ¹

4−0.02 ≈ ⁴7(1 − 0.01)⁻¹≈ 0.571429(1 + 0.01) ≈ 0.577。

想熟練約略速算, 有需要先記憶下列基本約略等式或等式: √

2 ≈ 1.414, √

3 ≈ 1.732,

√5 ≈ 2.236, √

10 ≈ 3.162, e ≈ 2.718, e⁻¹ ≈ 0.368, π ≈ 3.142, ln 2 ≈ 0.693, ln 10 ≈ 2.303; ¹7 = 0.142857, ²₇ = 0.285714, 等等。

—本文作者任教東吳大學—