标签:va vb 直线 C++ Yy Yx 交点 Vec2d
这里还是用直线的参数曲线方程:
A: f(t) = va*t+a
B: g(s) = vb*s +b
va是直线A的二维的单位方向向量,a是直线A的原点
计算直线A和B交点, f(t) = g(s) =>
va.x *t + a.x = vb.x*s +b.x
va.y *t + a.y = vb.y*s +b.y
=>
va.x*t - vb.x*s = b.x-a.x
va.y*t - vb.y*s = b.y-a.y
这就是 M*X = Y
其中M是如下的两行两列矩阵
| va.x -vb.x|
|va.y -vb.y|
X是向量 |t s| 转置
Y是列向量
|b.x-a.x|
|b.y-a.y|
其中Yx = b.x - a.x
Yy = b.y-a.y
这个方程组用克拉姆法则求解,不了解的看下线性代数
先M的行列式 m = -va.x*vb.y + va.y*vb.x
如果m 等于0(实际上会用fabs(m)<e判断小到什么限度, e是一个小的浮点值,比如1e-5)
说明两条直线几乎平行,就没有交点,直接返回false
否则,Ms =
| va.x Yx|
|va.y Yy|
ms = va.x*Yy - va.y * Yx
则 s= ms / m
Mt =
| Yx -vb.x|
|Yy -vb.y|
mt = -vb.y*Yx + vb.x* Yy
则 t = mt / m
计算出t和s的值以后,代入f(t)就可以得到交点c的值了,
这个可以不用什么3d向量库,按照公式推导用C++写公式代码就行了,这个一般
在matlab的m文件就叫算例,用C++表达式写出来就行了,或者包个函数
bool lineALineBCross(Vec2d a, Vec2d va, Vec2d b, Vec2d vb, double& t, double& s)
{...}
标签:va,vb,直线,C++,Yy,Yx,交点,Vec2d 来源: https://www.cnblogs.com/abcstar/p/15164493.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。