このページは[url=https://www.geogebra.org/m/ddm5798j]マス旅[/url]の一部です。[br][br]今回は「曲率の視覚化」がテーマ
[b][size=150]<角には平角もマイナスもある>[br][/size][/b][br]角を始めて学んだとき、[br]90度が角であることはよいとして、[br]「[b]平角、180度が尖がってないのに角[/b]」というんだ[br]というのに違和感を持った人もいるでしょう。[br][br]また、高校に入って、回転の向きで正の角だけでなく「[b]負の角[/b]」があることを知り、[br]角を「[b]弧の長さで表せる[/b]」ことも知ることで、「[b]角がただの1次元の数量[/b]」である[br]という落としどころに触れて、[br]「角の本質」に納得した人もいれば、[br]「角が抽象的」に感じられて数学嫌いになった人もいるかもしれません。[br][br][b][size=150]<多角形には凸も凹もある>[br][/size][/b]多角形を学んだときに、わざわざ「[b]凸多角形[/b]」というコトバに出会うと、[br]多角形にへこみがあってもいいんだとわかることで、[br]数学の「[b]ゆるい定義[/b]」に親しみを覚えた人もいるでしょう。[br][br][b][size=150]<曲面の曲がり方>[br][/size][/b][br]山があれば谷もある。[br][br]前回、水面の波の干渉をやりましたが、[br]水面という曲面がぶつかり合うと、[br][b]腹(山山、谷谷)と節(平地)[/b]ができました。[br]まさに1つの面なのに、いろんな見どころがあったね。[br][br]表面が一様な作りの立体もあれば、[br]ところ変われば、見どころも変わるという立体もあるでしょう。
[size=150][b]<2つの主曲率の積がガウス曲率>[/b][br][/size][br]曲面上の、ある1点(点P)の曲がり具合を調べるために、[br]その点で曲面をスパッと切った「切り口の曲線」を考えます。[br]切る方向によって曲がり方は変わりますが、[br]その中で「[b]最もきつく曲がっている方向(最大)[/b]」と[br]「[b]最も緩やかに曲がっている方向(最小)[/b]」の2つを取り出します。[br]この2つの曲がり具合の数値を[color=#0000ff][b]主曲率(κ_1,κ_2)[/b][/color]と呼びます。[br]この2つの掛け算[b]κ_1κ_2[/b]こそが、曲面のキャラクターを決定づける[b]ガウス曲率K[/b]です。[br][color=#0000ff][b][size=200]K=κ_1κ_2[/size][/b][/color][br][br]少し具体的な例を上げましょう。[br]みかん(球面)はどちらの方向にも同じ側に曲がる (κ_1>0,κ_2>0)だから[b]K > 0[/b][br][br]馬の鞍(サドル)は前後は凸で左右は凹のように逆に曲がる (κ_1>0,κ_2<0)だから[b]K < 0[/b][br][br]円柱側面は片方は曲がっているが、もう片方は真っ直ぐ (κ_1>0,κ_2=0)だから[b]K = 0[/b]。[br][br][b][size=150]<ガウス曲率の計算式>[br][/size][/b][br]関数fのガウス曲率は次の計算で求められることがわかっています。[br]f,fx,fy,fxx,fyy,fxyから求められます。[br][br][color=#0000ff][b][size=200]K = {fxx fyy - fxy^2}/{(1 + fx^2 + fy^2)^2}[br][/size][/b][/color][br]分子:fxx fyy - fxy^2は、[br][b]Kの符号決定[/b]。[br]同じ向きの曲がりは+、逆向きに曲がりは-です。[br][br]分母:(1 + fx^2 + fy^2)^2は、[br]曲面の傾きによる空間的な引き伸ばしを補正しています。[br][br][b][size=150]<6種類の2次曲面>[br][/size][/b][br]次は、毎度おなじみの6種の2次曲面で[br]ガウス曲率Kチェックをしよう。[br][br]楕円面はどこも凸でK>正です。[br]2葉双曲面、どこでもK>正ですね。[br]楕円放物面、どこも凸でK>0。[br]2次の錐面、ふくらむ楕円と直線が1点で直交しK=0。[br]1葉双曲面、ふくらむ楕円に凹双曲面の直交してK<0。[br]双曲放物面は鞍に乗るあたりは凸と凹が共存してK<0。[br][br][b][size=150]<ガウスの定理>[br][/size][/b][br]「曲面をを引きちぎったり、伸ばしたり(伸縮)しない限り[br]ガウス曲率を変えことはできない」[br]というのがガウスの定理です。[br]たとえば、みかんの皮(ゴムボール型)はK>0ですが、[br]これを切らずに伸ばさずに机の上においてK=0にはできません。[br][br]つまり、ガウス曲率とは、曲面を外から眺めたときの形(3次元的な見た目)ではなく、[br][color=#0000ff][b]「その曲面自身が持っている、伸縮させない限り絶対に変わらない固有のDNA(内在的幾何)」[/b][br][/color]を表していると言えますね。
3.コード化[br][br][b][size=150]<geogebra>[br][/size][/b][br]snum=Slider(1,3,1)[br]f(x, y) = 0.1 x^2 - 0.1 y^2 #サドル[br]g(x, y) = 0.1 x^2 + 0.1 y^2 #放物楕円面[br]h(x, y) = 0.1 x^2 #円柱[br]Fs={f,g,h}[br]#マーカーの色を際立たせるために、Fの属性である色は標準75を50%にさげましょう。[br]F= Element(Fs,snum)[br]P = Point(F)[br]Px=x(P)[br]Py=y(P)[br]fx(x,y) = Derivative(F,x)[br]fy(x,y) = Derivative(F,y)[br]fxx(x,y) = Derivative(F,x,2)[br]fyy(x,y) = Derivative(F,y,2)[br]fxy(x,y) = Derivative(Derivative(F,x),y)[br]K = (fxx(Px,Py) * fyy(Px,Py) - fxy(Px,Py)^2) / (1 + (fx(Px,Py))^2 + (fy(Px,Py))^2)^2[br][br]Ps=Sphere(P, 0.3) [br]#ガウス曲率マーカーです。[br]#設定>動的な色[br]#シグモイド関数的に0付近を漸近するカーブにしています。[br]#RGB 動的カラー計算用の式(0〜1 にスケール)[br]#赤 (R): K > 0 (正の曲率)放物楕円面[br]If(K > 0, K / (K + 0.05), 0)[br][br]#緑 (G): K < 0 (負の曲率)サドル[br]If(K < 0, -K / (-K + 0.05), 0)[br][br]#青 (B): K == 0 (平坦・円柱)[br]If(K == 0, 1, 0)
[br][b][size=150]<振り返り>[br][/size][/b][br]. 曲面切り替えスライダー [code]snum[/code] を動かすと[br]載せているマーカー [code]Ps[/code] の色[br]「 緑(サドル) ➔ 赤(楕円放物面) ➔ 青(円柱) 」[br]と鮮やかに変身しますね。[br]まさに、曲率は曲面のDNAです。[br][br]円柱では違和感がないですか。[br]曲がっているのに曲率は青(ゼロ)のままです。[br]見た目はまがってていても、いつもまっすぐで平らな道があるからです。[br]だから、いつも青(ゼロ)です。[br][br]楕円放物面はくぼんでいるから負になるという勘違いをしやすいです。[br]楕円放物面は外からみると凸に見えて、内からみると凹んで見えます。[br]でもガウス曲率は[b]主曲率の積だから、正×正でも、負×負でも正[/b]になりますね。[br]だから、いつも赤(正)です。[br]