不同於托勒密的弦表製作方法
戴尉庭
本文起源來自於簡單的建築學問題, 此方法在工程延伸討論中偶然的情況下發現, 此一發 現主要在於托勒密 [1] 的弦表製作方法有了另一解, 並且以此解進一步經由轉換而能快速建立 起三角函數表。 起源問題如下:
長度為 l2 與 l3 的兩根柱子分別於距離 l1 的 A, B 兩點並令其相互傾倒, 求其交點 G 的位置?
定義兩柱子的數學式為
x2+ y2 = l22, (x − l1)2+ y2 = l23, 求解後獲得
x= l12+ l22− l32 2l1
, y =
q
l22− x2 = 1 2l1
q
4l12l22− (l12+ l22− l32)2,
(1)
其中 y 又可改寫成類似海龍公式的形式 y= 1
2l1
p(l1+ l2+ l3)(l1+ l2− l3)(l1 − l2+ l3)(−l1+ l2+ l3).
而上式在經過移項後, 可獲得一個常數 ϕ, 如下式 ϕ = l1y= 1
2p(l1+ l2+ l3)(l1+ l2− l3)(l1− l2+ l3)(−l1+ l2+ l3). (2)
78
式 (2) 說明任意三角形只要知其三個邊長就可知道各對應的高。 舉例來說, 邊長 3、 4 與 5 的 直角三角形, 其 ϕ 常數經式 (2) 計算為 12, 則 ϕ 除以各邊長即為高, 也就是 4、 3 與 2.4。 另 一邊長 1、 1 與 √
2 直角三角形的例子也一樣, ϕ 等於 1, 各高為 1、 1 與 0.7071。 這常數就如 同身分證一樣, 適用於任意三角形, 下列手繪圖例即簡單闡明 ϕ 常數的存在, 未填的數值可由 讀者進行驗證。
l1 l2 l3 ϕ l1 的高 l2 的高 l3 的高 3 2.8 4.7 7.9793 2.6598
4.2 2.3 3.5 8.0498 1.9166 2 5 4 7.5993 3.7997
2.5 1.5 3 3.7417 1.4967 2.4944 1.2472 3.3 2.7 1.9 5.1300 1.5545 1.9000 2.7000 1 1 1 0.8660 0.8660 0.8660 0.8660
此 ϕ 常數在特殊三角形有不同的形式, 例如在隨後的弦長計算上會用到的等腰三角形在式 (2) 即可簡化成如下
ϕ = L1
2 q
4L2 − L21. (3) 在有了這個常數的概念後, 接下來就是弦長的計算, 而此法建立在弧長切割越小越接近弦長的 概念下進行。 首先描繪一個半徑為 1 的單位圓, 並在此圓上先畫一個以基準角 θ 構成的等腰三 角形, 此後的疊加都是以這三角形為基準, 這樣每次新構成等腰三角形所對應的圓心角就是以 θ 的倍數增加, 而等腰三角形的底邊就是弦長 C, 隨著基準三角形的疊加, 弦長 C 也隨之變化, 如下圖
基準疊加角 二次疊加角
的弦長變化
三次疊加角 的弦長變化
而此方法的展現, 則是以四邊形的疊加組合來達成, 也就是不論是疊加幾次的基準角都看成最 近一次與前幾次總合的基準角疊加而成的四邊形, 如圖
不管幾次疊加都看成 兩三角形疊加而成的 四邊形
當逐次疊加構成四邊形的概念建立後, 就有構成了前後次弦長的關係, 這裡主要仍是 ϕ 常數的 相關應用, 目的在於求取四邊形的 「對角線」 長度, 如上圖中的虛線部分, 由該圖中可知, 原本的 對角線在新一次的疊加時就變成四邊形的邊長, 也可說是拆解成其中一個三角形的邊長, 如圖:
GD = C1, BG= Cj, BD= Cj+1,
j = 1, 2, 3, 4, . . . ,
AB = AG = AD = 1.
由此圖, 我們的目的是求取弦長 BD = Cj+1, 這裡的 j + 1 是想建立起前後關係, 如 j = 1, 2, 3, 4, . . ., 這樣就能進行疊代計算, 先將焦點放置 △ABG 與 △AGD 上, 這兩個三角形 的共同邊是 AG, 所以兩者連結由此建立, 利用式 (1) 的 x 計算法先計算 E、 F 兩點距離 A 點的長度
E點 → AE = AG2+ AD2− GD2
2AG = 2 − C12 2 ,
F 點 → AF = AG − GF = AG − AG2+ BG2− AB2
2AG = 2 − Cj2
2 . 而此兩三角形的高資訊
DE = ϕAGD
AG = ϕ1, BF = ϕABG
AG = ϕj;
這裡的 ϕABG = ϕj, 也是想建立起前後關係, 每疊加一次 ϕ 常數就改變, 在得知這幾個資訊 後, 再來就是關鍵的連結動作, 我們將對角線部分與兩三角形的高提取分離出來觀察, 如圖
因為兩三角形的高分別垂直於 AG 而使得這樣的操作讓對角線 BD 能夠用畢氏定理計算出來, 並且能與上一次的弦長做連結。
BH = AE − AF = 2 − C12
2 − 2 − Cj2
2 = Cj2− C12 2 , HD= DE + BF = ϕ1+ ϕj,
⇒ BD = Cj+1 = s
Cj2− C12 2
2
+ (ϕ1+ ϕj)2
⇒ Cj+1 = 1 2
q
(Cj2 − C12)2+ 4(ϕ1+ ϕj)2 j = 1, 2, 3, . . . (4) 其中由式 (3) 可得
ϕj = Cj
2 q
4 − Cj2.
在獲得了這個關係式後, 就宣告了只要有了基準角的弦長, 其餘的弦長都是可計算的。 舉例來 說, 若需要繪製 37π 這個角度, 僅需將其切割得更小例如 3700π 並作為初始弦長再疊加回去即可, 切割的多寡端看要求的精度而定, 此一方法的適用範圍是 0 ∼ π 之間的任一角度均可用此法計 算得出弦長。 以 0.5 度
π 360
等分點舉例 :
先計算基準角 0.5 度, 將其分成精度較低的 3 等分 (初始弦長 = 1080π ) 後利用式 (4) 再疊加 三次可獲得 0.5 度的弦長 C ∼= 0.00872, 此時有了 0.5 度的基準弦長則可再利用式 (4) 繼續 疊加至 180 度即可製作出半圓周的弦表。 式 (4) 實現了不需使用托勒密定理與三角函數計算 的弦表製作方法。
除此之外, 這裡也存在另一形式使應用更為廣泛, 首先將 C1 與 Cj 分別視為一個三角形, 則可用式 (1) 與式 (2) 來描述等分點的座標
x1 = 2 − C12 2 , y1 = ϕ1,
xj = 2 − Cj2
2 , yj = ϕj.
將座標點資訊代回式 (4)
Cj+12 = (x1− xj)2+ (y1+ yj)2, 又依據畢氏定理 x2+ y2 = 12, 則又可改寫成
Cj+12 = 2 − 2x1xj + 2y1yj. 而 Cj+1 所包含的等分點資訊利用式 (1)
xj+1 = 2 − Cj+12
2 . 這些關係又可進一步簡化成
xj+1 = x1xj − y1yj, j = 1, 2, 3, . . . . (5) 其中由畢氏定理
yj =q
12− x2j.
式 (5) 若仔細觀察可發現其為三角函數複角公式的變化形 「cos(α + β) = cos α cos β − sin α sin β」, 也因為這變化形揭露了等分點的訊息, 所以經由簡單的代數運算即可獲得各個等 分點的座標, 而在單位圓下, x 與 y 分別代表 cos 與 sin, 這樣一來只要知道基準角的資訊則 半圓周任意等分角的三角函數表很快就能建立出來。 以前述製弦表的例子:
利用疊代 3 次而得的 0.5 度基準弦長 C ∼= 0.00872 來計算基準等分點座標資訊, x1 = 2 − C2
2 ∼= 0.99996 y1= ϕ1 = C 2
√4 − C2 = 0.00872
由此利用式 (5) 就能快速建立起每 0.5 度間隔增加的三角函數表。 同樣地, 此表的精度取決於 基準角的精度。
利用本文所描述的方法來求取任意角度的正餘弦值或是製作三角函數表, 方法其實相當精 巧簡單, 表示式也簡潔漂亮, 本以為前人已討論過, 但在資料搜尋的情況下, 卻無發現有人闡述 這方法, 因此簡易描述此方法的過程, 在此分享給大家。
文末附帶一提, 在觀察到 ϕ 常數的特性後, 便思索將此思路延伸至三維的四面體上, 過程 中也發現到一個類似 ϕ 常數的常數存在, 思路由下圖開始
AB = l1 AE = L1
AD = l2 BE = L2
BD= l3 DE = L3
則
x2 + y2+ z2 =L21, (6) (x − l1)2 + y2+ z2 =L22, (7) (x − a)2 + (y − b)2 + z2 =L23, (8) 其中 a, b 由式 (1) 獲得
a=l21+ l22 − l23 2l1
, (9)
b=ϕABD
l1
= 1 2l1
q
4l12l22− (l12+ l22− l32)2. (10) 解聯立求得 x、 y 與 z 後, 從 z 的表示式中可獲得一個常數 v
zϕ = v. (11)
在經過數個實例代值運算後發現, 這個 v 常數如同 ϕ 常數的性質一樣, 任一四面體均有, 只要 v 除以此四面體任一平面的 ϕ 即可獲得相對應的立體高, 相當方便, 可惜的是, 至本文定稿為 止尚未找到如同 ϕ 的簡潔表示式。 除此之外, 由此思路若繼續延伸至更多維度, 則各維度間的 構型是否仍然具有特定常數, 則仍待後續探討。
參考文獻
1. 蔡聰明。 星空燦爛的數學 (I) — 托勒密如何編製弦表? 數學傳播季刊, 23(2), 57-67, 1999。
—本文作者為機械工程師, 曾服務於華泰電子 —
