matlab-光学薄膜资料讲解.docxVIP

安徽工业大学光信息科学与技术专业 《光学软件课程设计》课程报告 年级： 光082 姓名： 指导教师： 黄仙山（2900） 日期： 2011年6月28日 报告正文: 一、设计题目 一维周期性光学薄膜结构中光学电磁传输特性的理论研究 二、设计理论 薄膜特征矩阵 在界面1 E1 Ei1 Er1 Et1 Er2 HiH ii cos iiH ri cos iiH ti cos i2 Hi H ii cos ii H ri cos ii H ti cos i2 H r2 cos i2 由《电磁场与电磁波》 H与E的关系H ez E 得: H 0则Hi 0 (Eii Eri) n0 cos ii -(Eti Er2)ni cos i2 在界面2 1 1 E2 Ei2 Er2 Et2 H2 :0(Ei2 Er2)n1cos i2 0Et2nGcos t2 在不考虑薄膜对光能的吸收时， Ei2和Et1， E「2和E「2的关系如下: Ei2 Et1exp(i 1) I Er2 Er2exp(i 1) 平面波通过薄膜一次 A,B两点的位相变化 —n 1h1 cos i2 E2 Et1 exp(i 1) E；2exp( i i) H2 [ Et1 exp(i 1) E；2exp( i 1] R| cos i2 Et1 E；2 exp( i J(E —(E2 exp(i )(E ^(E2 氏） 1 生） 2 o -n-i cos o i2 Et1 E2 COS 1 H1 E2 1i sin 战叫) 1 H 2 cos(i 1) 写成矩阵形式： E1 H1 COS 1 i 1 sin —riihi cos i .Sin 1 1 cos E2 H2 i2 （ （对S波） o -nj cos i 2 o n1 o-——1 (对P波) 0 COS i 2 当膜系包含 N层膜时，则有 Ei Hi M1M2...M n En 1 H N 1 整个膜系的特征矩阵M MiM2...Mn 膜系反射率的计算 令M 令M的矩阵元为 A,B,C,D 即 M A EN 1 H N 1 E r1 膜系反射系数 膜系反射系数2 膜系反射系数 2 o A 0 B 0 G C D G 反射率：R r 投射率：T t 由两种不同介电常数（A 由两种不同介电常数（A, b ）和厚度（dA,dB ）的电磁材料所组成的一维周期 性光学薄膜结构。利用传输矩阵计算光在介质薄膜。 A E A B A 0 A B A E 本文中所选取的色散材料，如 [11,12]，LiF （模型中的 A层）的折射率在 1.25 ~ 2.35 m范围内的关系为 nA 1.38761 2 0.001796 L 0.00041 L 0.0023045 0Si （模型中的B层）的折射率在同一波段为 nB 3.41696 0.138497 L0.013924 L nB 3.41696 0.138497 L 0.013924 L2 0.0000209 0.000000148 其中 L 1/ 其中 L 1/ 2 0.028。 在本课程设计中取 1.5 m，贝U na 1.3831, nb 3.4770 1.386率射折介质A折射率随波长变化关系1.3841.3821.381.3781.3761.4 1.61.8 1.386 率射折 介质A折射率随波长变化关系 1.384 1.382 1.38 1.378 1.376 1.4 1.6 1.8 2 波长d/ m 2.2 2.4 2.6 3.52介质B折射率随波长变化关系nbn率射折539438 7 6 3.52 介质B折射率随波长变化关系 nb n率射折 5 3 9 4 3 8 7 6 4 4 4 3 3 3 5 4 4 4 3 3 1.4 1.6 1.8 2 波长d/ m 2.2 2.4 2.6 三、计算程序 （运行环境 Matlab R2009a ） clear all; n0=1;na=1.3831;nb=3.4770;% 折射率 n=i nput（输入光学薄膜的周期数:）; theta仁input（ 输入入射角:）; da=input（输入介质A的实际厚度:）; db=input（输入介质B的实际厚度:）; for lambda=1250:2350; theta2=asin(n0*sin(theta1)/na);% 介质 A 中的折射角 theta3=asin(n0*sin(theta1)/nb);% 介质 B 中的折射角 deltaa=2*pi*na*da*cos(theta2)/lambda*1e9; deltab=2*pi*nb*db*cos(theta3)/lambda*1e9; %由于波长取值不是 SI 单位制，此处化为 SI 单