0001 function direction_cosine = T_ZYX_rot(vector)
0002
0003
0004
0005 sizeV = size(vector);
0006 if (sum(sizeV)==4 && ( ( sizeV(1)==1 && sizeV(2)==3 ) || ( sizeV(2)==1 && sizeV(1)==3 ) ) )
0007 alpha = vector(1); beta = vector(2); gamma = vector(3);
0008
0009 direction_cosine = [ cos(beta)*cos(gamma) cos(beta)*sin(gamma) -sin(beta) ;...
0010 sin(alpha)*sin(beta)*cos(gamma)-cos(alpha)*sin(gamma) sin(alpha)*sin(beta)*sin(gamma)+cos(alpha)*cos(gamma) sin(alpha)*cos(beta);...
0011 cos(alpha)*sin(beta)*cos(gamma)+sin(alpha)*sin(gamma) cos(alpha)*sin(beta)*sin(gamma)-sin(alpha)*cos(gamma) cos(alpha)*cos(beta)];
0012
0013 else
0014 display('the input vector for convertion is:');
0015 display(Vector);
0016 error('Vector for convertion is not a 3*1 nor a 1*3 matrix (see line above)!');
0017 end
0018 end