void ratateAxis(Eigen::Vector3d& axis, double angle, Eigen::MatrixXd& rm) {
double n1 = axis(0);
double n2 = axis(1);
double n3 = axis(2);
double cosAngle = std::cos(angle);
double sinAngle = std::sin(angle);
rm << n1 * n1 * (1 - cosAngle) + cosAngle, n1* n2* (1 - cosAngle) - n3 * sinAngle, n1* n3* (1 - cosAngle) + n2 * sinAngle,
n1* n2* (1 - cosAngle) + n3 * sinAngle, n2* n2* (1 - cosAngle) + cosAngle, n2* n3* (1 - cosAngle) - n1 * sinAngle,
n1* n3* (1 - cosAngle) - n2 * sinAngle, n2* n3* (1 - cosAngle) + n1 * sinAngle, n3* n3* (1 - cosAngle) + cosAngle;
}
1.axis
是轴方向的单位向量
eg. 若绕x轴旋转 axis
为(1, 0, 0)
的转置
2.rm
为旋转矩阵
3.angle
是旋转的角度