《数学实验教程_实验19(格林公式及其应用)》由会员分享，可在线阅读，更多相关《数学实验教程_实验19(格林公式及其应用)（4页珍藏版）》请在金锄头文库上搜索。
1、实验 19 格林公式与形心实验目的作为计算二重积分的一种手段，格林公式把函数在一个域上的二重积分和一个相关的向量函数在围绕这个区域的边界的线积分联系起来，本实验中我们将要考虑如何求一些比较复杂的平面区域的形心。预备知识 曲线积分、格林公式、形心实验内容计算复杂区域的形心【原 理】【格林公式格林公式】设区域D由分段光滑的曲线L围成，函数及在D上具有一阶连P x y( , )Q x y( , )续偏导数，则有DLQPxydxdyPdxQdy() 其中L是D的区正向的边界曲线。格林公式的三种特殊情况：格林公式的三种特殊情况：（）若，则；Py Qx, DLdxdyxdyydx1 2() （）若，则；x
2、PQ220 , DLxdxdyx dy21 2 （）若，则。PyQ21 2,0 DLydxdyy dx21 2 （）可以用来计算区域D的面积， （2）和（）可以用来计算区域D关于y轴和轴的矩。【步 骤】：众所周知，平面区域D的形心的坐标由公式x y( , )DDDDxdxdyydxdy xydxdydxdy, 给出，其中分子的积分是关于区域的坐标轴的矩，分母的积分是区域的面积。对于形状比较好的区域，这些几分可以通过累次积分算出，但是对于比较复杂的不规则的区域，如一个国家的州（省）的几何中心，用一定的方式作积分计算几乎是不可能的。为了计算形心坐标，一般可采用如下两步骤：实验 19 格林公式及其应
3、用 - 107 -【Step1】：用一条简单曲线（如多边形）路径近似代替区域的边界；【Step2】：利用格林公式将多边形区域上的二重积分化为围绕边界的线积分。【公式推导】：下面推导区域D为多边形，其中nnnnxyxyxyxy112211(,),(,),(,),(,) L L时的公式。我们只需要分别计算出分子和分母即可。nnxxyy1111, 设连接点和的直线段的参数方程为iixy(,)iixy11(,) iLiiiix tt xt xty tt yt y11( )(1),0,1( )(1) 则1 2 111 11112011 112 1213211 12231221(1)()(1)()(),i
4、nDLiniiiiiiii iniiii innnndxdyxdyydxt xt xyyt yt yxxdtx yxyyyyyxxxxxxyy L LL LMMMM12211 1122011221 1116 121322222221 112222331161(1) ()()(),inniiiiDLiiniiiiii innnnnnxdxdyx dyt xt xyy dtxx xxyyyyyyxx xxxx xxxx xxyy L LMM- 108 - 第二章 专题实验12211 1122011221 1116 121322222221 112222331161(1) ()()(),inniii
5、iDLiiniiiiii innnnnnydxdyy dxt yt yxx dtyy yyxxxxxxyy yyyy yyyy yyxx L LMM有了上面的公式，我们可以利用它们来计算顶点为(0,0)，(1,0)，(2,1)，(3,2)，(1,4)的五边形区域的形心。利用所附程序计算得到。x y164 39( , )( ,) 实验练习实验练习选择一个复杂的平面图形，然后利用上述方法计算其形心。参考程序：参考程序：(Matlab(Matlab程序程序) )clearclear n=5;n=5; x=0x=0 1 1 2 2 3 3 1 1 0;0; y=0y=0 0 0 1 1 2 2 4 4
6、 0;0; S=(x(1:n)*y(2:n+1)-x(2:n+1)*y(1:n)/2;S=(x(1:n)*y(2:n+1)-x(2:n+1)*y(1:n)/2; Mx=1/6*(x(1:n).2+x(1:n).*x(2:n+1)+x(2:n+1).2)*(y(2:n+1)-y(1:n);Mx=1/6*(x(1:n).2+x(1:n).*x(2:n+1)+x(2:n+1).2)*(y(2:n+1)-y(1:n); My=-1/6*(y(1:n).2+y(1:n).*y(2:n+1)+y(2:n+1).2)*(x(2:n+1)-x(1:n);My=-1/6*(y(1:n).2+y(1:n).*y(
7、2:n+1)+y(2:n+1).2)*(x(2:n+1)-x(1:n); x0=Mx/Sx0=Mx/S y0=My/Sy0=My/S（mathematicsmathematics程序）程序）AreaMomentvertices_:=Modulepts,AreaMomentvertices_:=Modulepts,pts=Appendvertices,vertices1;pts=Appendvertices,vertices1;Sum(ptsi,1Sum(ptsi,1 ptsi+1,2-ptsi+1,1ptsi+1,2-ptsi+1,1 ptsi,2),i,1,Lengthpts-ptsi,2)
8、,i,1,Lengthpts-1/21/2xMomentvertices_:=Modulepts,xMomentvertices_:=Modulepts,pts=Appendvertices,vertices1;pts=Appendvertices,vertices1;Sum(ptsi,12+ptsi,1Sum(ptsi,12+ptsi,1 ptsi+1,1+ptsi+1,12)(ptsi+1,2-ptsi+1,1+ptsi+1,12)(ptsi+1,2-ptsi,2),i,1,Lengthpts-1/6ptsi,2),i,1,Lengthpts-1/6实验 19 格林公式及其应用 - 109 -yMomentvertices_:=Modulepts,yMomentvertices_:=Modulepts,pts=Appendvertices,vertices1;pts=Appendvertices,vertices1;-Sum(ptsi,22+ptsi,2-Sum(ptsi,22+ptsi,2 ptsi+1,2+ptsi+1,22)(ptsi+1,1-ptsi,1),i,1,Lengthpts-ptsi+1,2+ptsi+1,22)(ptsi+1,1-ptsi,1),i,1,Lengthpts-1/61/6
