语言判断点面位置关系系统.docVIP

C语言判断点面位置关系系统 一：程序 /*此程序适用于任意复杂程度的单连通和复连通多边形区域*/ #includestdio.h #includestring.h #define M 36 /*给定多边形边数*/ float X[M]，Y[M]; float x，y，t，s，Ymin，Ymax，C，D，sn=O，sp=O,xp，yp,l; int i，j=O，k=l; char ch[2]; void input() /*定义输入坐标函数*/ { printf(”Please input points made of area: \n”); for(i=0;iM;i++) { scanf(%f，％r，X[i]，Y[i]); /*从键盘输入围成多边形的顶点坐标*/ } printf(”Please input one point:); scanf(n%f，％f，x，y); /* 从键盘任意输入一点 */ } float func(float n) /*定义求绝对值的函数*/ { if(n〉0) ( return n; } else { return (-n); } } float Min(float m,float n) /*定义求最小值的函数*/ { if(m=n) { return m; } else return n; float Max(float m,float n) /*定义求最大值的函数 { if(m〉=n) { return m; } else { return n; } } void simplefun() /*定义判断点是否在凸多边形上、内、外函数*/ { xp=(XL0]+Xl2J)/2; /*起点与第三点的巾点横坐.标5V yp=(Y[0]+Y[2])/2; /*起点与第三点的中点纵坐标*/ printf(xp=%f\typ=%f\n’’，xp，yp); for(i=0;iM-1 ;i++) { t=((X[i]-xp)*(Y[i+l ]-yp)-(Y[i]-yp)*(X[i+l ]-xp))/2; t 二 func(t); sn=sn+t; if(i==M-2) { sn=sn+func(((X[0]-xp)*(Y[M-l]-yp)-(Y[0]-yp)*(X[M-l]-xp))/2); /* 求出凸多 边形的总面积并输出*/ } } printf(sn=%f\n”，sn); for(i=0;iM-l;i++) { t=((X[i]-x)*(Y[i+l ]-y)-(Y[i]-y)*(X[i+l l-x))/2; t=func(t); if(t==O) { if((X[i]〉X[i+l]xX[i+l]xX[i])||(X[i]X[i+l]x〉X[i]xX[i+l])) { sp=0; /*如果给定点的横坐标在相邻两点横坐标之间， 则给定点与各个相邻顶点依次围成的面积之和sp赋值为0,并跳 出循环*/ break; } else if((X[i]==X[i+1 ])(( Y[i]Y[i+1 ]yY[i+1 ]y Y[i])||(Y[i] Y[i+1 ]yY[i]yY[i+1 ]))) sp=O; /*如果给定点的纵坐标在相邻两点纵坐标之间，并且相邻两点横 坐标相等 则sp=O，并跳出循环*/ break; } else if(x==X[i]y==Y[i]||x==X[i+l ]y==Y[i+l ]) { sp=O; /*如果给定点的坐标等于相邻两点任意一点的坐标，则sp=O，并 跳出循环*/ break; } } sp=sp+t; if(i=M-2) /*最后一点的情况必须单独讨论，同上*/ { l=func(((X[O]-x)*(Y[M-l ]-y)-(Y[O]-y)*(X[M-l ]-x))/2); if(l==O) if((X[O]X[M-1 1x〉X[M- 1 1xX[01)||(X[01X[M- 1 1x〉X[01xX[M- 1 ])) { sp=O; break; } else if((X[O]==X[M-1 ])((Y[O]Y[M-1 ]yY[M-l ]yY[O])||( Y[0] Y[M-1 ]yY[O]y Y[M-1]))) { sp=O; break; } else if(x==X[O]y==Y[O]||x==X[M-1 ]y==Y[M-1 ]) { sp=O; break; } } sp=sp+l; /*若前面的if都不执行，则sp就是给定点与各个相邻顶点依次围 成的真实的面积之和*/ } } 终得出结论*/ if(sp〉sn) printfCThis point (%f,%f) is out of the region\n,x,y); } else if(sp==sn) { print