squarepixelband_1d.asv

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

%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];
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=5;                    
%Maximal half width of spectrum
NM=2*M+1;

Nk0=30;
Nk=2*Nk0;
%Step number for calculation of k  

gz=(-M:M)*2*pi/cw;
Gz=diag(gz,0);
ang=60;
ang=ang*pi/180;
k1=pi/cw;
kx=linspace(0,0,Nk);
kz=linspace(0,k1,Nk);
kx=cat(2,kx,linspace(0,k1*sin(ang),Nk));
kz=cat(2,kz,linspace(k1,k1,Nk));
ON=5;
[OH,OE]=piecewise_band_1D(FCx,FCz,Gz,ON,kx,kz,k1,Nk);

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');