Iteratively concatenate matrices held in a structure with column vectors stored in an 3xn matrix

조회 수: 1 (최근 30일)
I have a series of 3x3 rotation matrices calculated from a 3xn matrix where columns are euler rotations. These are stored in a structure (euler) provided by:
%Calculates rotation matrices for each camera position, stores them
%inside a structure (euler) and create dynamics fieldnames
%(euler.PV1,..,PVn)
euler=struct();
for k=1:size(MPV,2)
Rxk = [1 0 0; 0 cosd(MPV(1,k)) -sind(MPV(1,k)); 0 sind(MPV(1,k)) cosd(MPV(1,k))];
Ryk = [cosd(MPV(2,k)) 0 sind(MPV(2,k)); 0 1 0; -sind(MPV(2,k)) 0 cosd(MPV(2,k))];
Rzk = [cosd(MPV(3,k)) -sind(MPV(3,k)) 0; sind(MPV(3,k)) cosd(MPV(3,k)) 0; 0 0 1];
RMk = Rzk*Ryk*Rxk;
euler.(['PV' num2str(k)])=RMk;
end
The resultant structure is:
PV1: [3x3 double]
PV2: [3x3 double]
PV3: [3x3 double]
PV4: [3x3 double]
PV5: [3x3 double]
PV6: [3x3 double]
PV7: [3x3 double]
Where individual matrices have the form: >> euler.PV1
ans =
0.9994 -0.0251 0.0256
0.0254 0.9996 -0.0082
-0.0254 0.0088 0.9996
I also have a 3xn matrix (MPV where columns contain x, y, z components of column vectors:
MPV =
-0.9741 -1.7308 -1.6931 -1.7115 ...
0.9437 -0.0551 -0.0210 -0.0079 ...
0.9437 -0.0551 -0.0210 -0.0079 ...
My problem is that I want to iteratively concatenate the 3x3 matrices stored in the euler structure with the x, y, z components stored in MPV to give a matrix of the form:
r11 r12 r13 x
r21 r22 r23 y
r31 r32 r33 z
I thought that a for loop of the form:
for j=1:size(euler.PV)
MM = vertcat(PV(j),MPV(i));
Mat.(['M' num2str(j)])=MMj;
end
might do the trick, but gives the error:
Undefined function or variable 'M'.
Error in matcon (line 42)
for j=1:length(M)
Note that I need to perform a further horizontal concatenation for each matrix and a single row array of the form ID = [0 0 0 1](i.e. horzcat(mat.M1,..Mn,ID) to create 4 x 4 matrices (homogenous coordinates), so having an output as a structure may be desirable at this stage. Apologies in advance for this noob-like question. Any help with the code would be much appreciated.

채택된 답변

Andrei Bobrov
Andrei Bobrov 2012년 10월 11일
편집: Andrei Bobrov 2012년 10월 11일
outcl = cellfun(@(x,y)[x,y],struct2cell(euler),num2cell(MPV,1)','un',0);
eulernew = cell2struct(outcl,fields(euler),1);
  댓글 수: 1
Thomas
Thomas 2012년 10월 11일
Wow! fantastic. Code works beautifully. It may take me a while to work out what's going on here. Thanks again Andrei. Thomas

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Creating and Concatenating Matrices에 대해 자세히 알아보기

제품

Community Treasure Hunt

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

Start Hunting!

Translated by