资源描述：
《体心立方晶胞的matlab仿真》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、functionmain()clc,clearall;single(0,0,0),holdonaxisoff;title('体心立方结构');view(55,20);%spread();%cube_4();functionballs(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([001]);functionsingle(x0,y0,z0)globalxyz;x=[1-1-1111-1-1]+x0;y=[1111-1-1-1-1]+y0;z=[11-1-1-111-1]+z0;fori
2、=1:8balls(.1,x(i),y(i),z(i)),holdon;endcolormap([111])balls(0.2,x0,y0,z0),holdonfori=1:7line(x(i:i+1),y(i:i+1),z(i:i+1),'LineWidth',2),holdonendline([-1,1]+x0,[-1,-1]+y0,[-1,-1]+z0,'LineWidth',2),holdonline([-1,-1]+x0,[1,-1]+y0,[-1,-1]+z0,'LineWidth',2),holdonline([-1,-1]+x0,[1,-1]+y0,[1,1]+
3、z0,'LineWidth',2),holdonline([1,1]+x0,[1,-1]+y0,[1,1]+z0,'LineWidth',2),holdonline([1,1]+x0,[1,1]+y0,[-1,1]+z0,'LineWidth',2)functionspread()single(0,0,0),holdonF(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;
4、spread_y(-i);F(1)=getFrame;spread_z(-i);F(1)=getFrame;movie(F);functionspread_x(e)globalxyz;x_front=[11110]*sign(e)+2*e;y_front=[11-1-10];z_front=[1-1-110];forj=1:4balls(.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(
5、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],
6、[-1,-1],'LineWidth',2)functionspread_y(e)globalxyz;x_right=[x(1:4),0];y_right=[y(1:4),0]*sign(e)+2*e;z_right=[z(1:4),0];forj=1:4balls(.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)
7、,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
