0001
0002
0003
0004
0005 global noCheckSum;
0006 noCheckSum=true;
0007
0008
0009
0010 nNmax=50;
0011
0012 nNmaxESA=50;
0013
0014
0015 a=25;
0016
0017 d=1;
0018
0019
0020 lambda=transpose(300:1:600);
0021
0022
0023 epsilonM=1.0;
0024
0025
0026 scrsz = get(0,'ScreenSize');
0027 figAspectRatio=1;
0028 figHeight=scrsz(4)-150;
0029 figWidth=figAspectRatio*figHeight;
0030 figure('Name','Fig. 4.3', ...
0031 'Position',[(scrsz(3)-figWidth)/2 50 figWidth figHeight]);
0032
0033
0034 for gg=1:2
0035
0036
0037
0038 if gg==1
0039 epsilonIn=2.25;
0040 sTitle='Glass sphere in air';
0041 sYscale='linear';
0042 axScale=[275 625 0.3 2.8];
0043 sLegLoc='East';
0044 else
0045 epsilonIn=epsAg(lambda);
0046 sTitle='Silver sphere in air';
0047 sYscale='log';
0048 axScale=[275 625 0.09 600];
0049 sLegLoc='NorthEast';
0050 end
0051
0052
0053 stMdip=DipSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,d,nNmaxESA,'noplot','coeff');
0054
0055
0056
0057
0058 nNbTheta=361;
0059 stEfarAllP=DipFarFieldAllP(stMdip,nNbTheta);
0060 theta=stEfarAllP.theta;
0061
0062 indNegTheta=length(theta):-1:2;
0063 indAllTheta=[indNegTheta,1:length(theta)];
0064
0065 thetaAll=[-theta(indNegTheta), theta] - pi/2;
0066
0067
0068
0069
0070
0071
0072 subplot(3,2,gg);
0073 plot(lambda,stMdip.MRadPerp,lambda,stMdip.MRadPara);
0074 set(gca,'YScale',sYscale);
0075 axis(axScale);
0076 title(sTitle);
0077 legend({'M_{Rad}^{perp}','M_{Rad}^{//}'},'Location',sLegLoc);
0078 xlabel('Wavelength [nm]');
0079 ylabel('M_{Rad}');
0080
0081
0082
0083 indlambda=find(lambda==360);
0084 subplot(3,2,2+gg);
0085 polar(thetaAll,stEfarAllP.MdRadPerpAllPhi(indlambda,indAllTheta));
0086 ylabel('M^{d-Perp}_{Rad}');
0087 title([sTitle '- \lambda=360nm - M^{Perp}_{Rad}=' ...
0088 num2str(stMdip.MRadPerp(indlambda))]);
0089
0090
0091
0092 subplot(3,2,4+gg);
0093 polar(thetaAll,stEfarAllP.MdRadParaPhi0(indlambda,indAllTheta));
0094 ylabel('M^{d-Para}_{Rad}');
0095 title([sTitle '- \lambda=360nm - M^{Para}_{Rad}=' ...
0096 num2str(stMdip.MRadPara(indlambda))]);
0097
0098 end