Weekend Mathematicsコロキウム室/NO.156

コロキウム室



NO.1317 2002.11.1.物理君誤差関数

はじめまして。僕は東京都内の大学で物理学を学んでいる者です。
誤差関数(エラーファンクション)とは何か、教えていただけませんか?



NO.1318 2002.11.5.DDTガウスの発散定理の応用

最近、落石防護擁壁の転倒計算に関わりました。落石防護擁壁とは、図-1に示すようなコンクリートの壁のことで、斜面を転げ落ちて来た落石に衝突されても、倒れないだけの自重を持たなければなりません。

どれほどの自重(重さ)があれば良いのでしょうか?。斜面を転げ落ちて来る落石は、運動エネルギーを持ちます。もし落石の運動エネルギーがそのまま全部、擁壁に伝わるなら、面倒はありません。擁壁重心の位置エネルギーを落石の運動エネルギーに等しいとおいて、擁壁重心の変位を計算し、それが回転中心の位置を越えなければOKです。しかし現実のエネルギー伝達には、常にエネルギー損失が憑き物(?)で、全てが擁壁に伝わるわけではありません。エネルギー損失の定量理論は現在に到っても難しい問題で、一般的な傾向を述べる概論めいたものはあっても、個々の特殊事情までも考慮して使用できるような、完璧な理論化はなされていません。いきおい実験ということになります。ところがエネルギーは直接測定できる量ではありません。

ではいったい何が測定できるのかというと、例えば、コンクリート擁壁に鉄球をぶつけるような実験では、鉄球に加速度計を仕込んでおきます。加速度計は、鉄球が擁壁に衝突する際に生じる、鉄球の加速度を記録します。その衝突加速度波形に鉄球の質量をかければ、衝突の際に擁壁から鉄球に作用した衝撃力波形が得られます。これは鉄球から擁壁に作用した衝撃力でもあるはずです。要するに直接測れるのは、力の時間変動の様子だけです。衝撃力波形の記録を数値的に時間積分すれば、衝突時に鉄球から擁壁に伝達された運動量(力積)が出ます。この運動量によって擁壁は動き出すのですから、伝達運動量に見合う擁壁の運動速度を定めれば、運動エネルギーの形で鉄球から擁壁へ伝達されるエネルギーを計算できます。この伝達エネルギーは、衝突の際に生じる力の実測から得られたものなので、衝突時に起きたエネルギー損失を全て含む、正しい伝達エネルギーであるはずです。

すいません。やっと本題です。愚痴っぽく言うと、伝達運動量から擁壁の運動速度を算出するために、剛体回転の運動方程式を解く破目に陥りました。剛体回転の運動方程式には、次の3つの定数が必要です。

  質量M,重心位置(x0,y0,z0),固有慣性モーメントI0(重心回りの慣性モーメント)

で、密度が一様な2次元剛体の場合は、それぞれ、

  面積A,重心位置(x0,y0),固有極断面2次モーメントIp

の計算に帰着されます。計算式は以下となります。


対象とする積分領域Rは、右図にように比較的簡単なものです。まず(1)の面積は台形公式で簡単に計算できますが、(2)になると、もう単純な公式はすぐには出てきません。それでも構造力学公式集やデザインデータブックを検索すると、何とか使う気になる程度の公式は出てきました。でも(3)については公式すらありません。
しかたないので、右図のような座標系を入れ、x,yに関する2重積分を使って計算しましたが、かなり面倒な微積分計算となり(高校レベルです)、久しぶりにやったので3時間もかかってしまい、その上結果に自信が持てないという有様でした。
じつは最初に述べた実験では、せっかくだからということで、加速度計のほかに歪み計(擁壁の変形を測定する)やレーザー変位計(擁壁の変位を測定する)を沢山盛り込み、実験は相当に高額なものでした。もちろん実験資金は、発注者の金です。「こりゃ、ちゃんとした結果を出せないと首がとぶぞ」と皆で薄笑いを浮かべる始末。(3)の結果の検算に迫られました。

(1)〜(3)は、ガウスの発散定理を利用して計算できます。これが思いのほか便利だったので、報告します。しかも、どんなにいびつな領域に対しても、それを詳しく折れ線近似しておけば、十分な精度で答えを出せます。任意多角形に対してはもちろん正解を与え、Cadの利用を前提に考えれば、折れ線近似をCad上でスケッチしただけで、折れ点の座標データをファイル出力できるので、プログラム化(自動化)にはもってこいです。


Rは積分領域,Sはその境界,dsは境界の接線線素ベクトルdnと右回りに直行する外法線線素ベクトル, 一周積分の方向は、もちろん左回り.

式(4)において、
f=(x/2,y/2)TでAが,f=(x2/4,xy/2)Tでx0が,f=(xy/2,y2/4)Tでy0が,f=((x−x0)3/3,(y−y0)3/3)TでIpが出ます。



NO.1319 2002.11.10.水の流れ最大・最小となる点

第108回数学的な応募問題

太郎さんは、数学Aの平面幾何を眺めていたところ、次のような問題をみつけました。皆さん、チャレンジください





NO.1320 2002.11.11.DDTガウスの発散定理の応用(2)

ガウスの発散定理について説明します。

1.2重積分について

2重積分とは、図-1のように平面領域Rがあり、Rに含まれる点(x,y)において2変数関数f(x,y)が定義されているとき、

  

のように表わされるもので、関数f(x,y)をxで積分した後、yで積分せよ、という意味です。 最も簡単な例として、f(x,y)が変数(x,y)について定数1である場合を考えます。

  

なので、2重積分(1)は、

  

となります。式(3)は、領域Rの面積を与えます。どうしてかというと、まずyを一つ固定して、定数1をxで積分するので、

  

となりますが、g1(y)とg2(y)とは図-1に示したように、領域Rの境界Sをx方向の高さに関して上下に分けて、縦線@〜Aまでの境界Sの部分をx=g2(y),縦線A〜@の部分をx=g1(y)と表わしています。縦線@とAに直交する点線は、x軸に平行な境界Sの接線を表わします。g2(y)−g1(y)は 明らかに、固定したyに応じた領域Rの横幅です。 次にg1(y)とg2(y)の不定積分をG1(y),G2(y)とした場合、 (4)の最右辺では最初に固定したyを、今度は動かして積分するので、


  

となり、確かに領域Rの面積です。何故なら式(5)は、各yに対して変化する領域Rの横幅を、yの下限から上限まで順番に積算していったものだからです(図-2)。同じ結果は、yで積分した後、xで積分しても得られるはずです。図-3は、各xに対して変化する領域Rの縦幅を、xの下限から上限まで順番に積算していったことを意味します。 よって、



  

です。うるさい事さえ言わなければ、もっと一般に任意の関数f(x,y)に関して、

  

も成り立ちます。横に積算してから縦に積算するか、それともその逆かの違いしかないからです。式(7)の左辺を、図-2方式で2重積分すれば、

  

となります。

2.ガウスの曲面論[1]

多様体の概念の基礎を与えたのは、恐らくリーマンです。こんにち多様体の一分野となった、微分幾何学や微分多様体を出発させたのは、しかしガウスです。ガウスの曲面論は疑いもなく、リーマンの多様体の概念の出発点にもなりました。幾重にも広がった量として、その就任公演で多様体のイメージを初めて述べたリーマンですが、就任公演は当時、大学への就職試験のようなものでした。その面接官であったガウスは、リーマンに完全な満点評価を付けました。多様体の概念は、ガウスの曲面論をより一般的に具体化するものだったからです。
多様体の重大な理論的成果の一つに、図形はその内部性質のみによって、完全に自己を記述できる、というのがあります。これはある図形を、その図形を含むより高次元の空間から、外から見なくても、図形に貼りついた視点で図形を詳細に調べれば、外から見た図形の姿を再現できるという意味です。考えようによっては、ものすごい結果だと思います。例をあげれば一般相対性理論です。
我々は3次元空間に住んでいますが、その3次元空間を詳細に調べれば、一般相対性理論が言うように宇宙は曲がった4次元時空だとわかる、という意味です。実際、4次元の中の3次元宇宙の曲率を測る実験の数々が、現在進行形で考案されている現実があります。これを3次元の中の2次元曲面論の形で最初に提出したのは、ガウスでした。



曲面を外から見た姿を曲面内部から知る方法は、私の狭い理解範囲内で言えば、曲面に沿った座標系を曲面に入れることによって実現されます。曲面に沿った座標系とは、図-4の右図のようなものです。対比として、平面に対する平面に沿った座標系を書いてあります。図-4においては、平面の格子点が、そのまま曲面の格子点に対応していると想像して下さい。余りにも面倒なので図-4の座標格子では省略してしまいましたが、曲面に沿った座標格子の格子点間隔は一定でなくても良いのです。曲面の格子点間隔を平面格子に対応付けて、どのように自然に与えれば良いかは、ガウスの理論に述べられています。そして、曲面の格子点間隔の変化を丹念に追跡することで、曲面の曲がり方がわかり、ひいては曲面を外から見た姿も再現でき、かつ格子点間隔の測定は曲面内部の性質だけから決定できるはずだ、というのがガウスの曲面論です。だって格子点間隔は、曲面に沿って測るだけですから。逆に言えば、曲面と仮定して測った格子点間隔がどこでも一定なら、図-5に示したように、どんなに傾いていても、どんなに遠い場所にあろうと、どんなに大きさが違っても、平面は平面で平らだということが内的性質のみを使って、もってまわったやり方ではありますが、ともかく一目瞭然となります。ここで、図形の傾きや位置やスケールは、外的性質であって、図形の形とは無関係です。



現在ではガウスの曲面論のアイデアは、3次元CadのRuled Surface(線織面)や、Meshed Surfaceの計算法として、コンピュータ内で実現されています。もしガウスがいなかったら、CGも造れなかったかもしれません。

3.内部座標系の威力

もう一度、式(8)と図-1を、式(9)と図-5として並べて書きます。

  



式(9)の右辺において、F(x,y)のxは、g1(y)とg2(y)におきかわっています。g1(y),g2(y)は、二つ合わせて式(9)左辺の積分領域Rの境界Sを構成します。ということは、式(9)の右辺のFは、一見2変数関数の形をしていますが、実質は境界Sに束縛された、1次元の(1変数の)関数です。それを、FS(y)と書きます。
FS(y)のyに関する積分値は、関数x=g1(y)とx=g2(y)の形で決まるので、じつは境界Sの形で決まります。
よって式(9)は本来、

  

のように書かれるべきです。ここで右辺の積分範囲pr2(S)は、境界曲線Sのy軸への正射影で、区間[y1,y2]を表わします。ところで曲線は曲面の一種と考えることができます。境界曲線Sの形とは、まさに曲線Sに固有な内部的性質です。式(9)の値が、ということは式(10)の値が境界Sの形で決まるということは、言い直せば、積分値が境界Sの内部的性質から定まることを意味します。領域Rは境界Sの形をその定義とするので、領域Rに固有な、形という内部性質だけから、積分値が定まるといっても同じです。これが式(10)の意味です。この見方は事の本質を捉えていると思います。それは面積の例を考えれば明らかです。
事の本質が上記の通りなら、その通りに式を変更すべきです。2.で述べたように、曲面に沿った座標系とは、座標系自身が載る曲面の形すら再現するほど強力なものでした。そして式(10)の積分値は、積分領域の形から定まるものです。積分を実行する座標系として、境界に沿った座標系sを採用します。以後内部座標と呼びます。



次式は正しい式ではありません。 話を簡潔にするための模式図と考えて下さい。式(9)より、

  

と大変簡潔な形になります。式(11)と(9)を比べてみると、何故あんなに式(9)が不細工になったのかが理解できます。境界Sを外部座標x-yで表わしていたからです。まさに境界Sの形を外から観察し、x=g1(y),x=g2(y)と境界の形を表わす情報を外から与える必要に迫られました。式(10)においても、積分区間pr2(S)=[y1,y2]は、境界Sに固有な表現ではなく、やはり不満足なものです。[y1,y2]をy軸への正射影とする曲線は無数にあるからです。その責任は、外部座標yを積分パラメータとし、それを偏重したからです。式(11)では、それらの情報の全てが、内部座標系sの中に含まれています。よって式を不細工にした責任は外部座標にある!。お前が悪い!!。
とはいえ計算を実行する人間は、境界内部には住んでいません。問題をスケッチするために、最初は何らかの外部座標に頼らざるえません。理屈の上では、領域Rも境界Sも最初からその内部座標で表わしておけば、外部座標不要で済ませられますが、上の理由からそれは実用的ではありません。だから式(11)のように、外部座標を介在させる書き方が、最もバランスの良いやり方と思えますが、そうすると外部座標と内部座標とを結ぶ変換則が必要になります。Frenet-Serretの公式がそれですが、それを用いてx(s)やy(s)を求めようとすると、初等関数で積分できなくなったりしてなかなか辛い目に会います。だから式(9)のタイプも多用されるわけです。

4.ガウスの発散定理

長い前ふりでした。すいません。理屈の上では、2重積分の実行は、境界の内部座標で実行すべきだということになりましたが、その方法にはほとんど実用性がないのも事実です。しかし次の形の2重積分、

  

に関しては、内部座標での積分の実行を行うと、結果をきれいなベクトル形式にまとめられます。外部座標と内部座標の変換もFrenet-Serretの公式を使わずに、その関係をスケッチするだけで済みます。式(12)は、ガウスの発散公式(発散定理)の一部で、その名の示すとおり、彼はこのタイプの積分を考える機会が多かったはずです。そしてガウスは、内部座標の威力に、恐らく世界で初めて気づいた人です。

ガウスの発散定理

  

Rは積分領域,Sはその境界,dsはSの外法線線素ベクトル,・は内積を表わし、線積分の方向は左回り.
[説明]
まず問題を定式化するために、表式がいかに不細工になろうとも、外部座標を利用します。問題の定式化に関しては、それが唯一の現実的手段だからです。そのためにRの境界Sを、4つの部分に分けます(図-7)。
図-7で、境界Sの@→Aの部分をy=h1(x),A→@の部分をy=h2(x)で、B→Cの部分をx=g1(y),C→Bの部分をx=g2(y)で表わします。



次に式(13)の、

  

の部分ですが、ガウスでなくてもというか、ガウスでさえも、とりあえずは第1項はxで、第2項はyで積分したくなったに違いない思われます。1.で述べたように、x,yの順序で積分しても、y,xの順序で積分しても、2重積分の積分結果に違いはありません。f1=f1(x,y),f2=f2(x,y)と堂々と書くことにして、

 

です。ただし後で便利に書き換えられるように、積分区間の上下限の順序を適宜変えて、積分記号の前の符号を揃えています。内部座標sでの積分を実行したいので、式(14)をsで書き換えることを考えます。
上式右辺上段のf1(g2(y),y)は、f1(x,y)の(x,y)が曲線Sに束縛されたために、(x,y)=(g2(y),y)とパラメータ表示された状態だと読めますが、Sは、(x,y)=(g1(y),y)または(x,y)=(g2(y),y)と表わされる点全体のことなので、右辺上段の積分を実行するために、yをy1からy2まで動かせば、曲線Sに束縛された点(x,y)は、自動的に曲線S上の全ての点を通ります。しかも同じ点は2度と通りません。よって曲線Sの内部座標sと点(x,y)との間に関数関係が成立し、(x,y)=(x(s),y(s))という記述が可能です。下段のf2の(x,y)=(x,h1(x)),(x,h2(x))についても同様で、(x,y)=(x(s),y(s))の記述が可能です。そこで内部座標sによるパラメータ表現を用いれば、式(14)は、

  

と書けてしまいます。このときsの原点は、S上であればどこでもOKです。どこに原点をおいても結果は変わらないのは、積分の一般的性質です。この段階で、外部座標(x,y)と内部座標sとの関係をスケッチします。



式(15)において、sへの書き換え作業として残っているのは、dx,dyをsで表わすことだけです。 図-8に示すように、S上の一点(x(s),y(s))での単位接線ベクトルを(cosθ,sinθ)とすれば、接線ベクトルは線素をdsとして、

   (ds cosθ,ds sinθ)      (16)

となります。これのx,y軸への正射影が−dxとdyなので、

   dx=−ds cosθ
   dy= ds sinθ        (17)

です。dxに負符号がつくのは、図-8のdsがsに載って動いたときの正射影の方向と、式(14)の積分方向を比較するとわかります。
式(17)を(15)右辺に代入すれば、

  

を得ます。図-8に示したように、単位外法線ベクトルは(sinθ,−cosθ)なので外法線線素ベクトルを、

  

で導入し、f1とf2をまとめて、

  

と書くことにすれば、式(18)の右辺の積分記号内は、f・dsとなり、

  

が得られます。

5.発散公式の実用性

式(13)または(21)の計算を内部座標で行う必要は全くありません。まず領域Rが任意多角形の場合は、隣り合った角点の(x,y)座標から、辺上の座標値をたやすく補間できるので、f1(x,y)やf2(x,y)の関数値は、その座標値を代入することで即座に出ます。線素ベクトルdsの方向についても、隣り合った角点を結ぶ矢線ベクトルの方向から簡単に計算できます。そして線素の大きさdsは、多くの場合、矢線ベクトルの長さで代用すれば十分です。Sが任意の曲線の場合は、解析的計算は不可能なケースがほとんどなので、Sを細かく折れ線近似し、任意多角形の場合と同様に扱います。

参考文献
[1]多様体,萩上紘一,共立講座 21世紀の数学6,共立出版1997



NO.1321 2002.11.16.中川 幸一素数定理(2)

素数定理がすでに証明されているということをことを聞いたのですが、いつ頃誰によって証明されたのか分かりますか? また、証明するのに用いた数学の背景はどういうものですか? (証明方法の大まかな手順でも構いません。)



NO.1322 2002.11.16.yokodon最大・最小となる点(2)

やってみました。手っ取り早さのため、反則?(高校数学の範囲外の知識)を使っ ていますが、範囲内知識での答案が出てくることを期待して、露払いということにさ せて下さい。(^^;

(1)BC/PD + CA/PE + AB/PF を最小にする点P

辺BC、CA、ABの長さを各々 a, b, c 、線分PD、PE、PFの長さを各々 x, y, z とし、三角形ABCの内接円の半径を r とおきます。題意の点Pの関数を F(P) とすると、当然

   F(P) = a/x + b/y + c/z

です。ところで、三角形PAB、三角形PBC、三角形PCAの面積の和が三角形A BCの面積に等しいことから、

   ax + by + cz = (a + b + c)・r  …[1]

という条件式が成り立ちますね。そこで、多変数関数の条件付き極値問題に持ち込み ます(#嗚呼、反則哉)。
以下の関数を考えます。K は未定乗数です。

  G(P) = F(P) + K・{(ax + by + cz) - (a + b +c)・r}

偏導関数を計算すると、各々以下の通りです。

  ∂G/∂x = -a/x2 + Ka
  ∂G/∂y = -b/y2 + Kb
  ∂G/∂z = -c/z2 + Kc

従って、∂G/∂x = 0 、∂G/∂y = 0 、∂G/∂z = 0 から、

  K > 0 、x = y = z = 1/√K   …[2]

を得ます。これを [1] 式に代入して、K = 1/r2 を得ます。
念のため、ヘシアンを計算すると(一変数関数の微分による増減調べのときの、2 次導関数に相当する判定基準の一つ;定義から、x, y, z の変域が“> 0”であるこ とに注意して)、

  Hes[G(P)] = 8・abc/(x3・y3・z3) > 0

なので(ヘシアンが正である以外の条件を満たすことも議論には必要ですが、実は満 たしています)、確かにこのとき G(P) は最小で、従って F(P) も最小です。
さて、K = 1/r2 から x = y = z = r ですから、結論として、求める点Pの 位置は三角形ABCの内接円の中心、即ち内心です。…(答)

(2)PD・PE・PF を最大にする点P
同様に、多変数関数の微分で考えます。でも、こちらの方だと一変数消去、一変数 固定という“高校数学での多変数関数問題の定法(?)”にそってやる分にはやりや すいでしょうね(後で説明します)。
H(P) = PD・PE・PF とおき、以下の関数 J(P) を考えます。K は、やはりラグラ ンジュの未定乗数です。

J(P) = H(P) - K・{(ax + by + cz) - (a + b + c)r}

同様に偏導関数を計算すると、H(P) = xyz ですから、

  ∂J/∂x = yz - Ka
  ∂J/∂y = zx - Kb
  ∂J/∂z = xy - Kc

各偏導関数の値が 0 となる条件から、以下が得られます。

  K > 0 、(x , y, z) = √K(√(bc/a), √(ca/b), √(ab/c))  …[3]

この [3] 式を [1] 式に代入して、以下のような“厳つい”式を得ます。

  √K = (a + b + c)・r/{3・√(abc)}  …[4]

これを [3] 式に戻して、以下のように x, y, z が定まります。

  (x, y, z) =(a + b + c)・r/3 ×(1/a, 1/b, 1/c)  …[5]

長くなるので、ヘシアンとかの説明は省略しますね(笑)。
このとき、確かに J(P) 、従って H(P) は最大(のはず)です。
このときの、[5] 式の意味するところを考えてみましょう。ここで、三角形ABC の面積を S0 とします。 内接円の半径を用いて、S0 = (a + b + c)・r/2 である ことはすぐに分かると思います。これを、[4] に適用してみると、

  (x, y, z) =2/3・S0 ×(1/a, 1/b, 1/c)

例えば、これの x 成分を考えてみましょう。

  x = 2/3・S0 × 1/a

ですから、単純な変形で、

  S0 = 1/2×(3x)・a

となります。これは、底辺の長さが a 、高さが 3x の三角形の面積を表しますね。 つまり、点Aから辺BCに下ろした垂線の長さが 3x というわけです(#我ながら、 ちょっと無理矢理だな)。同様のことが、y, z でも言えます。ところで、点D、E 、Fは各々、点Pから辺BC、CA、ABに下ろした垂線の足でした。従って、題意 で求める点の位置は、三角形ABCの垂心です。…(答)
#本当は、点Pが線分AD上にあり、且つAD⊥BCであることなどを確かめなけれ ばならないのですが、それは宿題にさせて下さい。

=おまけ=

以上の議論では、高校数学の範囲を超えた議論になっちゃったので、せめて高校の 範囲の“一変数の微積分”で解決できる方針を示しておこうと思います。基本的な方 針は、3変数の間の関係式 [1] から変数を一つ消去し、2変数関数の問題にすると いうものです。そして、一変数を固定して、そのもとでの最大・最小値を固定した変 数の関数として求め、最後にその変数を動かすというものです。ただ、そのままでは やりにくいので、s = ax 、t = by 、u = cz などと変数変換すると、式の形が綺 麗になって見通しがよいでしょう。なお、3m = (a + b + c)・r とおきます(実は 、この値は 2・S0 に等しい)。

(1)では、s + t + u = 3m のもとで、

  F(P) = a2/s + b2/t + c2/u

から例えば u を消去し、s と t の2変数関数にします。あとは、例えば s を固定 して t を動かすことで、“暫定の最小値”F1(s) が得られます。あとは、これを s の関数として増減を調べれば結論にたどり着けますね。

(2)では、もっと簡単に

  H1(P) = stu (= abc×H(P))

を考えます。これは、u を消去して s を固定すると t の2次関数であり、その“暫 定の最大値”H2(s) は s の3次関数なので、(1)の場合よりも更に議論しやすい ですね。蛇足ながら、この場合は、形式的には1点に集う3辺の長さの和が一定の直 方体の体積の最大値を求める問題と同じ形です。

尤も、これでも数学Aの範囲を超えて(かな?;最近の科目区分には疎いので、正 確に把握していないのですが)、数学3(一変数関数の微積分)の範囲になっちゃう のですが(笑)、問(2)に限れば数学2(かな?;僕の時代の科目区分では、「基 礎解析」の範囲)まで落とせますね。

#初等幾何的な洗練された解法は、僕には思いつかないです(笑)。



NO.1323 2002.11.16.yokodon誤差関数(2)

ガウスの誤差関数のことです。正規分布の累積分布関数ですね。

  

詳しくは適当な統計学の教科書、または  こちら をご覧下さい。

#データソースによって、微妙に定義式が違いますが、変数変換すれば等価です。







E-mail 戻る