Solves the problem of dipole emission close to a sphere multilayer with Mie theory and produces 4 figures to illustrate the results This function can also be called as MulDipFullMonty(stMP,lambda0,nNbTheta) where stMP contains the parameter fields. See MulDipSolveMultiSphere for details of the parameters. The results are returned in 2 structures (for further processing if needed) See MulDipSolveMultiSphere (stMdip) and DipFarFieldAllP (stEfarAllP) for a description of these structures. This file is part of the SPlaC v1.0 package (copyright 2008) Check the README file for further information
0001 % rearrange function parameters 0002 if (nargin<7) 0003 nNmax=stMP.nNmax; 0004 Ca=stMP.Ca; 0005 lambda=stMP.lambda; 0006 Cepsilon=stMP.Cepsilon; 0007 d=stMP.d; 0008 nNmaxESA=stMP.nNmaxESA; 0009 end 0010 0011 % Figure 1 0012 % solve and plot decay rate EFs for perpendicular and parallel dipoles 0013 stMdip=MulDipSolveMultiSphere(nNmax,Ca,lambda,Cepsilon,d,nNmaxESA,'plot','coeff'); 0014 0015 % Figure 2 0016 % Plot of decay rates EFs for other dipole orientation 0017 px=[0 0.25 0.5 0.75 1]; 0018 pz=[1 0.75 0.5 0.25 0]; 0019 DipPlotM(stMdip,px,pz); 0020 0021 % Compute Radiation profiles 0022 stEfarAllP=DipFarFieldAllP(stMdip,nNbTheta); 0023 0024 % Figure 3 0025 % Plot of radiation profiles for selected wavelengths for standard cases 0026 DipPlotEfarAllP(stEfarAllP,lambda0); 0027 0028 % Figure 4 0029 % Example of plot of radiation profiles for selected wavelengths for special cases 0030 DipPlotEfarForP(stEfarAllP,lambda0,-1,1,0,pi/4);