【SPONSORED LINK】

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

3次元ユークリッド空間\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の公式]を導出します.

【SPONSORED LINK】

 準備

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

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

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

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

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

  • 次の\anb{a,b}ab内積 (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 bab外積 (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*}

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

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

  • aともbとも直交する.
  • abのなす角を\thetaとするとき,\|a\times b\|=\|a\|\|b\|\sin\thetaが成り立つ.

微分の定義

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

[定義2] パラメータ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*}

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

いくつかの補題

ここでは準備として[補題3],[補題4]を示します.

[補題3] パラメータ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*}

が従う.

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

[補題4] パラメータ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が定数であることと,[補題3]から

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

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

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

Frenet-Serretの公式

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

  1. rC^\infty級である.
  2. 任意のt\in(\alpha,\beta)に対して,r'(t)\neq0である.
  3. 任意のt\in(\alpha,\beta)に対して,r''(t)\times r'(t)\neq0である.

条件1は「rtに関して何回でも微分可能」,条件2は「rの速度は0にならない」,条件3は「曲線Cの軌跡は曲がっている」ということを述べています.

3つのベクトル

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

[定義5] 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の「ねじれる向き」

に相当します.

ここで,次の[命題6]を示します.

[命題6] {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)の両方と直交する.

Rendered by QuickLaTeX.com

\|v_1(t)\|=1は定数だから[補題4]より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は定数だから[補題4]より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によらず平行である.

[命題6]より,次の[定義7]ができる.

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

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

したがって,同一平面上を動く曲線の捩率は0ですね.

[命題8] 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*}

は[命題6]の証明中で示した.

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*}

[命題8]と同じく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)は定数なので,[補題4]より\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*}

である.また,[命題8]より\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*}

である.[命題8]より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の公式]が従う.

参考文献

シェアする

  • このエントリーをはてなブックマークに追加

フォローする

最後までありがとうございました!

以下の関連記事もいかがですか?

SPONSORED LINK
関連記事

記事一覧はこちらからどうぞ!

記事一覧は

こちら

Twitterを

フォロー

大学院入試

解答例

大学受験

解説ブログ