基于MATLAB实现的RCWA计算方案的MATLAB代码+使用说明文档

classdef RCWA < handle properties Name Polarization LightSource Period ModeNumber IndexProfileResolution LayerNumber Layers FieldDetectorNumber FieldDetectors ComputingResults Reflection Transmission end methods function mr=RCWA(varargin) mr.Name='rcwa'; mr.ModeNumber=50; mr.Period=800; mr.ComputingResults=Results; for icount=1:nargin if ischar(varargin{icount}) switch lower(varargin{icount}) case {'nm','name'} mr.Name=varargin{icount+1}; case {'period','lambda'} mr.Period=varargin{icount+1}; case {'nmode','modenumber'} mr.ModeNumber=varargin{icount+1}; otherwise disp('Error: Unknown input.'); return; end end end mr.IndexProfileResolution=2^nextpow2(4*mr.ModeNumber+1); mr.FieldDetectorNumber=0; mr.LayerNumber=2; mr.Layers=cell(mr.LayerNumber,1); mr.Layers{1}=Layer('period',mr.Period,... 'indresol',mr.IndexProfileResolution); mr.Layers{1}.SetThickness(0); mr.Layers{1}.SetLayerThickness(0); mr.Layers{2}=Layer('period',mr.Period,... 'indresol',mr.IndexProfileResolution); mr.Layers{2}.SetThickness(0); mr.Layers{2}.SetLayerThickness(0); end function SetModeNumber(mr,nmode) mr.ModeNumber=nmode; mr.IndexProfileResolution=2^nextpow2(4*mr.ModeNumber+1); for icount=1:mr.LayerNumber mr.Layers{icount}.SetIndexProfileResolution(... mr.IndexProfileResolution); end end function SetPeriod(mr,period) mr.Period=period; for icount=1:mr.LayerNumber mr.Layers{icount}.SetPeriod(... mr.Period); end end function AddLayer(mr,varargin) mr.Layers{mr.LayerNumber+1}=mr.Layers{mr.LayerNumber}; mr.Layers{mr.LayerNumber}=Layer('period',mr.Period,... 'IndexProfileResolution',mr.IndexProfileResolution,... 'lyt',0); if nargin>1 for icount=1:2:nargin-1 if ischar(varargin{icount}) switch lower(varargin{icount}) case {'layerthick','layerthickness','lyt'} mr.Layers{mr.LayerNumber}... .SetLayerThickness(... varargin{icount+1}); case {'nm','name'} mr.Layers{mr.LayerNumber}.SetName(... varargin{icount+1}); case {'index','profile',... 'indexprofile',... 'profileparameters',... 'indexprofileparameters'} mr.Layers{mr.LayerNumber}... .SetIndexProfile(... varargin{icount+1}); otherwise disp('Error: Unknown input.'); return; end end end end mr.Layers{mr.LayerNumber}.Thickness... =mr.Layers{mr.LayerNumber+1}.Thickness... +mr.Layers{mr.LayerNumber}.LayerThickness; mr.Layers{mr.LayerNumber+1}.Thickness... =mr.Layers{mr.LayerNumber+1}.Thickness... +mr.Layers{mr.LayerNumber}.LayerThickness; mr.LayerNumber=mr.LayerNumber+1; end function SetLayer(mr,varargin) if nargin==1 disp('Error: Please specify the layer by index or name.'); elseif mod(nargin,2) disp('Error: Input argument number incorrect.'); else if ischar(varargin{1}) for icount=1:mr.LayerNumber if strcmpi(mr.Layers{icount}.Name,varargin{1}) for jcount=2:2:nargin-1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% switch lower(varargin{jcount}) case {'nm','name'} mr.Layers{icount}.Name=varargin{jcount+1}; case {'layerthick','layerthickness','lyt'} delthick=varargin{jcount+1}-mr.Layers{icount}.LayerThickness; mr.Layers{icount}.LayerThickness=varargin{jcount+1}; for kcount=icount+1:length(mr.Layers) mr.Layers{kcount}.Thickness... =mr.Layers{kcount}.Thickness+delthick; end case {'index','profile',... 'indexprofile',... 'profileparameters',... 'indexprofileparameters'} mr.Layers{icount}.SetIndexProfile(varargin{jcount+1}); otherwise disp('Error: Unknown input.'); return; end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% end end end else icount=varargin{1}; for jcount=2:2:nargin-1 switch lower(varargin{jcount}) case {'nm','name'} mr.Layers{icount}.Name=varargin{jcount+1}; case {'layerthick','layerthickness','lyt'} delthick=varargin{jcount+1}-... mr.Layers{icount}.LayerThickness; mr.Layers{icount}.LayerThickness=... varargin{jcount+1}; for kcount=icount+1:length(mr.Layers) mr.Layers{kcount}.Thickness... =mr.Layers{kcount}.Thickness... +delthick; end case {'index','profile',... 'indexprofile',... 'profileparameters',... 'indexprofileparameters'} mr.Layers{icount}... .SetIndexProfile(varargin{jcount+1}); otherwise disp('Error: Unknown input.'); return; end end end end end function ApplyLightSource(mr,varargin) mr.LightSource=Source; mr.Polarization=mr.LightSource.Polarization; mr.SetLightSource(varargin{:}); end function SetLightSource(mr,varargin) if nargin==1 return; else for icount=1:2:nargin-1 switch lower(varargin{icount}) case {'pol','polarization'} mr.LightSource... .SetPolarization(varargin{icount+1}); mr.Polarization=upper(varargin{icount+1}); case {'wv','wavelength'} mr.LightSource... .SetWavelength(varargin{icount+1});