% 《工程線性代數（MATLAB版）》第七章例題7.7程序ea707
%
clear
% 飛行器的外形數據集
Gw=[-4,-3,0;4,-3,0;0,7,0;-4,-3,0]'; % 主翼的頂點坐標
Gt=[0,-3,0;0,-3,3;0,2,0;0,-3,0]';	   % 尾翼的頂點坐標
G=[Gw,Gt]		 		 % 整個飛行器外形的數據集
plot3(Gw(1,:),Gw(2,:),Gw(3,:),'r'),hold on
plot3(Gt(1,:),Gt(2,:),Gt(3,:),'g'),
axis equal,grid on
% 飛行器轉動時的變換矩陣生成
syms u w v
Y=[cos(u),sin(u),0;-sin(u),cos(u),0;0,0,1]
R=[1,0,0;0,cos(w),-sin(w);0,sin(w),cos(w)]
P=[cos(v),0,-sin(v);0,1,0;sin(v),0,cos(v)]
Q=Y*P*R
A=subs(Q,{u,v,w},{pi/4,0,pi/3})
% 轉動后的飛行器外形數據集G1