0001
0002
0003 scrsz = get(0,'ScreenSize');
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
0011 lambda=transpose(300:1:1500);
0012 eAg=epsAg(lambda);
0013 eAu=epsAu(lambda);
0014
0015
0016 for nn=1:4
0017 if (nn==1 || nn==2)
0018 hArray=[1.000001 3 7];
0019 else
0020 hArray=[2 5 10];
0021 end
0022 if (nn==1 || nn==3)
0023 stEllG=EllProlateGeomFact(hArray,1,'L');
0024 sTitle='Prolate spheroids';
0025 sLegLoc='SouthEast';
0026 lambdaMax=1500;
0027 else
0028 stEllG=EllOblateGeomFact(hArray,1,'L');
0029 sTitle='Oblate spheroids';
0030 sLegLoc='NorthEast';
0031 lambdaMax=1000;
0032 end
0033
0034
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
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