MATLAB常用代码片段整理分享

N=10;

w=1;

v=4;

state=4;

num=51;

T0=[0,w,w,0;w,0,0,w;w,0,0,w;0,w,w,0];

Tx=zeros(state);

Tx(1,3)=v;

Tx(2,4)=v;

Ty=zeros(state);

Ty(1,2)=v;

Ty(3,4)=v;

H00=kron(eye(N),T0)+kron(diag(ones(1,N-1),1),Tx)+kron(diag(ones(1,N-1),-1),Tx');

H01=kron(eye(N),Ty);

E=zeros(state*N,num);

k=linspace(-1,1,num);

for i=1:num

H=H00+H01*exp(1j*k(i)*pi)+H01'*exp(-1j*k(i)*pi);

eigenvalue=eig(H);

E(:,i)=real(eigenvalue);

end

plot(k,E,'black')

引用自

https://zhuanlan.zhihu.com/p/269595149?share_code=59GL9pv98Sec&utm_psn=2015370794268664118

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删