一般方程与参数方程求直线交点

一般方程与参数方程求直线交点

一、             一个例子:


如上图,有两条直线,设L1L2L1上有两点(0, 0)(10,10)L2上有两点(0,10)(10,0),它们的交点是(5,5)。求解交点有两种效率较高的常用方法,一般方程法与参数方程法,以下将分别描述其原理及实现。

二、             一般方程法:

直线的一般方程为F(x) = ax + by + c = 0。既然我们已经知道直线的两个点,假设为(x0,y0), (x1, y1),那么可以得到a = y0 – y1, b = x1 – x0, c = x0y1 – x1y0

因此我们可以将两条直线分别表示为

F0(x) = a0*x + b0*y + c0 = 0, F1(x) = a1*x + b1*y + c1 = 0

那么两条直线的交点应该满足

a0*x + b0*y +c0 = a1*x + b1*y + c1

由此可推出

x = (b0*c1 – b1*c0)/D

y = (a1*c0 – a0*c1)/D

D = a0*b1 – a1*b0 (D0时,表示两直线重合)

源代码:


复制代码
 1 #include <iostream>
 2 using namespace std;
 3 
 4 typedef struct
 5 {
 6     int x, y;
 7 } Point;
 8 int main()
 9 {
10     //一般方程法
11     Point line1[2], line2[2];
12     int a[2], b[2], c[2], x, y, D;
13     cout << "Frist Line(x0 y0 x1 y1):";
14     cin >> line1[0].x >> line1[0].y >> line1[1].x >> line1[1].y;
15     cout << "Second Line(x0 y0 x1 y1):";
16     cin >> line2[0].x >> line2[0].y >> line2[1].x >> line2[1].y;
17     
18     a[0= line1[0].y - line1[1].y;b[0= line1[1].x - line1[0].x;
19     c[0= line1[0].x * line1[1].y - line1[1].x * line1[0].y;
20     a[1= line2[0].y - line2[1].y;b[1= line2[1].x - line2[0].x;
21     c[1= line2[0].x * line2[1].y - line2[1].x * line2[0].y;
22     D = a[0* b[1- a[1* b[0];
23     if (D != 0)
24     {
25         x = (b[0* c[1- b[1* c[0]) / D; y = (a[1* c[0- a[0* c[1]) / D;    
26         cout << "一般方程求解的交点:" << x << "," << y << endl;
27     }
28     else
29     {
30         cout << "两直线重合" << endl;
31     }
32     return 0;
33 }
34 
复制代码


三、             参数方程法:

设直线上的两个点为A0(x0, y0), B0(x1, y1),那么线段可用向量

 =  –=(x1 – x0, y1 – y0)

表示,的方向即为直线上的方向,那么直线上的任意点便可表示为

s0(t0) =  +  t0 * 

同理,另一条直线(A1(x2, y2), B1(x3, y3)为其上两个点)可类似表示为

s1(t1) =  +  t1 * 

因此,满足两条直线的交点必满足以下条件:

 + t0 *  + t1 * 

可求出

t0 = (x0(y3 – y2) + x2(y0 – y3) + x3(y2 – y0))/D

t1 = -(x0(y2 – y1) + x1(y0 – y2) + x2(y1 – y0))/D

D = x0(y3 – y2) + x1(y2 – y3) + x3(y1 – y0) + x2(y0 – y1),(D0时,表示两直线重合)

则交点为:

       x = x0 + t0 * (x1 – x0); y = y0 + t0 * (y1 – y0);

       

       x = x2 + t1 * (x3 – x2); y = y2 + t1 * (y3 – y2);

如果将t0,t1限定在[0,1]内,则变为求线段的求点

源代码:


复制代码
 1 #include <iostream>
 2 using namespace std;
 3 
 4 typedef struct
 5 {
 6     int x, y;
 7 } Point;
 8 int main()
 9 {
10     //参数方程法
11     Point pt[4];
12     int t1, t2, dx, dy, D, x, y;
13     cout << "Frist Line(x0 y0 x1 y1):";
14     cin >> pt[0].x >> pt[0].y >> pt[1].x >> pt[1].y;
15     cout << "Second Line(x0 y0 x1 y1):";
16     cin >> pt[2].x >> pt[2].y >> pt[3].x >> pt[3].y;    
17     t1 = pt[0].x * (pt[3].y - pt[2].y) + pt[2].x * (pt[0].y - pt[3].y) + pt[3].x * (pt[2].y - pt[0].y);
18     t2 = - (pt[0].x * (pt[2].y - pt[1].y) + pt[1].x * (pt[0].y - pt[2].y) + pt[2].x * (pt[1].y - pt[0].y));
19     D = pt[0].x * (pt[3].y - pt[2].y) + pt[1].x * (pt[2].y - pt[3].y) + pt[3].x * (pt[1].y - pt[0].y) + pt[2].x * (pt[0].y - pt[1].y);
20     if (D != 0)
21     {
22         dx = pt[1].x - pt[0].x; dy = pt[1].y - pt[0].y;
23         x = pt[0].x + t1 * dx / D; y = pt[0].y + t1 * dy / D;
24         cout <<"参数法求交点:" << x << "," << y <<  endl;
25     }
26     else
27     {
28         cout << "两直线重合" << endl;
29     }
30     return 0;
31 }
32 
复制代码

展开阅读全文

没有更多推荐了,返回首页