Code covered by the BSD License
-
[KX,KY]=bz_irr2(M)
%%% 2D discretization of irreducible Brillouin zone (IBZ); here, example for triangular lattice;
-
[P, beta]=oblic_eigs(omega,kG...
%%% this function solves the eigenvalue problem beta=f(omega,epsi,kx,ky) for
-
epsi=epsgg(r, na, nb, b1, b2,...
%%% calculation of the 'epsgg' matrix for circular holes using
-
kvect2(kx,ky,b1,b2,N1,N2)
-
pwem2Db.m
-
View all files
Wave-vector diagram for a 2D photonic crystal
by Cazimir-Gabriel Bostan
25 Jan 2009
Calculates and plots the photonic bands at constant frequency for a 2D photonic crystal.
|
Watch this File
|
| File Information |
| Description |
This prog calculates and plots the wave-vector diagram (i.e. photonic bands at constant frequency) for a 2D photonic crystal consisting of cylinders with circular cross-section and
infinite height, arranged in a triangular lattice; oblique propagation is implicit, so the polarization states cannot be separated in E-pol and H-pol; 'omega'is taken as input;
Fourier coefficients for the expansion of dielectric constant are calculated analytically;
the materials considered here are dielectric and dispersionless, but dispersion can be added as a function e(omega)
The package contains the following programs:
pwem2Db.m - main program
epsgg.m - routine for calculating the matrix of Fourier coefficients of dielectric function
bz_irr2.m - routine for 2D discretization of irreducible Brillouin zone polygon;
kvect2.m - routine for calculating diagonal matrices with elements (kx+Gx) and (ky+Gy), where G=(Gx,Gy) is a reciprocal lattice vector
oblic_eigs.m - routine for solving the eigenvalue problem for
H-field |
| MATLAB release |
MATLAB 7.2 (R2006a)
|
|
Tags for This File
|
| Everyone's Tags |
|
| Tags I've Applied |
|
| Add New Tags |
Please login to tag files.
|
|
Contact us at files@mathworks.com
