资源描述：
《电磁场fdtd算法以及仿真图》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、function[output_args]=Untitled2(input_args)%UNTITLED2Summaryofthisfunctiongoeshere%Detailedexplanationgoeshere%***********************************************************************%3-DFDTDcodewithPECboundaries%***********************************************************************%%Prog
2、ramauthor:SusanC.Hagness%DepartmentofElectricalandComputerEngineering%UniversityofWisconsin-Madison%1415EngineeringDrive%Madison,WI53706-1691%608-265-5739%hagness@engr.wisc.edu%%Dateofthisversion:February2000%%ThisMATLABM-fileimplementsthefinite-differencetime-domain%solutionofMaxwell'scu
3、rlequationsoverathree-dimensional%Cartesianspacelatticecomprisedofuniformcubicgridcells.%%Toillustratethealgorithm,anair-filledrectangularcavity%resonatorismodeled.Thelength,width,andheightofthe%cavityare10.0cm(x-direction),4.8cm(y-direction),and%2.0cm(z-direction),respectively.%%Thecompu
4、tationaldomainistruncatedusingPECboundary%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.%%Thecavityisexcitedbyanadditivecurrentsourceorient
5、ed%alongthez-direction.Thesourcewaveformisadifferentiated%Gaussianpulsegivenby%J(t)=-J0*(t-t0)*exp(-(t-t0)^2/tau^2),%wheretau=50ps.TheFWHMspectralbandwidthofthiszero-dc-%contentpulseisapproximately7GHz.Thegridresolution%(dx=2mm)waschosentoprovideatleast10samplesper%wavelengthupthrough15GH
6、z.%%ToexecutethisM-file,type"fdtd3D"attheMATLABprompt.%ThisM-filedisplaystheFDTD-computedEzfieldsateveryother%timestep,andrecordsthoseframesinamoviematrix,M,which%isplayedattheendofthesimulationusingthe"movie"command.%这个MATLAB的m文件实现了麦克斯韦旋度方程的有限差分方法在三维笛卡尔空间点阵组成的统一的立方网格细胞。为了说明算法,建模一个充气矩形空腔谐
7、振器。腔的长度、宽度和高度是10.0厘米(x方向)(y方向)4.8厘米,2.0厘米(z三个方向),分别为。计算域边界截断使用压电陶瓷这些压电陶瓷外无损耗腔壁边界形式。腔是由面向添加剂电流源的兴奋沿着z三个方向。源波形是一个差异化的高斯脉冲的半最大值光谱带宽zero-dc-内容脉冲大约是7GHz。网格分辨率(dx=2毫米)被选为至少10样品每波长通过15GHz。执行这个m文件,输入“fdtd3D”MATLAB提示。这个m文件显示在每个其他FDTD-computedEz字段时间步长,并记录这些框架在电影矩阵,M播放结束时模拟使用“
