《積分幾何学》のソースを表示

'''【せきぶんきかがく (integral geometry)】'''

　“ビュッフォンの針”のように図形に関係した確率を幾何確率といい, これらの理論的な部分は, 積分幾何学(integral geometry)を基礎にしている. 積分幾何学は幾何確率に応用可能であるばかりでなく, もっと広く図形に関する様々な局面で基礎となり得ると考えられる. そこで積分幾何学の基礎概念と直線に関する主公式にかぎって述べることにする. 詳しくは[1], [2]を参照されたい. 積分幾何学の基礎概念とは合同変換によって不変な測度を求めることにあるが, 例を直観では分かりづらい直線の集合に採り, 以下で議論を展開する. 

        
　図1のように, 直線 <math>g\, </math> があって, 原点からこれに下ろした垂線の長さを <math>p\, </math> , 垂線と <math> x\, </math> 軸との角度を <math>\theta\, </math> とすると, この直線は


<center><math>
x\cos \theta+y\sin \theta=p
\, </math>　　　　　<math>(1)\, </math></center>


と表される. そこで, この<math>p\,</math>と<math>\theta \, </math>で, 直線の集合(たとえばある領域と交わる直線の集合など)の測度を


<center>
<math>
m(X)=\int_{X}f(p,\theta)\; \mbox{d}p \; \mbox{d}\theta 
\, </math>　　　　　<math>(2)\, </math>
</center>


とおき, 以下, 合同変換で不変の条件を考えよう. 直線 (1)で表されるような, 直線の集合 <math>X\, </math> が, 


<center><table><tr><td align=center>[[画像:sk-0207-c-e-01-1.png]]</td></tr>
<td align=center>図1：直線に下ろした垂線</td></table></center>

<!--
%\begin{figure}[h] %\begin{center} %%\vspace*{40mm} %%%%%%%%%%%%%%%%%%%%%%%%%% %\includegraphics{C-E-01f+2.eps} %\caption{直線に下ろした垂線} \label{fig:C-E-01+2} %\end{center} %\end{figure} 

\begin{figure}[h] \begin{center} \begin{picture}(300,170)(-50,-20) %\put(-50,-20){\framebox(300,170)} \put(0,0){\vector(1,0){220}} \put(220,0){\makebox(0,0)[l]{\large $ x $ }} \put(0,0){\vector(0,1){130}} \put(0,130){\makebox(0,0)[b]{\large $ y $ }} \put(0,0){\line(1,2){50}} \put(-10,100){\line(2,-1){220}} \put(140,30){\makebox(0,0)[bl]{\large $ g $ }} \put(22,71){\line(2,-1){11}} \put(22,71){\line(1,2){5}} \put(20,50){\makebox(0,0)[tr]{\large $ p $ }} \put(0,0){\makebox(0,0)[tr]{O}} \put(5,0){\oval(15,20)[tr]} \put(15,15){\makebox(0,0)[l]{\large $ \theta $ }} \end{picture} \vspace{-10mm} \caption{直線に下ろした垂線} \label{fig:C-E-01+2} \end{center} \end{figure} 
-->


合同変換


<center>
<math>\begin{array}{lll}
x' &=& x\cos \alpha-\sin \alpha+a \\
y' &=& x\sin \alpha+y \cos \alpha+b
\end{array}\, </math>　　　　　<math>(3)\, </math>
</center> 


によって, 


<center>
<math>
x'\cos\theta'+y'\sin\theta'=p' 
\, </math>
</center>


なる直線の集合 <math>X' \, </math> に変換されたものとする. すると, 上式より


<center>
<math>\begin{array}{lll}
\theta &=& \theta'-\alpha, \\
p &=& p'-a\cos \theta'-b\sin \theta' 
\end{array}\, </math>　　　　　<math>(4)\, </math>
</center>


が得られる. ここで, 一様な直線の分布を考えているので,  <math>X\, </math> と <math> X' \, </math> の測度は等しく, すなわ
ち


<center>
<math>\int_{X}f(p,\theta)\; \mbox{d}p \; 
\mbox{d}\theta=\int_{X}f(p',\theta') \; \mbox{d}p' \; 
\mbox{d}\theta'\, </math>
</center>


でなければならない. ところが, 式 (4)から <math>(p',\theta')\, </math>  より <math>(p,\theta)\, </math> への変数変換のヤコビアンは1なので, 上式の右辺を変形して


<center>
<math>
\int_{X}f(p,\theta)\; \mbox{d}p \;
\mbox{d}\theta=\int_{X}f(p',\theta')\; \mbox{d}p \; 
\mbox{d}\theta 
</math>
</center>


と表すことができ, これは <math> X\, </math> をどのようにとっても成立しなければならないので, 


<center><math>f(p,\theta)=c\, </math> <math>( \, </math>定数<math>) \, </math></center>


が得られる. そこで <math> c=1 \, </math> とし, <math> \mbox{d}p \; \mbox{d} \theta \, </math> を <math> \mbox{d}G\, </math>  とおいて, これを直線の集合 <math> G \, </math> の密度を表わすとすれば, 直線の集合の測度が式 (2)より


<center><math>m(X)=\int_X \mbox{d}G\, </math> 　　　　　<math>(5)\, </math></center>


と導かれる. ここで, 図2のような滑らかな曲線ABがあり, その長さを <math> L \, </math> とする. AB上の点Pでこの曲線と交わる直線を  <math>g\, </math>, 曲線ABのPにおける接線と <math> g\, </math> とのなす角度を <math> \varphi \, </math>, AからPまでの曲線上の長さを <math> s \, </math> とおく. すると,  <math> (p,\theta) \, </math> から <math> (s,\varphi) \, </math> への変数変換におけるヤコビアンは <math> |\sin\varphi| \, </math> となり, 曲線ABと交わる直線のすべては,  <math> 0\leq s\leq L, 0\leq\varphi \leq \pi \, </math> の範囲を考えればよい. ゆえに


<center>
<math>\int_{0}^{L}\! 
\int_{0}^{\pi}\sin\varphi \; \mbox{d}s \;\mbox{d}\varphi=2L\, </math>
</center> 


となるので, この <math> 2L \, </math> を曲線 <math> AB \, </math> と交わる直線の集合の測度 <math>\textstyle \int\mbox{d}p \; \mbox{d}\theta \, </math> としてよいかにみえる. ところが上記の積分は曲線 <math> AB\, </math> に沿って行なわれたので, この曲線と二つ以上の点で交わる直線については, それぞれの交点で計算されている. したがって, 上記の <math> 2L \, </math> という値は求めたい測度よりも大きい数値であり, <math> n(p,\theta) \, </math> でこの重複を表現すれば


<center><math>\int{n}(p,\theta) \; \mbox{d}p \; \mbox{d}\theta=2L\, </math></center>


が得られる. ところで, 曲線ABを滑らかなものとして議論してきたが, 滑らかな曲線を有限個つなげたものや, 不連続なものについても, ほとんどいたるところで接線は存在するので, 上式は成立する. そして対象となる曲線と交わらない直線については,  <math> n(p,\theta)=0 \, </math> と解釈し, 


<center><math>\int{n}\mbox{d}G=2L\, </math>　　　　　<math>(6)\, </math></center>


が導かれる. さて, これをもとに, 簡単な図形と交わる直線の集合の測度を求めよう. まず上記の曲線ABが直線分であれば, ほとんどいたるところ <math> n=1 \, </math> なので, 長さ <math> L \, </math> の線分と交わる直線の集合 <math> X \, </math> の測度は


<center><table><tr><td align=center>[[画像:sk-0207-c-e-01-2.png]]</td></tr>
<td align=center>図２：曲線と直線</td></table></center>

<!--
*{40mm} %\includegraphics{C-E-01+4.eps} %\caption{曲線と直線} \label{fig:C-E-01+4} %\end{center} %\end{figure} 

\begin{figure}[h] \begin{center} \begin{picture}(300,170)(-50,-20) %\put(-50,-20){\framebox(300,170)} \put(0,0){\vector(1,0){220}} \put(220,0){\makebox(0,0)[l]{\large $ x $ }} \put(0,0){\vector(0,1){130}} \put(0,130){\makebox(0,0)[b]{\large $ y $ }} \put(0,0){\line(2,1){150}} \put(100,-10){\line(-1,2){70}} \put(90,10){\makebox(0,0)[bl]{\large $ g $ }} \put(71,22){\line(-1,2){5}} \put(71,22){\line(2,1){11}} \put(40,30){\makebox(0,0)[tr]{\large $ p $ }} \put(0,0){\makebox(0,0)[tr]{O}} \put(15,0){\oval(11,15)[tr]} \put(30,10){\makebox(0,0)[tl]{\large $ \theta $ }} \put(30,110){\oval(60,80)[br]} \put(56.2,77.6){\vector(2,3){40}} \put(56.6,85.5){\oval(9,8)[t]} \put(56,89.5){\makebox(0,0)[b]{\large $ \varphi $ }} \put(55.5,72){\makebox(0,0)[t]{\large P}} \put(30,70){\makebox(0,0)[r]{\large A}} \put(60,115){\makebox(0,0)[b]{\large B}} \end{picture} \vspace{-5mm} \caption{曲線と直線} \label{fig:C-E-01+4}
-->


<center><math>m(X)=\int\mbox{d}G=2L\, </math>　　　　　<math>(7)\, </math></center>


となる. 次に曲線が凸閉曲線の場合, ほとんどいたるところ <math> n=2 \, </math> が成立するので, 長さ <math> L \, </math> の凸閉曲線と交わる直線の集合 <math> X \, </math> の測度は


<center><table><tr><td align=center>[[画像:sk-0207-c-e-01-3.png]]</td></tr>
<td align=center>図３：<math>C\, </math> と交わる確率</td></table></center>

<!--
%\begin{figure}[p] %\begin{center} %\vspace*{40mm} %\includegraphics{C-E-01+5.eps} %\caption{ $ C $ と交わる確率} \label{fig:C-E-01+5} %\end{center} %\end{figure} 

\begin{figure}[h] \begin{center} \begin{picture}(300,140)(-50,-10) % \put(-50,-10){\framebox(300,140)} \put(-10,130){\line(2,-1){220}} \put(-30,130){\line(1,-1){130}} \put(40,50){\circle{80}} \put(90,60){\oval(200,100)} \put(150,90){\makebox(0,0)[bl]{\large $ C_0 $ }} \put(60,60){\makebox(0,0)[l]{\large $ C $ }} \end{picture} \vspace{-5mm} \caption{ $ C $ と交わる確率} \label{fig:C-E-01+5} \end{center} \end{figure} 
-->


<center><math>m(X)=\int\mbox{d}G=L\, </math>　　　　　<math>(8)\, </math></center>


となる. これを用いると, 以下のような確率を簡単に求めることができる. 図3のように凸閉曲線 <math> C_0,C \, </math> があって, 長さをそれぞれ,  <math> L_0,L \, </math> とする. 曲線 <math> C_0,C \, </math> に交わる直線の集合を <math>X_0,X \, </math> とすれば,  <math> C_0 \, </math> に交わる直線が,  <math> C\, </math>  にも交わる確率 <math> P \, </math> は, 式(8)より
 
 
<center><math>P=m(X)/m(X_0)=L/L_0\, </math>　　　　　<math>(9)\, </math></center>



となる.

----
'''参考文献'''

[1] L.A.Santal&oacute;, ''Integral Geometry and Geometric Probability,'' Addison-Wesley, Reading, Massachusetts, 1976.

[2] 腰塚武志, 「積分幾何学について(1)～(4)」, 『オペレーションズ・リサーチ』, '''21''' (1976), 524-529, 591-596, 654-659, 711-717.

[[category:都市システム|せきぶんきかがく]]