clc;clear;
load ener1.data
ener=ener1;
len=length(ener)^(1/3);
len1=51;
max(ener(:,5))
min(ener(:,5))
cc=reshape(ener(:,5),len1,len1,len1);
% image(cc);figure(gcf);
ccccc=0:(len1-1);
[X,Y,Z] = meshgrid(ccccc,ccccc,ccccc);
% slice(X,Y,Z,cc,[25],25,[25])
%%
count=0;
for i=-62:5:175
cla;
p = patch(isosurface(X,Y,Z,cc,i));
isonormals(X,Y,Z,cc,p)
daspect([1 1 1])
isocolors(X,Y,Z,flipdim(cc,2),p)
shading interp
grid off
axis(volumebounds(X,Y,Z,cc))
axis tight
view(3)
camlight
lighting phong
count=count+1;
name=['U=',num2str(i)];
legend(name,'North');
name=[num2str(count),'.jpg'];
saveas(gcf,name)
end
转载本文请联系原作者获取授权,同时请注明本文来自杜波科学网博客。
链接地址:http://blog.sciencenet.cn/blog-464042-753219.html
上一篇:关闭WINDOWS的 watch dog 功能(cuda可以在windows不限时运行)
下一篇:简历模板