plotbands.m

来自「外国人开发的电磁时域有限差分方法工具包 Electromagnetic Fin」· M 代码 · 共 30 行

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