www.gusucode.com > FDTD计算二维光子带隙的程序 > FDTD计算二维光子带隙的程序/247787/InitialField.m
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