linearte_g.m

来自「可以用来求复合半导体能带随着偏压的变化情况」· M 代码 · 共 25 行

function linearTE_g

global ks ind kn Ny Nx E dx dy ind Hr

handle=findobj('Tag','vid');
nr=str2num(get(handle,'String'));

for n=1:(Ny+1)
  Hm(n,:) = Hr(1+(n-1)*(Nx+1):n*(Nx+1),ind(nr))';
end

for n=2:Nx
  for m=2:Ny
    Ey(m,n) = 0.5*(Hm(m,n+1)-Hm(m,n-1))/dx(n);
    Ex(m,n) = 0.5*(Hm(m+1,n)-Hm(m-1,n))/dy(m);
  end
end
Emax = max( max( sqrt( abs(Ex*Ey'))));

figure;  surf( Hm/max(max(abs(Hm))));
title(['Hz, kn = ',num2str( ks( nr))]);
figure;  quiver( Ex/Emax, Ey/Emax);
title(['Et, kn = ',num2str( ks( nr))]);
figure;  quiver( Ey/Emax, Ex/Emax);
title(['Ht, kn = ',num2str( ks( nr))]);