求算两条十字交叉直线交点的程序

Question

求算两条十字交叉直线交点的程序

Answer 1

// VC 6.0 编译通过
#include
// 求两直线交点
// 直线 1 上的两点：(ax1,ay1),(ax2,ay2)
// 直线 2 上的两点：(bx1,by1),(bx2,by2)
// 返回值：true:有交点,保存在 (ix,iy) 中
// ：false:两直线平行
bool Intersection(double ax1,double ay1,double ax2,double ay2,double bx1,double by1,double bx2,double by2,double *ix,double *iy)
{
double d = (ay2-ay1) * (bx2-bx1) - (by2-by1) * (ax2-ax1);
if (d==0)
return false; // 两直线平行就退出,避免除数为 0 的情况
*ix = (
(ax2-ax1)*(bx2-bx1)*(by1-ay1)+(ay2-ay1)*(bx2-bx1)*ax1-(by2-by1)*(ax2-ax1)*bx1
) / d;
*iy = (
(ay2-ay1)*(by2-by1)*(bx1-ax1)+(ax2-ax1)*(by2-by1)*ay1-(bx2-bx1)*(ay2-ay1)*by1
) / -d;
return true;
}
void main()
{
double ix,iy;
// 简化示例,没有写 cin 语句.
// 示例,求 (7,3)-(5,8) 和 (3,5)-(10,6) 的交点
if ( Intersection(7,3,5,8,3,5,10,6,&ix,&iy) )
cout