Matlab scripts examples

These figures are in relation with the computation of a photonic crystal lens using the " Combined Fictitious Sources - Scattering Matrix Method"

Photonic crystal lens. The source location is  x = 0 and y = 5.3 µm. Left: modulus of the complex field. Right: real field.

The method is described in the paper
[1] G. Tayeb, S. Enoch, "Combined Fictitious Sources - Scattering Matrix method", J. Opt. Soc. Am. A 21, No.8, p.1417-1423, 2004. PDF
and used in the paper
[2] T. Decoopman, G. Tayeb, S. Enoch, D. Maystre, B. Gralak, "Photonic crystal lens: from negative refraction and negative index to negative permittivity and permeability", Phys. Rev. Lett. 97 , 073905 (2006). PDF

by using the following files (replace XXX by the string used in the first line of the input file. In this case, XXX = fig11_JOSAA_2004):
cartechamp.XXX.dat :  gives the modulus and the phase of the field for the various positions
XXX.xy.dat :  gives the positions of the centers of the rods (columns 2 and 3)
prof_norm.XXX :  gives the jacket profile (columns 1 and 2)

If you have Matlab, you can automatically compute maps and movies using the scripts map_script.m and movie_script.m.
You will need the file colormap_03 in order to run movie_script.m
Examples:

• Map created with Matlab
• Movie created with Matlab (15 MB, plays fine in Media Player Classic HC. The movie duration is one time period, and if you activate Options > Playback > Repeat forever, il will give something like the movie above). Matlab creates non compressed movies that are huge, and quite difficult to compress without loosing too much quality:
• Same movie recompressed in animated gif with little degradation (3 MB, should play fine in your browser).
• Same movie recompressed in DivX with SUPER (300 kB).
• Same movie recompressed in Quicktime MOV format using MPEG Streamclip (6 MB).