传输矩阵法计算薄膜的反射率和透射率

传输矩阵法，一般用来研究周期性分层介质中的传输特性。对于分层媒质的微分方程在特定条件下的解，最简单便捷的方式便是使用传输矩阵。传输矩阵法就是用矩阵的形式来描述电磁波在多层介质中的传播情况，在每一层介质中传播过程中的运动规律满足的Maxwell 方程组。

clear
clc
 pi=3.1415926;
 cc=2.99792458e8;             
 muz=4.0*pi*1.0e-7;          
 epsz=1.0/（cc*cc*muz）;      
 
 nm=1e-9;
       
eQ=1.60217733e-19;
h_b=1.05457266e-34;

eV_omegap=7.75;
eV_tao=0.08267;
omegaP=eV_omegap*eQ/h_b;
omegaV=eV_tao*eQ/h_b;
 
n0=1.0;
 ng=1.0;

mlayer=2;
 
 Rplot=[];
Romega=[];
Rapha=[];
lambda0=0.5e-6;
omega0=2*pi*cc/lambda0;
for omega=0*omega0:0.001*omega0:6*omega0
nSiO2L=1.378;
 epsSiO2L=nSiO2L^2;
nAgL=2.356;
  
epsAgL=nAgL^2;
 
dSiO2=lambda0/4/nSiO2L;
 
dAg=lambda0/4/nAgL;
 
n=[epsAgL^（1/2）epsSiO2L^（1/2）epsAgL^（1/2）epsSiO2L^（1/2）epsAgL^（1/2）...
     epsSiO2L^（1/2）epsAgL^（1/2）epsSiO2L^（1/2）epsAgL^（1/2）epsSiO2L^（1/2）...
     epsAgL^（1/2）epsSiO2L^（1/2）epsAgL^（1/2）...
     epsSiO2L^（1/2）epsAgL^（1/2）epsSiO2L^（1/