■地図と三角法
ちょうど今,地図上の座標(x,y,z)をジオイドと呼ばれる仮想的回転楕円体面上に投影する問題を考えている.この問題では,標高の高い地点の距離は楕円体面上の距離に直されると短くなるので,密に投影されることになる.三角法なので使用する式は難しくはないが,相当複雑になることは避けられないし,予備知識なしに理解できるほど単純な問題ではなさそうである.
ところで,三角法といえば,影の長さから棒の高さを知るのに三角比の表が必要となったことが始まりであることは容易に想像されるだろう.このことに関連して,
エリ・マオール「素晴らしい三角法の世界」青土社
が,実に多くのことを教えてくれた.今回のコラムでは,この本を参考にして,地図測量と三角法に関連する話題を取り上げることにする.
===================================
【1】メルカトール図法
われわれが最も慣れ親しんでいる世界地図は,メルカトールの地図と呼ばれるものである.前回のコラムでは等角写像を取り上げたが,メルカトールの地図も等角写像,すなわち,任意の点から任意の他の点までの正しい方向(方位角)を示す地図であって,海図として航海のナビゲーション用に広く使用されている.また,確かめたわけではないのだが,メルカトールの地図は宇宙衛星の軌道制御にも使われているらしい.
地図はその目的によって,地球上の2点の正確な距離,形(面積),方向を保存するように工夫されているが,地球は丸いので,これらのどれかを保存することは他のものを犠牲にすることになる.
地図投影の中で,もっとも単純なのは円筒を使って平面にする円筒図法であろう.円筒投影では赤道で接触している円筒に対して,地球の中心から投影するもので,円筒が開けられたとき平面地図が得られる.
x=Rλ,y=Rtanφ
円筒図法では経線(子午線)・緯線(卯酉線)は直交するが,緯線間隔は緯度とともに増加し,高緯度での歪みが大きくなる.これはtanφがもたらす結果である.
次に,南極に平面が触れるようにして,北極から投影する図法を考えてみよう.最初から平面を使ったこの投影法を仮に平面投影と呼ぶことにするが,この場合,北緯φの点と南緯φの2つの地点は,単位円に関して反転した点に投影される.
y =2Rtan(π/4+φ/2)
y’=2Rtan(π/4−φ/2)=1/y
このことから,平面投影法は方向が保存される正方位図法であることが示される.ただし,正距図法ではないので,地図上で2つの点を結んでいる線分は最短距離(測地線)を表すものではない.
平面投影では,経線は南極から放射する直線として,緯線は同心円として写されることになる.それを受けて,メルカトールが
(1)経線・緯線は直交し,
(2)任意の点から任意の他の点までの正しい方向(方位角)を示す
地図の製作にとりかかったのは,1568年のことといわれている.
すなわち,円筒投影と平面投影の両方の特徴を併せもつ地図製作であるが,この計画が実現されるためには,緯線(平行線)の間隔を決めてやる必要がある.メルカトールがどのようにしてこれを計算したかは知られていないため,正確に述べることは非常に難しいが,以下のように推測されている.
地球上の緯線は2πRcosφの長さであるが,地図上では2πRに引き延ばされる.その際の拡張比は
2πR/2πRcosφ=secφ
である.そして,写像が等角であるためには,相似条件
Δy=RsecφΔφ
を満たさなければならない.
これを現代的に表記すると,微分方程式:
dy/dφ=Rsecφ
の解:
y=R∫(0,φ)sectdt=Rlntan(π/4+φ/2)
となるが,もちろん,メルカトールの時代は,ニュートンとライプニッツが微積分を発明するよりずっと前のことである.
この解は平面投影のtan(π/4+φ/2)を引き継いでいるが,これは偶然の一致ではない.複素関数論によれば,複素関数の写像は等角写像となり,複素対数関数
w=f(z)=lnz
は同心円と放射線よりなる極座標系を直交座標系に変換してくれるのである.
メルカトールの地図には平面図法の性質が遺伝するので,正距図法とはならない.すなわち,地図上の任意の2つの点を結ぶ線分は測地線ではない.また,logがついている分,高緯度における歪みは円筒図法ほど誇張されにくいと思われるが,それでもグリーンランドは南アメリカより大きく表される(実際は1/9の大きさである).
===================================
メルカトール図法と円筒図法はどちらも直交座標系であって,表面的には似ているため,しばしば混同される.しかし,これまでの説明から,両者の投影原理は全く異なっていることがおわかり頂けたであろうか?
すなわち,メルカトール図法は解析的・計算的に得られたものであって,幾何的・投影的なものではない.今日でさえ,多くの地理の教科書で誤った記事が掲載されていることはまことに残念なことである.
===================================
【2】レギオモンタヌスの問題
レギオモンタヌスの問題とは「空中に垂直に支えられた棒が,棒の延長線と地面との交点から水平に測って,どの地点で最大に見えるか」という問題である.近すぎる場合,棒はひどくつぶれて見えるだろうし,遠すぎる場合は単に小さく見えることになる.
最大視角を張る地点を求める問題であるから,交点から測った棒の上端をa,下端をb,水平方向の位置をxとして,三角法を使って解いてみよう.
tanθ=tan(α−β)
=(tanα−tanβ)/(1+tanαtanβ)
=(a/x−b/x)/(1+a/x・b/x)
=(a−b)x/(ab+x^2)
これを最大にするxは,微分によりx=√(ab)と計算され,求める点は棒の上部と下部の端点の高さの幾何平均の距離に位置することがわかる.
この解は,幾何学的には方ベキの定理(OA・OB=OP^2)によっているが,目の高さを考慮すると,この位置は棒の上端と下端を通る円が見ている人の眼の位置に水平に接するときの地点となる.円周角(同じ弧の上の角)を考えてみれば容易に理解できるであろう.
レギオモンタヌスの問題は,ラグビーでトライが決まったときに,どこからゴールを狙えばよいかなど,いろいろな場面で応用できる問題となっている.
===================================
【3】役に立たない三角関数公式集
レギオモンタヌスの問題の証明では,正接の加法公式が用いられているが,正接のn倍角公式は,パスカルの三角形を用いて,次のように書くことができる.
tannα=(nC1tanα−nC3tan^3α+nC5tan^5α−・・・)/(nC0−nC2tan^2α+nC4tan^4α−・・・)
分母・分子の係数は,2項係数の符号が対で交代するパスカルの正接三角形
1
1 1
1 2 −1
1 3 −3 −1
1 4 −6 −4 1
1 5 −10 −10 5 1
の形に並べることができる.ほとんどの教科書から消えてしまったが,美しい公式である.以下に,このような公式を収集してみた.
===================================
[1]三角形,四角形についての公式
任意の三角形に対して
tanα+tanβ+tanγ=tanαtanβtanγ
が成り立つ.
この式は
γ=π−(α+β)
として,tanの加法公式を用いることにより容易に証明される.役に立つかどうかは別として,私にとってこの公式は対称性のある美しい公式と感じられる.もちろん,美しく感じるかどうかは主観的であり,強制すべきものではないが,きっと多くの人の美意識にも訴えるに違いない(希望).
同様に,任意の三角形に対して
sinα+sinβ+sinγ=4cosα/2cosβ/2cosγ/2
sin2α+sin2β+sin2γ=4sinαsinβsinγ
sin3α+sin3β+sin3γ=−4cos3α/2cos3β/2cos3γ/2
cosα+cosβ+cosγ=1+4sinα/2sinβ/2sinγ/2
===================================
等式の世界も面白いが,不等式の世界だって奥深いものがある.
鋭角三角形ならば,算術平均≧幾何平均より
tanα+tanβ+tanγ≧33√tanαtanβtanγ
前項より,
tanαtanβtanγ≧33√tanαtanβtanγ
したがって,
tanαtanβtanγ≧√27=3√3
であるから,
tanα+tanβ+tanγ≧3√3 (等号は正三角形のとき)
を容易に証明することができる.
少し気分を変えて,次の不等式はどうだろうか?
(問題)
sinαsinβsinγ≦3√3/8
(証明)
2sinβsinγ=cos(β−γ)−cos(β+γ)
=cos(β−γ)+cosα
sinαsinβsinγ
=1/2sinα(cos(β−γ)+cosα)
≦1/2sinα(1+cosα)
これより極大値を計算すると,3√3/8が得られる.なお,この不等式は三角形の外接円,内接円および面積をR,r,△とすれば,
abc=4R△,(a+b+c)r=2△
また,正弦法則
a/sinα=b/sinβ=c/sinγ=2R
より,
abc≦3√3R^3
と同値である.
===================================
(問題)△ABC内の1点をP,面積を△とすると,
AP+BP+CP≧2√((√3)△)
等号は△ABCが正三角形で,Pが重心のときに限る.
(証明)
n角形の周の長さが与えられているとき,面積の最大のものは正n角形であるから,
L^2≧4nStan(π/n)
等号は正n角形に対してのみ成り立つ.
そこで,Pの辺BC,CA,ABに関する対称点をA’,B’,C’とし,六角形AC’BA’CB’に不等式
L^2≧4nStan(π/n)
を使えばよい.ここで,L=2(AP+BP+CP),S=2△
なお,△≧√(27)r^2が成り立つので,
AP+BP+CP≧6r
これは,エルデシュの定理の特別な場合になっていて,シュライバーの定理とも呼ばれる.
一方,四角形については,プトレマイオス(トレミー)の定理「円に内接する四角形の対角線の積は,対辺の積の和に等しい」がある.
AC・BD=AB・CD+BC・DA
この定理において,もし四角形が長方形ならば
AC^2=AB^2+BC^2
となり,ピタゴラスの定理に帰着する.
また,4点が同一円周上にないとき,不等式
AC・BD<AB・CD+BC・DA
が成り立つ.
===================================
[2]正弦・余弦の和公式
等差級数
S=1+2+3+・・・+n
の値を求めるのに,逆順にして
S=n+(n−1)+(n−2)+・・・+1
辺同士を加えると
2S=(n+1)+(n+1)+(n+1)+・・・+(n+1)
より,
S=n(n+1)/2
これが等差級数の和公式で,これを使うと,たとえば1から100まで数の合計が5050であることが瞬時に計算できることはご存知であろう.
この取り扱いと似た方法で,正弦の和公式
sinα+sin2α+sin3α+・・・+sinnα
=sinnα/2sin(n+1)α/2/sinα/2
を証明してみよう.
(証明)
T=sinα+sin2α+sin3α+・・・+sinnα
T=sinnα+sin(n−1)α+sin(n−2)α+・・・+sinα
ここで,和から積への式
sinα+sinβ=2sin(α+β)/2cos(α−β)/2
を用いると
2T=2sin(n+1)α/2{cos(1−n)α/2+cos(3−n)α/2+・・・+cos(n−3)α/2+cos(n−1)α/2}
両辺にsinα/2を掛けて,積から和への公式
sinαcosβ=1/2{sin(α+β)+sin(α−β)}
を用いると
2Tsinα/2
=sin(n+1)α/2{sinnα/2+sin(1−n/2)α+・・・+sin(−1+n/2)α+sinnα/2}
=2sin(n+1)α/2sinnα/2
同様に,余弦の和公式
cosα+cos2α+cos3α+・・・+cosnα
=sinnα/2cos(n+1)α/2/sinα/2
も証明できる.これらの式において,α=π/nとおくと
Σsinkπ/n=cotπ/2n
Σcoskπ/n=1
さらに,
sinα+sin3α+sin5α+・・・+sin(2n−1)α=sin^2nα/sinα
cosα+cos3α+cos5α+・・・+cos(2n−1)α=sin2nα/2sinα
α=π/(2n+1)とおくと,
Σcos(2k−1)π/(2n+1)=1/2
===================================
[3]正弦・余弦の積公式
正弦・余弦の和公式はフーリエ級数との関連で研究された歴史がある.一方,和公式ほどよく知られていないが,正弦・余弦の積公式としていろいろな公式が登場してくる.ここでは証明は省いたが,複素数を使って証明するのが一番の近道であろう.
Πsinkπ/n=sinπ/n・・・sin(n−1)π/n
=n/2^(n-1)
Πsin(θ+kπ/n)
=sin(θ+π/n)・・・sin(θ+(n−1)π/n)
=sinnθ/2^(n-1)sinθ
ここで,θ→θ−π/2nと置き換えれば
Πsin(θ+(2k−1)π/n)=cosnθ/2^(n-1)
θ=0とおけば
Πsin((2k−1)π/n)=1/2^(n-1)
また,θ=π/2とおけば
Πcoskπ/n=sin(nπ/2)/2^(n-1)
などを導き出すことができる.
===================================
[4]漸近挙動の世界(無限的公式)
微積の学び初めに,x→0としたとき,
sinx/x→1
に出会う.この結果は
(sinx)’=cosx,(cosx)’=−sinx
を示すのに用いられる.
その後,sinxのテイラー展開によって,無限級数
sinx=x−x^3/3!+x^5/5!−x^7/7!+・・・
sinx/x=1−x^2/3!+x^4/5!−x^6/7!+・・・
が示される.
それでは,任意のxに対して,無限積公式
sinx/x=cosx/2cosx/4cosx/8・・・
も示しておこう.
(証明)
sinx=2sinx/2cosx/2
=4sinx/4cosx/4cosx/2
=8sinx/8cosx/8cosx/4cosx/2
・・・・・
=2^nsinx/2^ncosx/2^n・・・cosx/2
書き直すと
sinx=x[sin(x/2^n)/(x/2^n)]cosx/2・・・cosx/2^n
ここで,n→∞のとき,
sin(x/2^n)/(x/2^n)→1
であるから,sinxの無限積表示
sinx=xΠcosx/2^n
=x(1−x^2/π^2)(1−x^2/4π^2)(1−x^2/9π^2)・・・
が得られる.この結果は,sinxがx=0,±π,±2π,±3π,・・・のとき,0になることに一致している.
[補]正弦積分とは,
Si(x)=∫(0,t)sint/tdt
=x−x^3/3・3!+x^5/5・5!−・・・
として定義される特殊関数(初等関数によって表し得ない関数)である.また,その特殊値
Si(∞)=∫(0,∞)sint/tdt=π/2
は,ディリクレ積分と呼ばれる.
===================================
前項では,cosxの無限積表示
cosx=(1−4x^2/π^2)(1−4x^2/9π^2)(1−4x^2/25π^2)・・・
を用いているが,tanxに対しては,部分分数の無限級数表示
tanx=8x[1/(π^2−4x^2)+1/(9π^2−4x^2)+1/(25π^2−4x^2)+・・・]
が成り立つ.
x=π/4とすると,
1=4/π(1−1/3+1/5−1/7+・・・)
であるから,グレゴリー・ライプニッツ級数
π/4=1−1/3+1/5−1/7+・・・
が導かれる.
グレゴリー・ライプニッツ級数はπを含んでいる無限級数として最初のものなのだが,オリジナルは
arctanx=x−x^3/3+x^5/5−x^7/7+・・・
から発見されたものである.
また,x→0としたときのtanx/xの漸近挙動から,
π^2/8=1/1^2+1/3^2+1/5^2+・・・
さらに,
S=1/1^2+1/2^2+1/3^2+・・・
=1/1^2+1/3^2+1/5^2+・・・+1/2^2+1/4^2+1/6^2
=1/1^2+1/3^2+1/5^2+・・・+1/4[1/1^2+1/2^2+1/3^2+・・・]
=π^2/8+S/4
したがって,
S=1/1^2+1/2^2+1/3^2+・・・=π^2/6
となるが,これはオイラーにより発見された有名な級数である.
===================================
[5]無限級数とフーリエ展開
テイラー展開(ベキ級数展開)は連続な関数にしか適用できないが,フーリエ展開(三角級数展開)はたとえ不連続な関数でも存在するため,19世紀の解析学に大きな貢献をしるした.
−π<x<πの周期関数と見なされる不連続関数f(x)=xのフーリエ展開は
f(x)=2[sinx/1−sin2x/2+sin3x/3−・・・]
x=π/2とおくと,
π/4=1−1/3+1/5−1/7+−・・・(グレゴリー・ライプニッツ級数)
また,x=π/4とおくと
π√2/4=1+1/3−1/5−1/7++・・・
この式で,符号は2項毎に交代する.
同様に,−π<x<πの周期関数と見なすと,f(x)=x^2のフーリエ展開は
f(x)=π^2/3−4[cosx/1^2−cos2x/2^2+cos3x/3^2−・・・]
x=πとおくと,
π^2/6=1/1^2+1/2^2+1/3^2+・・・(オイラー級数)
===================================
[6]虚の三角法
オイラーの公式(1748年)
exp(iθ)=cosθ+isinθ
ド・モアブルの定理(1722年)
(cosθ+isinθ)^n=cosnθ+isinnθ
などは虚(あるいは複素数)の三角法の公式と捉えることがができるであろう.
複素数の世界では,z=x+iyとしたとき,
sinz=sinxcoshy+icosxsinhy
cosz=cosxcoshy−isinxsinhy
expz=exp(x)cosy+iexp(x)siny
などのように書くことができる.
sinz,coszが周期2πをもつこと
sin(z+2π)=sinz
は複素数の世界でも同じであるが,一方,expzは虚の周期2πiをもつ.
exp(z+2πi)=expz
複素関数は等角写像になるのであるが,前述のごとく,複素関数lnzも等角写像となる.このことから,メルカトール図法を極座標から直交座標上に等角変換できるのである.
なお,平面正弦定理では,
sinα:sinβ:sinγ=a/R:b/R:c/R
であるが,球面正弦定理は
sinα:sinβ:sinγ=sin(a/R):sin(b/R):sin(c/R)
で表される.
それに対して,双曲的三角法では,球面正弦定理のRをiRに置き換えることによって,
sinα:sinβ:sinγ=sinh(a/R):sinh(b/R):sinh(c/R)
が得られる.
===================================
[7]エピサイクロイドとヒポサイクロイド
いま,5歳の娘は歯車の穴に鉛筆を差し込んで回転させると花びら模様が描かれるおもちゃに夢中である.これをスパイログラフというらしい.この装置がヒポサイクロイドの応用であることはすぐに理解される.固定円と回転円の半径比R/rが無理数なら曲線は決して閉じないから,有理数倍になっているのであろう.
ヒポサイクロイドで,R/r=4の場合がアステロイドである.
x^2/3+y2/3=R^2/3
アステロイドはx,y軸上に端点のある長さRの線分により作られる包絡線であり,また,長半径と短半径の和がRである楕円
x^2/a^2+y^2/(R−a)^2=1
の包絡線ともなっている.
アステロイドが囲む周長は6Rであるが,サイクロイドと同様,定数πに依存していない.また,面積は3πR^2/8で,固定円の面積の3/8(回転円の6倍)である.
一方,エピサイクロイドは地球から見たときの惑星の逆行運動の説明に用いられた曲線で,古代ギリシア人は,惑星の動きを表現するために周転円(円の周りをまわる円)を考えていたことが知られている.
エピサイクロイドで,R/r=1の場合がカージオイドである.カージオイドは定円上に中心があり,カスプを通るような円の包絡線でもある(周長16R,面積6πR^2).
===================================