Home > SPlaC v1_0 > BookFigures > MakeFig6_17.m

MakeFig6_17

PURPOSE ^

Function to reproduce Fig. 6.17

SYNOPSIS ^

function MakeFig6_17()

DESCRIPTION ^

 Function to reproduce Fig. 6.17
 Wavelength dependence of A_1^{perp}, maximum LFIEF on the surface
 (i.e. at the tip) for excitation along the long axis (x-axis)
 for silver or gold prolate (left) or oblate (right) spheroids of
 various aspect ratios.

 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.17', ...
0008     'Position',[(scrsz(3)-figWidth)/2 50 figWidth figHeight]);
0009 
0010 % wavelength in nm (column vector)
0011 lambda=transpose(300:1:1500);
0012 eAg=epsAg(lambda);
0013 eAu=epsAu(lambda);
0014 
0015 % loop 4 times for each subplot
0016 for nn=1:4
0017     if (nn==1 || nn==2) % h=1,3,7
0018         hArray=[1.000001 3 7]; % 1.000001 to model sphere without having to call EllSphereGeomFact
0019     else % h=2,5,10
0020         hArray=[2 5 10];
0021     end
0022     if (nn==1 || nn==3) % prolate
0023         stEllG=EllProlateGeomFact(hArray,1,'L');
0024         sTitle='Prolate spheroids';
0025         sLegLoc='SouthEast';
0026         lambdaMax=1500;
0027     else % oblate
0028         stEllG=EllOblateGeomFact(hArray,1,'L');
0029         sTitle='Oblate spheroids';
0030         sLegLoc='NorthEast';
0031         lambdaMax=1000;
0032     end
0033 
0034     % computes EM properties (A's only)
0035     stEllEMAgWat=EllEMproperties(stEllG,1.77,eAg,'L');
0036     stEllEMAuWat=EllEMproperties(stEllG,1.77,eAu,'L');
0037     stEllEMAuAir=EllEMproperties(stEllG,1.0,eAu,'L');
0038     
0039     % plot results
0040     subplot(2,2,nn);
0041     h=semilogy(lambda,stEllEMAgWat.Aperp1, '-');
0042     set(h,'LineWidth',1.5);
0043     legend({['h=' num2str(hArray(1))], ...
0044         ['h=' num2str(hArray(2))], ...
0045         ['h=' num2str(hArray(3))]},'Location',sLegLoc);
0046     hold on;
0047     h=semilogy(lambda,stEllEMAuWat.Aperp1, '--');
0048     set(h,'LineWidth',1.5);
0049     h=semilogy(lambda,stEllEMAuAir.Aperp1, ':');
0050     set(h,'LineWidth',1.5);
0051     text(350,2.5e5,0,{'Ag/Water (solid)','Au/Water (dashed)','Au/Air (dotted)'});
0052     
0053     axis([300 lambdaMax 10 8e5]);
0054     title(sTitle);
0055     xlabel('Wavelength [nm]');
0056     ylabel('Maximum LFIEF (A_1^{perp})');
0057     set(gca,'YGrid','on','YMinorGrid','off');
0058     hold off;
0059 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)