0001
0002 global noCheckSum;
0003 noCheckSum=true;
0004
0005
0006
0007 lambda=transpose(300:1:800);
0008
0009
0010 nNmax=50;
0011
0012
0013
0014 epsilonIn=epsAu(lambda);
0015
0016
0017
0018 scrsz = get(0,'ScreenSize');
0019 figAspectRatio=1;
0020 figHeight=scrsz(4)-150;
0021 figWidth=figAspectRatio*figHeight;
0022 figure('Name','Fig. 6.6', ...
0023 'Position',[(scrsz(3)-figWidth)/2 50 figWidth figHeight]);
0024
0025
0026 epsilonMArray=[1.0, 1.77, 1.0, 1.77, 1.0, 1.77];
0027
0028
0029 aArray=[10, 10, 25, 25, 50, 50];
0030
0031 sLegendCell={'air', 'water', 'air', 'water', 'air', 'water'};
0032
0033
0034 for nn=1:6
0035
0036 epsilonM=epsilonMArray(nn);
0037 a=aArray(nn);
0038
0039
0040 stM=PweSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,'noplot');
0041
0042
0043
0044 axtmp=subplot(3,2,nn);
0045 postmp=get(axtmp,'Position');
0046 delete(axtmp);
0047
0048 ax1=axes('Position',postmp,...
0049 'XAxisLocation','bottom','XColor','k', ...
0050 'YAxisLocation','left','Ycolor','b');
0051
0052 if mod(nn,2)==0
0053 xlim1=[300 800];
0054 set(ax1,'XLim',xlim1,'YLim',[0 8]);
0055 else
0056 xlim1=[300 800];
0057 set(ax1,'XLim',xlim1,'YLim',[0 4.5]);
0058 end
0059
0060 if mod(nn,3)==0
0061 set(get(ax1,'XLabel'),'String','Wavelength [nm]');
0062 end
0063 set(get(ax1,'YLabel'),'String','Q_{ext},Q_{sca},Q_{abs}');
0064
0065 line(lambda,stM.Qext,'Color','b','LineWidth',1.5);
0066 line(lambda,stM.Qsca,'Color',[0 0.5 0]);
0067 line(lambda,stM.Qabs,'Color','k','LineStyle','--');
0068
0069 if nn==1
0070 legend({'Q_{ext}','Q_{sca}','Q_{abs}'});
0071 end
0072 title(['Au in ' sLegendCell{nn} ' - a=' num2str(stM.a) 'nm']);
0073 ax2=axes('Position',get(ax1,'Position'),...
0074 'XAxisLocation','top','XColor','k', ...
0075 'XTickLabel',[], ...
0076 'YAxisLocation','right','Ycolor','r', ...
0077 'Color','none', ...
0078 'YScale','linear', ...
0079 'YLim',[0 1],'XLim',xlim1);
0080 set(get(ax2,'Ylabel'),'String','Q_{abs}/Q_{ext}');
0081
0082 line(lambda,stM.Qabs./stM.Qext,'Color','r','LineStyle',':');
0083 end