| 查看: 376 | 回复: 1 | |||
[交流]
【求助】请教下高手，是程序有错误，还是思路就不对呢？ 已有1人参与
|
|
昨天想起之前考虑过标准正交基和旋转变换矩阵的问题，所以编了一个小程序，如下： BRP = eul2r(pi/6,pi/4,pi/3); L1P = [-0.866;0;-0.5]; L2P = [0;0;-1]; L1B = BRP*L1P; L2B = BRP*L2P; B = Ortho_basis(L1B,L2B); P = Ortho_basis(L1P,L2P); R = B1*inv(P); function A = Ortho_basis(a1,a2) beta1 = a1; beta2 = a2-(dot(a2,beta1)/dot(beta1,beta1))*beta1; a1 = beta1/norm(beta1); a2 = beta2/norm(beta2); a3 = cross(a1,a2); A = [a1';a2';a3']; end 本意是想给定一个欧拉旋转矩阵BRP和P坐标系中两个单位向量L1P和L2P 将这两个单位向量经BRP旋转变换后得到在B坐标系中相应的单位向量L1B=BRP*L1P和L2B=BRP*L2P， 然后想通过L1P与L1B，L2P与L2B，这两对向量之间的关系反求出其欧拉旋转变换矩阵R， 反求思路是由L1P和L2P得出P坐标系中的标准正交基P，由L1B和L2B得出B坐标系中的标准正交基B，然后由R*P=B R=B*inv(P)算出R 如果这个思路没问题的话，应该是R=BRP的吧，但是为什么我得出的R不等于BRP呢？ 请教下高手，是我的程序有错误，还是利用标准正交基来求欧拉旋转变换矩阵的思路就不对呢？ |
回复此楼» 猜你喜欢
288求调剂，一志愿华南理工大学071005
已经有5人回复
-
一志愿中国海洋大学，生物学，301分，求调剂
已经有5人回复
-
294求调剂材料与化工专硕
已经有12人回复
-
0856调剂，是学校就去
已经有5人回复
-
0817 化学工程 299分求调剂 有科研经历 有二区文章
已经有13人回复
-
梁成伟老师课题组欢迎你的加入
已经有11人回复
-
复试调剂
已经有6人回复
-
0703化学调剂 ，六级已过，有科研经历
已经有11人回复
-
321求调剂
已经有10人回复
-
东华理工大学化材专业26届硕士博士申请
已经有8人回复
2楼2011-02-09 08:53:57
