差分法matlab程序,FDTD(时域有限差分法)算法Matlab源程序.doc

%*********************************************************************** [color=#FFFFFF]本文来自:[/color] ?%? ? 3-D FDTD code with PEC boundaries [color=#FFFFFF]本文来自:[/color] ?%*********************************************************************** [color=#FFFFFF]本文来自:[/color] ?% [color=#FFFFFF]本文来自:[/color] ?%? ? Program author: Susan C. Hagness [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ? ? ? ? ? Department of Electrical and Computer Engineering [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ? ? ? ? ? University of Wisconsin-Madison [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ? ? ? ? ? 1415 Engineering Drive [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ? ? ? ? ? Madison, WI 53706-1691 [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ? ? ? ? ? 608-265-5739 [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ? ? ? ? ? hagness@ [color=#FFFFFF]本文来自:[/color] ?% [color=#FFFFFF]本文来自:[/color] ?%? ? Date of this version:? February 2000 [color=#FFFFFF]本文来自:[/color] ?% [color=#FFFFFF]本文来自:[/color] ?%? ? This MATLAB M-file implements the finite-difference time-domain [color=#FFFFFF]本文来自:[/color] ?%? ? solution of Maxwell's curl equations over a three-dimensional [color=#FFFFFF]本文来自:[/color] ?%? ? Cartesian space lattice comprised of uniform cubic grid cells. [color=#FFFFFF]本文来自:[/color] ?% [color=#FFFFFF]本文来自:[/color] ?%? ? To illustrate the algorithm, an air-filled rectangular cavity [color=#FFFFFF]本文来自:[/color] ?%? ? resonator is modeled.? The length, width, and height of the [color=#FFFFFF]本文来自:[/color] ?%? ? cavity are 10.0 cm (x-direction), 4.8 cm (y-direction), and [color=#FFFFFF]本文来自:[/color] ?%? ? 2.0 cm (z-direction), respectively. [color=#FFFFFF]本文来自:[/color] ?% [color=#FFFFFF]本文来自:[/color] ?%? ? The computational domain is truncated using PEC boundary [color=#FFFFFF]本文来自:[/color] ?%? ? conditions: [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ex(i,j,k)=0 on the j=1, j=jb, k=1, and k=kb planes [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ey(i,j,k)=0 on the i=1, i=ib, k=1, and k=kb plane