Solves the problem of dipole emission close to a single sphere with Mie theory and produces 4 figures to illustrate the results This function can also be called as DipFullMonty(stMP,lambda0,nNbTheta) where stMP contains the parameter fields. See DipSolveSingleSphere for details of the parameters. The results are returned in 2 structures (for further processing if needed) See DipSolveSingleSphere (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<8) 0003 nNmax=stMP.nNmax; 0004 a=stMP.a; 0005 lambda=stMP.lambda; 0006 epsilonM=stMP.epsilonM; 0007 epsilonIn=stMP.epsilonIn; 0008 d=stMP.d; 0009 nNmaxESA=stMP.nNmaxESA; 0010 end 0011 0012 % Figure 1 0013 % solve and plot decay rate EFs for perpendicular and parallel dipoles 0014 stMdip=DipSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,d,nNmaxESA,'plot','coeff'); 0015 0016 % Figure 2 0017 % Plot of decay rates EFs for other dipole orientation 0018 px=[0 0.25 0.5 0.75 1]; 0019 pz=[1 0.75 0.5 0.25 0]; 0020 DipPlotM(stMdip,px,pz); 0021 0022 % Compute Radiation profiles 0023 stEfarAllP=DipFarFieldAllP(stMdip,nNbTheta); 0024 0025 % Figure 3 0026 % Plot of radiation profiles for selected wavelengths for standard cases 0027 DipPlotEfarAllP(stEfarAllP,lambda0); 0028 0029 % Figure 4 0030 % Example of plot of radiation profiles for selected wavelengths for special cases 0031 DipPlotEfarForP(stEfarAllP,lambda0,-1,1,0,pi/4);