squarepixelband_1d.m

通过计算晶体结构来确定其光谱分布

%SquarePixelBand_1D.m
%Light propagates in the plane of x-z.
%Only periodicity of premittivity is at z-axis.
clear;
clear global;
close all;

pw=10e-9;
%Pixel width      

eb=1;
ea=11.4;
ebxn=[eb,eb,eb,ea,ea,eb,eb,ea];
t=size(ebxn);
cw=t(2)*2*pw;
ebxn0=cat(2,fliplr(ebxn),ebxn);
hN=0;
hN=hN+1;
figure(hN);
plot(ebxn0);
title('Original ebxn');


M=12;                    
%Maximal half width of spectrum
NM=2*M+1;

Nk0=60;
Nk=3*Nk0;
%Step number for calculation of k  

gz=(-M:M)*2*pi/cw;
Gz=diag(gz,0);
ang=80;
ang=ang*pi/180;
k1=pi/cw;
k0=2*pi/cw;
kx=linspace(0,0,Nk0);
kz=linspace(0,k1,Nk0);
kx=cat(2,kx,linspace(0,k1*sin(ang),Nk0));
kz=cat(2,kz,linspace(k1,k1,Nk0));
kx=cat(2,kx,linspace(k1*sin(ang),0,Nk0));
kz=cat(2,kz,linspace(k1,0,Nk0));

ON=8;
[OH,OE]=piecewise_band_1D(ebxn,pw,M,ON,Nk,kx,kz,Gz,k0);

hN=hN+1;
figure(hN);
plot(OH);
xlabel('k');
ylabel('omega');
title('H-polarization');
hN=hN+1;
figure(hN);
plot(OE);
xlabel('k');
ylabel('omega');
title('E-polarization');
disp('SquarePixelBand ends');