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

C语言判断点面位置关系系统 一：程序 /*此程序适用于任意复杂程度的单连通和复连通多边形区域*/ #includestdio.h #includestring.h #define M 36 /*给定多边形边数*/ float X[M],Y[M]; float x,y,t,s,Ymin,Ymax,C,D,sn=0,sp=0,xp,yp,l; int i,j=0,k=1; char ch[2]; void input() /*定义输入坐标函数*/ { printf(Please input points made of area: \n); for(i=0;iM;i++) { scanf(%f,%f,X[i],Y[i]); /*从键盘输入围成多边形的顶点坐标*/ } printf(Please input one point: ); scanf(%f,%f,x,y); /*从键盘任意输入一点*/ } float func(float n) /*定义求绝对值的函数*/ { if(n0) { 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=(X[0]+X[2])/2; /*起点与第三点的中点横坐标*/ 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+1]-yp)-(Y[i]-yp)*(X[i+1]-xp))/2; t=func(t); sn=sn+t; if(i==M-2) { sn=sn+func(((X[0]-xp)*(Y[M-1]-yp)-(Y[0]-yp)*(X[M-1]-xp))/2); /*求出凸多边形的总面积并输出*/ } } printf(sn=%f\n,sn); for(i=0;iM-1;i++) { t=((X[i]-x)*(Y[i+1]-y)-(Y[i]-y)*(X[i+1]-x))/2; t=func(t); if(t==0) { if((X[i]X[i+1]xX[i+1]xX[i])||(X[i]X[i+1]xX[i]xX[i+1])) { sp=0; /*如果给定点的横坐标在相邻两点横坐标之间， 则给定点与各个相邻顶点依次围成的面积之和sp赋值为0，并跳出循环*/ break; } else if((X[i]==X[i+1])((Y[i]Y[i+1]yY[i+1]yY[i])||(Y[i]Y[i+1]yY[i]yY[i+1]))) { sp=0; /*如果给定点的纵坐标在相邻两点纵坐标之间，并且相邻两点横坐标相等 则sp=0，并跳出循环*/ break; } else if(x==X[i]y==Y[i]||x==X[i+1]y==Y[i+1]) { sp=0; /*如果给定点的坐标等于相邻两点任意一点的坐标，则sp=0，并跳出循环*/ break; } } sp=sp+t; if(i==M-2) /*最后一点的情况必须单独讨论，同上*/ { l=func(((X[0]-x)*(Y[M-1]-y)-(Y[0]-y)*(X[M-1]-x))/2); if(l==0) { if((X[0]X[M-1]xX[M-1]xX[0])||(X[0]X[M-1]xX[0]xX[M-1])) { sp=0; break; } else if((X[0