function [x,y,z] = rotation(x,y,z,az,el); % function [x,y,z] = rotation(x,y,z,az,el) % modifie les coodronnées en x,y et z d'un angle d'élévation et d'azimuth (en degrés) % rotation selon axe y d'angle el (élévation, vertical) a = el/180*pi; roty = [cos(a) 0 -sin(a);0 1 0;sin(a) 0 cos(a)]; % rotation selon axe z d'angle az (azimut, horizontal) a = az/180*pi; rotz = [cos(a) -sin(a) 0;sin(a) cos(a) 0;0 0 1]; for i = 1:size(x,1) M = [x(i,:);y(i,:);z(i,:)]; M = rotz * roty * M; x(i,:) = M(1,:); y(i,:) = M(2,:); z(i,:) = M(3,:); end