plotbands.m

一组很好的FDTD计算二维光子带隙的程序

function PlotBands(fname,dk)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Electromagnetic Finite-Difference Time-Domain %
% Version 1.20, Release 1                       %
%                                               %
%   (C) Copyright 2005                          %
%   Sharif University of Technology             %
%   School of Electrical Engineering            %
%   All Rights Reserved                         %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

global Counter Phase NoDiv AbsolutePhase 
global Peaks

Counter=0;

mMax=NoDiv+1;


for m=1:mMax
    MakeFileName(fname)
    RunFFT
    plot(AbsolutePhase,Peaks,'o')
    hold on
    AbsolutePhase=AbsolutePhase+dk/mMax;
    Counter=Counter+1;
    fprintf('.')
end