%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%ljNo
%%%%%%%%%%%%%%%%%%开始计算M@ZI\
forT=1:
TimeTL_s:
l9!
r
%%%保存前一时间的电磁场#o2[hibq
Ex_Pre=Ex;o$.fhD
Hy_Pre=Hy;bYPKh
%%%%中间差分计算Dx8sCv]|cn
fori=2:
KEO|hp
XkV
Dx（i）=Dx（i）-（dt/zdelta）*（Hy（i）-Hy（i-1））;cFWc<55aX6
endI`p;F!
s
%%%%%%%%加入源!
Rt>xD
Dx（kc）=cos（2*pi*f*T*dt）*exp（-4*pi*（（T*dt-t0）/d）^2）;1!
gbTeVlY
>dG[G>
%%%计算电场Exw+{LAS
fori=1:
kpstart-109Cez\0
Ex（i）=Dx（i）/epsz;tNX
|U:
Y*
endDDH:
）=;z
fori=kpstop+1:
KEU`m54f@U
Ex（i）=Dx（i）/epsz;_f:
W$\ho
endJ9[r|`gJ（
Dx3=Dx2;C6AUNRpl
Dx2=Dx1;w*JGUk
Dx1=Dx;>"=>3
fori=kpstart:
kpstopFG*r'tC~r
Ex（i）=（1/b0）*（a0*Dx1（i）+a1*Dx2（i）+a2*Dx3（i）-b1*Ex1（i）-b2*Ex2（i））;）TH@#1
end`^Em&6!
!
Ex2=Ex1;'CkIz"Wd
Ex1=Ex;vOpKNp
Sx3=Sx2;kq,ucU%>p
Sx2=Sx1;}d}Ke_Q0
Ex
（1）=ex_low_m2;^
ex_low_m2=ex_low_m1;4u5-7[TZ
ex_low_m1=Ex
（2）;HqT#$}rv
DG:
Z=LuJr
Ex（KE）=ex_high_m2;76h,]xi
ex_high_m2=ex_high_m1;SmSH2m-
ex_high_m1=Ex（KE-1）;X=fYWj[H,
%%%%%%%%%%%%%%%%%%计算磁场O*
）Vhw'pK
fori=1:
KE-1XBu"-（
Hy（i）=Hy（i）-（dt/（mu*zdelta））*（Ex（i+1）-Ex（i））;GMf`A,>
end*kDCliL
plot（Ex）;）g#T9tx2D
gridon;CxOob1@
pause;:
hk5.
[
end
%FDTDMainFunctionJobstoWorkers%
%***********************************************************************
%3-DFDTDcodewithPECboundaries
%***********************************************************************
%
%ThisMATLABM-fileimplementsthefinite-differencetime-domain
%solutionofMaxwell'scurlequationsoverathree-dimensional
%Cartesianspacelatticecomprisedofuniformcubicgridcells.
%
%Toillustratethealgorithm,anair-filledrectangularcavity
%resonatorismodeled.Thelength,width,andheightofthe
%cavityareXcm（x-direction）,Ycm（y-direction）,and
%Zcm（z-direction）,respectively.
%
%ThecomputationaldomainistruncatedusingPECboundary
%conditions:
%ex（i,j,k）=0onthej=1,j=jb,k=1,andk=kbplanes
%ey（i,j,k）=0onthei=1,i=ib,k=1,andk=kbplanes
%ez（i,j,k）=0onthei=1,i=ib,j=1,andj=jbplanes
%ThesePECboundariesformtheouterlosslesswallsofthecavity.
%
%Thecavityisexcitedbyanadditivecurrentsourceoriented
%alongthez-direction.Thesourcewaveformisadifferentiated
%Gaussianpulsegivenby
%J（t）=-J0*（t-t0）*exp（-（t-t0）^2/tau^2）,
%wheretau=50ps.TheFWHMspectralbandwidthofthiszero-dc-
%contentpulseisapproximately7GHz.Thegridresolution
%（dx=2mm）waschosentoprovideatleast10samplesper
%wavelengthupthrough15GHz.
%
%ToexecutethisM-file,type"fdtd3D"attheMATLABprompt.
%ThisM-filedisplaystheFDTD-computedEzfieldsateveryother
%timestep,andrecordsthoseframesinamoviematrix,M,which
%isplayedattheendofthesimulationusingthe"movie"command.
%
%***********************************************************************
function[Ex,Ey,Ez]=FDTD3D_Main（handles）
globalSimRunStop
%if~isdir（'C:
\MATLAB7\work\cavity\figures'）
%mkdir'C:
\MATLAB7\work\cavity\figures'
%end
%***********************************************************************
%GridPartition
%***********************************************************************
=get,'Value'）;
=get,'Value'）;
=get,'Value'）;
%***********************************************************************
%GridDimensons
%***********************************************************************
ie=get,'Value'）;%numberofgridcellsinx-direction
je=get,'Value'）;%numberofgridcellsiny-direction
ke=get,'Value'）;%numberofgridcellsinz-direction
ib=ie+1;
jb=je+1;
kb=ke+1;
%***********************************************************************
%AllDomainsFieldsIni.
%***********************************************************************
Ex=zeros（ie,jb,kb）;
Ey=zeros（ib,je,kb）;
Ez=zeros（ib,jb,ke）;
Hx=zeros（ib,je,ke）;
Hy=zeros（ie,jb,ke）;
Hz=zeros（ie,je,kb）;
%***********************************************************************
%Fundamentalconstants
%***********************************************************************
=;%speedoflightinfreespace
=*pi*;%permeabilityoffreespace
=**;%permittivityoffreespace
%***********************************************************************
%Gridparameters
%***********************************************************************
=get,'Value'）;%locationofz-directedcurrentsource
=get,'Value'）;%locationofz-directedcurrentsource
=floor（ke/2）;%Surfaceofobservation
=get,'Value'）;;%spaceincrementofcubiclattice
=*;%timestep
=get,'Value'）;;%totalnumberoftimesteps
%***********************************************************************
%DifferentiatedGaussianpulseexcitation
%***********************************************************************
=get,'Value'）*100e-12;
=;
=3*;
=get,'Value'）*10e11;
=*;
%***********************************************************************
%Materialparameters
%***********************************************************************
=get,'Value'）;
=get,'Value'）;
%***********************************************************************
%Updatingcoefficients
%***********************************************************************
=*/**）/+*/**）;
=/+*/**）;
=;
=;
%***********************************************************************
%CallingFDTDAlgorithm
%***********************************************************************
ex=zeros（ib,jb,kb）;
ey=zeros（ib,jb,kb）;
ez=zeros（ib,jb,kb）;
hx=zeros（ib,jb,kb）;
hy=zeros（ib,jb,kb）;
hz=zeros（ib,jb,kb）;
[X,Y,Z]=meshgrid（1:
ib,1:
jb,1:
kb）;%Gridcoordinates
Psim=zeros,1）;
Panl=zeros,1）;
if（==1）&==0））
x=ceil（ie/
=[1:
x-1:
ie-x];
=[x+1:
x-1:
ie];
=[1,1];
=[je,je];
m2=1;
forn=1:
1:
form1=1:
1:
[ex,ey,ez]=Efields（param,handles,ex,ey,ez,hx,hy,hz,ie,je,ke,ib,jb,kb,n,m1,m2,p）;
[hx,hy,hz]=Hfields（param,hx,hy,hz,ex,ey,ez,ie,je,ke,ib,jb,kb,n,m1,m2,p）;
end
[Psim（n）,Panl（n）]=Cavity_Power（param,handles,ex,ey,ez,n）;
field_viz（param,handles,ex,ey,ez,X,Y,Z,n,Psim,Panl,p）;
Dyn_FFT
pause;
end
elseif（==0）&==1））
y=ceil（je/;
=[1:
y-1:
je-y];
=[y+1:
y-1:
je];
=[1,1];
=[ie,ie];
m1=1;
forn=1:
1:
form2=1:
1:
[ex,ey,ez]=Efields（param,handles,ex,ey,ez,hx,hy,hz,ie,je,ke,ib,jb,kb,n,m1,m2,p）;
[hx,hy,hz]=Hfields（param,hx,hy,hz,ex,ey,ez,ie,je,ke,ib,jb,kb,n,m1,m2,p）;
end
[Psim（n）,Panl（n）]=Cavity_Power（param,handles,ex,ey,ez,n）;
field_viz（param,handles,ex,ey,ez,X,Y,Z,n,Psim,Panl,p）;
pause;
end
elseif（==1）&==1））
x=ceil（ie/;
=[1:
x-1:
ie-x];
=[x+1:
x-1:
ie];
y=ceil（je/;
=[1:
y-1:
je-y];
=[y+1:
y-1:
je];
forn=1:
1:
form2=1:
1:
form1=1:
1:
[ex,ey,ez]=Efields（param,handles,ex,ey,ez,hx,hy,hz,ie,je,ke,ib,jb,kb,n,m1,m2,p）;
[hx,hy,hz]=Hfields（param,hx,hy,hz,ex,ey,ez,ie,je,ke,ib,jb,kb,n,m1,m2,p）;
end
end
[Psim（n）,Panl（n）]=Cavity_Power（param,handles,ex,ey,ez,n）;
field_viz（param,handles,ex,ey,ez,X,Y,Z,n,Psim,Panl,p）;
pause;
end
else
=1;
=ie;
=1;
=je;
m1=1;m2=1;
forn=1:
1:
[ex,ey,ez]=Efields（param,handles,ex,ey,ez,hx,hy,hz,ie,je,ke,ib,jb,kb,n,m1,m2,p）;
[hx,hy,hz]=Hfields（param,hx,hy,hz,ex,ey,ez,ie,je,ke,ib,jb,kb,n,m1,m2,p）;
SimRunStop=get,'value'）;
ifSimRunStop==1
h=warndlg（'SimulationRunisStoppedbyUser!
','!
!
Warning!
!
'）;
waitfor（h）;
break;
end
[Psim（n）,Panl（n）]=Cavity_Power（param,handles,ex,ey,ez,n
升级会员 