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.16', ...
0008 'Position',[(scrsz(3)-figWidth)/2 50 figWidth figHeight]);
0009
0010
0011 h=3;
0012
0013 stEllGobl=EllOblateGeomFact(h,1,'L');
0014 stEllGpro=EllProlateGeomFact(h,1,'L');
0015
0016
0017 lambda=transpose(200:1:1000);
0018 eAg=epsAg(lambda);
0019 eAu=epsAu(lambda);
0020
0021
0022 stEllGsph=EllSphereGeomFact(1,'L');
0023 stEllEMsphAg=EllEMproperties(stEllGsph,1.77,eAg,'L');
0024 stEllEMsphAu=EllEMproperties(stEllGsph,1.0,eAu,'L');
0025
0026
0027 for nn=1:4
0028 if (nn==1 || nn==3)
0029 stEllG=stEllGpro;
0030 sTitle='Prolate spheroid, h=3';
0031 else
0032 stEllG=stEllGobl;
0033 sTitle='Oblate spheroid, h=3';
0034 end
0035 if (nn==1 || nn==2)
0036 epsilonM=1.77;
0037 epsilonIn=eAg;
0038 stEllEMsph=stEllEMsphAg;
0039 sTitle=[sTitle ' - Ag in water'];
0040 else
0041 epsilonM=1.0;
0042 epsilonIn=eAu;
0043 stEllEMsph=stEllEMsphAu;
0044 sTitle=[sTitle ' - Au in air'];
0045 end
0046
0047
0048 stEllEM=EllEMproperties(stEllG,epsilonM,epsilonIn,'L');
0049
0050
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