rotate acceleration vector using rotation matrix
9 次查看(过去 30 天)
显示 更早的评论
Hello, first time question asker here.
I am trying to orient a smartphone whlist it is in a persons pocket. I collected accelerometer, gyroscope, and magnetometer data from the phone. Next i used the ahrsfilter to output quaternions relative to the NED reference frame.
Next i created a rotation matrix using rotationMatrix = quat2rotm(orientation). this gave me my rotation matrix.
Now here is the problem. My acceleration data which i want to rotate is in a 940x3 (X Y Z) array and my rotation matrix is in a 3x3x940 multidimensional array. So, how can i multiply each of my acceleration points with the appropriate rotation matrix. I suspect i will have to use a loop, however i am unsure about how to do this. How would you go about doing this?
Thank you
0 个评论
采纳的回答
James Tursa
2020-3-17
编辑:James Tursa
2020-3-17
You can use a loop, e.g.
acc = your 940x3 matrix
r = your 3x3x940 array
result = zeros(size(acc));
for k=1:size(acc,1)
result(k,:) = r(:,:,k) * acc(k,:)';
end
Depending on how r is defined (you didn't specify) you might need this instead:
result(k,:) = acc(k,:) * r(:,:,k);
Or you can explore using one of these FEX submissions (some of which require a C/C++ compiler to create mex routines):
2 个评论
James Tursa
2020-3-17
The zeros( ) stuff is siimply pre-allocating the memory for the result before getting into the loop.
更多回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!