光子晶体传输矩阵PPT课件_第1页
光子晶体传输矩阵PPT课件_第2页
光子晶体传输矩阵PPT课件_第3页
光子晶体传输矩阵PPT课件_第4页
光子晶体传输矩阵PPT课件_第5页
已阅读5页,还剩6页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

.,1,Matlab仿真一维光子晶体禁带和缺陷模,理论模型:,传输矩阵,每一层的特征矩阵:,ni,di分别第i层介质的折射率和厚度,是光线在第i层介质中的传播方向与界面法线的夹角,是入射光的波长.,空气中,由单层模推广到多层模:,=n1cos(),.,2,a,b,(ab)N二元结构,一个ab单元的特征矩阵:,任意一层的特征矩阵,整个一维光子晶体的特征矩阵:,M=(M0)N,.,3,光子晶体周围材料的折射率设为n0,在这里我们只考虑TE波0=ncos,光在光子晶体中传播时的反射系数,光在光子晶体中传播时的透射系数,反射率,投射率,因此:从以上方程可以看出,只要知道了一维光子晶体的特征矩阵中的每个元素,就可以求出反射率和透射率。,.,4,M=(M0)N,Matlab程序执行,clearall,%清除变量,na=2.35,%a介质层折射率,nb=1.38,%b介质层折射率,ha=63.8e-9,%a介质层厚度,hb=108.7e-9,%b介质层厚度,yeta1=2.35,yeta2=1.38,bo=400:1:900,%b入射波长范围,derta1=2*pi*na*ha./(bo*1e-9),%a介质中相位变化,正入射情况,二元周期结构,derta2=2*pi*nb*hb./(bo*1e-9),%b介质中相位变化,num=length(bo);,.,5,M=(M0)N,forj=1:num,Ma=cos(derta1(j),-i*sin(derta1(j)./yeta1;-i.*yeta1.*sin(derta1(j),cos(derta1(j);,%a介质特征矩阵,Mb=cos(derta2(j),-i*sin(derta2(j)./yeta2;-i.*yeta2.*sin(derta2(j),cos(derta2(j);,%b介质特征矩阵,Mab=Ma*Mb;,%ab单元特征矩阵,N=10;,%周期数,M=MabN;,%一维光子晶体特征矩阵特征矩阵,Rfan(j)=abs(M(1,1)*(yeta0)+M(1,2)*(yeta0)*(yeta0)-M(2,1)-M(2,2)*(yeta0)./(M(1,1)*(yeta0)+M(1,2)*(yeta0)*sqrt(yeta0)+M(2,1)+M(2,2)*(yeta0)2;,%反射率,end,figure(1),holdon,p

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论