求任意多边形的重心

原本是一道题,从中分离出来的:点击打开链接

计算多边形重心方法:

(1)划分多边形为三角形:
以多边形的一个顶点V为源点(V可取输入的第一个顶点),作连结V与所有非相邻顶点的线段,即将原N边形或分为(N-2)个三角形;
(2)求每个三角形的重心和面积:
设某个三角形的重心为G(cx,cy),顶点坐标分别为A1(x1,y1),A2(x2,y2),A3(x3,y3),则有cx = (x1 + x2 + x3)/3.同理求得cy。求面积的方法是s =  ( (x2 - x1) * (y3 - y1) - (x3 - x1) * (y2 - y1) ) / 2,当A1,A2,A3顺时针排列时取-,否则取正(此定理不证)。事实上,在求每个三角形时不需要辨别正负,之后有方法抵消负号,见下述。
(3)求原多边形的重心:
公式:cx = (∑ cx[i]*s[i]) / ∑s[i];  cy = (∑ cy[i]*s[i] ) / ∑s[i];其中(cx[i], cy[i]), s[i]分别是所划分的第i个三角形的重心坐标和面积。由题“ connect the points in the given order”知每个s[i]的正负号相同,故而∑ cx[i]*s[i]能与∑s[i]消号,所以根本不需要在第(2)步判断每个s[i]的正负。另外,在(2)中求每个重心坐标时要除以3,实际上不需要在求每个三角形坐标时都除以3,只需要求出∑ cx[i]*s[i]后一次性除以3即可。即是多边形重心坐标变为:cx = (∑ cx[i]*s[i]) / (3*∑s[i]);  cy = (∑ cy[i]*s[i] ) / (3*∑s[i]);


总结:
每个三角形重心:cx = x1 + x2 + x3;cy同理。
每个三角形面积:s =  ( (x2 - x1) * (y3 - y1) - (x3 - x1) * (y2 - y1) ) / 2;

多边形重心:cx = (∑ cx[i]*s[i]) / (3*∑s[i]);  cy = (∑ cy[i]*s[i] ) / (3*∑s[i]);


#include<stdio.h>
int main()
{
    int n,i;
    int x1,y1,x2,y2,x3,y3;
    double sum_x=0,sum_y=0,sum_s=0;
    scanf("%d",&n);
    scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
    for(i=1; i<=n-2; i++)
    {
        scanf("%d%d",&x3,&y3);
        double s=((x2-x1)*(y3-y1)-(x3-x1)*(y2-y1))/2.0;
        sum_x+=(x1+x2+x3)*s;
        sum_y+=(y1+y2+y3)*s;
        sum_s+=s;
        x2=x3;
        y2=y3;
    }
    printf("%.2lf %.2lf\n",sum_x/sum_s/3.0,sum_y/sum_s/3.0);
    return 0;
}



计算任意多边形的几何中心和重心需要用到多边形的顶点坐标。以下是使用MATLAB计算多边形几何中心和重心的方法: 假设多边形顶点坐标为 $P_1(x_1,y_1), P_2(x_2,y_2),...,P_n(x_n,y_n)$。 1. 计算几何中心 几何中心是所有顶点坐标的平均值,即: $$x_c=\frac{1}{n}\sum_{i=1}^{n} x_i$$ $$y_c=\frac{1}{n}\sum_{i=1}^{n} y_i$$ MATLAB代码如下: ```matlab % 输入顶点坐标 x = [x1, x2, ..., xn]; y = [y1, y2, ..., yn]; % 计算几何中心 xc = mean(x); yc = mean(y); ``` 2. 计算重心 重心是所有三角形重心的加权平均值,其中每个三角形的重心是三个顶点坐标的平均值。具体而言,对于第 $i$ 个三角形,其重心坐标为: $$x_i=\frac{1}{3}(x_{i1}+x_{i2}+x_{i3})$$ $$y_i=\frac{1}{3}(y_{i1}+y_{i2}+y_{i3})$$ 其中 $(x_{i1},y_{i1})$、$(x_{i2},y_{i2})$ 和 $(x_{i3},y_{i3})$ 是第 $i$ 个三角形的三个顶点坐标。 重心的坐标为所有三角形重心坐标的加权平均值,其中每个三角形的面积为其三个顶点构成的三角形面积的和,即: $$x_g=\frac{\sum_{i=1}^{n} S_i x_i}{\sum_{i=1}^{n} S_i}$$ $$y_g=\frac{\sum_{i=1}^{n} S_i y_i}{\sum_{i=1}^{n} S_i}$$ 其中 $S_i$ 是第 $i$ 个三角形的面积。 MATLAB代码如下: ```matlab % 输入顶点坐标 x = [x1, x2, ..., xn]; y = [y1, y2, ..., yn]; % 计算三角形面积和重心 S = 0; xg = 0; yg = 0; for i = 1:n-2 for j = i+1:n-1 for k = j+1:n % 计算第 i,j,k 个顶点组成的三角形面积 Si = abs((x(j)-x(i))*(y(k)-y(i))-(x(k)-x(i))*(y(j)-y(i)))/2; % 计算三角形重心坐标 xi = (x(i) + x(j) + x(k))/3; yi = (y(i) + y(j) + y(k))/3; % 更新重心坐标和三角形面积 xg = xg + Si*xi; yg = yg + Si*yi; S = S + Si; end end end % 计算重心坐标 xg = xg/S; yg = yg/S; ``` 注意,如果多边形是凸多边形,则几何中心和重心重合。如果多边形是凹多边形,则重心会偏离几何中心。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值