环形光波导耦合器件程序

环形光波导耦合器件程序

ID:20829481

大小:101.50 KB

页数:15页

时间:2018-10-16

环形光波导耦合器件程序_第1页
环形光波导耦合器件程序_第2页
环形光波导耦合器件程序_第3页
环形光波导耦合器件程序_第4页
环形光波导耦合器件程序_第5页
资源描述:

《环形光波导耦合器件程序》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、10度:%Thisprogramcalculatethewaveguidebyaphotoniccrystal.%Fortwodimensioncase.%OnlyforTMcase(Ez,Hx,Hy),inthisversion试验版本.clear;%tic%Initialparametersandotherthings.W=0.36;%Normalizedfrequency就是能带图中的纵坐标,单位为(a/lamda)%Thefollowingparametersarecontrolparameters.WaveGuide=1;%IfthisprogramisforWaveGuide

2、?Ifso,pleasespecify1.IsMovie=0;%Ifyouwanttoplaymovie,pleaseuse1.IsFigure=1;%Ifitwillplotthefigures?Ifso,pleasespecify1.WantToSeeEp=1;%DoyouwanttoseethedistrubutionofEp(TE极化,S极化)?Ifso,pleasespecify1.%EndofdefiningcontrolparametersMLatx=11;%HowmanyLatticecellinxdirection.栅格MLaty=11;%HowmanyLatticec

3、ellinydirection.NMlat=21;%ThegirdnumberineachLatticeCell.每个栅格中包含的fdtd网格个数%SHOULDBEODDINTEGER!!!!奇数整数个ifmod(NMlat,2)==0NMlat=NMlat+1;end%Forceittobeaoddinteger!NTx=MLatx*NMlat+1;%ItisthenumberoftheGridalongxaxis.总的fdtd网格坐标点(所以要加1)NTy=MLaty*NMlat+1;%ItisthenumberoftheGridalongyaxis.ifWaveGuide==1Nr

4、ow=5;%TherownumberofcoloumnsbetweenthePMLboundaryandthewaveguide.(此处定义,现在还不知所措)endNPML=12;%HowmanyPMLlayerswillbeusedinourcomputation.PML的层数NTimeSteps=5000;%TotalnumberofTimeStepsMeach=20;%DefinetheintervalforplotfiguresifIsFigure==1.%Thisalsoworksforsavingintervals.每隔20个时间步显示一次。R=0.2;%Theradiuso

5、fdielectriccolumns,ea=11.4;%Thedielectricconstantofthesecolumns.Zmax=0.8;%Themaximumvalueforzaxiswhenplottingfigures.Colormax=2.8;%Themaximumvalueforcolormapwhenplottingfigures.%Someconstantsmu0=4*pi*1.0e-7;%EpsilonZero,ifusingGaussUnit,itequalsto1.e0=8.85*1e-12;%MuZero,ifusingGaussUnit,itequalst

6、o1.c=1/sqrt(mu0*e0);%Thelightspeed.factor=mu0/e0;%Thefactorbetweenconductivityandpermeability.%Permeability=Conductivity*factor,inPML.PML中的匹配条件a=1;%e-6;%Thelatticeconstant.栅格常数。W=W*(2*pi*c/a);%frequencyDx=a/NMlat;%Deltax.dx=1/(一个栅格中的网格数)Dy=Dx;%Deltay.Dt=1/sqrt(1/(Dx*Dx)+1/(Dy*Dy))/c;%Timeinterval

7、%tic%Inthefollowingpartmwedefinethedielectric%constants:,定义FDTD中平面波导与空气的介电常数Ep=ones(NTx-1,NTy-1)*e0;%FDTD总的网格中的介质常数,真空。Ep_cell=ones(NMlat,NMlat)*e0;%每个栅格中的fdtd介电常数,真空。x=-(NMlat-1)/2*Dx:Dx:(NMlat-1)/2*Dx;%先对栅格操作,然后进行rep

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。