All All 'Enter Axis:' 'S' 'Enter Angle in Degrees:': Function If
All All 'Enter Axis:' 'S' 'Enter Angle in Degrees:': Function If
clc
clear all
close all
ax=input('Enter Axis :','s');
ang=input('Enter Angle in degrees :');
ang=ang*pi/180.0;
ax
ang
r=rot(ax,ang);
r
function r=rot(ax,ang)
if(ax=='x')
r(1,1)=cos(ang); r(1,2)=(-1)*sin(ang); r(1,3)=0;
r(2,1)=sin(ang); r(2,2)=cos(ang); r(2,3)=0;
r(3,1)=0; r(3,2)=0; r(3,3)=1;
elseif(ax=='y')
r(1,1)=cos(ang); r(1,2)=0; r(1,3)=sin(ang);
r(2,1)=0; r(2,2)=1; r(2,3)=0;
r(3,1)=(-1)*sin(ang); r(3,2)=0; r(3,3)=cos(ang);
else
r(1,1)=1; r(1,2)=0; r(1,3)=0;
r(2,1)=0; r(2,2)=cos(ang); r(2,3)=(-1)*sin(ang);
r(3,1)=0; r(3,2)=sin(ang); r(3,3)=cos(ang);
end
end
using program for question 1 and 2:
clc
clear all
close all
elseif(ax=='y')
r(1,1)=cos(ang); r(1,2)=0; r(1,3)=sin(ang);
r(2,1)=0; r(2,2)=1; r(2,3)=0;
r(3,1)=(-1)*sin(ang); r(3,2)=0; r(3,3)=cos(ang);
else
r(1,1)=1; r(1,2)=0; r(1,3)=0;
r(2,1)=0; r(2,2)=cos(ang); r(2,3)=(-1)*sin(ang);
r(3,1)=0; r(3,2)=sin(ang); r(3,3)=cos(ang);
end
end
Question 6:
clc
close all
clear all
a=input('Enter angle alpha:');
b=input('Enter angle beta:');
c=input('Enter angle gamma:');
px=input('Enter px:');
py=input('Enter py:');
pz=input('Enter pz:');
t=trans(a,b,c,px,py,pz);
fprintf('The transform is :');
t
function t=trans(a,b,c,px,py,pz)
a=a*pi/180;b=b*pi/180;c=c*pi/180;
t(1,1)=cos(a)*cos(b); t(1,2)=cos(a)*sin(b)*sin(c)-sin(a)*cos(c);
t(1,3)=cos(a)*sin(b)*cos(c)+sin(a)*sin(c);t(1,4)=px;
t(2,1)=sin(a)*cos(b); t(2,2)=sin(a)*sin(b)*sin(c)+cos(a)*cos(c);
t(2,3)=sin(a)*sin(b)*cos(c)-cos(a)*sin(c);t(2,4)=py;
t(3,1)=(-1)*sin(b); t(3,2)=cos(b)*sin(c); t(3,3)=cos(b)*cos(c);t(3,4)=pz;
t(4,1)=0;t(4,2)=0;t(4,3)=0;t(4,4)=1;
end