用格林定理求Jordan曲线包围的封闭区域的面积

发现很多教材讲微积分中的格林定理忽略其引申，显得粗糙。看了不同版本教材比对之后，这种感受更深了。 Green′stheorem 联系着二重积分和第二类平面曲线积分,是个漂亮的结果.

对原始定理稍作引申，不仅加深理解，在计算几何的某些算法实现中灵活应用起来也很方便。不但格林定理，散度定理也有类似的应用，让人惊讶。

定理 ( Green's theorem, also Jordan curve theorem ): 向量场 (x,y)−→F⟨X,Y⟩ , ∂X∂y 和 ∂Y∂x 在有界单连通区间 D 上连续；D 的边界为Jordan曲线 C ：r(s)=⟨x(s),y(s)⟩。则环量形式 circulation form:

适当更换 F ，⟨x,0⟩, ⟨0,y⟩ , 12⟨x,y⟩ ，可以类似得到flux form的等价的引申结果。

用户评价：使用了上面公式之后发现，

Jordan曲线作为适用范围只是定理适用的充分但非必要条件。边界为封闭曲线，区域内任何点关于边界曲线的 卷绕数 都为 1 似乎更理想, 尤其是如果面积也允许有符号时; 因为这样得到的面积与曲线方向有关, 否则也有反例存在： Lemniscates of Bernoulli, 跟曲线参数方程特定的形式有关( 因为这影响曲线的走向),此时极坐标积分反而可能更好。

http://math.uga.edu/~pete/handouteight.pdf

国外教科书中多讲椭圆面积的计算。大概因为这个最具有代表性。用其它方法计算都没有这个省力。 比如 C=⟨acost,bsint⟩,t∈[0,2π] 所包围的区域的面积。

格林定理的引申： 令向量场 F=12⟨−y,x⟩ , 利用定理的引申:

如果用极坐标呢？

因为参数方程中的参数 t 并不是极坐标下的角度 θ, 所以变极坐标形式还是有一点点繁琐的。先写成隐函数形式：

如果用一个形状更特殊, 常规方法繁琐的例子则容易给人更深的印象。

比如：

这个例子四个花瓣有公共点(0,0)应该就不算Jordan曲线了,但封闭的开区域内 卷绕数 都为 1 且曲线方向一致(如果方向相反可能有负卷绕数), 所以仍然适用格林定理；此外，四个象限中,曲线都正好等分单位正方形,这让人也感觉很好。

此外, 这个案例 ⟨sin2t,sint⟩ 也不错。

前面提到这类曲线能否适用Green’s theorem实际上跟曲线参数方程特定的形式有关(因为这影响曲线的走向), 这里举两个例子, 同样都是 Lemniscates of Bernoulli , 但参数方程决定的曲线走向差异, 一个适用, 一个太难适用。

先说适用的:

参数方程

而这条曲线:

调整参数方程的形式之后用前面方法计算它的面积其实还是可以得到解析形式 751−−√−π2+cot−1(751−−√)≈48.43920 的。这类曲线可以看成可以应用格林定理的不同封闭区域之间单连通。

参数方程为:

如果直接用 格林定理 则得到 87π , 跟半径 87−−√ 的圆一比较, 发现实际面积并没有圆大，显然是错误的。中心卷绕数为2、近似于正五边形的部分实际被重复计算了。

它的面积的计算是这样实现的，也用到了Green定理!

(James Stewart, Calculus, 7th Ed. p1111) 用于计算不规则图形区域的面积, 比如生物学里测量植物叶片的面积或鸟类翅膀的面积, 医学上器官或肿瘤切面面积, 以及地图学或遥感林业研究中地图或遥感影像图上特定地理区域的面积.(数字化时代可能已经有更简单快捷的替代方法)

仪器原理 即格林定理,介绍性材料(Oliver Knill, 12/9/2000)

The planimeter is a mechanical device for measuring areas in the plane. It has the shape of a ruler with two legs. One leg of length 1 connects the fixed origin (0,0) to (a,b) . A second leg of length 1 connects (a,b) with the end point (x,y) . The point (x,y) determines (a,b) as the intersection of two unit circles centered at (0,0) and (x,y) . The intersection is unique if the angle between the two planimeter legs is smaller than 180 degrees.

The planimeter vector field Let F(x,y)=⟨P(x,y),Q(x,y)⟩ be the Planimeter vector field. It is defined by attaching a unit vector orthogonal to the vector (x−a,y−b) at (x,y) . The wheel rotation is the line integral of F along the boundary of R. Green’s theorem tells that this integral is the double integral of curl(F) over the region R . The planimeter vector field is explicitely given by F(x,y)=(P(x,y),Q(x,y))=(−(y−b(x,y)),(x−a(x,y))). Furthermore, curl(F)=Qx−Py is equal to 2+(−ax−by) which is 2 plus the curl of the vector field G(x,y)=⟨b(x,y),−a(x,y)⟩. A direct verification shows that curl(G)=−1 . The Planimeter lineintegral is therefore the area of the enclosed region.

上面这段原理,前半部分可以看看; 后面讲到向量场的旋量时似乎糊涂了, 这个向量场的旋量实际上绝对值就是 1 , 只需要把两个函数a(x,y),b(x,y)解析形式的两种可能显式表示出来一计算就有了,讲得太罗嗦了.

最早的原理介绍1947年就有了:

Bibliography. Green’s theorem is the classic way to explain the planimeter. The explanation of the planimeter through Green’s theorem seems have been given first by G. Ascoli in 1947:

更新已完成