體心立方晶胞的matlab仿真_第1頁
體心立方晶胞的matlab仿真_第2頁
體心立方晶胞的matlab仿真_第3頁
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)

文檔簡介

1、精選優(yōu)質(zhì)文檔-傾情為你奉上function main()clc,clear all;single(0,0,0),hold onaxis off;title(體心立方結(jié)構(gòu));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

2、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

3、+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)

4、; 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),Li

5、neWidth,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

6、_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)l

7、ine(-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)

8、,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,2*i);

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論