風行電磁波模擬 (FonSinEM) - Free ElectroMagnetic Simulator

內建結構

中文/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==========================

Page updated

Google Sites

Report abuse