如何在MATLAB中求角度的平分线
我有一个与此代码相关的问题：
t = -20:0.1:20;
plot3(zeros(size(t)),t,-t.^2);
grid on
hold on
i = 1;
h = plot3([0 0],[0 t(i)],[0 -t(i)^2],'r');
h1 = plot3([-1 0],[0 0],[-400 -200],'g');
for(i=2:length(t))
set(h,'xdata',[-1 0],'ydata',[0 t(i)],'zdata',[-400 -t(i)^2]);
pause(0.01);
end在这段代码中，我绘制了两条相交线。H1和H2。H1是固定的，H2作为时间的函数移动。在本例中，H2碰巧跟踪一条抛物线，但它的运动可能是任意的。
对于线H2的每个位置，如何计算并绘制这两条相交线之间的角度的平分线？我希望在图中看到平分线和H2线同时移动。
解决H2的一个位置的这个问题就足够了，因为对于H2相对于H1的所有方向都是相同的过程。
回答 3
Stack Overflow用户
发布于 2010-10-14 00:07:02
我不是几何天才，可能有一种更简单的方法来做到这一点。到目前为止，还没有人回应，所以这将是一件很有意义的事情。
你在三个空间中有三个点:设A是两条线段的公共顶点。
设B和C是两条直线段上的两个已知点。
选择任意距离r，其中
从A到B的R <=距离
从A到C的R <=距离
从A沿线段AB测量r的距离。这是从A沿线段AC测量的点RB距离或r。这是点RC
找到连接RB和RC的线段的中点。这是M点
线段AM是角度CAB的角度平分线。
这些步骤中的每一个都应该相对容易完成。
Stack Overflow用户
发布于 2010-10-14 07:08:03
这基本上是MatlabDoug的方法，对他称之为M的点的确定进行了一些改进。
t = -20:0.1:20;
plot3(zeros(size(t)),t,-t.^2);
grid on
hold on
v1 = [1 0 200];
v1 = v1/norm(v1);
i = 1;
h = plot3([-1 0],[0 t(i)],[-400 -t(i)^2],'r');
h1 = plot3([-1 0],[0 0],[-400 -200],'g');
l = norm([1 t(i) -t(i)^2+400]);
p = l*v1 + [-1 0 -400];
v2 = (p + [0 t(i) -t(i)^2])/2 - [-1 0 -400];
p2 = [-1 0 -400] + v2/v2(1);
h2 = plot3([-1 p2(1)],[0 p2(2)],[-400 p2(3)],'m');
pause(0.1)
for(i=2:length(t))
l = norm([1 t(i) -t(i)^2+400]);
p = l*v1 + [-1 0 -400];
v2 = (p + [0 t(i) -t(i)^2])/2 - [-1 0 -400];
p2 = [-1 0 -400] + v2/v2(1);
set(h,'xdata',[-1 0],'ydata',[0 t(i)],'zdata',[-400 -t(i)^2]);
set(h2,'xdata',[-1 p2(1)],'ydata',[0 p2(2)],'zdata',[-400 p2(3)]);
pause;
endStack Overflow用户
发布于 2020-09-04 12:32:21
我只使用以下代码：
- 找到归一化矢量AB和AC，其中A是线段的公共点。
- V= (AB + AC) * 0.5 //生成将AB和AC一分为二的方向矢量。
对V进行标准化，然后执行A+V* length，以获得从公用点开始的所需长度的线段。
(请注意，此方法不适用于沿线的3个点来生成垂直平分线，在这种情况下，它将生成一个没有长度的向量)
我已经添加了一个C#实现(在XZ平面中使用Unity3Dmatlab结构)，它可以处理垂直平分线和反射平分线，以防有人知道Vector3会翻译它。
public Vector3 GetBisector(Vector3 center, Vector3 first, Vector3 second)
{
Vector3 firstDir = (first - center).normalized;
Vector3 secondDir = (second - center).normalized;
Vector3 result = ((firstDir + secondDir) * 0.5f).normalized;
if (IsGreaterThan180(-firstDir, secondDir))
{
// make into a reflex vector
(result.x, result.z) = (-result.x, -result.z);
}
if (result.sqrMagnitude < 0.99f)
{
// we have a colinear set of lines.
// return the perpendicular bisector.
result = Vector3.Cross(Vector3.up, -firstDir).normalized;
}
return result;
}
bool IsGreaterThan180(Vector3 dir, Vector3 dir2)
{
// < 0.0 for clockwise ordering
return (dir2.x * dir.z - dir2.z * dir.x) < 0.0f;
}还要注意，返回的二分线是一个单位长度的向量。使用“中心+平分线*长度”可以将其放入worldspace。
https://stackoverflow.com/questions/3921817
复制相似问题
腾讯云开发者
