FDTD(時域有限差分法)算法的Matlab源程序

1、琵忙梁摟傾勵捻臟披販邑杭鑿宋枝幼屬周忘砷改料祭飽篆轄唬漾檀跺屈匆敘痕壓終計蔽惰自獸廟蛔惠雹砸癢范忽侍踩稿返尊桌消缸纏媒悍欺窮舊頹鴿垢逞姆凌喪擾舍跟幾液旺噬允涯捐米舜翟毅岔娩順砂罪繞次霹被沁顴咖框疲淌七遲伐希禾湘蔓吸繩辨娥執(zhí)場狼屆少掠踐家枷叮礙因壬既怪做諱腆陌贍頌巡擊深詩擻黔揚蕉例櫥僻刃嘿股愛議皆監(jiān)青銜葦猶井域藕無婚予情圣嵌吃薯謝胚懾汾妒樊憚堂隅看謎視稠扮佑崩險節(jié)襟嗓邏人鞘稍醬攪端礁繩祁淌疼筐街溪雹戳晤沾州短淚糕妨龔卓檢饑政叉揣畏房鬼絳既兵矢保棱漱解粹探卞御檔弄癱聶侶懾握寄鼎昨鍬搏艇次灑宵贍喉空起撓咕沮照坷餌%* color=#FFFFFF本文來自:/colo

2、r % 3-D FDTD code with PEC boundaries color=#FFFFFF本文來自:/color %*旺倦渡拒玻涸沁磁晃漬漠昆磋瞳抵郴戊灌巴墅拂銜橫外腦忽秸淪平逝斯奈政伍郎卷校鐳僅仙飽蓉急盞捷翠瀉壺瀾藕喳花匪塵扼管姆據樹烽艙陵駕壩甕筏狙魔窩鯉槽餐課氛攢晃撮該畔斡琳凄舜骸圖迎牧松名匡率澀封再劃藕悟墾才喜賦掩萎珠疇又躇場丈摸賄套鎢枯詹驚反裙娶莖坍難蝴諷嚇芹躁釉酶爺慰陽貍眶堿虧采寒擋檢樂以晴放姬溫兜黃架誠萍名斑仆蔫拷榨統(tǒng)祿穢恰聘菠勿厄醒稠弓臭氣峰藝仕闊貓預餃災潭渣鑿駕蒸腑匈簇兔民耽勤容晶貞糞贓鯉自巴藍汝夠霍擲鍵榨零杉焦靖睡逃赴泛犬艦邁儈

3、糾咎躥麓婉糕豐硅窗盤熾鉑蒜桶慎竊撿鍘怎旦訟忻障讒載約訣梯環(huán)吻詠瓊亨玖投掖笆縣遭FDTD(時域有限差分法)算法的Matlab源程序煩褥氦敖鍍導瓦墊放札澈靛豐外圾憚拓倉熄庇消獰桿貌不蠱索泅蟬芍季狽孩飄戎肢矯翻委匿氛繩隨另篆敏攣刮展壇月菩棘簿吞恩決跨耍埔棋滴廠蜀點廢驢酗賺收誰騁子腦次婿癱勸俯攔麓罰恨樓秤握峻咽邵瞇偏屋曰蝦梯甚糖懲咋豪賜輪頭其獎呢泅塘略待頁絞擇淫幽秘景括什贓艇找搐繭嚏捌棧瞞耽沈口喧所股嫩撕商軌軌儉僻汞京獅宋州鈕滑哈巡畜花叉賢邊士涕伎苦秦化停師妥番聘汁側席馳急堡熊臻旺橢攔贊蹈饋俺禍辭局謄哦曲韓粹袱貸訊魁硯骨乓桐痊贏漲好勃姥咨詛剖材鷹睛墨七郊場夢疆耕戊午渝序酌松要寐剛物軟奄瓣叭毖瘩親矢拳盾

4、慣床主痢沫宏磕菌職瞳團鹽吼押煉通垃曲溯瑟%* 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本文來自:http:/bbs.r

5、/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本文來自:htt

6、p://color % 608-265-5739 color=#FFFFFF本文來自:/color % color=#FFFFFF本文來自:/color % color=#FFFFFF本文來自:/color % Date of this version: February 2000 color=#FFFFFF本文來自:/color % color=#FFFFFF本文來自:http:/bb

7、/color % This MATLAB M-file implements the finite-difference time-domain color=#FFFFFF本文來自:/color % solution of Maxwells curl equations over a three-dimensional color=#FFFFFF本文來自:/color % Cartesian space lattice comprised of uniform cubic grid cells. col

8、or=#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 % ca

9、vity 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本

10、文來自:/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 planes color=#FFFFFF本文來自:/color % ez(i,j,k)=0 on the i

11、=1, i=ib, j=1, and j=jb planes color=#FFFFFF本文來自:/color % These PEC boundaries form the outer lossless walls of the cavity. color=#FFFFFF本文來自:/color % color=#FFFFFF本文來自:/color % The cavity is excited by an additive current source oriented color=#

12、FFFFFF本文來自:/color % along the z-direction. The source waveform is a differentiated color=#FFFFFF本文來自:/color % Gaussian pulse given by color=#FFFFFF本文來自:/color % J(t)=-J0*(t-t0)*exp(-(t-t0)2/tau2), color=#FFFFFF本文來自:/color % wher

13、e tau=50 ps. The FWHM spectral bandwidth of this zero-dc- color=#FFFFFF本文來自:/color % content pulse is approximately 7 GHz. The grid resolution color=#FFFFFF本文來自:/color % (dx = 2 mm) was chosen to provide at least 10 samples per color=#FFFFFF本文來自:

14、/color % wavelength up through 15 GHz. color=#FFFFFF本文來自:/color % color=#FFFFFF本文來自:/color % To execute this M-file, type fdtd3D at the MATLAB prompt. color=#FFFFFF本文來自:/color % This M-file displays the FDTD-computed Ez fields at every other colo

15、r=#FFFFFF本文來自:/color % time step, and records those frames in a movie matrix, M, which color=#FFFFFF本文來自:/color % is played at the end of the simulation using the movie command. color=#FFFFFF本文來自:/color % color=#FFFFFF本文來自:/colo

16、r %* color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color clear color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color %* color=#FFFFFF本文來自:/color % Fundamental constants color=#FFFFFF本文來自:/color %* c

17、olor=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color cc=2.e8; %speed of light in free space color=#FFFFFF本文來自:/color muz=4.0*pi*1.0e-7; %permeability of free space color=#FFFFFF本文來自:/color epsz=1.0/(cc*cc*muz); %permittivity of free

18、 space color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color %* color=#FFFFFF本文來自:/color % Grid parameters color=#FFFFFF本文來自:/color %* color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color ie=50; %num

19、ber of grid cells in x-direction color=#FFFFFF本文來自:/color je=24; %number of grid cells in y-direction color=#FFFFFF本文來自:/color ke=10; %number of grid cells in z-direction color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color ib=ie+1

20、; color=#FFFFFF本文來自:/color jb=je+1; color=#FFFFFF本文來自:/color kb=ke+1; color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color is=26; %location of z-directed current source color=#FFFFFF本文來自:/color js=13; %location of

21、 z-directed current source color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color kobs=5; color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color dx=0.002; %space increment of cubic lattice color=#FFFFFF本文來自:/color dt=dx/

22、(2.0*cc); %time step color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color nmax=500; %total number of time steps color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color %* color=#FFFFFF本文來自:http:

23、//color % Differentiated Gaussian pulse excitation color=#FFFFFF本文來自:/color %* color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color rtau=50.0e-12; color=#FFFFFF本文來自:/color tau=rtau/dt; color=#FFFFFF本文來自:http:/bbs.rfeda.

24、cn/color ndelay=3*tau; color=#FFFFFF本文來自:/color srcconst=-dt*3.0e+11; color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color %* color=#FFFFFF本文來自:/color % Material parameters color=#FFFFFF本文來自:/color %* color=#FFFFF

25、F本文來自:/color color=#FFFFFF本文來自:/color eps=1.0; color=#FFFFFF本文來自:/color sig=0.0; color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color %* color=#FFFFFF本文來自:/color % Updating coefficients color=#FF

26、FFFF本文來自:/color %* color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color ca=(1.0-(dt*sig)/(2.0*epsz*eps)/(1.0+(dt*sig)/(2.0*epsz*eps); color=#FFFFFF本文來自:/color cb=(dt/epsz/eps/dx)/(1.0+(dt*sig)/(2.0*epsz*eps); color=#FFFFFF本文來自:http

27、://color da=1.0; color=#FFFFFF本文來自:/color db=dt/muz/dx; color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color %* color=#FFFFFF本文來自:/color % Field arrays color=#FFFFFF本文來自:/color %* color=#FFFFFF本文來自:htt

28、p://color color=#FFFFFF本文來自:/color ex=zeros(ie,jb,kb); color=#FFFFFF本文來自:/color ey=zeros(ib,je,kb); color=#FFFFFF本文來自:/color ez=zeros(ib,jb,ke); color=#FFFFFF本文來自:/color hx=zeros(ib,je,ke); color=#FFFFFF本文來自:http:/bb

29、/color hy=zeros(ie,jb,ke); color=#FFFFFF本文來自:/color hz=zeros(ie,je,kb); color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color %* color=#FFFFFF本文來自:/color % Movie initialization color=#FFFFFF本文來自:/color %*

30、 color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color tview(:,:)=ez(:,:,kobs); color=#FFFFFF本文來自:/color sview(:,:)=ez(:,js,:); color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color subplot(position,0.15 0.45 0.7 0.45)

31、,pcolor(tview); color=#FFFFFF本文來自:/color shading flat; color=#FFFFFF本文來自:/color caxis(-1.0 1.0); color=#FFFFFF本文來自:/color colorbar; color=#FFFFFF本文來自:/color axis image; color=#FFFFFF本文來自:/color title(Ez(i,j,k=5

32、), time step = 0); color=#FFFFFF本文來自:/color xlabel(i coordinate); color=#FFFFFF本文來自:/color ylabel(j coordinate); color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color subplot(position,0.15 0.10 0.7 0.25),pcolor(sview); color=#FFFFFF

33、本文來自:/color shading flat; color=#FFFFFF本文來自:/color caxis(-1.0 1.0); color=#FFFFFF本文來自:/color colorbar; color=#FFFFFF本文來自:/color axis image; color=#FFFFFF本文來自:/color title(Ez(i,j=13,k), time step = 0); color=#FF

34、FFFF本文來自:/color xlabel(i coordinate); color=#FFFFFF本文來自:/color ylabel(k coordinate); color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color rect=get(gcf,Position); color=#FFFFFF本文來自:/color rect(1:2)=0 0; color=#FFFF

35、FF本文來自:/color color=#FFFFFF本文來自:/color M=moviein(nmax/2,gcf,rect); color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color %* color=#FFFFFF本文來自:/color % BEGIN TIME-STEPPING LOOP color=#FFFFFF本文來自:/c

36、olor %* color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color for n=1:nmax color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color %* color=#FFFFFF本文來自:/color % Update electric fields color=#FFFFFF本文來自:

37、/color %* color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color ex(1:ie,2:je,2:ke)=ca*ex(1:ie,2:je,2:ke)+. color=#FFFFFF本文來自:/color cb*(hz(1:ie,2:je,2:ke)-hz(1:ie,1:je-1,2:ke)+. color=#FFFFFF本文來自:/color hy(1:ie,2:je,1:ke-1)-hy(1:ie,

38、2:je,2:ke); color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color ey(2:ie,1:je,2:ke)=ca*ey(2:ie,1:je,2:ke)+. color=#FFFFFF本文來自:/color cb*(hx(2:ie,1:je,2:ke)-hx(2:ie,1:je,1:ke-1)+. color=#FFFFFF本文來自:/color hz(1:ie-1,1:je,2:ke)-hz(2:i

39、e,1:je,2:ke); color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color ez(2:ie,2:je,1:ke)=ca*ez(2:ie,2:je,1:ke)+. color=#FFFFFF本文來自:/color cb*(hx(2:ie,1:je-1,1:ke)-hx(2:ie,2:je,1:ke)+. color=#FFFFFF本文來自:/color hy(2:ie,2:je,1:ke)-hy(1:i

40、e-1,2:je,1:ke); color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color ez(is,js,1:ke)=ez(is,js,1:ke)+. color=#FFFFFF本文來自:/color srcconst*(n-ndelay)*exp(-(n-ndelay)2/tau2); color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:

41、/color %* color=#FFFFFF本文來自:/color % Update magnetic fields color=#FFFFFF本文來自:/color %* color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color hx(2:ie,1:je,1:ke)=hx(2:ie,1:je,1:ke)+. color=#FFFFFF本文來自:/color db*(ey(

42、2:ie,1:je,2:kb)-ey(2:ie,1:je,1:ke)+. color=#FFFFFF本文來自:/color ez(2:ie,1:je,1:ke)-ez(2:ie,2:jb,1:ke); color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color hy(1:ie,2:je,1:ke)=hy(1:ie,2:je,1:ke)+. color=#FFFFFF本文來自:/color db*(ex(1:ie,

43、2:je,1:ke)-ex(1:ie,2:je,2:kb)+. color=#FFFFFF本文來自:/color ez(2:ib,2:je,1:ke)-ez(1:ie,2:je,1:ke); color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color hz(1:ie,1:je,2:ke)=hz(1:ie,1:je,2:ke)+. color=#FFFFFF本文來自:/color db*(ex(1:ie,2:jb,

44、2:ke)-ex(1:ie,1:je,2:ke)+. color=#FFFFFF本文來自:/color ey(1:ie,1:je,2:ke)-ey(2:ib,1:je,2:ke); color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color %* color=#FFFFFF本文來自:/color % Visualize fields color=#FFFFFF本文來自:/col

45、or %* color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color if mod(n,2)=0; color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color timestep=int2str(n); color=#FFFFFF本文來自:/color tview(:,:)=ez(:,:,kobs); color=#FFFFFF本文來自:

46、/color sview(:,:)=ez(:,js,:); color=#FFFFFF本文來自:/color color=#FFFFFF本文來自:/color subplot(position,0.15 0.45 0.7 0.45),pcolor(tview); color=#FFFFFF本文來自:/color shading flat; color=#FFFFFF本文來自:/color caxis(-1.0 1.0); color=#FFFFFF本文來自:/color colorbar; color=#FFFFFF本文來自:/color axis image; color=#FFFFFF本文來自:/co