*****************************************************************
[x,y,z]=meshgrid(-2:.02:2,-2:.02:2,-2:.02:2);%建立需要的填充实体,其实是一个正方体
[t,p,r]=cart2sph(x,y,z);%将笛卡尔坐标系换算成球体坐标系
v=NaN*zeros(size(x));%乘上NaN表明实体外我们不考虑
v(r<=1.5)=1;%半径小于1.5的球体
v(r<=1)=2;%半径小于1的球体 O(∩_∩)O哈哈~ 其实是同心圆 不过质量(数值)不一样
xslice=[-2:0.5:2];%X轴切片向量
yslice=[-2:0.5:2];%Y...
zslice=[-2:0.5:2];%Z...
h=slice(x,y,z,v,[],[],zslice);%将球体按需要的方向切
set(h,'LineStyle','none');%填充
view(-60,30);%观看视角
axis equal;%坐标轴等长显示
alpha(0.5);%设置透明度
********************************************************************
自己可以中间加循环 让切片动态显示