传输矩阵法计算薄膜的反射率和透射率_matlab迭代法求多层介质反射率, matlab介质膜传输矩阵透射特性-C/C文档类资源

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/2),epsAgL^(1/2)]; h=[dAg,dSiO2,dAg,dSiO2,dAg,dSiO2,dAg,dSiO2,dAg,dSiO2,dAg,dSiO2,dAg,dSiO2,dAg,dSiO2,dAg]; eta0=(epsz/muz)^(1/2)*n0; etag=(epsz/muz)^(1/2)*ng; delt=zeros(1,mlayer); eta=zeros(1,mlayer); MM=zeros(2,2,mlayer); Mtotal=[1 0;0 1]; for m=mlayer:-1:1 delt(m)=(omega/cc)*n(m)*h(m); eta(m)=(epsz/muz)^(1/2)*n(m); MM(:,:,m)=[cos(delt(m)),-i*sin(delt(m))/eta(m);-i*eta(m)*sin(delt(m)),cos(delt(m))]; Mtotal=MM(:,:,m)*Mtotal; end A=Mtotal(1,1); B=Mtotal(1,2); C=Mtotal(2,1); D=Mtotal(2,2); r=(A*eta0+B*eta0*etag-C-D*etag)/(A*eta0+B*eta0*etag+C+D*etag); R=(abs(r))^2; apha=atan(imag(r)/real(r)); Rapha=[Rapha,apha]; Rplot=[Rplot,R]; Romega=[Romega,omega]; end plot(Romega./omega0,Rplot,'g') xlabel('g');ylabel('R')