高对称点扫描
L= 12 ; %% ;
H= 10 ; %% ;
NN= 1;
reso= 0.2;
int1= 5;
int2= 11;
KX= 4*pi/L/3/sqrt(3);
MX =pi/L/sqrt(3);
MY= pi/L/3;
HZ= pi/H;
cinx= 4 ; %% ;
ciny= 3; %% ;
coutx= 1; %% ;
couty= 0.8; %% ;
cout= 1; %% ;
a= L*sqrt(3);
AA=1
for NN=0:0.02:11
AA=AA+1;
kx1=NN*2*pi/3/a;
kx2=(2-NN)*2*pi/3/a;
kx3=0;
kx4=(NN-4)*2*pi/3/a;
kx5=(6-NN)*2*pi/3/a;
kx6=0;
kx7=(NN-8)*2*pi/3/a;
kx8=(10-NN)*2*pi/3/a;
kx9=0;
ky1=NN*2*pi/sqrt(3)/a;
ky2=2*pi/sqrt(3)/a;
ky3=(3-NN)*2*pi/sqrt(3)/a;
ky4=0;
ky5=(NN-4)*2*pi/sqrt(3)/a;
ky6=2*pi/sqrt(3)/a;
ky7=(7-NN)*2*pi/sqrt(3)/a;
ky8=0;
ky9=(NN-8)*2*pi/sqrt(3)/a;
ky10=2*pi/sqrt(3)/a;
ky11=(11-NN)*2*pi/sqrt(3)/a;
kz1=0;
kz2=(NN-3)*pi/2/H;
kz3=pi/2/H;
kz4=(NN-7)*pi/2/H+pi/2/H;
kz5=pi/H;
kx=(NN>=0&NN<1)*kx1+(NN>=1&NN<2)*kx2+(NN>=2&NN<4)*kx3+(NN>=4&NN<5)*kx4+(NN>=5&NN<6)*kx5+...
(NN>=6&NN<8)*kx6+(NN>=8&NN<9)*kx7+(NN>=9&NN<10)*kx8+(NN>=10&NN<11)*kx9;
ky=(NN>=0&NN<1)*ky1+(NN>=1&NN<2)*ky2+(NN>=2&NN<3)*ky3+(NN>=3&NN<4)*ky4+(NN>=4&NN<5)*ky5+...
(NN>=5&NN<6)*ky6+(NN>=6&NN<7)*ky7+(NN>=7&NN<8)*ky8+(NN>=8&NN<9)*ky9+(NN>=9&NN<10)*ky10;
kz=(NN>=0&NN<3)*kz1+(NN>=3&NN<4)*kz2+(NN>=4&NN<7)*kz3+(NN>=7&NN<8)*kz4+(NN>=8&NN<10)*kz5;
plot3(kx,ky,kz,'*'); hold on;
title(NN)
pause(0.01)
end
页:
[1]