%exm07422_1.m

%以下指令用切面上的色彩表現射流速度
figure(1),clf
[X,Y,Z,V]=flow;	%取4個 的射流數據矩陣，V是射流速度。
x1=min(min(min(X)));x2=max(max(max(X)));	%取x坐標上下限
y1=min(min(min(Y)));y2=max(max(max(Y)));	%取y坐標上下限
z1=min(min(min(Z)));z2=max(max(max(Z)));	%取z坐標上下限
sx=linspace(x1+1.2,x2,5);	%確定5個垂直x軸的切面坐標
sy=0;								%在y=0處，取垂直y軸的切面
sz=0;								%在z=0處，取垂直z軸的切面
slice(X,Y,Z,V,sx,sy,sz);	%畫切片圖
view([-12,30]);shading interp;colormap jet;axis off;colorbar  

disp('按任意鍵，表演動態切片圖！')
pause
clf,
sx=linspace(x1+1.2,x2,20);
nsx=length(sx);
slice(X,Y,Z,V,sx(1),0,0);	%畫切片圖
view([-34,28]);shading interp;colormap jet;axis off;colorbar
title('四維動態切片圖')
for k=1:nsx
   slice(X,Y,Z,V,sx(k),0,0)
   view([-34,28]);shading interp;colormap jet;axis off;colorbar
   title('四維動態切片圖')
   MM(:,k)=getframe;
end
movie(MM,10,40)