initialfield.m

来自「fdtd计算二维光子晶体带隙」· M 代码 · 共 36 行

function InitialField

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 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 Alpha Power PowerPlanes
global xCenter yCenter Source 
global U X Y Tx Ty Ux Uy

Power=zeros(length(PowerPlanes),1);

Tx=zeros(X,Y);
Ty=zeros(X,Y);
Ux=zeros(X,Y);
Uy=zeros(X,Y);

if Source>0
    U=zeros(X,Y);
else
    xCoordinate=(1:X).'*ones([1,Y]);
    yCoordinate=ones([X,1])*(1:Y);

    x=xCoordinate-xCenter;
    y=yCoordinate-yCenter;
    r=sqrt(x.^2+y.^2);

    U=exp(-Alpha*r.^2);
end