huitucycj.m

来自「包括MATLAB7.0实用指南（1）书中所有的例子的源代码。」· M 代码 · 共 18 行

load wind
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)
set(hpatch,'facecolor','r','edgecolor','none');
[f,vt]=reducepatch(isosurface(x,y,z,wind_speed,45),0.05);
daspect([1,1,1]);
hcone=coneplot(x,y,z,u,v,w,vt(:,1),vt(:,2),vt(:,3),2);
set(hcone,'facecolor','b','edgecolor','none');
camproj perspective
camva(5);
%hlight=camlight('headlight');
set(hpatch,'ambientstrength',1...
'spectularstrength',1...
'diffusestrength',1);