內建結構
中文/English
其中有些結構有
x=; 對應角度 alpha
y=; 對應角度 gama
z=; 對應角度 beta
這三個是輸入尤拉轉角(Euler angle),可以將結構旋轉,如上圖
尤拉轉角請參考
http://zh.wikipedia.org/wiki/%E6%AC%A7%E6%8B%89%E8%A7%92
這裡的x對應於上述網址的alpha, y對應於gama, z對應beta
這是程式裡我所使用的轉角程式碼
cosx=cos(x*pi/180);sinx=sin(x*pi/180);
cosy=cos(y*pi/180);siny=sin(y*pi/180);
cosz=cos(z*pi/180);sinz=sin(z*pi/180);
e1=cosx*cosz-cosy*sinx*sinz;
e2=cosy*cosz*sinx+cosx*sinz;
e3=sinx*siny;
e4=-cosz*sinx-cosx*cosy*sinz;
e5=cosx*cosy*cosz-sinx*sinz;
e6=cosx*siny;
e7=siny*sinz;
e8=-cosz*siny;
e9=cosy;
下方nindex是指er或mr
程式裡,結構寫法的起始點為xposition,yposition,zposition
高度h則是在結構中間往上下擴張
例如:
其中
icenter =結構尺吋 X/2=200/2= 100
jcenter =結構尺吋 Y/2=200/2= 100
kcenter =結構尺吋 Z/2=200/2= 100
=================================
內建結構
=================================
%=================橢圓孔洞 Hole(Elliptical)==========================
radius_x=600e-9; %diameter=2*radius
radius_y=800e-9; %diameter=2*radius
h=400e-9;
nindex=1^2;
sigma=0;
choice='E_Iso'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=kcenter*gdz;
x=0; %角度
y=0;
z=0;
Iso_EllipticalHole(choice,gridtype,nindex,sigma,xposition,yposition,zposition,radius_x,radius_y,h,x,y,z)
%=================橢圓孔洞 Hole(Elliptical)==========================
%=================橢圓柱形 Rod(Elliptical)==========================
radius_x=600e-9; %diameter=2*radius
radius_y=800e-9; %diameter=2*radius
h=400e-9;
nindex=1^2;
sigma=0;
choice='PEC'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=kcenter*gdz;
x=0; %角度
y=0;
z=0;
Iso_EllipticalRod(choice,gridtype,nindex,sigma,xposition,yposition,zposition,radius_x,radius_y,h,x,y,z)
%=================橢圓柱形 Rod(Elliptical)==========================
%=================橢圓環狀 Ring(Elliptical)==========================
radius_x=500e-9; %diameter=2*radius
radius_y=600e-9; %diameter=2*radius
radius_x2=800e-9; %diameter=2*radius
radius_y2=900e-9; %diameter=2*radius
h=400e-9;
nindex=1^2;
sigma=0;
choice='PEC'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=kcenter*gdz;
x=0; %角度
y=0;
z=0;
Iso_EllipticalRing(choice,gridtype,nindex,sigma,xposition,yposition,zposition,radius_x,radius_y,radius_x2,radius_y2,h,x,y,z)
%=================橢圓環狀 Ring(Elliptical)==========================
%=================橢圓環狀孔洞 Ring Hole(Elliptical)==========================
radius_x=500e-9; %diameter=2*radius
radius_y=600e-9; %diameter=2*radius
radius_x2=800e-9; %diameter=2*radius
radius_y2=900e-9; %diameter=2*radius
h=400e-9;
nindex=1^2;
sigma=0;
choice='PEC'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=kcenter*gdz;
x=0; %角度
y=0;
z=0;
Iso_EllipticalRingHole(choice,gridtype,nindex,sigma,xposition,yposition,zposition,radius_x,radius_y,radius_x2,radius_y2,h,x,y,z)
%=================橢圓環狀孔洞 Ring Hole(Elliptical)==========================
%=================雙球 Two_sphere===========================
radius1=600e-9; %diameter=2*radius
xposition1=icenter*gdx;
yposition1=jcenter*gdy;
zposition1=kcenter*gdz-300e-9;
radius2=600e-9; %diameter=2*radius
xposition2=icenter*gdx;
yposition2=jcenter*gdy;
zposition2=kcenter*gdz+300e-9;
nindex=1^2;
sigma=0;
choice='PEC'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %electric material -1; magnetic material 1; other 0
Iso_Two_Sphere(choice,gridtype,nindex,sigma,xposition1,yposition1,zposition1,radius1,xposition2,yposition2,zposition2,radius2)
%=================雙球 sphere===========================
%=================矩形孔洞 SquareHole===========================
lengthx=200e-9;
lengthy=150e-9;
lengthz=100e-9;
nindex=1^2;
sigma=0;
choice='E_Model1'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=kcenter*gdz;
Iso_SquareHole(choice,gridtype,nindex,sigma,xposition,yposition,zposition,lengthx,lengthy,lengthz)
%=================矩形孔洞 SquareHole===========================
%=================球 sphere===========================
radius=150e-9; %diameter=2*radius
nindex=1^2;
sigma=0;
choice='E_Model1'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %electric material -1; magnetic material 1; other 0
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=kcenter*gdz;
Iso_Sphere(choice,gridtype,nindex,sigma,xposition,yposition,zposition,radius)
%=================球 sphere===========================
%=================狹縫 Slit===========================
lengthx=100e-9;
lengthz=200e-9;
nindex=1.5^2;
sigma=0;
choice='E_Iso'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=kcenter*gdz;
Iso_Slit(choice,gridtype,nindex,sigma,xposition,zposition,lengthx,lengthz)
%=================狹縫 Slit===========================
%=================橢圓 oval===========================
radiusx=50e-9; %diameter=2*radius
radiusy=100e-9; %diameter=2*radius
radiusz=150e-9; %diameter=2*radius
nindex=1^2;
sigma=0;
choice='PMC'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=kcenter*gdz;
Iso_oval(choice,gridtype,nindex,sigma,xposition,yposition,zposition,radiusx,radiusy,radiusz)
%=================橢圓 oval===========================
%=================矩形 Brick===========================
xstart=50e-9;
xend=100e-9;
ystart=50e-9;
yend=200e-9;
zstart=50e-9;
zend=300e-9;
nindex=1^2;
sigma=0;
choice='E_Iso'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
Iso_Brick(choice,gridtype,nindex,sigma,xstart,xend,ystart,yend,zstart,zend)
%=================矩形 Brick===========================
%=================橢圓錐形 EllipticalCone==========================
R1=100e-9;
h=300e-9;
nindex=1^2;
sigma=0;
choice='E_Model1'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=kcenter*gdz;
x=0; %角度
y=0;
z=0;
Iso_EllipticalCone(choice,gridtype,nindex,sigma,xposition,yposition,zposition,R1,h,x,y,z)
%=================橢圓錐形 EllipticalCone==========================
%=================錐形 Cone==========================
R1=100e-9;
R2=10e-9;
h=150e-9;
nindex=1^2;
sigma=0;
choice='E_Model1'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=10*gdz;
x=0; %角度
y=0;
z=0;
Iso_Cone(choice,gridtype,nindex,sigma,xposition,yposition,zposition,R1,R2,h,x,y,z)
%=================錐形 Cone==========================
%=================柱形(3邊) rod(3-side)==========================
R1=100e-9;
R2=150e-9;
h=200e-9;
nindex=1^2;
sigma=0;
choice='E_Model1'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=10*gdz;
x=0; %角度
y=0;
z=0;
Iso_3sideRod(choice,gridtype,nindex,sigma,xposition,yposition,zposition,R1,R2,h,x,y,z)
%=================柱形(3邊) rod(3-side)==========================
%=================柱形(4邊) rod(4-side)==========================
R1=150e-9;
R2=100e-9;
h=200e-9;
nindex=1^2;
sigma=0;
choice='E_Model1'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=10*gdz;
x=0; %角度
y=0;
z=0;
Iso_4sideRod(choice,gridtype,nindex,sigma,xposition,yposition,zposition,R1,R2,h,x,y,z)
%=================柱形(4邊) rod(4-side)==========================
%=================柱形(5邊) rod(5-side)==========================
R1=150e-9;
R2=100e-9;
h=200e-9;
nindex=1^2;
sigma=0;
choice='E_Model1'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=10*gdz;
x=0; %角度
y=0;
z=0;
Iso_5sideRod(choice,gridtype,nindex,sigma,xposition,yposition,zposition,R1,R2,h,x,y,z)
%=================柱形(5邊) rod(5-side)==========================
%=================柱形(6邊) rod(6-side)==========================
R1=150e-9;
R2=100e-9;
h=200e-9;
nindex=1^2;
sigma=0;
choice='E_Model1'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=10*gdz;
x=0; %角度
y=0;
z=0;
Iso_6sideRod(choice,gridtype,nindex,sigma,xposition,yposition,zposition,R1,R2,h,x,y,z)
%=================柱形(6邊) rod(6-side)==========================
%=================柱形(圓底) rod(circle)==========================
radius=150e-9; %diameter=2*radius
h=200e-9;
nindex=1^2;
sigma=0;
choice='E_Model1'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=kcenter*gdz;
x=0; %角度
y=0;
z=0;
Iso_CylinderRod(choice,gridtype,nindex,sigma,xposition,yposition,zposition,radius,h,x,y,z)
%=================柱形(圓底) rod(circle)==========================
%=================孔洞(圓底) Hole(cylinder)==========================
radius=100e-9; %diameter=2*radius
h=250e-9;
nindex=1^2;
sigma=0;
choice='E_Iso'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=kcenter*gdz;
x=33; %角度
y=22;
z=11;
Iso_CylinderHole(choice,gridtype,nindex,sigma,xposition,yposition,zposition,radius,h,x,y,z)
%=================孔洞(圓底) Hole(cylinder)==========================
綜合的範例 (加轉角)
%=================柱形(3邊) rod(3-side)==========================
R1=20e-9;
R2=70e-9;
h=150e-9;
nindex=1.5^2;
sigma=0;
choice='E_Iso'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=kcenter*gdz;
x=33; %角度
y=66;
z=0;
Iso_3sideRod(choice,gridtype,nindex,sigma,xposition,yposition,zposition,R1,R2,h,x,y,z)
%=================柱形(3邊) rod(3-side)==========================
%=================柱形(6邊) rod(6-side)==========================
R1=20e-9;
R2=70e-9;
h=140e-9;
nindex=1^2;
sigma=0;
choice='E_Model1'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=kcenter*gdz;
x=35; %角度
y=135;
z=0;
Iso_6sideRod(choice,gridtype,nindex,sigma,xposition,yposition,zposition,R1,R2,h,x,y,z)
%=================柱形(6邊) rod(6-side)==========================
%=================錐形 EllipticalCone==========================
R1=60e-9;
h=250e-9;
nindex=1^2;
sigma=0;
choice='E_Model2'; %E_Iso,PEC,M_Iso,PMC,E_Model1,M_Model1,EM_Model1
gridtype=-1; %
xposition=icenter*gdx;
yposition=jcenter*gdy;
zposition=kcenter*gdz;
x=55; %角度
y=60;
z=49;
Iso_EllipticalCone(choice,gridtype,nindex,sigma,xposition,yposition,zposition,R1,h,x,y,z)
%=================錐形 EllipticalCone==========================