function HCPlattice
function HCPlattice
a=1;  %½«a¹éÒ»»¯¡£
x0=a/2*[0, 2, 1, -1, -2, -1, 1, 0, -1, 1, 0, 2, 1, -1, -2, -1, 1];
y0=sqrt(3)*a/6*[0, 0, 3, 3, 0, -3, -3, 2, -1, -1, 0, 0, 3, 3, 0, -3, -3];
z0=sqrt(6)*a/3*[0,0,0,0,0,0,0,1,1,1,2,2,2,2,2,2,2];
% x0=[0,a,a/2,-a/2,-a,-a/2,a/2,0,-a/2,a/2,0,a,a/2,-a/2,-a,-a/2,a/2];
%
% y0=[0,0,sqrt(3)*a/2,sqrt(3)*a/2,0,-sqrt(3)*a/2,-sqrt(3)*a/2,sqrt(3)*a/3,-sqrt(3)*a/6,-sqrt(3)*a/6,0,0,sqrt(3)*a/2,sqrt(3)*a/2,0,-sqrt(3)*a/2,-sqrt(3)*a/2];
%
% z0=[0,0,0,0,0,0,0,sqrt(6)*a/3,sqrt(6)*a/3,sqrt(6)*a/3,2*sqrt(6)*a/3,2*sqrt(6)*a/3,2*sqrt(6)*a/3,2*sqrt(6)*a/3,2*sqrt(6)*a/3,2*sqrt(6)*a/3,2*sqrt(6)*a/3];

%(×¢Ò⣬ÔÚmatlabÖÐÔËÐÐʱ£¬Òª½«x0,y0,z0ºóÃæµÄÈý¸ö±í´ïʽÀ­µ½1ÅÅÀï¡£)


%ÏÂÃæÕâ¾ä¿ÉÒÔÓÃÀ´³õ²½²é¿´×ø±êλÖà    

figure, plot3(x0,y0,z0,'o')


%ÏÂÃæ»­Ô­×Óͼ

r=a/2;   %¿É×ÔÐе÷Õû°ë¾¶  ²»Í¬°ë¾¶µÄ»­Í¼Ð§¹û²î±ðÂù´óµÄ

rho=0:0.05*pi:2*pi;

phi=-pi:0.05*pi:pi;

figure

for i=1:17

    x=x0(i)+r*(cos(phi))'*(cos(rho));

    y=y0(i)+r*(cos(phi))'*(sin(rho));

    z=z0(i)+r*(sin(phi))'*(ones(size(rho)));

    surf (x,y,z);

    shading interp   %ʹͼÐι⻬£¬Îޱ߽çÒõÓ°

    hold on

end

axis equal