wby19851228
2008-11-25, 19:26
我想画单晶铜的面心立方结构，把各个铜原子的坐标算出来了，可是怎么用matlab画出来不知道？如果哪位大侠知道请帮帮我啊，谢谢了！下面是铜原子坐标程序
%建立面心立方单晶铜的模型
%设单晶铜最近邻原子的距离为d=0.256nm,则单晶铜的边长l=2*d*sin(45*pi/180)
d=0.256;
l=2*d*sin(45*pi/180);
for i=1:1
x(1)=0;
y(1)=0;
z(1)=0;
x(i+1)=x(i);
y(i+1)=y(i)+l;
z(i+1)=z(i);
end
for i=3:3
x(3)=l/2;
y(3)=l/2;
z(3)=0;
end
for i=4:4
x(4)=l;
y(4)=0;
z(4)=0;
x(i+1)=x(i);
y(i+1)=y(i)+l;
z(i+1)=z(i);
end
for i=6:6
x(6)=0;
y(6)=l/2;
z(6)=l/2;
end
for i=7:7
x(7)=l/2;
y(7)=0;
z(7)=l/2;
x(i+1)=x(i);
y(i+1)=y(i)+l;
z(i+1)=z(i);
end
for i=9:9
x(9)=l;
y(9)=l/2;
z(9)=l/2;
end
for i=10:10
x(10)=0;
y(10)=0;
z(10)=l;
x(i+1)=x(i);
y(i+1)=y(i)+l;
z(i+1)=z(i);
end
for i=12:12
x(12)=l/2;
y(12)=l/2;
z(12)=l;
end
for i=13:13
x(13)=l;
y(13)=0;
z(13)=l;
x(i+1)=x(i);
y(i+1)=y(i)+l;
z(i+1)=z(i);
end
%建立面心立方单晶铜的模型
%设单晶铜最近邻原子的距离为d=0.256nm,则单晶铜的边长l=2*d*sin(45*pi/180)
d=0.256;
l=2*d*sin(45*pi/180);
for i=1:1
x(1)=0;
y(1)=0;
z(1)=0;
x(i+1)=x(i);
y(i+1)=y(i)+l;
z(i+1)=z(i);
end
for i=3:3
x(3)=l/2;
y(3)=l/2;
z(3)=0;
end
for i=4:4
x(4)=l;
y(4)=0;
z(4)=0;
x(i+1)=x(i);
y(i+1)=y(i)+l;
z(i+1)=z(i);
end
for i=6:6
x(6)=0;
y(6)=l/2;
z(6)=l/2;
end
for i=7:7
x(7)=l/2;
y(7)=0;
z(7)=l/2;
x(i+1)=x(i);
y(i+1)=y(i)+l;
z(i+1)=z(i);
end
for i=9:9
x(9)=l;
y(9)=l/2;
z(9)=l/2;
end
for i=10:10
x(10)=0;
y(10)=0;
z(10)=l;
x(i+1)=x(i);
y(i+1)=y(i)+l;
z(i+1)=z(i);
end
for i=12:12
x(12)=l/2;
y(12)=l/2;
z(12)=l;
end
for i=13:13
x(13)=l;
y(13)=0;
z(13)=l;
x(i+1)=x(i);
y(i+1)=y(i)+l;
z(i+1)=z(i);
end