フルネ-セレの公式の導出｜曲線に関するベクトルの関係式

３次元ユークリッド空間$\R^3$上の$t$をパラメータとする滑らかな曲線$C:r=r(t)$に対して

「進む向き」を表す接ベクトル$v_1(t)$
「曲がる向き」を表す法線ベクトル$v_2(t)$
「ねじれる向き」を表す従法線ベクトル$v_3(t)$

を考えることができます．

このときの$[v_1(t),v_2(t),v_3(t)]$と$[{v_1}'(t),{v_2}'(t),{v_3}'(t)]$との関係を[Frenet(フルネ)-Serret(セレ)の公式]といいます．

[Frenet-Serretの公式]は

1847年にジャン・フレデリック・フルネ(Jean Frédéric Frenet)によって
1851年にジョセフ・アルフレッド・セレ(Joseph Alfred Serret)によって

それぞれ独立に発見されました．

この記事では，[Frenet-Serretの公式]を導出します．

1 準備
2 Frenet-Serretの公式
- 2.1 ３つのベクトル
- 2.2 Frenet-Serretの公式
3 参考文献

準備

パラメータ$t$はある開区間$(\alpha,\beta)$ ($\alpha$, $\beta$は実数)上を動くとします．

$\R^3$の元は太字で書くこともありますが，この記事では太字にはしていません．また，ゼロも零ベクトルもどちらも０で表しますが，混乱のおそれはないでしょう．

内積，ノルム，外積の定義

以下で$\R^3$上の内積，ノルム，外積を定義します．

[定義１]　$a=\bmat{a_1\\a_2\\a_3},b=\bmat{b_1\\b_2\\b_3}\in\R^3$に対し，次を定義する．

次の$\anb{a,b}$を$a$と$b$の内積 (inner product)という：
$\begin{align*} \anb{a,b}:=a_1b_1+a_2b_2+a_3b_3. \end{align*}$
次の$\|a\|$を$a$のノルム (norm)という：
$\begin{align*} \|a\|:=\sqrt{\anb{a,a}}=\sqrt{{a_1}^2+{a_2}^2+{a_3}^2}. \end{align*}$
次の$a\times b$を$a$と$b$の外積 (outer product)という：
$\begin{align*} a\times b:=\bmat{a_2b_3-a_3b_2\\a_3b_1-a_1b_3\\a_1b_2-a_2b_1}. \end{align*}$

定義から，内積と外積に関して

$\begin{align*} \anb{a,b}=\anb{b,a},\quad a\times b=-b\times a \end{align*}$

が成り立つことが分かります．また，$a$と$b$が直交するとは，$\anb{a,b}=0$となることをいいます．

さらに，外積$a\times b$は以下を満たします(証明略)．

$a$とも$b$とも直交する．
$a$と$b$のなす角を$\theta$とするとき，$\|a\times b\|=\|a\|\|b\|\sin\theta$が成り立つ．

微分の定義

次に，以下で$\R^3$上の導関数を定義します．

[定義２]　パラメータ$t$に関する$\R^3$上の元$p(t)=\bmat{p_1(t)\\p_2(t)\\p_3(t)}$に関して，各成分$p_1(t),p_2(t),p_3(t)$が$t$に関して微分可能であるとする．このとき，$p(t)$は微分可能であるといい，$p(t)$の$t$に関する導関数$p'(t)$を次で定義する：

$\begin{align*} p'(t)=\bmat{{p_1}'(t)\\{p_2}'(t)\\{p_3}'(t)} \end{align*}$

すなわち，ベクトルの微分は各成分で微分したものと定める．

いくつかの補題

ここでは準備として[補題３]，[補題４]を示します．

[補題３]　パラメータ$t\in(\alpha,\beta)$に関する$\R^3$上の元$p(t),q(t)$が微分可能であるとする．このとき，次が成り立つ：

$\begin{align*} &\frac{d}{dt}\anb{p(t),q(t)}=\anb{p'(t),q(t)}+\anb{p(t),q'(t)}, \\&\frac{d}{dt}(p(t)\times q(t))=p'(t)\times q(t)+p(t)\times q'(t) \end{align*}$

$p(t)=\bmat{p_1(t)\\p_2(t)\\p_3(t)}$, $q(t)=\bmat{q_1(t)\\q_2(t)\\q_3(t)}$とする．計算により

$\begin{align*} &\od{}{t}\anb{p(t),q(t)} \\=&\od{}{t}\{p_1(t)q_1(t)+p_2(t)q_2(t)+p_3(t)q_3(t)\} \\=&\sum_{i=1}^{3}\bra{{p_i}'(t)q_i(t)+p_i(t){q_i}'(t)} \\=&\sum_{i=1}^{3}{p_i}'(t)q_i(t)+\sum_{i=1}^{3}p_i(t){q_i}'(t) \\=&\anb{p'(t),q(t)}+\anb{p(t),q'(t)} \end{align*}$

と，

$\begin{align*} &\od{}{t}(p(t)\times q(t)) \\=&\od{}{t}\bmat{p_2(t)q_3(t)-p_3(t)q_2(t)\\p_3(t)q_1(t)-p_1(t)q_3(t)\\p_1(t)q_2(t)-p_2(t)q_1(t)} \\=&\bmat{{p_2}'(t)q_3(t)+p_2(t){q_3}'(t)-{p_3}'(t)q_2(t)-p_3(t){q_2}'(t)\\ {p_3}'(t)q_1(t)+p_3(t){q_1}'(t)-{p_1}'(t)q_3(t)-p_1(t){q_3}'(t)\\ {p_1}'(t)q_2(t)+p_1(t){q_2}'(t)-{p_2}'(t)q_1(t)-p_2(t){q_1}'(t)} \\=&\bmat{{p_2}'(t)q_3(t)-{p_3}'(t)q_2(t)\\{p_3}'(t)q_1(t)-{p_1}'(t)q_3(t)\\{p_1}'(t)q_2(t)-{p_2}'(t)q_1(t)} +\bmat{p_2(t){q_3}'(t)-p_3(t){q_2}'(t)\\p_3(t){q_1}'(t)-p_1(t){q_3}'(t)\\p_1(t){q_2}'(t)-p_2(t){q_1}'(t)} \\=&p'(t)\times q(t)+p(t)\times q'(t) \end{align*}$

が従う．

この[補題３]から，内積，外積の導関数に関して積の微分公式(のようなもの)が成り立つことが分かりますね．

[補題４]　パラメータ$t$に関する$\R^3$の元$p(t)$が微分可能で$\|p(t)\|$が$t$によらず常に一定であるとする．このとき，$p(t)$と$p'(t)$は$t$によらず常に直交する．すなわち，次が成り立つ：

$\begin{align*} \anb{p(t),\ p'(t)}=0 \end{align*}$

$t$によらない定数$C>0$を

$\begin{align*} C^2:=\|p(t)\|^2=\anb{p(t),p(t)} \end{align*}$

で定める．$C^2=\anb{p(t),p(t)}$の両辺を$t$で微分すると，左辺$C^2$が定数であることと，[補題３]から

$\begin{align*} 0=&\anb{p'(t),p(t)}+\anb{p(t),p'(t)} \\=&2\anb{p(t),p'(t)} \end{align*}$

が分かる．よって，両辺を２で割って$0=\anb{p(t),p'(t)}$が従う．

この[補題４]は[Frenet-Serretの公式]の導出でキーとなる補題です．

Frenet-Serretの公式

$t\in(\alpha,\beta)$をパラメータとする$\R^3$内の曲線$C:r=r(t)$を考えます．ただし，$r$は次の条件を満たしているとします．

$r$は$C^\infty$級である．
任意の$t\in(\alpha,\beta)$に対して，$r'(t)\neq0$である．
任意の$t\in(\alpha,\beta)$に対して，$r”(t)\times r'(t)\neq0$である．

条件１は「$r$は$t$に関して何回でも微分可能」，条件２は「$r$の速度は0にならない」，条件３は「曲線$C$の軌跡は曲がっている」ということを述べています．

３つのベクトル

$r$に対し，接線ベクトル，主法線ベクトル，従法線ベクトルを次のように定義します．

[定義５]　$v_1(t),v_2(t),v_3(t)\in\R^3$を

$\begin{align*} &v_1(t):=\frac{r'(t)}{\|r'(t)\|}, \\&v_2(t):=\frac{{v_1}'(t)}{\|{v_1}'(t)\|}, \\&v_3(t):=v_1(t)\times v_2(t) \end{align*}$

で定義し，$v_1(t)$を$C$の接線ベクトル，$v_2(t)$を$C$の主法線ベクトル，$v_3(t)$を$C$の従法線ベクトルという．また，$\kappa(t):=\|{v_1}'(t)\|$を$C$の曲率という．

なお，冒頭で書いた通り

接ベクトル$v_1(t)$は曲線$C$の「進む向き」
法線ベクトル$v_2(t)$は曲線$C$の「曲がる向き」
従法線ベクトル$v_3(t)$曲線$C$の「ねじれる向き」

に相当します．

ここで，次の[命題６]を示します．

[命題６]　${v_3}'(t)$と$v_2(t)$は$t$によらず平行である．

$v_3(t)=v_1(t)\times v_2(t)$より，両辺$t$で微分して

$\begin{align*} {v_3}'(t)={v_1}'(t)\times v_2(t)+v_1(t)\times{v_2}'(t) \end{align*}$

が成り立つ．また，${v_1}'(t)$と$v_2(t)$は平行だから${v_1}'(t)\times v_2(t)=0$なので，

$\begin{align*} {v_3}'(t)=v_1(t)\times{v_2}'(t) \end{align*}$

となって，${v_3}'(t)$は$t$によらず$v_1(t)$, ${v_2}'(t)$の両方と直交する．

$\|v_1(t)\|=1$は定数だから[補題４]より$v_1(t)$と${v_1}'(t)$は$t$によらず直交し，さらに$v_2(t)=\dfrac{{v_1}'(t)}{\|{v_1}'(t)\|}$だから，$v_1(t)$と$v_2(t)$は$t$によらず直交する．

また，$\|v_2(t)\|=1$は定数だから[補題４]より$v_2(t)$と${v_2}'(t)$は$t$によらず直交する．

したがって，${v_2}'(t)$は$t$によらず$v_1(t)$, $v_2(t)$の両方と直交する．

いまは$R^3$で考えているので，${v_3}'(t)$と$v_2(t)$は$t$によらず平行である．

[命題６]より，次の[定義７]ができる．

[定義７]　${v_3}'(t)=-\tau(t)v_2(t)$によって定まる$\tau(t)$を$C$の捩率れいりつという．

捩率は「接線ベクトル，主法線ベクトルの張る平面からどれくらいの勢いではみ出そうとするのか」ということを表します．

したがって，同一平面上を動く曲線の捩率は０ですね．

[命題８]　$t$によらず，$A(t):=[v_1(t),v_2(t),v_3(t)]$は直交行列である．

$\anb{v_i(t),v_j(t)}=\delta_{i,j}$ $(i,j=1,2,3)$を示せば良いが，内積は可換だから$i\ge j$の場合を示せば十分である．

$v_1(t),v_2(t)$の定義から

$\begin{align*} \|v_1(t)\|=\|v_2(t)\|=1 \end{align*}$

である．また，定義$v_3(t)=v_1(t)\times v_2(t)$から

$\begin{align*} &\|v_3(t)\|=\|v_1(t)\|\|v_2(t)\|\sin\frac{\pi}{2}=1, \\&\anb{v_3(t),v_1(t)}=\anb{v_3(t),v_2(t)}=0 \end{align*}$

である．さらに，

$\begin{align*} \anb{v_2(t),v_1(t)}=0 \end{align*}$

は[命題６]の証明中で示した．

Frenet-Serretの公式

準備が整ったので，[Frenet-Serretの公式]を示します．

[Frenet-Serretの公式]　次の等式が成り立つ．

$\begin{align*} [{v_1}'(t),{v_2}'(t),{v_3}'(t)] =[v_1(t),v_2(t),v_3(t)]\bmat{0&-\kappa(t)&0\\\kappa(t)&0&-\tau(t)\\0&\tau(t)&0} \end{align*}$

[命題８]と同じく$A(t):=[v_1(t),v_2(t),v_3(t)]$とする．

$\begin{align*} {}^{t}A(t)A'(t) =&\bmat{{}^{t}v_1(t)\\{}^{t}v_2(t)\\{}^{t}v_3(t)}[{v_1}'(t),\ {v_2}'(t),\ {v_3}'(t)] \\=&\bmat{\anb{v_1(t),{v_1}'(t)}&\anb{v_1(t),{v_2}'(t)}&\anb{v_1(t),{v_3}'(t)}\\ \anb{v_2(t),{v_1}'(t)}&\anb{v_2(t),{v_2}'(t)}&\anb{v_2(t),{v_3}'(t)}\\ \anb{v_3(t),{v_1}'(t)}&\anb{v_3(t),{v_2}'(t)}&\anb{v_3(t),{v_3}'(t)}} \end{align*}$

である．$\|v_i\|=1$ ($i=1,2,3$)は定数なので，[補題４]より$\anb{v_i(t),{v_i}'(t)}=0$である．また，${v_3}'(t):={v_1}'(t)\times v_2(t)+v_1(t)\times {v_2}'(t)$であり

$\begin{align*} {v_1}'(t)\times v_2(t)={v_1}'(t)\times \frac{{v_1}'(t)}{\|{v_1}'(t)\|}=0 \end{align*}$

であることと，

$\begin{align*} \anb{v_1(t)\times {v_2}'(t),v_1(t)}=0 \end{align*}$

であることから$\anb{v_1(t),{v_3}'(t)}=0$が成り立つ．さらに，

$\begin{align*} v_3(t):=v_1(t)\times v_2(t)=v_1(t)\times \frac{{v_1}'(t)}{\|{v_1}'(t)\|} \end{align*}$

だから，$\anb{v_3(t),{v_1}'(t)}=0$である．次に，

$\begin{align*} &\anb{v_2(t),{v_1}'(t)}=\anb{\dfrac{{v_1}'(t)}{\|{v_1}'(t)\|},{v_1}'(t)}=\|{v_1}'(t)\|=\kappa(t), \\&\anb{v_2(t),{v_3}'(t)}=\anb{v_2(t),-\tau(t)v_2(t)}=-\tau(t) \end{align*}$

である．また，[命題８]より$\anb{v_1(t),v_2(t)}=0,\anb{v_2(t),v_3(t)}=0$だったから，これらの両辺を$t$で微分して整理すると，

$\begin{align*} &\anb{v_1(t),{v_2}'(t)}=-\anb{{v_1}'(t),v_2(t)}=-\kappa(t), \\&\anb{{v_2}'(t),v_3(t)}=-\anb{v_2(t),{v_3}'(t)}=\tau(t) \end{align*}$

である．[命題８]より$A(t)$は直交行列なので，${A(t)}^{-1}={}^{t}A(t)$である．よって，

$\begin{align*} A(t)^{-1}A'(t)=\bmat{0&-\kappa(t)&0\\\kappa(t)&0&-\tau(t)\\0&\tau(t)&0} \end{align*}$

が成り立つから，

$\begin{align*} A'(t)=A(t)\bmat{0&-\kappa(t)&0\\\kappa(t)&0&-\tau(t)\\0&\tau(t)&0} \end{align*}$

となって，[Frenet-Serretの公式]が従う．

参考文献

「曲線と曲面の微分幾何」(小林昭七著，裳華房)
- Amazon
- 楽天市場
「線型代数入門」(斎藤正彦著，東京大学出版会)
- Amazon
- 楽天市場