Home > SPlaC v1_0 > BookFigures > MakeFig6_16.m

MakeFig6_16

PURPOSE ^

Function to reproduce Fig. 6.16

SYNOPSIS ^

function MakeFig6_16()

DESCRIPTION ^

 Function to reproduce Fig. 6.16
 Wavelength dependence of the maximum and minimum LFIEF
 on a spheroid surface for oblate and prolate (aspect ratio of h=3)
 silver spheroid in water or gold spheroid in air.

 This file is part of the SPlaC v1.0 package (copyright 2008)
 Check the README file for further information

EXAMPLE OF OUTPUT ^

Example figure output

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 
0002 % create, position, and resize figure
0003 scrsz = get(0,'ScreenSize'); % scrsz(4) contains screen height
0004 figAspectRatio=1.3;
0005 figHeight=scrsz(4)-150;
0006 figWidth=figAspectRatio*figHeight;
0007 figure('Name','Fig. 6.16', ...
0008     'Position',[(scrsz(3)-figWidth)/2 50 figWidth figHeight]);
0009 
0010 % aspect ratio
0011 h=3;
0012 % get geometrical factors for oblate and prolate spheroids
0013 stEllGobl=EllOblateGeomFact(h,1,'L');
0014 stEllGpro=EllProlateGeomFact(h,1,'L');
0015 
0016 % wavelength in nm (column vector)
0017 lambda=transpose(200:1:1000);
0018 eAg=epsAg(lambda);
0019 eAu=epsAu(lambda);
0020 
0021 % sphere case:
0022 stEllGsph=EllSphereGeomFact(1,'L');
0023 stEllEMsphAg=EllEMproperties(stEllGsph,1.77,eAg,'L');
0024 stEllEMsphAu=EllEMproperties(stEllGsph,1.0,eAu,'L');
0025 
0026 % loop 4 times for each subplot
0027 for nn=1:4
0028     if (nn==1 || nn==3) % prolate
0029         stEllG=stEllGpro;
0030         sTitle='Prolate spheroid, h=3';
0031     else % oblate
0032         stEllG=stEllGobl;
0033         sTitle='Oblate spheroid, h=3';
0034     end
0035     if (nn==1 || nn==2) % silver in water
0036         epsilonM=1.77;
0037         epsilonIn=eAg;
0038         stEllEMsph=stEllEMsphAg;
0039         sTitle=[sTitle ' - Ag in water'];
0040     else % gold in air
0041         epsilonM=1.0;
0042         epsilonIn=eAu;
0043         stEllEMsph=stEllEMsphAu;
0044         sTitle=[sTitle ' - Au in air'];
0045     end
0046 
0047     % computes EM properties (A's only)
0048     stEllEM=EllEMproperties(stEllG,epsilonM,epsilonIn,'L');
0049     
0050     % plot results
0051     subplot(2,2,nn);
0052     h=semilogy(lambda,stEllEM.Aperp1, '-', ...
0053         lambda,stEllEM.Aperp3, '-', ...
0054         lambda,stEllEMsph.Aperp1, '-');
0055     set(h,'LineWidth',1.5);
0056     line(lambda,stEllEM.Apara1,'LineStyle',':','LineWidth',1.5,'Color',get(h(1),'Color'));
0057     line(lambda,stEllEM.Apara3,'LineStyle',':','LineWidth',1.5,'Color',get(h(2),'Color'));
0058     line(lambda,stEllEMsph.Apara1,'LineStyle',':','LineWidth',1.5,'Color',get(h(3),'Color'));
0059     
0060     axis([200 1000 0.1 5e4]);
0061     title(sTitle);
0062     xlabel('Aspect ratio, h');
0063     ylabel('LFIEF');
0064     legend({'A_1^{perp}', 'A_3^{perp}', 'A_{sphere}^{perp}', ...
0065         'A_1^{//}'},'Location','NorthEast');
0066 end

This web page is part of the SPlaC package © 2008. Contact: Eric Le Ru
Generated on Wed 03-Dec-2008 11:10:14 by m2html © 2003 (adapted)