• 検索結果がありません。

第 1 章曲線と曲面

N/A
N/A
Protected

Academic year: 2021

シェア "第 1 章曲線と曲面"

Copied!
44
0
0

読み込み中.... (全文を見る)

全文

(1)

1章 曲線と曲面

1.1 平面曲線

この節では、平面曲線の基本事項を復習し、平面曲線の曲率が曲線の形を決め ていることと曲率から曲線を構成できることを示す。

一つのパラメータによって平面の点の軌跡として平面曲線を定める。Rの開区 I上で定義されR2に値を持つ写像p:I R2について考える。pの接線の変化 を調べることによって、pの曲り方をとらえる。Iの元tに対して微分p(t)0 なければ、p(t)におけるpの接線を引くことができる。そこで、任意のtIに対 してp(t)̸= 0となるときに、像p(I)を平面曲線と呼び、写像p:I R2を平面 曲線のパラメータ表示 という。

平面曲線のパラメータ表示が

p(t) = (x(t), y(t))

によって与えられているとする。

p(t) = (x(t), y(t))

を速度ベクトルとも呼ぶことにする。運動する点の時刻tのときの位置がp(t) あると考えると、p(t)は平面の点の軌跡を表し、p(t)はその軌跡の速度ベクトル とみなせる。速度ベクトルの長さは

p(t)=

x(t)2+y(t)2

で与えられる。パラメータがaからb (a b)まで動くときの曲線Cの長さを L(C) =

b a

p(t)dt

で定める。この曲線の長さの定義がパラメータの変更によって変わらないことは、

積分の変数変換の公式よりわかる。始点t=aからパラメータtまでの曲線の長さ sで表すことにすると、stの関数であり

s=

t a

p(t)dt, ds dt = d

dt

t a

p(t)dt =p(t)>0

(2)

が成り立つ。最後の不等式はp(t)̸= 0という仮定からわかる。微分が正なので、

stに関して単調増加になり、逆関数が存在する。よって、tsの関数として t=t(s)を考えることもできる。これを元の曲線に代入し(x(t(s)), y(t(s)))とする と、パラメータtによる表示からパラメータsによる表示を得る。このパラメータ sを曲線の弧長パラメータと呼ぶ。幾何学的な意味を考えると弧長パラメータs 関する速度ベクトルは単位ベクトルになることがわかる。曲線に対して弧長パラ メータは平行移動と向きを逆にすることを除けば一意的に定まる。そのため曲線 の一般論を展開する際には、弧長パラメータは使いやすいパラメータである。そこ で、R2の曲線を弧長パラメータsを使ってp(s)と表す。sが弧長パラメータであ ることからp(s)の速度ベクトルp(s)は単位ベクトルになり、R2 の内積を⟨ ·, · ⟩ で表すと、p(s),p(s)= 1が成り立つ。この両辺をsで微分すると

0 =p(s),p(s) =p′′(s),p(s)+p(s),p′′(s)= 2p(s),p′′(s)

となり、⟨p(s),p′′(s) = 0を得る。これはp(s)p′′(s)が直交することを意味す る。e1(s) = p(s)とおくとe1(s)e1(s)が直交することになる。e1(s)を反時計 回りに90度回転したベクトルをe2(s)で表すと、e1(s)e2(s)は線形従属になり

e1(s) =κ(s)e2(s)

を満たす関数κ(s)が定まる。κ(s)が曲線の曲率である。ここで、e1(s),e2(s)= 0 の両辺を微分すると

0 = e1(s),e2(s)+e1(s),e2(s) となり、

e1(s),e2(s)=−⟨e1(s),e2(s)=κ(s) がわかる。⟨e2(s),e2(s)= 1 の両辺を微分すると

e2(s),e2(s)= 0 がわかる。e1(s),e2(s)R2の正規直交基底なので、

e2(s) =e1(s),e2(s)e1(s) +e2(s),e2(s)e2(s) = κ(s)e1(s).

以上より次の定理を得る。

定理 1.1.1 (フレネの公式) 弧長sをパラメータに持つ平面曲線p(s)の曲率をκ(s) とする。p(s)の速度ベクトルe1(s) = p(s)e1(s)を半時計回りに90度回転した ベクトルe2(s)は次を満たす。

e1(s) =κ(s)e2(s), e2(s) = κ(s)e1(s).

(3)

R2の元を縦ベクトルとみなして、二つ並べたものを2次正方行列とみると、フ レネの公式を

d

ds[e1(s)e2(s)] = [e1(s)e2(s)]

[

0 κ(s) κ(s) 0

]

(1.1) と表せる。e1(s),e2(s)は曲線の動標構と呼ばれている。曲線の曲率は単位接ベク トルe1(s)の変化を記述するものとして定めたが、曲率は動標構の変化も記述して いることをフレネの公式は示している。さらにフレネの公式から曲率が曲線の形 を決めていることや、曲率からもとの曲線を構成できることを示す。これらを示 すためと後で曲面の場合に話を進めるために、Rnの運動群を考えておく。

n次実正方行列の全体をM(n,R)で表し、

SO(n) ={g M(n,R)|tgg =In, detg = 1}

によってn次回転群SO(n)を定義する。ここで,tggの転置行列であり、In n次単位行列である。回転群SO(n)が行列の積に関して群になることを示してお こう。g, hSO(n)に対して

t(gh)(gh) = thtggh=thInh=thh=In, det(gh) = det(g) det(h) = 1 となるので、ghSO(n)である。

t(g1)g1 =gg1 =In, det(g1) = (detg)1 = 1

となるので、g1 SO(n)である。以上より、SO(n)が行列の積に関して群にな ることがわかる。

平面曲線の話に戻って、動標構の定め方より[e1(s) e2(s)]SO(2)の元である ことを以下で証明する。

e1(s) = [

x(s) y(s) ]

とおく。これは単位ベクトルなので

x(s)2+y(s)2 = 1

が成り立つ。e2(s)e1(s)を反時計回りに90度回転したベクトルだから e2(s) =

[y(s) x(s)

]

(4)

となる。これらより

[e1(s)e2(s)] = [

x(s) y(s) y(s) x(s)

] ,

t[e1(s)e2(s)][e1(s)e2(s)] = [

x(s) y(s)

y(s) x(s) ] [

x(s) y(s) y(s) x(s)

]

= 1, det[e1(s)e2(s)] =

x(s) y(s) y(s) x(s)

=x(s)2+y(s)2 = 1 となり、[e1(s)e2(s)]SO(2)がわかる。

再びRnの話に戻り、Rnの運動群を定める。Rnの向きを変えない等長変換は RSO(n)vRnによって次のように表せる。

Rn Rn; x7→Rx+v この作用を行列の積で実現しようとすると、

[ R v

0 1 ] [

x 1 ]

= [

Rx+v 1

]

とできる。この表現を使って、Rnの運動群M(Rn) M(Rn) =

{[

R v 0 1

]R SO(n),v Rn }

によって定める。RSO(n)v Rnの定める等長変換とR SO(n)v Rn の定める等長変換の合成は

x7→R(Rx+v) +v=RRx+Rv+v

となるので、RRRv+vの定める等長変換になっている。他方、これら等長 変換に対応するn+ 1次正方行列の積は

[ R v

0 1 ] [

R v

0 1

]

= [

RR Rv+v

0 1

]

となり、上記のRRRv+vの定める等長変換に対応するn+ 1次正方行列に 一致する。したがって、M(Rn)の積は等長変換の合成とみなしても、n+ 1次正 方行列の積とみなしてもよい。

M(Rn)の元のRn成分を対応させる写像を π0 :M(Rn)Rn;

[ R v

0 1 ]

7→v

(5)

で表すことにする。

平面曲線の議論に戻る。平面曲線p(s)とその動標構e1(s),e2(s)を合わせたもの

˜ p(s) =

[

e1(s) e2(s) p(s)

0 0 1

]

(sI)

は定義域IからM(R2)への写像になりp=π0p˜を満たす。このような性質を持 p˜pの持ち上げと呼ぶ。

平面曲線の持ち上げを考える理由は,曲率と運動群の作用との関係を表しやす くするためである。p(s)˜ を微分すると

˜ p(s) =

[

e1(s) e2(s) p(s)

0 0 0

]

= [

κ(s)e2(s) κ(s)e1(s) e1(s)

0 0 0

]

= [

e1(s) e2(s) p(s)

0 0 1

] 

0 κ(s) 1

κ(s) 0 0

0 0 0

= ˜p(s)

0 κ(s) 1

κ(s) 0 0

0 0 0

となり、

˜

p(s) = ˜p(s)

0 κ(s) 1

κ(s) 0 0

0 0 0

(1.2)

を得る。そこで

K(s) =

0 κ(s) 1

κ(s) 0 0

0 0 0

(1.3)

と書くことにすると、(1.2)

˜

p(s) = ˜p(s)K(s) (1.4)

となり、この両辺にp(s)˜ 1を左から掛けると、

˜

p(s)1p˜(s) = K(s) (1.5)

を得る。

もし、同じ開区間Iで定義された平面曲線p1(s)p2(s)が,RSO(2)v R2 によって

p2(s) =Rp1(s) +v (s I) を満たすとき、R,vから定まるM(R2)の元を

A= [

R v 0 1 ]

(1.6)

(6)

で表すと、˜p2(s) = p1(s)が成り立つ。これは次のようにわかる。まず、p1(s) p2(s)の関係式の両辺をsで微分すると、p1(s)の速度ベクトルにRを作用させる p2(s)の速度ベクトルに一致する。これらを反時計回りに90度回転したベクト ルについても同様である。よって、˜p2(s) = p1(s)が成り立つことがわかる。こ の等式から

˜

p2(s)1p˜2(s) = (A˜p1(s))1(A˜p1(s)) = ˜p1(s)1A1p1(s) = ˜p1(s)1p˜1(s) がわかり、(1.5)よりp1(s)p2(s)の曲率が等しいと結論できる。

逆に,p1(s)p2(s)が同じ曲率κ(s)を持てば、

˜

p1(s) = ˜p1(s)K(s), p˜2(s) = ˜p2(s)K(s)

が成り立つ。この後、n次正則行列に値を持つ関数g(s)に対して、g(s)1の微分 を計算する必要があるので、これを求めておく。

補題 1.1.2 n次正則行列に値を持つ関数g(s)に対して次が成り立つ。

d

dsg(s)1 =g(s)1g(s)g(s)1. (1.7) 証明 g(s)1g(s) = Inの両辺をsで微分すると、積の微分の公式より

( d

dsg(s)1 )

g(s) +g(s)1g(s) = 0

となり、左辺の第二項を移項して両辺に右からg(s)1をかけると(1.7)を得る。

(1.7)を使って計算すると

d ds

(p˜2(s)˜p1(s)1)

= ˜p2(s)˜p1(s)1+ ˜p2(s) d

dsp˜1(s)1

= ˜p2(s)K(s)˜p1(s)1p˜2(s)˜p1(s)1p˜1(s)˜p1(s)1

= ˜p2(s)K(s)˜p1(s)1p˜2(s)˜p1(s)1p˜1(s)K(s)˜p1(s)1

= ˜p2(s)K(s)˜p1(s)1p˜2(s)K(s)˜p1(s)1 = 0

がわかり、˜p2(s)˜p1(s)1 sに依存しない一定のAM(R2)に一致する。すなわ ち、˜p2(s)˜p1(s)1 =Aである。この両辺に右からp˜1(s)をかけるとp˜2(s) = p1(s) となり、(1.6)と同様にAを記述すると、π0p˜2(s) =π0p1(s)より

p2(s) =Rp1(s) +v が成り立つ。以上の議論をまとめると次の定理を得る。

定理 1.1.3 二つの平面曲線が同じ曲率を持つ必要十分条件は、二つの曲線が運動

群の作用で写り合うことである。

(7)

開区間I上の関数κ(s)に対して、κ(s)を曲率として持つ曲線を構成する。(1.1) を未知関数[e1(s)e2(s)]に関する常微分方程式とみなすと、(1.1)は線形常微分方 程式である。線形常微分方程式の一般論から、任意の初期条件に対して(1.1)の解 Iにおいて一意的に存在することが知られているが、この場合は以下のように 直接解を記述できる。その記述に必要になる回転行列の性質をまず示しておく。

R(θ) = [

cosθ sinθ sinθ cosθ

]

によって回転行列R(θ)SO(2)を定める。θsの関数の場合には d

dsR(θ) =

[sinθds cosθds cosθds sinθds

]

= [

cosθ sinθ sinθ cosθ

] [

0 ds

ds 0

]

=R(θ) [

0 ds

ds 0

]

が成り立つ。したがって、常微分方程式(1.1)は、常微分方程式

ds =κ(s)

に帰着する。これは最も簡単な常微分方程式であり、初期条件θ(s0) = θ0, s0 I に対して

θ(s) =

s s0

κ(u)du+θ0

が一意的な解である。これより、R(θ(s))は初期条件R(θ(s0)) = R(θ0)に対する (1.1)の一意的な解になる。R(θ(s)) = [e1(s) e2(s)]によってe1(s), e2(s)を定め ると、

e1(s) = [

cosθ(s) sinθ(s) ]

=

cos(∫s

s0κ(u)du+θ0 ) sin(∫s

s0κ(u)du+θ0 )

となる。v0 R2に対して

p(s) =

s

s0cos(∫s

s0κ(u)du+θ0 )

s ds

s0sin(∫s

s0κ(u)du+θ0 )

ds

+v0

とおくと、p(s)は曲率κ(s)の平面曲線であることがわかる。

上で構成したp(s)からp(s)˜ を定めると、˜p(s)(1.4)を満たすこともわかる。

すなわち、(1.4)を未知関数p(s)˜ に関する常微分方程式とみなしたとき、任意の初 期条件に対して(1.4)の解を構成したことになる。

(8)

1.2 線形群

前節ではM(R2)への写像を利用して、平面曲線の曲率とM(R2)の作用との関 係を明らかにし、曲率が平面曲線の一意性と存在を支配していることを示した。こ の方法を曲面の場合に適用するために、行列の成す群への写像の一意性と存在の 条件を考える。

n次実正方行列の全体M(n,R)は自然にRn2 と同一視でき、この同一視によっ Rn2 の標準的な位相からM(n,R)の位相を定める。

GL(n,R) ={g M(n,R)|detg ̸= 0}

によって一般線形群GL(n,R)を定める。行列式は行列の成分の多項式になり、特 に連続関数になる。よってGL(n,R)M(n,R)の開集合である。GL(n,R)の部 分群が閉集合になっているとき、線形群と呼ぶ。特に、GL(n,R)GL(n,R)の部 分群であり閉集合なので、GL(n,R)は線形群である。

SO(n)は線形群であることを示す。SO(n)の定義より

SO(n) ={g GL(n,R)|tgg =In, detg = 1}

である。そこで、φ(g) = tgg(g GL(n,R))によって写像φ:GL(n,R)GL(n,R) を定めると、φ(g)gの成分の二次式になり、特に連続写像である。したがって、

SO(n) =φ1(In)det1(1)GL(n,R)の閉集合になり、SO(n)は線形群である。

M(Rn)も線形群であることを示す。

A(Rn) = {[

R v 0 1

]R GL(n,R), v Rn }

とおくと、GL(n+ 1,R)の元の第n+ 1行が[0 1]に固定されたもの全体になって いるので、A(Rn)GL(n+ 1,R)の閉集合である。M(Rn)の場合と同様に、Rn のアフィン変換x7→Rx+v の合成は対応するA(Rn)の行列としての積に対応す る。特に、行列の積によって群になる。さらに上の形よりA(Rn)は位相空間とし GL(n,R)×Rnと同相になる。M(Rn)SO(n)×Rnと同相になり、SO(n) GL(n,R)の閉集合なので、M(Rn)A(Rn)の閉集合になる。したがって、M(Rn) GL(n+ 1,R)の閉集合になり、M(Rn)は線形群である。

定理 1.2.1 線形群GGL(n,R)の単位元Inを通る曲線 c: (ε, ε)G (c(0) = In) 0における速度ベクトルc(0)の全体をgで表す。

[X, Y] =XY Y X (X, Y M(n,R))

によってM(n,R)の二項演算[X, Y]を定める。すると、gM(n,R)の部分ベクト ル空間になり、X, Y gに対して[X, Y]gが成り立つ。gGのリー環と呼ぶ。

(9)

証明 c(t) =Intで微分すると0になり、0gである。X gに対して、In

を通る曲線

c: (ε, ε)G (c(0) =In)

が存在してc(0) =Xとなる。r Rに対してd :t 7→c(rt)Inを通るGの曲線 になり、d(0) =rc(0) =rXgに含まれる。X, Y gに対してInを通るGの曲 x(t), y(t)であってx(0) =X, y(0) =Y を満たすものをとる。GGL(n,R) 部分群なのでx(t)y(t)Inを通るGの曲線になる。さらに、

d

dtx(t)y(t) t=0

= d dtx(t)

t=0

y(0) +x(0) d dty(t)

t=0

=x(0) +y(0) =X+Y となるので、X+Y gが成り立つ。したがって、gM(n,R)の部分ベクトル空 間である。

g GX gに対してgXg1 gが成り立つことを次に示す。Inを通る曲線 c: (ε, ε)G (c(0) =In, c(0) =X)

をとる。gc(t)g1Gの曲線になり、gc(0)g1 =Inが成り立つ。よって d

dtgc(t)g1 t=0

=gXg1 g

が成り立つ。これを使ってX, Y gに対して[X, Y] gが成り立つことを示す。

c(t) G, Y gより、c(t)Y c(t)1 gが成り立つ。積の微分の公式と(1.7)を使 うと

d

dtc(t)Y c(t)1 t=0

= d dtc(t)

t=0

Y c(0)1+c(0)Y d

dtc(t)1 t=0

=XY Y X = [X, Y]

を得る。c(t)Y c(t)1M(n,R)内の部分ベクトル空間gの曲線なので、[X, Y]g が成り立つ。

実は線形群G GL(n,R)は、GL(n,R)の部分多様体になることが知られてい る。この結果は認めることにする。上で定義したGのリー環は部分多様体として Gの単位元Inにおける接ベクトル空間TIn(G)に他ならない。

定理 1.2.2 SO(n)のリー環so(n)

so(n) = {XM(n,R)|tX+X = 0}.

証明 SO(n)の単位元を通る曲線c : (ε, ε) SO(n) に対してtc(t)c(t) = In が成り立つ。この両辺をtで微分すると

d dt

tc(t) t=0

c(0) +tc(0) d dtc(t)

t=0

= 0

(10)

を得る。行列に対して微分する操作と転置行列をとる操作は可換なので、

t( d dtc(t)

t=0

) + d

dtc(t) t=0

= 0 となり、

so(n)⊂ {X M(n,R)|tX+X = 0} を得る。逆にtX+X = 0を満たすX M(n,R)をとる。

c(t) = (

In+ t 2X

) ( In t

2X )1

とおく。c(t)0の近傍で定義されたSO(n)の曲線になり、

c(0) =In, c(0) = X

を満たすことを以下で示す。c(0) = Inc(t)の定め方よりわかり、c(t)0の近 傍で定義できるGL(n,R)の曲線になる。

tc(t)c(t) =

t( In t

2X

)1t( In+ t

2X ) (

In+ t 2X

) ( In t

2X )1

= (

In+ t 2X

)1( In t

2X ) (

In+ t 2X

) ( In t

2X )1

((

In t 2X

)

( In+ t

2X )

は可換なので )

= (

In+ t 2X

)1( In+ t

2X ) (

In t 2X

) ( In t

2X )1

=In

となりc(t)は直交行列になる。よって、detc(t) = ±1であるが、detc(t)tに関 して連続でありdetc(0) = 1なので、detc(t) = 1が成り立つ。これらより、c(t) Inを通るSO(n)の曲線になる。積の微分の公式と(1.7)を使うと

d dtc(t)

t=0

= d dt

( In+ t

2X) t=0

+ d dt

( In t

2X )1

t=0

= 1 2X+1

2X =X.

以上より

so(n)⊃ {X M(n,R)|tX+X = 0} を得る。したがって、次の等式を得る。

so(n) ={X M(n,R)|tX+X = 0}.

参照

関連したドキュメント

俗曲藝》第 70 期,1991)、曹琳《覡教傳人-南通童子胡錫蘋》(《民俗曲 藝》第

(石川県立松任農業高校教諭)

第I章 文献曲二研究目的       2)妊娠第4月末期婦人原尿注射成種

Internal finishing of cooling channel in molding die with free abrasive grains Effects of face protuberance on internal face in curvature channels.. Tatsuaki FURUMOTO, Daiki

本報 で は,布 の曲 げ振 動特 性 を支配 す る摩擦 機構 *金 沢大学教育学部,金 沢市 角間町, Faculty Member,同... Models of yarn contacting at the cross-—over point

Stress-strain curves of Spandex 46.7 tex yarn A : Tensile from natural length B : Tensile from 157mN in initial loading aNominal stress-nominal strain curves bActual

(Robertson, Sanders, Seymour, Thomas,

Van de Ven, Compact Complex Surfaces (second enlarged edition), Ergebnisse der Mathematik und ihrer Grenzgebiete (3), 4, Springer-Verlag, 2004..