一维光子晶体MATLAB仿真代码吸收率折射率透射率参数可修改

clear all

tic
eps0=8.85e-12;
mu0=4*pi*1e-7;
ll=sqrt（eps0/mu0）;
c=3e+8;
lamd0=800e-9;
i=sqrt（-1）;
w0=2*pi*c/lamd0;

ct0=0*pi/180;

c0=cos（ct0）;
NN=5;

na=1.7;
a=142e-9;
nb=3;%（3:5.5）硅光电池
%nb=1.66;
b=200e-9;
%4*lamb
n0=2.86;
% k0=0.15;
k0=0;
nd=n0+k0*i;
d=500e-9;
%l=（a+b+d）*NN*2;
l=（a+b）*2*NN+d;

scta=sin（ct0）/na;
sctb=scta*na/nb;
sctd=scta*nb/nd;

ca=sqrt（1-scta*scta）;
cb=sqrt（1-sctb*sctb）;
cd=sqrt（1-sctd*sctd）;

deltaA=（ca*na*a/c）;
deltaB=（cb*nb*b/c）;
deltaD=（cd*nd*d/c）;

    etaA=（ll*na*ca*ca）;
    etaB=（ll*nb*cb*cb）;
    etaD=（ll*nd*cd*cd）;
%     etaA=（ll*na/ca/ca）;
%     etaB=（ll*nb/cb/cb）;
%     etaD=（ll*nd/cd/cd）;


n=1;nw（1）=nb;m=-1;k=0;
dx=0.0001;

shui=[];shuo=[];heni=[];heno=[];dw=[];ww0=[];
p=0;q=0;
for x=8/11:dx:8/3
    Ma=[cos（deltaA*x*w0） -i*sin（deltaA*x*w0）/etaA;-i*etaA*sin（deltaA*x*w0） cos（deltaA*x*w0）];  %% x=w/w0
    Mb=[cos（deltaB*x*w0） -i*sin（deltaB*x*w0）/etaB;-i*etaB*sin（deltaB*x*w0） cos（deltaB*x*w0）];
    Md=[cos（deltaD*x*w0） -i*sin（deltaD*x*w0）/etaD;-i*etaD*sin（deltaD*x*w0） cos（deltaD*x*w0）];
    
M=（Mb*Ma）^NN*Md*（Mb*Ma）^NN;

    t=2*ll*c0/（M（11）*ll*c0  +  M（12）*ll*ll*c0*c0  +  M（21）+  ll*M（22））;
    r=（M（11）*ll*c0  +  M（12）*ll*ll*c0  -  M（21）  -  ll*M（22））/（ll*M（11）*c0+  ll*ll*M（12）*c0+M（21）+ll*M（22））;

    %t=2*ll/（M（11）*ll+M（12）*ll*ll+M（21）+ll*M（22））;
    %r=（M（11）*ll+M（12）*ll*ll-M（21）-