这个问题在这里也被问到了，但是它是在非编程的背景下提出的。一个建议的解决方案是使用 { y, -x, 0 }。这适用于所有具有x或y分量的向量，但如果向量等于+或-{0, 0, 1}则失败。在这种情况下，我们将得到{0, 0, 0}。
我的当前解决方案（c++）：
// floating point comparison utilizing epsilon
bool is_equal(float, float);
// ...
vec3 v = /* some unit length vector */
// ...
// Set as a non-parallel vector which we will use to find the
// orthogonal vector. Here we choose either the x or y axis.
vec3 orthog;
if( is_equal(v.x, 1.0f) )
orthog.set(1.0f, 0.0f, 0.0f);
else
orthog.set(0.0f, 1.0f, 0.0f);
// Find orthogonal vector
orthog = cross(v, orthog);
orthog.normalize();
这种方法可行，但我觉得可能有更好的方法，而我的搜索结果并没有更多信息。
[编辑]
仅供娱乐，我快速编写了每个建议答案的天真实现（c++），并验证它们都有效（尽管有些不总是自然返回一个单位向量，我在需要时添加了normalize（）调用）。
我的原始想法：
vec3 orthog_peter(vec3 const& v)
{
vec3 arbitrary_non_parallel_vec = v.x != 1.0f ? vec3(1.0, 0.0f, 0.0f) : vec3(0.0f, 1.0f, 0.0f);
vec3 orthog = cross(v, arbitrary_non_parallel_vec);
return normalize( orthog );
}
https://dev59.com/JnjZa4cB1Zd3GeqPe4Eb#19650362
vec3 orthog_robert(vec3 const& v)
{
vec3 orthog;
if(v.x == 0.0f && v.y == 0.0f)
orthog = vec3(1.0f, 1.0f, 0.0f);
else if(v.x == 0.0f)
orthog = vec3(1.0f, v.z / v.y, 1.0f);
else if(v.y == 0.0f)
orthog = vec3(-v.z / v.x, 1.0f, 1.0f);
else
orthog = vec3(-(v.z + v.y) / v.x, 1.0f, 1.0f);
return normalize(orthog);
}
https://dev59.com/JnjZa4cB1Zd3GeqPe4Eb#19651668
// NOTE: u and v variable names are swapped from author's example
vec3 orthog_abhishek(vec3 const& v)
{
vec3 u(1.0f, 0.0f, 0.0f);
float u_dot_v = dot(u, v);
if(abs(u_dot_v) != 1.0f)
return normalize(u + (v * -u_dot_v));
else
return vec3(0.0f, 1.0f, 0.0f);
}
https://dev59.com/JnjZa4cB1Zd3GeqPe4Eb#19658055
vec3 orthog_dmuir(vec3 const& v)
{
float length = hypotf( v.x, hypotf(v.y, v.z));
float dir_scalar = (v.x > 0.0) ? length : -length;
float xt = v.x + dir_scalar;
float dot = -v.y / (dir_scalar * xt);
return vec3(
dot * xt,
1.0f + dot * v.y,
dot * v.z);
};
