体心立方晶胞的matlab仿真.docx

function main()clc,clear all;single(0,0,0),hold onaxis off;title(体心立方结构);view(55,20);%spread();%cube_4();function balls(r,x,y,z)X,Y,Z=sphere();X=x+r*X;Y=y+r*Y;Z=z+r*Z;surf(X,Y,Z);colormap(0 0 1);function single(x0,y0,z0)global x y z;x=1 -1 -1 1 1 1 -1 -1+x0;y=1 1 1 1 -1 -1 -1 -1+y0;z=1 1 -1 -1 -1 1 1 -1+z0;for i=1:8 balls(.1,x(i),y(i),z(i),hold on;endcolormap(1 1 1)balls(0.2,x0,y0,z0),hold onfor i=1:7 line(x(i:i+1),y(i:i+1),z(i:i+1),LineWidth,2),hold onendline(-1,1+x0,-1,-1+y0,-1,-1+z0,LineWidth,2),hold online(-1,-1+x0,1,-1+y0,-1,-1+z0,LineWidth,2),hold online(-1,-1+x0,1,-1+y0,1,1+z0,LineWidth,2),hold online(1,1+x0,1,-1+y0,1,1+z0,LineWidth,2),hold online(1,1+x0,1,1+y0,-1,1+z0,LineWidth,2)function spread()single(0,0,0),hold onF(1)=getFrame;i=-1; spread_x(-i); F(2)=getFrame; spread_y(-i); F(3)=getFrame; spread_z(-i); F(4)=getFrame;i=1; spread_x(-i); F(1)=getFrame; spread_y(-i); F(1)=getFrame; spread_z(-i); F(1)=getFrame;movie(F);function spread_x(e)global x y z;x_front=1 1 1 1 0*sign(e)+2*e;y_front=1 1 -1 -1 0;z_front=1 -1 -1 1 0;for j=1:4 balls(.1,x_front(j),y_front(j),z_front(j);endballs(.2,x_front(5),y_front(5),z_front(5);line(x_front(1:4),y_front(1:4),z_front(1:4),LineWidth,2) line(x_front(1),x_front(4),y_front(1),y_front(4),z_front(1),z_front(4),LineWidth,2)line(1*sign(e)+2*e,1,1,1,1,1,LineWidth,2)line(1*sign(e)+2*e,1,1,1,-1,-1,LineWidth,2)line(1*sign(e)+2*e,1,-1,-1,1,1,LineWidth,2)line(1*sign(e)+2*e,1,-1,-1,-1,-1,LineWidth,2)function spread_y(e)global x y z;x_right=x(1:4),0;y_right=y(1:4),0*sign(e)+2*e;z_right=z(1:4),0;for j=1:4 balls(.1,x_right(j),y_right(j),z_right(j);endballs(.2,x_right(5),y_right(5),z_right(5);line(x_right(1:4),y_right(1:4),z_right(1:4),LineWidth,2);line(x_right(1),x_right(4),y_right(1),y_right(4),z_right(1),z_right(4),LineWidth,2)line(-1,-1,1*sign(e)+2*e,1,-1,-1,LineWidth,2)line(1,1,1*sign(e)+2*e,1,-1,-1,LineWidth,2)line(1,1,1*sign(e)+2*e,1,1,1,LineWidth,2)line(-1,-1,1*sign(e)+2*e,1,1,1,LineWidth,2);function spread_z(e)global x y z;x_up=1 1 -1 -1 0;y_up=1 -1 -1 1 0;z_up=1 1 1 1 0*sign(e)+2*e;for j=1:4 balls(.1,x_up(j),y_up(j),z_up(j);endballs(.2,x_up(5),y_up(5),z_up(5); line(x_up(1:4),y_up(1:4),z_up(1:4),LineWidth,2);line(x_up(1),x_up(4),y_up(1),y_up(4),z_up(1),z_up(4),LineWidth,2);line(1,1,1,1,1*sign(e)+2*e,1,LineWidth,2)line(-1,-1,1,1,1*sign(e)+2*e,1,LineWidth,2)line(-1,-1,-1,-1,1*sign(e)+2*e,1,LineWidth,2)line(1,1,-1,-1,1*sign(e)+2*e,1,LineWidth,2)function cube_4()single(0,0,0),hold onF(1)=getFrame;i=1; single(2*i,0,0); F(2)=getFrame; single(0,2*i,0); F(3)=getFrame; single(2*i,2*i,0); F(4)=getFrame; single(0,0