matlab 石墨烯代码_石墨烯能带_

x=-pi:pi/100:pi; y=-pi:pi/100:pi; t=pi/3; m=1; [X,Y]=meshgrid(x,y); e1=[0,1]; e2=[-sqrt(3)/2,-0.5]; e3=[sqrt(3)/2,-0.5]; %紧束缚近似中最相近的三个原子的坐标 v1=[sqrt(3),0]; v2=[-sqrt(3)/2,3/2]; v3=[-sqrt(3)/2,-3/2]; Hx=(-cos(X*e1(1,1)+Y*e1(1,2))-cos(X*e2(1,1)+Y*e2(1,2))-cos(X*e3(1,1)+Y*e3(1,2)))*10; Hy=(-sin(X*e1(1,1)+Y*e1(1,2))-sin(X*e2(1,1)+Y*e2(1,2))-sin(X*e3(1,1)+Y*e3(1,2)))*10; Hz=(-sin(X*v1(1,1)+Y*v1(1,2))-sin(X*v2(1,1)+Y*v2(1,2))-sin(X*v3(1,1)+Y*v3(1,2)))*sin(t)+m; H0=(-cos(X*v1(1,1)+Y*v1(1,2))-cos(X*v2(1,1)+Y*v2(1,2))-cos(X*v3(1,1)+Y*v3(1,2)))*cos(t); Z=H0+(Hx.^2+Hy.^2+Hz.^2).^(1/2); c=zeros(size(x)); subplot(1,3,1) surf(X,Y,Z); hold on Z=H0-(Hx.^2+Hy.^2+Hz.^2).^(1/2); surf(X,Y,Z)