■Divergence の定理の副産物■
Divergence の定理の実用的な応用として、面積の計算があります。
Heat flow をqx=xとqy=y とすると(qx=2xとqy=0 or qx=0とqy=2y でもOK) 、前ページに表示されているDivergence の定理の領域積分は、divq=2ですから、2∬dA になります。
つまり、2倍の領域の面積(Area)です。また、面積を積分で表すと、下式になりますよね。
ですから、Greenの定理から次の式が出来上がります。
上の式を利用すると、任意の境界を持つ領域の面積を計算できます。 例えば、左下の図の様に境界が、直線のセグメントで書かれていると、上の式は、右下の式の様に離散化できます。 この離散化について詳しく知りたい方は、次をクリックして下さい。すると、手書きのノートPage one(52KB), Page two(56KB), Page three(40KB) が画面に表示されます。 表示された画像は、マウスのボタン操作で保存できます。
ここに、Δx=xi+1-xi、Δy=yi+1-yi、
図形の重心(図心)の位置は、断面1次モーメントを計算することにより、求めることができますよね。つまり、重心の位置(X0, Y0)は、次の様に定義されています。
これらの計算をGreenの定理またはDivergenceの定理の式で行う場合、qxとqyを次の様に設定するとX0, Y0が求められます。 X0の計算では、qx=(1/2)x2、qy=0 とします。そして、Y0の計算では、qy=(1/2)y2、qx=0 とします。計算式の離散化については、次のページのJavaScriptを見て理解するか、連絡下さい。