load wind                                           % 裝入數(shù)據(jù)
wind_speed = sqrt(u.^2 + v.^2 + w.^2);                  % 風速計算
hpatch = patch(isosurface(x,y,z,wind_speed,35));        % 等值面繪制
isonormals(x,y,z,wind_speed,hpatch)                   % 計算等值面的頂點范數(shù)  
set(hpatch,'FaceColor','red','EdgeColor','none');          % 設置等值面屬性
 [f vt] = reducepatch(isosurface(x,y,z,wind_speed,45),0.05);  % 減少面片數(shù)量
daspect([1,1,1]);                                        % 繪圖比例
hcone = coneplot(x,y,z,u,v,w,vt(:,1),vt(:,2),vt(:,3),2);          % 流錐圖
set(hcone,'FaceColor','blue','EdgeColor','none');             % 設置流錐對象屬性
camproj perspective                 % 相機透視投影
camva(25)                           % 設置相機
hlight = camlight('headlight');          % 設置相機光照類型
set(hpatch,'AmbientStrength',.1,...      % 設置等值面屬性 
      'SpecularStrength',1,...
      'DiffuseStrength',1);
set(hcone,'SpecularStrength',1);         % 設置流錐體屬性
set(gcf,'Color','k')                     % 設置圖形窗口屬性
lighting gouraud                   % 光照屬性
set(gcf,'Renderer','OpenGL')       % 設置渲染引擎
lighting phong                    % 光照屬性
set(gcf,'Renderer','zbuffer')        % 設置渲染引擎