Solving Complex Matrix Operations

2020 3 4
1 个回答
回答已采纳
更新时间：2020 3 4
3 次查看（30 天）

 采纳的回答

Fabio Freschi
Fabio Freschi 2020-3-4

0 个投票

You can do this in a for loop or using arrayfun
% your A vector
A = rand(1,1001); % dummy vector
% calculation of the matrix B
B = arrayfun(@(a)[cos(-a) sin(-a) 0; -sin(-a) cos(-a) 0; 0 0 1],A,'UniformOutput',false);
Than you can access each matrix using the {} indexing, e.g. the first matrix is
B{1}
etc

5 个评论
显示 3更早的评论 隐藏 3更早的评论

Missael Hernandez
编辑：Missael Hernandez 2020-3-4
u_body=[y(:,1)];
v_body=[y(:,2)];
w_body=[y(:,3)];
phi=[y(:,7)];
theta=[y(:,8)];
psi=[y(:,9)];
Rz = arrayfun(@(a)[cosd(-a) sind(-a) 0; -sind(-a) cosd(-a) 0; 0 0 1],psi,'UniformOutput',false);
Ry = arrayfun(@(a)[cosd(-a) 0 -sind(-a) ; 0 1 0 ; sind(-a) 0 cosd(-a)],theta,'UniformOutput',false);
Rx = arrayfun(@(a)[1 0 0; 0 cosd(-a) sind(-a) ; 0 -sind(-a) cosd(-a)],phi,'UniformOutput',false);
Ok so phi, theta, and psi are three 1001X1 column vectors. I want to take phi_(1,m), theta_(1,m), and psi_(1,m), and plug them into Rx, Ry, and Rz respectively. Then I want to multiply those results to obtain a 3X3 matrix. How do I do that?
Fabio Freschi
Fabio Freschi 2020-3-4
% choose an index
index = 1;
% get the rotation matrix
R = Rx{index}*Ry{index}*Rz{index};
is this what you want?
Missael Hernandez
编辑：Missael Hernandez 2020-3-4
Yes but I want to do it for all 1001 values at the same time
Fabio Freschi
Fabio Freschi 2020-3-4
Add
R = cellfun(@(Rx,Ry,Rz)Rx*Ry*Rz,Rx,Ry,Rz,'UniformOutput',false);
what if I want to multiply each individual R by another vector?

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Operating on Diagonal Matrices 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by