■Divergence の定理の副産物■
Divergence の定理の実用的な応用として、面積の計算があります。Heat flow をq_x=xとq_y=y とすると(q_x=2xとq_y=0 or q_x=0とq_y=2y でもOK) 、前ページに表示されているDivergence の定理の領域積分は、2∬dA になります。つまり、2倍の領域の面積(Area)です。また、normal vector n は、

上の式を利用すると、任意の境界を持つ領域の面積を計算できます。例えば、左下の図の様に境界が、直線のセグメントで書かれていると、上の式は、右下の式の様に離散化できます。

ここに、Δx=xi+1-xi、Δy=yi+1-yi、 xm=(xi+xi+1)/2.、ym=(yi+yi+1)/2. です。 n は、多角形のセグメントの数です。Index の i+1 が n+1 になる場合は、i+1 を 1 とします。

図形の重心(図心)の位置は、断面1次モーメントを計算することにより、求めることができます。つまり、重心の位置(X₀, Y₀)は、次の様に定義されています。

これらの計算をDivergenceの定理の式で行う場合、q_xとq_yを次の様に設定するとX₀, Y₀が求められます。 X₀の計算では、q_x=(1/2)x²、q_y=0 とします。そして、Y₀の計算では、q_y=(1/2)y²、q_x=0 とします。計算式の離散化については、次のページのJavaScriptをみて下さい。