包絡線
曲線 E が、媒介変数 m を含む曲線群 F={ F( x,y,m )=0 }の包絡線であるとは、
Eの各点が、どれかの曲線 F( x,y,m )=0 上にあり、かつ、その点で2曲線が
互いに接している
ときを言う。( → 参考 )
例 曲線群 F={ y=mx+m2 ( m は実数) }の包絡線を求めよ。
以前は、手作業でいくつかの曲線を描き、その包絡線らしきものを想像していたが、最
近は、パソコンで簡単に残像つきのグラフが描けるようになったので、包絡線そのものが
視覚的にはっきりと認識できるようになった。
左図においては、放物線(赤色!)
y=−(1/4)x2
が包絡線となる。
上記のように、グラフ(幾何)的に求めた包絡線を、代数的、解析的に求めようとするの
が、このページの主テーマである。
包絡線の理論を整理しておこう。
F( x , y , m )=0 が包絡線を持つとする。このとき、この曲線と包絡線との接点の座
標は、m の関数で、 x=φ(m)、y=ψ(m) と書ける。
このとき、曲線 x=φ(m)、y=ψ(m) が包絡線の方程式である。
曲線 F( x , y , m )=0 と包絡線 x=φ(m)、y=ψ(m) は接点で共通接線を持
つことから、その点で、Fx 、Fy は同時に0になることはない。
そこで、 Fy≠0 としても一般性を失うことはない。
曲線 F( x , y , m )=0 上の接点における微分係数は、m は定数として、
Fx +Fy (dy/dx)=0 より、 dy/dx=−Fx/Fy
また、包絡線 x=φ(m)、y=ψ(m) 上の接点における微分係数は、
dy/dx=(dy/dm)/(dx/dm)
この両者が一致するので、 −Fx/Fy=(dy/dm)/(dx/dm) から、
Fx (dx/dm)+Fy (dy/dm)=0
が成り立つ。一方、F( x , y , m )=0 を m で偏微分すると、
Fx (dx/dm)+Fy (dy/dm)+Fm=0
が得られる。
このとき、 Fx (dx/dm)+Fy (dy/dm)=0 なので、 Fm=0 が成り立つ。
以上から、
包絡線の方程式を求めるには、F( x , y , m )=0 、Fm( x , y , m )=0
から、 m を消去すればよい。
上記のような計算を理解すれば、偏微分が出ることの必然性は疑う余地はないが、初
学者にとっては突拍子もないことと感じるかもしれない。
別な視点で、このいわれをのぞいておこう。
いま、円群 (x−α)2+y2=r2 を考え、その包絡線を求めたい。
包絡線を求めるために、ごく近い2つの円を考える。
(x−α)2+y2=r2 、 (x−α+ε)2+y2=r2
これらの交点を求めるには、この連立方程式を解けばよい。
すなわち、 (x−α)2+y2=r2 、 2ε(x−α)+ε2=0
ε≠0 としてよいので、 (x−α)2+y2=r2 、 2(x−α)+ε=0
ここで、 εを限りなく0に近づけると、2つの円の交点は、包絡線上の点に近づく。
すなわち、包絡線上の点は、連立方程式
(x−α)2+y2=r2 、 2(x−α)=0
を満たす。
この式は、 F( x , y , α )=(x−α)2+y2−r2=0 として、
Fα( x , y , α )=2(x−α)=0
としたものに等しい。
このように考えれば、包絡線を求めるのに偏微分が用いられることが実感できることと
思われる。
因みに、円群の包絡線の方程式は、 第2式から、x=α で、これを第1式に代入して
y=±r で、円群を頭の中で「ググ〜ッ!」と動かして得られるイメージと一致する。
さて、包絡線の公式を用いて、先の例における包絡線を求めてみよう。
(解) F( x , y , m )=mx−y+m2 とおくと、 Fm( x , y , m )=x+2m なので、
m=−x/2 を、 mx−y+m2=0 に代入して、 (−x/2)x−y+(−x/2)2=0
よって、 包絡線の方程式は、 y=−(1/4)x2 となる。 (終)
ところで、曲線群 F={ y=mx+m2 ( m は実数) }の包絡線は、y=−(1/4)x2 であ
ることが示されたが、この式は、m の2次方程式 m2+mx−y=0 が重解を持つ、すな
わち、判別式 D=0 を計算しても得られることに気づかされる。
実際に、F が m についての多項式関数のとき、
「 F=0 、Fm=0 」が成り立つということは、
m の方程式として、「 F=0 」が重解を持つことに等しい
ので、特に、2次関数 F については、偏微分というものを表に出さずとも、高校数学風に判
別式を利用して包絡線の方程式を求める方が簡明だろう。
(平成20年12月18日付けで、S(H)さんから、「判別式を求めるのに偏微分を用いているのでは?」
というご指摘を頂きました。正にその通りで、「高校数学風に」と文言を修正させていただきました。)
例 曲線群 F={ (x−m)2+(y−m)2=m2 (mは実数) }の包絡線の方程式を
求めよ。
(解) m2−2(x+y)m+x2+y2=0 の判別式をDとすると、
D/4=(x+y)2−(x2+y2)=2xy=0
よって、包絡線の方程式は、 x=0 または y=0 (終)
例 曲線群 F={ y=mx−(1+m2)x2 (mは実数) }の包絡線の方程式を求めよ。
(解) x2m2−xm+x2+y=0 は、x=0 のとき重解を持たないので、x≠0 としてよい。
このとき、判別式をDとすると、 D=x2−4x2(x2+y)=0 より、 y=−x2+1/4
よって、包絡線の方程式は、 y=−x2+1/4
ただし、点( 0 , 1/4 )は省く。 (終)
曲線群 F={ F( x , y , α , β )=0 ( α 、 β は実数 ) }を表すパラメータが2個
あって、2つのパラメータ間に拘束条件 G( α , β )=0 がある場合は、1パラメータと
同様の公式が成り立つ。
すなわち、
曲線群が、2つの方程式 F(X,Y,α,β)=0 、G(α,β)=0 で与えられる
とき、包絡線上の全ての点は、次の方程式 :
F(X,Y,α,β)=0 、G(α,β)=0 、
FαGβ−GαFβ=0 (ただし、Fαは、Fのαに関する偏微分を表す。他も同様)
からパラメーターα,βを消去して得られる方程式を満たす。
実際に、 G( α , β )=0 より、 β を α の関数 β=g(α) と見て、
H( x , y , α )=F( x , y , α , g(α) )
と書ける。(必要があれば、α を β の関数と思って式を作ってもよい。)
このとき、包絡線上の点は、方程式 H( x , y , α )=0 、 Hα( x , y , α )=0
を満たす。ここで、 Fx(dx/dα)+Fy(dy/dα)=0 であることに注意して、
Hα=Fα+Fβ(dβ/dα)=0 ・・・ (*)
ところで、 G( α , β )=0 の両辺を α で偏微分すると、
Gα+Gβ(dβ/dα)=0 ・・・ (**)
このとき、 (*)、(**)から、(dβ/dα)を消去すれば、
FαGβ−GαFβ=0
が得られる。
例 α2+β2=L2 のとき、直線群 βx+αy=αβ の包絡線の方程式を求めよ。
(解) F( x , y , α , β )=βx+αy−αβ 、G(α,β)=α2+β2−L2 とおく。
FαGβ−GαFβ=(yーβ)(2β)−(2α)(x−α)=0 より、
β(y−β)−α(x−α)=0
以上から、連立方程式 αx−βy=α2−β2 、 βx+αy=αβ が得られる。
これを解いて、 L2x=α3 、L2y=β3 となる。
これらを、α2+β2=L2 に代入して、
(終)
( → 参考 : 「星芒形」、「直角に曲がる通路の通過問題」)
例 α2+β2=L2 のとき、円群 x2+y2+2αx+2βy=0 の包絡線の方程式を
求めよ。
(解) F( x , y , α , β )=x2+y2+2αx+2βy 、G(α,β)=α2+β2−L2 とおく。
FαGβ−GαFβ=(2x)(2β)−(2α)(2y)=0 より、 βx−αy=0
このとき、実数 m を用いて、 α=mx 、 β=my とおける。
α2+β2=L2 に代入して、 m2x2+m2y2=L2 より、
これらを、 x2+y2+2αx+2βy=0 に代入して整理すると、 x2+y2=4L2 (終)
(コメント) 上記の場合、原点も包絡線に含まれるとしていいんですよね?
当HPがいつもお世話になっているS(H)さんから次のような問題が提起された。
(平成20年12月13日付け)
例 放物線 y=x2 上の2点A( α , α2 )、B( β , β2) (ただし、α<β)が、
AB=L(=一定)という条件を満たしながら動く。このとき、線分ABの包絡線の方
程式を求めよ。
S(H)さんによれば、L=4 のときに、特異点を有する包絡線が出現するという。
L=4 の場合を考えてみよう。
直線 AB の方程式は、明らかに、 y=(α+β)x−αβ である。
このとき、拘束条件は、 (α−β)2+(α2−β2)2=16
すなわち、 (α−β)2{1+(α+β)2}=16 である。
ここで、 α+β=m 、 αβ=n とおくと、
直線 AB の方程式は、 y=mx−n
拘束条件は、 (m2−4n)(1+m2)=16
となる。このことから、
このとき、 x2=mx−n すなわち、 x2−mx+n=0 において、解の公式により、
であるので、
の範囲にある線分ABを変化させると、下図が得られる。(一部手書き補正!)
(コメント) あ〜、本当だ!特異点っぽい所がありますね。予想では、線分が滑らかに変化
していくので特異点は出そうにないと思っていましたが、予想が外れました...。
拘束条件が意外にきついのと放物線の傾きが滑らかに変化することの障害にな
っているんでしょうね?
S(H)さんによれば、 y=x4 に対しては、L=7 のときに、特異点を有する包絡
線が出現するという。
(追記)平成21年12月26日付け
HN「Boo」さんから、この話題についてご教示いただいた。
計算したり、グラフを描いてみたりしたところ、
・ y=x2 では、線分の長さLが2より大きい場合
・ y=x4 では、線分の長さLが2の場合
でも特異点が出そうな気がします。
とのこと。
(コメント) HN「Boo」さん、ありがとうございます。今後検証したいと思います。
上記では包絡線の様子を視覚的に眺めたが、その方程式はどのように記述されるのだろ
うか?特異点が出現しているので、このことも大いに興味があるところである。
S(H)さんからは、「描写された包絡線では見えがたい微妙な代数曲線の部分があるので、
その包絡線の代数方程式を明記し、具体的に、特異点の状況を調べた方がよい。」というア
ドバイスを頂いた。(平成20年12月15日付け)
まさにその通りで、上記の関数の真実を探りたいと思う。
上記の計算から、包絡線上の点は、次の方程式を満たす。
F( x , y , m , n )=mx−y−n=0
G( m , n )=m4+(1−4n)m2−4n−16=0
FmGn−GmFn=4m3−4m2x+2(1−4n)m−4x=0
n=mx−y を他の2式に代入して式を整理すると、
m4−4xm3+(4y+1)m2−4xm+4y−16=0
2m3−6xm2+(4y+1)m−2x=0
これらの方程式は、共通解 m を持っているので、シルベスターの消去法を用いて、m
を消すことにより、包絡線上の点が満たす方程式を求めることができる。
すなわち、
行列式の性質を用いて、左辺を計算(手計算 ... f(^^;) )すると、
−1024x6+(−2048y2+38912y−113792)x4
+(−1024y4+7168y3−51584y2+155584y−13252)x2
+1024y5−5120y4+12672y3−38464y2+33540y−67600=0
となる。
(コメント) 平成20年12月21日付けのS(H)さんからの正しい計算結果の提示がなかっ
たら、手計算ゆえに計算結果には全く自信がなく、係数の大きさから不安に煽
られ、ただ恐れおののくだけだったろう。S(H)さんからのアドバイスが神々しく
感じられた。S(H)さんに感謝します。
それにしても、上記の手計算は、とてもハードでした...!
確認のために、上記の方程式のグラフをグラフ描画ソフトに描かせてみると、
もっとも、上図は直線群 AB の包絡線である。問題では、線分群 AB なので、実際の包
絡線は下図のようになるであろうか?
当HPがいつもお世話になっているS(H)さんから包絡線に関わる新しい問題が提起され
た。(平成20年12月27日付け)
例 放物線 y=x2 において、放物線の上方の1点を、P( a , b ) とおく。点Pを
通る直線が放物線と交わる2点をA( α , α2 )、B( β , β2) とする。 (た
だし、α<β) このとき、次の問いに答えよ。
(1) 放物線と直線ABで囲まれた図形の面積の最小値 S( a , b
) を求めよ。
(2) 曲線 S( a , b )=k (k は実数の定数)は、ある直線群の包絡線とな
ることを示せ。
(1) 直線ABの方程式は、傾きを m として、 y=m(x−a)+b とおける。
(
→ 参考 : 「放物線で掛け算」)
放物線の方程式 y=x2 と連立して、 x2−mx+ma−b=0
解と係数の関係から、 α+β=m 、 αβ=ma−b
放物線と直線ABで囲まれた図形の面積は、 (β−α)3/6 で与えられる。
よって、
(β−α)2=(α+β)2−4αβ=m2−4(ma−b)=(m−2a)2−4a2+4b
より、求める面積は、m=2a のとき最小で、最小値 S( a , b )
は、
となる。
(2) k=(4/3)n3 を満たす実数 n はただ1つ存在して、曲線 S( a , b )=k は、放物
線 b=a2+n2 を表す。
(1)より、 m=2a のとき、面積は最小なので、そのときの直線ABの方程式は、
y=2ax−2a2+b となる。このとき、面積は、
(β−α)3/6=k=(4/3)n3 より、 β−α=2n となる。
このとき、(α+β)2−4αβ=4n2 に、α+β=2a 、αβ=2a2−b を代入して、
4a2−4(2a2−b)=4n2 すなわち、 −a2+b=n2
そこで、拘束条件 G(a,b)=−a2+b−n2=0 のもとに、直線群
F( x , y , a , b )=2ax−y−2a2+b=0
を考えると、
FaGb−GaFb=(2x−4a)・1−(−2a)・1=2x−2a=0
より、 a=x となる。このとき、 b=a2+n2=x2+n2
これらを、 F( x , y , a , b )=2ax−y−2a2+b=0 に代入して、
2x2−y−2x2+x2+n2=0 より、 y=x2+n2
よって、直線群の包絡線の方程式は、y=x2+n2 となる。
したがって、求める直線群は、 y=2ax−a2+n2 (ただし、a は実数) である。
(コメント) 拘束条件が包絡線そのものですが、これは勇み足?
k=1/6 のとき、 n=1/2 で、 直線群は、 y=2ax−a2+1/4 となる。
この場合を図示してみよう。
上記の計算で、
S=(β−α)3/6
の最小値を求めるために、 (β−α)2=(m−2a)2−4a2+4b という2次関数を利用し
たが、この種の問題の一般的解法のため次のような解法を確認しておくことも大切だろう。
(別解その1) 放物線 y=x2 の上方にある点( a , b ) を通る直線と放物線との交点の
x 座標を α、β (α<β) とすると、直線の方程式は、 y=(α+β)x−αβ
と書ける。このとき、 b=(α+β)a−αβ であるので、
を、Sの式に代入して、βを消去し、整理すると、
ここで、 a−α>0 、 b−a2>0 に注意して、相加平均と相乗平均の関係から、
求める最小値は、
となることが分かる。 (終)
(コメント) この(別解その1)の解法も、Sの特殊な形に依存しているので、「一般的」とは
言えないかも...?
(別解その2) 放物線 y=x2 の上方にある点( a , b ) を通る直線と放物線との交点の
x 座標を α、β (α<β) とすると、直線の方程式は、 y=(α+β)x−αβ
と書ける。このとき、 b=(α+β)a−αβ であるので、
を、Sの式に代入して、βを消去し、整理すると、
ここで、 x=α−a (<0)、m=b−a2 (>0) とおくと、
両辺を x で微分して、
S’=0 を解いて、 で、このとき、Sは極小かつ最小となり、
最小値は、
となる。 (終)
上記の計算で、陰関数の微分法を用いてもよいだろう。
(別解その3)
の分母を払って、x の降べきの順に式を整理すると、
x6+3mx4+6Sx3+3m2x2+m3=0 ・・・ (*)
両辺を x で微分して、
6x5+12mx3+6S’x3+18Sx2+6m2x=0
x<0 なので、 x4+2mx2+S’x2+3Sx+m2=0
ここで、S’=0 として、 x4+2mx2+3Sx+m2=0 ・・・ (**)
(*)(**) より、Sを消去すると、 x6+mx4−m2x2−m3=0
よって、 x2(x4−m2)+m(x4−m2)=0 より、 (x2+m)2(x2−m)=0
これを解いて、 より、Sの最小値
を得る。 (終)
上記では極値を与える候補(この場合、最小値を与えることは幾何学的に自明としていいだろう)
を求めてから最小値を求めたが、シルベスターの消去法を用いれば、そういった遠回りを
しなくても直接的に最小値は求められる。当HPがいつもお世話になっているS(H)さんは
その方法で最小値を求められた。(平成20年12月29日付け)
(S(H)さんの計算) S=(β−α)3/6 、 b−aα−aβ+αβ= 0 からβを消去し
−6Sa3+b3+18Sa2α−6ab2α−18Saα2+12a2bα2+3b2α2
+6Sα3−8a3α3−12abα3+12a2α4+3bα4−6aα5+α6=0
また、上式を α で微分して、S’=0 とおくと、
18Sa2−6ab2−36Saα+24a2bα+6b2α+18Sα2−24a3α2
−36abα2+48a2α3+12bα3−30aα4+6α5=0
シルベスターの消去法により、(← ここの部分は手計算だと大変そう...!)
−3779136S4(a2−b)6(9S2+16a6−48a4b+48a2b2−16b3)=0
−3779136S4(a2−b)6(9S2+16(a2−b)3)=0
よって、 9S2=16(b−a2)3 より、
を得る。 (終)
S(H)さんの解法を意識して、ここでは、互除法の考えを用いて手計算してみよう。
このような機会を与えていただいたS(H)さんに感謝します。
(別解その4)
の分母を払って、x の降べきの順に式を整理すると、
x6+3mx4+6Sx3+3m2x2+m3=0 ・・・ (*)
両辺を x で微分して、
6x5+12mx3+6S’x3+18Sx2+6m2x=0
x<0 なので、 x4+2mx2+S’x2+3Sx+m2=0
ここで、S’=0 として、 x4+2mx2+3Sx+m2=0 ・・・ (**)
極値を与える候補の x は、(*)(**) の2式を同時に満たす。
このとき、 x6+3mx4+6Sx3+3m2x2+m3
=(x4+2mx2+3Sx+m2)(x2+m)+3Sx3−3mSx
なので、極値を与える候補の x は、
x4+2mx2+3Sx+m2=0 、 x3−mx=0
の2式を同時に満たす。
このとき、 x4+2mx2+3Sx+m2=(x3−mx)x+3mx2+3Sx+m2
なので、極値を与える候補の x は、
x3−mx=0 、 3mx2+3Sx+m2=0
の2式を同時に満たす。シルベスターの消去法により、
すなわち、
これより、 m2{m4−9S2m+3m4+3m(3m3+m3)}=0 なので、
簡約すると、 m3(−9S2+16m3)=0 となる。
よって、求める面積の最小値は、
で与えられることが分かる。 (終)
シルベスターの消去法を用いると、どうしても高次の行列式の計算が必須になる。手計算
をする身にとっては、これは避けなければならない由々しき事態である。
現在当面している問題設定において、放物線との囲まれた面積が最小となる直線群の包
絡線が求める曲線であるという事実を知れば、次のような別解も可能となる。この解法につ
いては、平成20年12月31日付けでS(H)さんが取り組まれた。
(別解その5)
放物線と直線ABで囲まれた図形の面積は、 (β−α)3/6 で与えられる。そこで、
面積の最小値 S( a , b )=k とし、 とおくと、 β=α+K となる。
このとき、直線ABの方程式 y=(α+β)x−αβ が、点( a , b
) を通るので、
b=(2α+K)a−α(α+K) すなわち、 α2+(K−2a)α+b−Ka=0
そこで、 F( a , b , α )=α2+(K−2a)α+b−Ka=0 とすると、
Fα( a , b , α )=2α+(K−2a)=0
シルベスターの消去法により、
よって、 (K−2a)2+4(b−Ka)−2(K−2a)2=0 すなわち、 K2=4(b−a2)
したがって、
(終)
(コメント) 十分手計算の範疇まで落ちてきましたね!でも、Fα=0 が1次方程式なので、
α=−(K−2a)/2 を、F=0 に代入して、
(K−2a)2/4−(K−2a)2/2+b−Ka=0 より、 K2=4(b−a2)
としたい誘惑が....。
ちょうど、代入による消去とシルベスターの消去法の端境期ですね!
x 、y の多項式 F( x , y , m ) または F( x , y ,α,β) によって定義される方程
式 F( x , y , m )=0 または F( x , y ,α,β)=0 を満たす点(
x , y )全体の集
合は、幾何学的には「代数曲線」と呼ばれる。
高校で扱う曲線、例えば、直線、円、放物線、楕円、双曲線などは、もちろん代数曲線で
ある。
ここで、Fx =0 かつ Fy=0 であるような曲線 F=0 上の点 Pを
F の特異点という。
代数曲線は、特異点を全く持たないか、有限個の特異点(結節点、尖点、孤立点)を持
つの何れかである。
特異点の種類の判定には、ヘシアンの符号が用いられる。ヘシアンは、ヘッセの行列
式(ヤコビの行列式)とも言われる。
曲線 F( x , y )=0 が、点P( a , b )で特異点を持つとすると、その近くの曲線上
の点Q( a+h ,b+k )に対して、
F( a , b )=0 、 Fx( a , b )=0 、Fy(
a , b )=0
であることに注意して、Taylor の定理より、
F( a+h ,b+k )≒ (1/2){ Fxx(a,b)h2+2Fxy(a,b)hk+ Fyy(a,b)k2 }
F( a+h ,b+k )=0 なので、
Fxx(a,b)h2+2Fxy(a,b)hk+ Fyy(a,b)k2≒0
そこで、 A=Fxx(a,b) 、 B=Fxy(a,b) 、 C= Fyy(a,b) とおき、
ヘシアン(Hessian) : H=AC−B2 を考える。
ヘシアンは、関数 F( x , y ) に対して、( a , b ) におけるヘッセ行列
の行列式として得られる。ヘシアンは、行列 H の固有値の積に等しい。
特異点の判定公式
(1) H<0 のとき、 特異点Pは、結節点である。
実際に、H<0 だと、 Fxx(a,b)h2+2Fxy(a,b)hk+ Fyy(a,b)k2 は一次式の積
に分解できるので、点Pの近くでは曲線は近似的に2直線となる。
(2) H=0 のとき、 特異点Pは、尖点である。
実際に、H=0 だと、 Fxx(a,b)h2+2Fxy(a,b)hk+ Fyy(a,b)k2 は一次式の平
方に分解できるので、点Pの近くでは曲線は近似的に2接線が一致した1直線となる。
(3) H>0 のとき、 特異点Pは、孤立点である。
実際に、H>0 だと、 Fxx(a,b)h2+2Fxy(a,b)hk+ Fyy(a,b)k2 は定符号とな
り、点Pの近くには曲線上の点はない。
特異点の具体例をあげて、その感覚を磨いておこう。
例(結節点) r2=8cos2θ すなわち、 (x2+y2)2=8(x2−y2)
この曲線は、連珠形(レムニスケート)と言われる。( → 参考 : 「二葉線」 )
F( x , y )=(x2+y2)2−8(x2−y2) とおくと、
Fx=2(x2+y2)・2x−16x=4x(x2+y2−4)=0
Fy=2(x2+y2)・2y+16y=4y(x2+y2+4)=0
よって、原点( 0 , 0 )で、曲線 F( x , y )=0 は特異点を持つ。
ヘシアン H=FxxFyy−FxyFyx=(12x2+4y2−16)(4x2+12y2+16)−8xy・8xy
において、 H( 0 , 0 )=−256<0 なので、特異点は結節点である。
例(尖点) y2=x4
F( x , y )=x4−y2 とおくと、
Fx=4x3=0
Fy=−2y=0
よって、原点( 0 , 0 )で、曲線 F( x , y )=0
は特異点を持つ。
ヘシアン H=FxxFyy−FxyFyx=12x2・(−2)−0=−24x2
において、 H( 0 , 0 )=0 なので、特異点は尖点である。
例(孤立点) x4+y4−6x2−8y2=0
F( x , y )=x4+y4−6x2−8y2 とおくと、
Fx=4x3−12x=4x(x2−3)=0
Fy=4y3−16y=4y(y2−4)=0
よって、原点(
0 , 0 )で、曲線 F( x , y )=0
は特異点を持つ。
ヘシアン H=FxxFyy−FxyFyx=(12x2−12)(12y2−16)−0
=(12x2−12)(12y2−16)
において、 H( 0 , 0 )=192>0 なので、特異点は孤立点である。
包絡線の計算で注意しなければいけないことは、
F( x , y , m )=0 、Fm( x , y , m )=0
または、
F( x , y , α , β)=0 、G( α , β )=0 、 FαGβ−GαFβ=0
からパラメーター m または、α,βを消去して得られる方程式によって定義される曲線
には、包絡線の他に、包絡線上にはない他の点が含まれることもありうるという点だあろ
う。
すなわち、パラメータを消去して得られる曲線は、
包絡線 および 曲線の特異点の軌跡
を含む。
この事情は、先に述べた視点で考えれば当然のことと理解できる。すなわち、
ごく近い位置にある2つの曲線の交点の極限点として、包絡線上の点を捉えたが、
特異点、例えば、結節点Pの近くでも2つの曲線の交点Qが存在する。
このとき、 Q → P なので、パラメータを消去して得られる曲線に点Pが含まれる
ことになる。
したがって、特異点を1つも持たない曲線群の包絡線は、パラメータを消去して得られる
曲線そのものであると言える。
次に、包絡線の話題の応用として、縮閉線と伸開線の話題を取り上げたい。
曲線 y=F(x) において、2点P( a , F(a) )、Q( a+h , F(a+h) )を考える。
2点 P 、Q において接線を引き、2接線の x 軸の正の向きとのなす角をそれぞれ、
θ 、 θ+Δθ (Δθは、点P における接線から点Qにおける接線へ測った角)
とする。
このとき、点Pから歩き始めて道のりΔsだけ歩い
て点Qに達するとき、Δθは感覚的に曲線の曲がり
具合を示す。
そこで、単位長さ当たりの曲がり具合
を点Pにおける曲線の曲率(curvature)と言う。
曲率の絶対値の逆数は、曲率半径と言われる。
定 理 曲線 y=F(x) 上の点P( a , F(a) )における曲率 κ は、
(略証) tanθ=F’(a) 、 tan(θ+Δθ)=F’(a+h) なので、平均値の定理より、
Δθ=tan-1F’(a+h)−tan-1F’(a)=hF”(c)/(1+F’(c)2)
ただし、 c は、a と a+h の間の数
また、 Δs=h(1+F’(d)2)1/2 ただし、 d は、a と a+h の間の数
よって、 (Δθ)/(Δs)は、 h→0 のとき、κ に近づく。 (略証終)
例 原点中心、半径 r の円では、その曲率の絶対値は、 1/r に等しい。
実際に、 x2+y2=r2 より、 2x+2yy’=0 なので、 y’=−x/y 、y”=−r2/y3
このとき、 曲率κの絶対値は、 1/r となる。
例 曲線 y=log x の曲率半径の最小値を求めよ。 | |||
曲率半径を r とすると、 y=1/x 、 y”=−1/x2 より、 | |||
このとき、 | なので、 | ||
のとき最小で、最小値は、 | |||
特別な場合として、原点で x 軸に接する曲線の曲率 κ は、
で与えられる。
実際に、 曲線 y=F(x) が、原点で x 軸に接するので、 F’(0)=0
よって、原点における曲率 κ=F”(0) である。
ところで、ロピタルの定理より、
よって、上記の公式が成り立つ。
例 y=x2 は、原点において、x 軸に接する。
したがって、 y=x2 から、 y’=2x 、 y”=2 なので、 κ=2 となる。
このとき、曲率半径は、1/2 で、曲率円の中心の座標は、( 0 , 1/2 )となる。
上図からも気がつくように、曲率円は曲線の凹の側にとる。このときの曲率円の中心は
曲率中心と言われる。
上記の放物線 y=x2 の原点に対する曲率中心は、( 0 , 1/2)で、この点は、放物線
の原点における法線上にある。
曲率中心の軌跡のことを、縮閉線といい、もとの曲線のことを伸開線と言う。
曲線 y=F(x)上の点P( t , F(t) )に対して、縮閉線上の点Q( x , y )の座標は、
で与えられる。
実際に、点Pにおいて、法線と接線は垂直なので、 y’・{ y−F(t)}/(x−t)=−1
また、 PQは曲率半径なので、 (x−t)2+{ y−F(t)}2=(1+y’2)3/y”2
これより、 { y−F(t)}2=(1+y’2)2/y”2
曲率中心は曲線の凹の側にとるので、 y−F(t) と y” は同符号。
よって、 y−F(t)=(1+y’2)/y” より、 y=F(t)+(1+y’2)/y”
このとき、 x=t−y’{ y−F(t) }=t−y’(1+y’2)/y” である。
例 y=x2 上の点P( t ,t2 )に対して、 y’=2t 、 y”=2 なので、
縮閉線上の点Q(x , y )の座標は、
x=t−2t(1+4t2)/2=−4t3 y=t2+(1+4t2)/2=3t2+1/2
したがって、放物線 y=x2 の縮閉線のグラフを図示すると、
上図において、縮閉線は、点( 0 , 1/2)で特異点(尖点)を持つことが分かる。
平成21年1月9日付けで、S(H)さんは次の問題を提起された。
曲線 y=x4 の法線群に対する包絡線の方程式を多様な発想で求めよ。また、放
物線 y=x2 との相違点に留意して特異点を求め、特異点と曲線 y=x4 との最短
距離を求めよ。
この問題提起を意識しながら、曲線 y=x4 の縮閉線の方程式を求めたいと思う。
y=x4 上の点P( t ,
t4 )に対して、 y’=4t3 、 y”=12t2 なので、
縮閉線上の点Q(x , y )の座標は、
x=(2t−16t7)/3 y=(28t6+1)/(12t2)
したがって、放物線 y=x4 の縮閉線のグラフを図示すると、
上図において、縮閉線は、2個の特異点(尖点)を持つことが分かる。(t≒±0.5?)
S(H)さんのアドバイスと上図を鑑みると、上記で求めた縮閉線(曲率中心の軌跡)は、
与えられた曲線の法線群の包絡線と一致していることに気がつく。
実は、次の定理が成り立つ。
定 理 曲線 y=F(x) 上の点P( t , F(t) )における法線は、その曲率中心
Q(x , y )で、縮閉線に接する。さらに、縮閉線上の2点間の弧の長さは、
その両端における曲率半径の差に等しい。
(証明)
より、
よって、
このことから、曲率中心Qにおける接線の傾きと点Pにおける法線の傾きが一致する
ことが分かる。
また、縮閉線上の2点間の弧の長さを s とすると、
である。さらに、 (x−t)2+{ y−F(t)
}2=r2 ( ただし、r は曲率半径 ) の両辺を
t で微分して、
ここで、
を代入して、
また、法線の方程式から、
これらの2式より、
なので、結局
が成り立つ。そこで、曲率半径が増加するように弧の長さを測ることにすると、
すなわち、 s=r+(定数) である。
このとき、縮閉線上の2点間の弧の長さは、その両端における曲率半径の差に等しい
と言える。 (証終)
この定理から、縮閉線は、曲率中心の軌跡と言っても良いし、または、曲線の法線群の
包絡線と言っても良いことになる。さらに、もとの曲線が伸開線と呼ばれる理由も頷ける。
伸開線と縮閉線との関係を図示すると下図のようになるであろうか...。
左図において、縮閉線上の任意
の点で紐を画鋲で止め、縮閉線に
沿わせた紐の先端Oを鉛筆に固定
して紐をピンと張ったまま移動させ
る。このとき、鉛筆の先端Pは曲線
y=F(x) 上にあり、弧QQ’の長さ
は、ちょうど、 P’Q’−PQ に等し
い。
(まさしく、縮閉線に巻き付けた紐がちょ
うど解けていくように...。)
上記の定理で述べたことは、次のように直観的に考えた方が、す〜っと頭の中に刻まれ
るかもしれない。
曲線 y=F(x) 上の2点P、P’はごく近い場所にあるものとする。2本の法線PQ、P’Q’
の交点をRとおく。
このとき、△RPP’は、ほとんど2等辺三角形の形をしているので、 RP≒RP’ として
よい。
また、縮閉線上の2点Q、Q’を結ぶ弧の長さ s について、 s≒RQ+RQ’ としてよい。
よって、 s≒RQ+RQ’=(RP−PQ)+(P’Q’−RP’)≒P’Q’−PQ
(コメント) 上記の定理のスゴイところは、「2点P、P’はごく近い場所にある」と仮定しなく
ても、キッチリと「 s=P’Q’−PQ 」が成り立つことですね!
さらに、この定理から気がつくことは、曲線 y=F(x) に対して縮閉線はただ一つ存在す
るが、逆に、与えられた曲線を縮閉線として持つ伸開線は、無数に存在しうるということだ。
縮閉線の媒介変数表示を見れば、第1次導関数や第2次導関数が含まれるので、この
ことは自然に納得されるし、上図を見ても、紐の長さを変えれば、また違った曲線が得ら
れるだろうことは十分推測される。
種々の曲線に対して、興味深い縮閉線の幾つかを計算してみよう。
例 楕円
に対して、その縮閉線の方程式は、( α , β )を楕円上の点として、
α2+9β2=9 のとき、法線群 9βx−αy−8αβ=0 の包絡線の方程式
に一致する。
そこで、 F( x , y , α , β )=9βx−αy−8αβ 、
G(α,β)=α2+9β2−9
とおくと、
FαGβ−GαFβ=(−yー8β)(18β)−(2α)(9x−8α)=0 より、
9β(y+8β)+α(9x−8α)=0
以上から、
連立方程式 9αx+9βy=8α2−72β2 、 9βx−αy=8αβ
が得られる。 これを解いて、
81x=8α3 、y=−8β3 となる。
これらを、α2+9β2=9 に代入して、
このグラフを図示すると、下図のようになる。
(コメント) 何となく宇宙の神秘を眺めているような気分にさせてくれますね!
曲線 y=F(x) 上の点P( t , F(t) )における法線上の曲率中心Q(x
, y )の座標は、
により与えられた。さらに、このときの、曲率 κ は、
であった。
それでは、曲線が X=φ(t) 、 Y=ψ(t) と媒介変数表示されている場合の曲率中
心の座標と曲率は、どのように表されるだろうか?興味あることなので計算してみた。
結論から言えば、
が成り立つ。
これらは、
y’=ψ’/φ’ 、 y”=(φ’ψ”−φ”ψ’)/(φ’)3
などの媒介変数表示の場合の微分の公式から明らかだろう。
例 サイクロイド X=r( t−sin t ) 、 Y=r( 1−cos t ) に対して、
X’=r( 1−cos t
) 、 Y’=r・sin t 、 X”=r・sint 、 Y”=r・cos t
なので、 X’Y”−X”Y’=r( 1−cos t
)・r・cos t−r・sint・r・sin t
=r2( cos t −1
)
また、 X’2+Y’2=r2( 1−cos t
)2+( r・sin t )2=2r2( 1−cos t
)
よって、縮閉線上の点( x , y )は、
x=r( t+sin t ) 、 y=−r( 1−cos t )
ここで、 t=π+θ と置換すれば、
x=r( θ−sin θ )+πr 、 y=−r( 1+cos θ
)=r( 1−cos θ )−2r
となり、やはり、サイクロイドである。
r=2 として、そのグラフ達を眺めておこう。
(コメント) もとのサイクロイドのどの部分が新しいサイクロイドのどの部分に対応するのか
想像するのも楽しいですね!
曲線の方程式を扱う場合、媒介変数表示の場合を押さえておけば十分だろうが、極座標
による場合も時として必要になるので、公式として、まとめておこう。
極座標で曲線の方程式が、 R=F(θ) と表されるものとする。
このとき、直交座標系に変換して、曲線の方程式が、
X=F(θ)cosθ 、 Y=F(θ)sinθ
と媒介変数表示されることに気がつけば、後は容易だろう。
結論から言えば、曲線上の点P( F(θ)cosθ , F(θ)sinθ )における法線上の曲率
中心Q(x , y )の座標と曲率 κ について、
が成り立つ。
実際に、 X’=F’cosθ−Fsinθ 、 Y’=F’sinθ+Fcosθ で、さらに、
X”=F”cosθ−F’sinθ−F’sinθ−Fcosθ=F”cosθ−2F’sinθ−Fcosθ
Y”=F”sinθ+F’cosθ+F’cosθ−Fsinθ=F”sinθ+2F’cosθ−Fsinθ
このとき、
X’2+Y’2=(F’cosθ−Fsinθ)2+(F’sinθ+Fcosθ)2=F2+F’2
X’Y”−X”Y’=(F’cosθ−Fsinθ)(F”sinθ+2F’cosθ−Fsinθ)
−(F”cosθ−2F’sinθ−Fcosθ)(F’sinθ+Fcosθ)
=F2+2F’2−FF’
となり、これらの結果を媒介変数表示の場合の公式に代入して、上記の公式を得る。
例 心臓形 R=F(θ)=1+cosθ に対して、F’=−sinθ 、F”=−cosθ なので、
F2+F’2=(1+cosθ)2+sin2θ=2(1+cosθ)
F2+2F’2−FF”=(1+cosθ)2+2sin2θ−(1+cosθ)(−cosθ)
=3(1+cosθ)
よって、縮閉線上の点( x , y )は、
x=(1+cosθ)cosθ−(2/3){(1+cosθ)cosθ−sin2θ}
=(1/3){(1+cosθ)cosθ+2sin2θ}
=(1/3)(1+cosθ+sin2θ)
y=(1+cosθ)sinθ−(2/3){(1+cosθ)sinθ+sinθcosθ}
=(1/3){(1+cosθ)sinθ−2sinθcosθ}
=(1/3)(1−cosθ)sinθ
ここで、 (1−cosθ)cosθ=cosθ−cos2θ=−1+cosθ+sin2θ なので、
x−2/3=(1/3)(−1+cosθ+sin2θ)=(1/3)(1−cosθ)cosθ
このとき、 求める曲線は、やはり(?)心臓形 R=(1−cosθ)/3 を、
x 軸方向に
2/3 だけ平行移動したものになる。
因みに、心臓形 r=1+cosθ の代数曲線表示は、
x=r・cosθ 、 y=r・sinθ で、 r2=x2+y2 に注意して展開すると、
x4−2x3+2x2y2−2xy2+y4−y2=0
と与えられる。
この形で縮閉線を求めることは、あまり得策ではないように感じるが、S(H)さんは、この
代数曲線表示の形での縮閉線の計算に何か意味を見い出されているようだ。
(コメント) 上記の2例では、偶然にも、サイクロイドの縮閉線はサイクロイド、心臓形の縮
閉線も心臓形と同種の曲線となった。この話題について、当HPがいつもお世話
になっているS(H)さんも指摘されている。平成21年1月18日付けでS(H)さん
は、「縮閉線を反復して計算しても名前が不変な曲線の他の事例は?」と問いか
けておられる。
そのような例としては、アステロイドが一番面白いのではないだろうか?
もとのアステロイド曲線を、 X=cos3θ 、 Y=sin3θ とおくと、
縮閉線上の点( x , y )は、 x=cos3θ+3sin2θcosθ
y=sin3θ+3cos2θsinθ
この曲線は、もとのアステロイド曲線を、原点中心に45°回転し、2倍に相似
拡大した、やはり、アステロイド曲線である。
実際に、 x−y=cos3θ+3sin2θcosθ−sin3θ−3cos2θsinθ
=(cosθ−sinθ)3
x+y=cos3θ+3sin2θcosθ+sin3θ+3cos2θsinθ
=(cosθ+sinθ)3
よって、
から明らかだろう。
構成法から、縮閉線を求める操作を繰り返しても、常にアステロイドが得られるこ
とが見て取れる。(とても美しい結果ですね!)
上記では、陽関数表示、媒介変数表示、極座標表示の各場合について、
曲線上の点P( t , F(t) )における法線上の曲率中心Q(x
, y )の座標
を求めたが、それらの公式の構造を簡単化すれば、次のようになる。
曲線上の点Pにおける接線が x 軸の正の向きとなす角を α とし、曲率半径を
r とすると、
曲率中心Q(x , y )の座標は、
x = t−r・sinα 、 y = F(t)+r・cosα
と表される。
tanα=F’ なので、これより、 sinα 、cosα は容易に求められるだろう。
(コメント) 考えてみれば当たり前の公式だが、グチャグチャだった公式がスッキリして覚
えやすくなりましたね!
上記の公式の簡単化について、S(H)さんから、ご教示を頂いた。
(平成21年1月21日付け)
曲線 X=φ(t)、Y=ψ(t) 上の点P(φ(t),ψ(t))に対する曲率中心をQ(x , y )とす
ると、曲率中心は点Pにおける法線上にあり、しかも、その曲率中心の軌跡は、法線群の
包絡線として得られることを用いれば、
無理して公式を覚えなくても、次の連立方程式
ベクトルの内積 F(x,y,t)=(φ’,ψ’)・(x−φ,y−ψ)=0 および
t に関する偏導関数 Ft(x,y,t)=0
を解くという基本方針に従えば、曲率中心の媒介変数表示が「自然に」導かれる
という。
確かに、原理は上記の通りなのだが、曲線のいろいろな表示方法によって、その結果が
どのような形になるのかをまとめてみることも捨てがたいという気持ちもある。(ちょっと、受
験数学に毒されているかな?)
実際に問題に対処する場合は、S(H)さんの仰るような基本方針で臨んで解決に至ると
いう流れも数学的に美しいものがある。
例 放物線 y=x2 について、以前公式を用いて、その縮閉線を求めたが、上記の考えに
従えば、公式を特段意識することなく確かに自然に縮閉線の媒介変数表示が得られ、
感動的である。
y=x2 の媒介変数表示は、 x= t 、 y= t2 なので、
F( x , y , t )=( 1 , 2t )・( x−t , y−t2 )=x−t+2ty−2t3=0
このとき、 Ft( x , y , t )=−1+2y−6t2=0
これより、 y=3t2+1/2 で、 x=−4t3 となり、縮閉線の媒介変数表示が自然に
導かれた。
(追記) 平成21年1月25日付け
この包絡線の定義に関連して、平成21年1月24〜25日付けで、当HPがいつもお世話
になっているS(H)さんに、ご教示いただいた。S(H)さんに感謝します。
包絡線を定義する場合、当初当HPでは、
「媒介変数 m を含む曲線群 F={ F( x,y,m )=0 }のすべてに接して、
しかもその接点の軌跡となる曲線のことを、曲線群の包絡線」
としていたが、どうもこの定義は勇み足らしい。(冒頭の定義は修正済み)
このことを考える材料として、S(H)さんから、次の問題に取り組むことを勧められた。
問 題 曲線群 (x−t)2+y2=1−t2/4 の包絡線の方程式を求めよ。
(解) F( x , y , t )=(x−t)2+y2−1+t2/4=0 とおくと、
Ft( x , y , t )=−2(x−t)+t/2=0
よって、 x=5t/4 で、 y2=1−5t2/16=1−x2/5 より、
包絡線は、 楕円 x2/5+y2=1 となる。
(コメント) 上記の計算で、今まで気づかなかったことに気づかされた。
曲線群の中には、例えば、 t=−2 すなわち 点( −2 ,
0 )も含まれるので、
「曲線群のすべてに接する」という定義では、包絡線が存在しなくなってしまう。
この事情は、t=2 のときも同様である。
上図で、点( −2 , 0 )の周辺の様子をもう少し精査してみよう。
1−t2/4≧0 より、−2≦t≦2 なので、円の中心( t , 0 )は、点( −2 , 0 )から x 軸
に沿って、点( 2 , 0 )まで移動する。この際、円の半径は0より単調に増加し、t=0 のとき
最大値 1 のあと単調減少に転じ、0に至る。
このときの円と x 軸との交点のうち左側の交点の x 座標は、
であり、右側の交点の x 座標は、
である。この両者のグラフを描くと、下図のようになっている。
上図から、x の取り得る値の範囲は、 で、円達の大きくなっていく様子が
想像できる。
また、この範囲は、楕円 x2/5+y2=1 における x の変域にも一致する。(当然!)
このことから、曲線群 (x−t)2+y2=1−t2/4 のうち、
、
の範囲にある t の値で定まる円達は、包絡線とは全く接していない。
サッカーの試合でオフサイドの位置にあっても、プレーに全く関係していない場合はオフサイドになら
ないのと同じような感覚だ!
この件について、平成21年1月27日付けで、当HPがいつもお世話になっている、らすか
るさんは、「接しないような t の値を除外して、t の定義域を限定して『包絡線』と呼べば問題
ない」と述べておられる。
そうすると、曲線群を最初から制限された形で与えればいいわけだが、それは数学の問題
として美しさが損なわれることに等しい。
包絡線上の点は必ず、「ある」曲線上の点でもあり、かつ両者はその点で接するとすれば、
特別に曲線群を制限することもなくなると思うのだが...。
包絡線の「包」は「包む」、「絡」は「絡む」ということで、曲線群を包み絡める存在として「包
絡線」がイメージできる。このページの冒頭の例のように、包絡線は、曲線群の通過領域の
境界線という認識に立てば、包絡線に接しない曲線が存在してもいいような...そんな雰囲
気。こう考えた方が私としては、包絡線のイメージがスッキリする。
例 曲線群 x2+y2=1+t2 の包絡線を求めよ。
t を実数とすると、全ての曲線群に接するような曲線は存在しない。実際に、同心円の集
まりであることから明らかだろう。
しかし、包絡線は、現実のものとして確かに存在する。
曲線群 x2+y2=1+t2 の通過領域の境界線という立場に立てば、包絡線は、原点中心
で半径1の円 : x2+y2=1 である。
この曲線は、包絡線を求める公式からも得られる。実際に、
F( x , y , t )=x2+y2−1−t2=0 とおくと、Ft( x , y , t )=−2t=0 なので、
両者を連立すれば、包絡線の方程式 x2+y2=1 が得られる。
また、包絡線 x2+y2=1 上の任意の点( x , y )に対して、その点を共有する曲線群
の中の曲線(x2+y2=1)が存在し、しかもその点で両者が接することは明らかである。
したがって、曲線 x2+y2=1 は、曲線群 x2+y2=1+t2 の包絡線であると言える。
この話題に関連して、平成21年1月28日付けで、S(H)さんから次の問題を解いてみる
ことを勧められたので解いてみた。
問 題 曲線群 x2/t2+y2/(1−t)2=1 の包絡線の方程式を求めよ。
(解) F( x , y , t )=(1−t)2x2+t2y2−t2(1−t)2=0 とおくと、
Ft( x , y , t )=−2(1−t)x2+2ty2−2t+6t2−4t3=0
連立方程式を解くと、 x2=t3 、 y2=(1−t)3 となる。
この計算から、包絡線が存在するためには、必然的に、 0<t<1 と媒介変数
t の
値が制限される。
このとき、包絡線は、 アステロイド となる。
(コメント) t の制限がないと、曲線群はほとんど平面を埋め尽くしてしまって包絡線が見え
てこない。このように、t に制限をつける方が包絡線の定義としてはベストだろうと
思った。
(追記) 平成27年3月2日付け
平成27年度入試 東京大学 前期理系 で次のような問題が出題された。包絡線問題の
適当な練習となろう。
問題 正の実数aに対して、座標平面上で次の放物線を考える。
C : y=ax2+(1−4a2)/(4a)
aが正の実数全体を動くとき、Cの通過する領域を図示せよ。
(解) 両辺を4a倍して整理すると、 4(x2−1)a2−4ay+1=0
題意より、このaに関する2次方程式が少なくとも一つ正の解を持つ条件を求めればよい。
x2−1=0 すなわち、 x=±1 のとき、 4ay=1 が正の解を持つためには、y>0
x2−1≠0 のとき、 F(a)=4(x2−1)a2−4ay+1 とおくと、F(0)=1>0 なので、
題意を満たすためには、
x2−1>0 のとき、 放物線の軸について、 y/{2(x2−1)}>0 すなわち、 y>0
判別式について、 4y2−4(x2−1)≧0 すなわち、 x2−y2≦1
x2−1<0 のときは必ず正の解を持つ。
以上から、Cの通過する領域を図示すれば下図を得る。
ただし、境界線は実線部分は含み、点(±1,0)
および点線部分は含まない。
(コメント) 受験問題集に必ず載っているような典型問題で、合格のためには絶対落とせな
い問題ですね。それにしても河合塾や駿台予備校が発表している解答例は随分と
難しく解いているような雰囲気。その点が大いに気になりました。
以下、工事中