finding principal axis and moments of inertia tensor using eig function

Question

Melissa 2013-8-27

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/85683-finding-principal-axis-and-moments-of-inertia-tensor-using-eig-function

采纳的回答： Roger Stafford

Good Morning,

I am trying to calculate the principal moment of inertia as well as the principal axis of inertia.

If I am given a tensor matrix of the rotational moment of inertia:

I=[Jx -Jxy -Jxz; -Jxy Jy -Jyz; -Jxz -Jyz Jz];

And wish to find the principal moments and principal axis of inertia then I would simply use matlab eig.

[Axes,Moments]=eig(I);

And then diagonal of the Moments matrix is the principal moments. Now say I want to find the cosine of the angle of the principal axis of inertia with the three axes X,Y,Z that are center at the CG...how would I go about this?

Is it cos(Axes(:,1)), cos(Axes(:,2)), cos(Axes(:,3)??

Thank you.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Roger Stafford 2013-8-27

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/85683-finding-principal-axis-and-moments-of-inertia-tensor-using-eig-function#answer_95183

"Is it cos(Axes(:,1)), cos(Axes(:,2)), cos(Axes(:,3)??"

No, the components of the eigenvectors themselves, Axes(:,1)), Axes(:,2), Axes(:,3), are already the cosines of the angles between the three principal axes of inertia respectively and the x, y, and z axes, provided the eigenvectors are normalized.

Of course there is an ambiguity about the sense of each of these vectors depending on which of two opposite directions are to be considered as the positive directions of these axes of inertia.