1234

所属分类:matlab编程
开发工具:matlab
文件大小:48KB
下载次数:28
上传日期:2012-05-21 13:23:05
上 传 者liuzecheng
说明:  计算光子晶体能带结构 使用matlab,程序非常容易理解。还配有说明文档。
(Calculation of the photonic crystal band structure using the Matlab program is very easy to understand. Document comes with instructions.)

文件列表:
光子晶体能带 Matlab (0, 2012-05-21)
光子晶体能带 Matlab\2d_fdtd.htm (3785, 2002-06-27)
__MACOSX (0, 2012-05-21)
__MACOSX\光子晶体能带 Matlab (0, 2012-05-21)
__MACOSX\光子晶体能带 Matlab\._2d_fdtd.htm (240, 2002-06-27)
光子晶体能带 Matlab\pmltm (129844, 2002-06-27)
光子晶体能带 Matlab\pwsrc.prb (365, 2002-06-27)
__MACOSX\光子晶体能带 Matlab\._pwsrc.prb (240, 2002-06-27)
__MACOSX\._光子晶体能带 Matlab (240, 2012-05-21)

Here is the input data file format Build V1.0.1 Date: 5/20/12 By: Michael Liu liu.zecheng@h.mbox.nagoya-u.ac.jp 1-757-4231507 PML2DTM v1.01 PML2DTE v1.01 are compiled and run on Solaris, copy and run. tar -xvf *.tar to extract all files in tar Note: 1. follow format strictly, else the data will not be read correctly and program will not run correctly 2. Verification is written in echo2dtd.prb 3. Lower/upper case matters 4. no blank lines between data 5. no comments allowed between data !**sampletstep, integer, the time interval for writing into GIF file !**negative means no writing** -1 !**nb,nu,R0, !**nb: integer, number of PML layers** !**nu: double, the order of PML function !**R0: double, the normal reflection coefficient, affect a lot 10 2.0 1.0E-10 !**nx,dx,ny,dy*** !**nx,ny: integer, number of cells along x and y direction !**dx,dy: double, size of cell along x and y direction 101 0.01 101 0.01 !**dt,nt** !**dt: double, the time step, should follow stability check !**nt: integer, number of time steps to be calculated 1.6667E-11 400 !**x1,y1,x2,y2,srcshape,wave,amp,freq,t0,tw,srctype,tsrcstop** !**to describe the source !**x1,y1,x2,y2: integer, the 2 points of a rectangle. !**x1,y1: the coordinate of the lower-left point !**x2,y2: the coordinate of the upper-right point !**srcshape: how the source is distributed in the rectangle, !** UN: uniform !**wave: only support these 3 sine: a sine wave gauss_lp: gauss pulse gauss_bp: gauss pulse modulated with a sine wave !**amp: double, amplitude !**freq: double, frequency !**t0: time delay, any for sine,should put 0 !**tw: pulse width, any for sine,should put 0 !**srctype: only have 2 value: hard: hard source soft: treat as a soft transparent source !**tsrcstop: double, when to stop the source, (in unit of tw) !** 0: will not stop during simulation gauss_bp 1000.0 477460000.0 1.3333E-9 5.0E-10 soft !**numofmats: integer, number of materials** 3 !**mater,epsr,mur** !**mater: material number, 100 is for source, 200 is for PEC !**epsr: double, epsr !**mur: double, mur 1 1.0 1.0 100 1.0 1.0 200 1.0 1.0 !**numofregions: number of regions** 1 !**mater,shape,x1,y1,x2,y1,patchtype,radius** !**mater: material number, refer above !**shape: the shape of this region R: rectangle D: disk, including circle and ellipse !**x1,y1,x2,y2, define the rectangle region !**patchtype: the sub-type in the rectangle region NA: rectangle LL: low-left triangle UL: upper-left triangle LR: lower right triangle UR: upper-right triangle !**radius: only used for disk, I am now considering to remove it 1 R 1 1 101 101 NA 0.0E+0 !**number of source rect 1 !**srcpatch, x1,y1,x2,y2,patchtype** R 51 51 51 51 UN !**PCsupport: only 2 valuse support: PC:YES and PC:NO** PC:NO !**if photonic crystal is present, ie. PC:YES, then the following data !**is needed for input !**1. numofatoms: integer !**2. mater, llx,lly,urx,ury,type,x_ext,y_ext !**mater: material number of the atom !**llx,lly,urx,ury: defines the rectangle region for the first atom !** on the lower-left !**type: the shape of the atom,"R" and "D" is now supported !**x_ext,y_ext: how many same atoms are going to place in u1 and !** u2 direction (U1 and U2 are the unit axis of the photonic crystal) !**numofplotfiles** 11 !**plottype,field,x1,y1,x2,y2,tstep** !**plottype: T: time series R: field in a rectangle M: material distribution in a region P: power flux !**fieldtype: Ex, Ey, Ez, Hx, Hy,Hz !**x1,y1,x2,y2,define the region to output !**tstep: the frequency to record T Ez 1 51 1 1 1 T Ez 11 51 1 1 1 T Ez 21 51 1 1 1 T Ez 31 51 1 1 1 T Ez 41 51 1 1 1 T Ez 51 51 1 1 1 T Ez 61 51 1 1 1 T Ez 71 51 1 1 1 T Ez 81 51 1 1 1 T Ez 91 51 1 1 1 T Ez 101 51 1 1 1

近期下载者

相关文件


收藏者