TR=stlread('f-16.stl');
T=TR.ConnectivityList;
X=TR.Points(:,1);
Y=TR.Points(:,2);
Z=TR.Points(:,3);
i=1;
for beta=0:-30:-120
beta=beta/57.3;
if i==2
beta=beta-15/57.3;
end
alpha=0/57.3;
gamma=0/57.3;
Cx=cos(gamma);Sx=sin(gamma);% X轴 滚转
Cy=cos(beta);Sy=sin(beta); % Y轴 俯仰
Cz=cos(alpha);Sz=sin(alpha);% Z轴 偏航
R=[Cy*Cz Cy*Sz -Sy
Sx*Sy*Cz-Cx*Sz Sx*Sy*Sz+Cx*Cz Sx*Cy
Cx*Sy*Cz+Sx*Sz Cx*Sy*Sz-Sx*Cz Cx*Cy];
NP=R*[X';Y';Z'];
x=NP(1,:)'+25*(i-1);
y=NP(2,:)';
z=NP(3,:)';
if i==1
z=z-12;
end
TR = triangulation(T,x,y,z);
h=trisurf(TR);
axis equal;
box;
hold on
i=i+1;
end