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
近期下载者:
相关文件:
收藏者: