t=0:.1:2*pi; t0=118*pi/180; % 設置 t 向量，將角度變弧度
x=11*cos(t+t0); y=-0.6+0.1*t; z=-9+11*sin(t+t0);
plot3(x,z,y), grid, set(gca,'box','off')
set(gca,'xdir','reverse','ydir',reverse') % 常規坐標的 x,y 軸反向
view(-67.5,52) % 旋轉坐標系到一個更好的視角

myworld=vrworld('myvr1.wrl');
open(myworld); view(myworld);
plane=vrnode(myworld,'Plane'); % 獲得飛機對象句柄
set(plane)   % 顯示飛機對象的屬性

for i=1:length(x)
    set(plane,'center',[x(i),y(i),z(i)]); pause(0.05);
end

v(1:3)=[0;1;0]; % 標明按 y 軸進行旋轉
for i=1:length(x)-1
    a1=x(i+1)-x(i); c1=z(i+1)-z(i);
    v(4)=atan2(a1,c1);% 設置旋轉角度
    set(plane,'rotation',v,'center',[x(i),y(i),z(i)]),
    pause(0.05); % 同時設定位置和角度
end