计算几何初步 线 段 属 性

    第一个知识点：三点求圆心： http://www.cnblogs.com/hxsyl/p/3223000.html

  补充一点关于圆的知识：http://wenku.baidu.com/view/8ceabb8cb9d528ea81c779ac.html

最近为了研究最小覆盖圆，需要求圆心，所以随手download了几份资料。

若是平面点得到的是平面圆，若是三维空间点得到的是空间圆(不是球)。

一.第一种方法

这种方法当k1或者k2为0时需要分别讨论，不便于计算机编程实现。

二.第二种方法

       这种方法下只要三点不共线就有结果，三点共面的充要条件是混合积(a*(bXc))为0。

 

 

第二个知识点：求直线交点和三角形外心: http://www.cnblogs.com/hxsyl/p/3225698.html

 

一.求两直线交点

class Point {
    double x;
    double y;
    
    public Point() {
        this.x = 0;
        this.y = 0;
    }
}

class Line {
    Point a;
    Point b;
    
    public Line() {
        this.a = new Point();
        this.b = new Point();
    }
    //求两直线的交点，斜率相同的话res=u.a
    Point intersection(Line u,Line v){
        Point res = u.a;
        double t = ((u.a.x-v.a.x)*(v.b.y-v.a.y)-(u.a.y-v.a.y)*(v.b.x-v.a.x))
            /((u.a.x-u.b.x)*(v.b.y-v.a.y)-(u.a.y-u.b.y)*(v.b.x-v.a.x));
        res.x += (u.b.x-u.a.x)*t;
        res.y += (u.b.y-u.a.y)*t;
        return res;
    }

二.求三角形外心

        1. 垂心: 三角形三条边上的高相交于一点.这一点叫做三角形的垂心.

        2. 重心: 三角形三条边上的中线交于一点.这一点叫做三角形的重心.

        3. 外心: 三角形三边的中垂线交于一点.这一点为三角形外接圆的圆心.

        4. 内心三角形三内角平分线交于一点.这一点为三角形内切圆的圆心.

        已知圆的3点，先求出3边长,由海伦公式得出面积S=sqrt(p*(p-a)*(p-b)*(p-c)) p=(a+b+c)/2;由三角形面积公式S=1/2*a*b*sin(C)和正弦定理a/sin(A)=b/sin(B)=c/sin(C)=直径(根据相同弦长对应的圆周角相同可证正弦定理)可得直径=a*b*c/2/S。

        求圆心坐标。利用：G是⊿ABC外心的充要条件是(向量GA+向量GB)·向量AB= (向量GB+向量GC)·向量BC=(向量GC+向量GA)·向量CA=向量0.
这个性质的证明很容易的，只需要想到外心是中垂线交点即可，就可以证明这个性质了，利用向量可以避免求斜率，以及考虑斜率不存在等很多情况。

//三角形外接圆圆心(外心)
    Point center(Point a,Point b,Point c) {
        //加上这个才没有编译器提示未初始化，因为new所以也写了构造方法
        Line u = new Line(),v = new Line();
        u.a.x=(a.x+b.x)/2;
        u.a.y=(a.y+b.y)/2;
        u.b.x=u.a.x+(u.a.y-a.y);
        u.b.y=u.a.y-(u.a.x-a.x);
        v.a.x=(a.x+c.x)/2;
        v.a.y=(a.y+c.y)/2;
        v.b.x=v.a.x+(v.a.y-a.y);
        v.b.y=v.a.y-(v.a.x-a.x);
        return intersection(u,v);
    }

三.求三角形内心

        由于内心到各边距离就是半径r，可以把三角形分成三部分，再根据海伦公式得到半径r=2*S/(a+b+c)。

       内切圆心坐标(x,y): 三角形三个顶点的坐标:A(x1,y1),B(x2,y2),C(x3,y3)则圆心为x=(x1*BC+x2*CA+x3*AB)/(AB+BC+CA)、y=(y1*BC+y2*CA+y3*AB)/(AB+BC+CA)。

        证明：内心是角平分线的交点，到三边距离相等.
　　设:在三角形ABC中,三顶点的坐标为：A(x1,y1),B(x2,y2),C(x3,y3) BC=a,CA=b,AB=c，内心为M （X，Y）则有aMA+bMB+cMC=0（三个向量） ，MA=（X1-X，Y1-Y） ，MB=(X2-X,Y2-Y) ，MC=(X3-X,Y3-Y)
　　则：a(X1-X)+b(X2-X)+c(X3-X)=0，a(Y1-Y)+b(Y2-Y)+c(Y3-Y)=0
　　∴X=（aX1+bX2+cX3)/(a+b+c)，Y=（aY1+bY2+cY3)/(a+b+c)
　　∴M（（aX1+bX2+cX3)/(a+b+c)，（aY1+bY2+cY3)/(a+b+c)）。

  

已知O为三角形ABC的内心,a,b,c分别是A.B.C边所对边长. 则aOA+bOB+cOC=0(OA,OB,OC均指向量) 


证明：设三角形ABC，AD为BC边上的角平分线，内心为O。

|BC|=a，|AC|=b，|AB|=c

aOA+bOB+cOC

=aOA+b(AB+OA)+c(AC+OA)

=(a+b+c)OA+b(DB-DA)+c(DC-DA)

设BC的方向向量e,则DB=e|DB|,DC=-e|DC|

又由角平分线定理，|DB|/|DC|=c/b，所以bDB+cDC=0

(a+b+c)OA+b(DB-DA)+c(DC-DA)= (a+b+c)OA- b DA- c DA =aOA+(b+c)OD

又因为OA、OD反向，用角平分线定理和合比定理：

b/CD=c/BD=(b+c)/(CD+BD)=(b+c)/a, b/CD=OA/OD,

所以OA/OD=(b+c)/a , 又因为OA、OD反向，

故aOA+bOB+cOC=aOA+(b+c)OD =0.

 第三个知识点：

点到直线到线段距离：

http://www.cnblogs.com/hxsyl/p/3249763.html

一.点到直线距离

        已知一个点P(X0, Y0), 求点到直线Ax + By + C = 0的距离公式为：d = [AX0 + BY0 + C的绝对值]/[(A^2 + B^2)的算术平方根]，如求点P(-1, 2)到直线2X + Y - 10 = 0的距离：X0 = -1, Y0 = 2, A = 2, B = 1, C = -10 代入公式
d =[2 * (-1) + 1 * 2 - 10 的绝对值] / 根号[2 * 2 + 1 * 1] = 10 / 根号5。

        已知两点的坐便(x1, y1)，(x2, y2) ，另外一个点的坐标是(x0, y0); 求(x0, y0)到经过(x1, y1) (x2, y2)直线的距离。
        直线方程中 A = y2 - y1,B = x1- x2,C = x2 * y1 - x1 * y2(叉积);点的直线的距离公式为: double d = (fabs((y2 - y1) * x0 +(x1 - x2) * y0 + ((x2 * y1) -(x1* y2)))) / (sqrt(pow(y2 - y1, 2) + pow(x1 - x2, 2)))。

二.点到线段最短距离

private static double distance(Point p, Point p1) {
    return Math.hypot(p.x-p1.x, p.y-p1.y);
  }
  //点到线段的最短距离,x0,y0是圆心
  private static double pointToLine(Point p1,Point p2, Point p) {
    double ans = 0;
    double a, b, c;
    a = distance(p1, p2);
    b = distance(p1, p);
    c = distance(p2, p);
    if (c+b==a) {//点在线段上
      ans = 0;
      return ans;
    }
    if (a<=0.00001) {//不是线段，是一个点
      ans = b;
      return ans;
    }
    if (c*c >= a*a + b*b) { //组成直角三角形或钝角三角形，p1为直角或钝角
      ans = b;
      return ans;
    }
    if (b * b >= a * a + c * c) {// 组成直角三角形或钝角三角形，p2为直角或钝角
      ans = c;
      return ans;
    }
    // 组成锐角三角形，则求三角形的高
    double p0 = (a + b + c) / 2;// 半周长
    double s = Math.sqrt(p0 * (p0 - a) * (p0 - b) * (p0 - c));// 海伦公式求面积
    ans = 2*s / a;// 返回点到线的距离（利用三角形面积公式求高）
    return ans;
  }