main.m

用传输矩阵的方法来计算光子晶体带隙

function main
%% This is the main function of TMMfS.
% Copyright. yhwsoft 2005
% Created at 2005.5.11.
% Modifying list:

%% Read the parameters from the input file.
InputFileName='TMMfS.inp';
% Parameters list:
%  N -- Number of layers including scattering layer.
%  r -- Radius of every layer, from the scattering layer to the core layer,
%  where the radius of the scattering layer does not have the substaintial
%  meaning, bu should have a radius larger than that of the first layer.
%  m -- Refraction index of every layer, having the same order as r.
%  MaxOrder -- The highest order of the Bessel function used in calculation. 
%  WavlenRange -- The range of wavelength to be calculated, which is
%  assumed to be the wavelength in the vaccume.
[varmLayers varmFiles r m ...           \\Layers
    MaxOrder varMode WavlenRange ...    \\Control
    outputFile outputData outputPlot ...\\Output
    ] = readTMMfSArgument(InputFileName);

[Qext Qsca Qabs] = TMM(varmLayers,varmFiles,r,m,MaxOrder,WavlenRange);

%% Transform the independent var from wavelength to wave number
outputTMMfS(Qext, Qsca, Qabs, ...
            WavlenRange, varMode, ...
            outputFile, outputData, outputPlot)

%% Summary