基于matlab的光学薄膜特性分析,基于matlab的光学薄膜特性分析.doc

本文探讨了基于MATLAB的光学薄膜特性分析,包括单层和多层光学薄膜的反射和透射率计算。通过矩阵理论,详细阐述了如何利用MATLAB进行数值仿真,展示了不同入射角度和波长对薄膜反射率和透射率的影响。此外,还设计了一个用于模拟光学特性的图形用户界面,以直观地展示结果。MATLAB的运用简化了计算过程,为光学薄膜研究提供了便利的计算工具。
摘要由CSDN通过智能技术生成

目录

1 技术指标1

2 基本原理1

2.1 单层光学薄膜1

2.2 多层光学薄膜及矩阵理论2

3 基于matlab的光学薄膜特性分析4

4 结论9

5 参考文献9

6 附录:基于matlab的光学薄膜特性分析源程序代码10

16

武汉理工大学《能力拓展训练》报告书

基于matlab的光学薄膜特性分析

1 技术指标

对光学薄膜的理论进行分析,从电磁理论出发分析单层膜和多层膜的反射和透射率特性;并分析各参数对其光学性能的影响;运用Matlab对单层介质薄膜和多层介质薄膜的光学特性进行仿真;编写可用于一般多层膜的图形用户界面,并利用该程序对几种典型的薄膜系统进行验证。

2 基本原理

光射到任何一个光学零件的表面上都要发生反射和折射,不同的光学系统对反射率和透射率有不同的要求.对于透射成像的光学系统,反射不仅降低了透过光学系统的能量,而且造成杂散光,降低系统的成像清晰度.而对于一个有多个镜片的光学系统,其反射损失更为可观.为了降低反射损失,需要在镜片上镀光学薄膜,是每个光学面的反射损失都降到约等于百分之一,甚至更小。

薄膜光学发展至今,已经形成一套比较完整的实用理论。本文将利用多光束的干涉理论和矩阵理论,并运用matlab对薄膜光学系统的光学性质做简要分析。

2.1 单层光学薄膜

折射率为的光学基底上镀上一层厚度为,折射率为的透明薄膜,就构成典型的单层光学薄膜结构。薄膜可以看成一块平行薄板,光在这个薄板上下表面发射和透射的光叠加将会发生多光束干涉。利用已有的结果,可得反射光叠加后的合震动的复振幅

(2-1-1)

式中,是薄膜表面上的反射振幅比,是下表面的反射振幅比,是两相邻光束之间的相位差

(2-1-2)

而反射率

(2-1-3)

且透射率和反射率的关系:。

2.2 多层光学薄膜及矩阵理论

图(a) 两个界面的多膜系统

下面根据电磁理论分析:

在界面1

(2-2-1)

由《电磁场与电磁波》H与E的关系得:,

取则

(2-2-2)

在界面2

(2-2-3)

在不考虑薄膜对光能的吸收时,,的关系如下:

(2-2-4)

平面波通过薄膜一次A,B两点的位相变化。

(2-2-5)

(2-2-6)

(2-2-7)

写成矩阵形式:

(2-2-8)

当膜系包含N层膜时,则有

(2-2-9)

整个膜系的特征矩阵。

令M的矩阵元为A,B,C,D即,则

(2-2-10)

膜系反射系数

(2-2-11)

膜系透射系数

(2-2-12)

反射率:,透射率:。

由两种不同介电常数()和厚度()的电磁材料所组成的一维周期性光学薄膜结构。利用传输矩阵计算光在介质薄膜。

图(b)一维周期性光学薄膜结构

本文中所选取的色散材料,如[11,12],LiF(模型中的A层)的折射率在范围内的关系为:

Si(模型中的B层)的折射率在同一波段为:

其中。

在本课程设计中取,则3.4770。

3 基于matlab的光学薄膜特性分析

matlab是一种以矩阵运算为基础的交互式程序语言,广泛应用于自动控制、数学运算、信号分析、航天工业和绘图等各行各业.本文使用Matlab计算膜系特征矩阵,进行数值的运算仿真研究与验证。

仿真模拟的光学特性matlab程序见附录。

图(c) 光学特性分析的GUI交互界面

图(d) 单膜反射率和透射率与入射波长的关系(正入射)

图(e) 单膜反射率和透射率与入射波长的关系(斜入射)

图(f) 多膜反射率和透射率与入射波长的关系(正入射)

图(g) 单膜反射率和透射率与入射波长的关系(斜入射)

图中,纵坐标表示反射率和透射率,其中,红实线表示反射率,蓝虚线表示透射率,二者之和等于1。横坐标则表示入射光的波长(um)。

从图中分别对比分析知,对于一个介质折射率和厚度固定等参数固定的薄膜系统,其反射率或透射率达到最大值的时候所对应的这一固定波长称为中心波长,当入射光偏离中心波长的时候,其反射率必定会相应的下降,因此,每一种具体的膜系,只对一定波长范围的光波才能有反射率。以上各图反射率最大值(接近于1)时所对应的波长范围称为反射带的宽度。

图(h)8层膜(偶数)系反射率和透射率与波长的关系

图(i)9层膜(奇数)系反射率和透射率与波长的关系

从上面两个图对比可知,在同一波长情况下,奇数层膜的反射率要比偶数层

反射率略高。这是因为,膜系第一层的折射率高于第二层折射率,如果膜系为奇数层,那么最后一层的折射率也是高反射率层。层数为时,。

犹豫,故可以获得更高的反射率。反之,膜系为偶数层时,达不到增加反射率的目的。

4 结论

由上述的计算公式可知,光学薄膜的特性计算在引入矩阵理论后成为一个相对简单的问题,只要给定每一层膜系的参数,多层膜就可以利用矩阵的连乘方便地求出在每一特定入射角和波长下的透射率、反射率等参数.许多计算机软件的出现,能方便进行自动数值计算,如matlab软件的应用,不仅能计算并保存数据,而且更好地调用资源,作图、仿真模拟使得研究清晰明了,易于操作理解,为光学薄膜的研究提供了有利的计算工具.因此,计算机软件的应用促进了光学薄膜的发展,它与传统的理论研究相辅相成,是理论应用与实践的有力保障,使光学薄膜的发展迈上了新的台阶.光学薄膜已经在光学技术发展的历史中创造了光辉的成绩,如今广泛地渗透到各个新

function [R,T] = multilayer_reflectance_transmittance(n,d,k,theta_in,wavelength) % n: vector of refractive indices of each layer % d: vector of thicknesses of each layer (in nm) % k: vector of extinction coefficients of each layer % theta_in: incidence angle (in radians) % wavelength: wavelength of incident light (in nm) c = physconst('LightSpeed'); lambda = wavelength*1e-9; lambda_p = 2*pi*c./sqrt(eps0.*n.^2 - k.^2 - j*2*pi*eps0.*n.*k./(lambda/(2*pi*c)))); lambda_s = 2*pi*c./sqrt(eps0.*n.^2 - k.^2 + j*2*pi*eps0.*n.*k./(lambda/(2*pi*c))); theta_t = asin(n(1)*sin(theta_in)/n(end)); % transmission angle Z0 = 1/376.73031; for i = 1:length(n) [r_p(i),t_p(i)] = fresnel_am_tf(lambda_p(i),theta_in,n(i),n(i+1),k(i),k(i+1)); [r_s(i),t_s(i)] = fresnel_am_tf(lambda_s(i),theta_in,n(i),n(i+1),k(i),k(i+1)); Z(i) = n(i)*sqrt(eps0)/cos(theta_in); Gamma_p(i) = (Z(i) - n(i+1)*sqrt(eps0)*cos(theta_t))/(Z(i) + n(i+1)*sqrt(eps0)*cos(theta_t)); Gamma_s(i) = (n(i)*sqrt(eps0)*cos(theta_t) - n(i+1)*sqrt(eps0)*cos(theta_in))/(n(i)*sqrt(eps0)*cos(theta_t) + n(i+1)*sqrt(eps0)*cos(theta_in)); Delta_p(i) = exp(-2*j*lambda_p(i)*d(i)/cos(theta_in)); Delta_s(i) = exp(-2*j*lambda_s(i)*d(i)/cos(theta_in)); end r_pp = r_p(1); r_ss = r_s(1); t_pp = t_p(1); t_ss = t_s(1); for i = 2:length(n) r_pp = r_pp + (t_pp^2 * Gamma_p(i-1) * r_p(i) * Delta_p(i-1)); r_ss = r_ss + (t_ss^2 * Gamma_s(i-1) * r_s(i) * Delta_s(i-1)); t_pp = t_pp * t_p(i) * Delta_p(i-1) * (1 + Gamma_p(i-1) * r_p(i)); t_ss = t_ss * t_s(i) * Delta_s(i-1) * (1 + Gamma_s(i-1) * r_s(i)); end R_p = abs(r_pp)^2; R_s = abs(r_ss)^2; T_p = real((Z0 * cos(theta_t) / n(end)) * (abs(t_pp)^2)); T_s = real((Z0 * cos(theta_t) / n(end)) * (abs(t_ss)^2)); R = (R_p + R_s)/2; T = (T_p + T_s)/2; end function [r_p,t_p,r_s,t_s] = fresnel_am_tf(lambda,theta_in,n1,n2,k1,k2) % lambda: wavelength (in m) % theta_in: incidence angle (in radians) % n1, n2: refractive indices of incident and transmitted media, respectively % k1, k2: extinction coefficients of incident and transmitted media, respectively theta_t = asin(n1*sin(theta_in)/n2); % transmission angle Z0 = 1/376.73031; Gamma_p = (n1*cos(theta_in) - n2*cos(theta_t))/(n1*cos(theta_in) + n2*cos(theta_t)); Gamma_s = (n1*cos(theta_t) - n2*cos(theta_in))/(n1*cos(theta_t) + n2*cos(theta_in)); Delta_p = exp(-2*j*pi*n2*k2*lambda*cos(theta_t)/(lambda^2*n2^2*cos(theta_t)^2 + k2^2)); Delta_s = exp(-2*j*pi*n2*k2*lambda*cos(theta_t)/(lambda^2*n2^2*cos(theta_in)^2 + k2^2)); r_p = (Gamma_p + (1-Gamma_p)*Delta_p); r_s = (Gamma_s + (1-Gamma_s)*Delta_s); t_p = 2*n1*cos(theta_in)/(n1*cos(theta_in) + n2*cos(theta_t)) * (1-Gamma_p)*sqrt(Delta_p); t_s = 2*n1*cos(theta_in)/(n1*cos(theta_t) + n2*cos(theta_in)) * (1-Gamma_s)*sqrt(Delta_s); end
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值