0001
0002
0003
0004
0005
0006 lambda=633;
0007
0008
0009 aideg=34.5:0.005:37;
0010
0011
0012 Lgap=850;
0013
0014
0015 eAg=epsAg(lambda);
0016
0017 nP=1.766;
0018 nM=1.0;
0019
0020
0021
0022 nNbSurf=2;
0023 Cepsilon{1}=nP^2+0*lambda;
0024 Cepsilon{2}=nM^2+0*lambda;
0025 Cepsilon{3}=eAg;
0026 CL{1}=0;
0027 CL{2}=Lgap;
0028
0029
0030 stResTM=PlnMultiRef ('TM',nNbSurf, lambda, Cepsilon, CL, aideg);
0031
0032
0033
0034
0035 Rp=abs(stResTM.rP).^2;
0036
0037
0038 MLoc=stResTM.MoutPerp{2}+stResTM.MoutPara{2};
0039
0040
0041 MLocEq3_27=2*abs(eAg)/(nM^2 * imag(eAg)) * imag(eAg*sqrt(1/(eAg+nM^2)))*nP*cos(aideg*pi/180).*(1-Rp);
0042
0043
0044
0045
0046
0047
0048 scrsz = get(0,'ScreenSize');
0049 figAspectRatio=1.5;
0050 figWidth=3/4*scrsz(3);
0051 figure('Name','Fig. 3.15(a)', ...
0052 'Position',[(scrsz(3)-figWidth)/2 scrsz(4)-150-figWidth/figAspectRatio figWidth figWidth/figAspectRatio]);
0053
0054 line(aideg,Rp,'Color','b');
0055 ax1 = gca;
0056 set(ax1,'XColor','k','YColor','b')
0057 set(ax1,'XLimMode','Manual','YLimMode','Manual');
0058 set(ax1,'XLim',[34.5 37],'YLim',[0 1.15]);
0059 set(get(ax1,'XLabel'),'String','Angle of incidence [degrees]');
0060 set(get(ax1,'YLabel'),'String','Reflectivity');
0061 legend('R','Location','NorthWest');
0062 title('Air/Ag PSPP excitation in the Otto configuration - L_{gap}=850 nm, \lambda=633 nm');
0063
0064 ax2=axes('Position',get(ax1,'Position'),...
0065 'XAxisLocation','top','XColor','k', ...
0066 'YAxisLocation','right','Ycolor','r', ...
0067 'Color','none', ...
0068 'YScale','Log', ...
0069 'XTickLabel',{}, ...
0070 'XLimMode','Manual','YLimMode','Manual', ...
0071 'XLim',[34.5 37],'YLim',[0.2 800]);
0072 set(get(ax2,'Ylabel'),'String','Local field intensity EF');
0073
0074 line(aideg,MLoc,'Color','r');
0075 indMarker=10:20:length(aideg);
0076 line(aideg(indMarker),MLocEq3_27(indMarker),'Color','r','Marker','o','LineStyle','none');
0077 legend({'M_{Loc}','M_{Loc} from Eq. 3.27'},'Location','SouthEast','Color','w');
0078
0079