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