%计算 transmittance DOI:10.1364/OE.26.008634
clear;
clc;
%基本参数的设置
c = 3e8;%真空中的光速
e0=1;%空气介电常数与磁导率 入射介质
u0=1;
e3=1;%出射介质
u3=1;
n1=1.46;%Si折射率
n2=2.82;%SiO2折射率
e1=n1^2;
e2=n2^2;
u1=1;
u2=1;
%PhC1
d1=680e-9;%Si厚度
d2=815e-9;%SiO2厚度
%PhC2
d3=658e-9;%Si厚度
d4=1290e-9;%SiO2厚度
d5=1045e-9;%SiO2缺陷层
theta=0;%入射角
angle=theta*pi/180;
for j=1:5000
lambda(1)=1300e-9;
lambda(j+1)=lambda(j)+1200e-9/5000;
w(j)=2*pi*c/lambda(j);%w=2pi*c/lambda啊！！！！
k0(j)=w(j)/c;
%********PhC1********si/sio2/si/sio2/si/sio2/si/sio2(AB)4
%********Si********
q1(j)=sqrt(e1)/sqrt(u1)*sqrt(1-sin(angle)^2/(e1*u1));
pz1(j)=k0(j)*sqrt(e1)*sqrt(u1)*sqrt(1-sin(angle)^2/(e1*u1));
Ma = zeros(2,2);
Ma(1,1) = cos(pz1(j)*d1);
Ma(1,2) =-1i*sin(pz1(j)*d1)/q1(j);
Ma(2,1) =-1i*q1(j)*sin(pz1(j)*d1);
Ma(2,2) = cos(pz1(j)*d1);
%********SiO2********
pz2(j)=k0(j)*sqrt(e2)*sqrt(u2)*sqrt(1-sin(angle)^2/(e2*u2));
q2(j)=sqrt(e2)/sqrt(u2)*sqrt(1-sin(angle)^2/(e2*u2));
Mb = zeros(2,2);
Mb(1,1) = cos(pz2(j)*d2);
Mb(1,2) = -1i*sin(pz2(j)*d2)/q2(j);
Mb(2,1) = -1i*q2(j)*sin(pz2(j)*d2);
Mb(2,2) = cos(pz2(j)*d2);
%********PhC2********(dc)4
%********Si********
pz3(j)=k0(j)*sqrt(e1)*sqrt(u1)*sqrt(1-sin(angle)^2/(e1*u1));
q3(j)=sqrt(e1)/sqrt(u1)*sqrt(1-sin(angle)^2/(e1*u1));
Mc = zeros(2,2);
Mc(1,1) = cos(pz3(j)*d3);
Mc(1,2) = -1i*sin(pz3(j)*d3)/q3(j);
Mc(2,1) = -1i*q3(j)*sin(pz3(j)*d3);
Mc(2,2) = cos(pz3(j)*d3);
%********SiO2********
pz4(j)=k0(j)*sqrt(e2)*sqrt(u2)*sqrt(1-sin(angle)^2/(e2*u2));
q4(j)=sqrt(e2)/sqrt(u2)*sqrt(1-sin(angle)^2/(e2*u2));
Md = zeros(2,2);
Md(1,1) = cos(pz4(j)*d4);
Md(1,2) = -1i*sin(pz4(j)*d4)/q4(j);
Md(2,1) = -1i*q4(j)*sin(pz4(j)*d4);
Md(2,2) = cos(pz4(j)*d4);
%********SiO2 defect********
pz5(j)=k0(j)*sqrt(e2)*sqrt(u2)*sqrt(1-sin(angle)^2/(e2*u2));
q5(j)=sqrt(e2)/sqrt(u2)*sqrt(1-sin(angle)^2/(e2*u2));
Me = zeros(2,2);
Me(1,1) = cos(pz5(j)*d5);
Me(1,2) = -1i*sin(pz5(j)*d5)/q5(j);
Me(2,1) = -1i*q5(j)*sin(pz5(j)*d5);
Me(2,2) = cos(pz5(j)*d5);
% X=(Ma*Mb)^4;%PhC1
% X=(Md*Mc)^4;%PhC2
X=(Ma*Mb)^4*(Md*Mc)^4;
% X=(Ma*Mb)^4*(Md*Mc)^4*Me*(Ma*Mb)^4;
%e0 e3空气 e1 Si e2 SiO2
p0(j)=sqrt(e1)/sqrt(u1)*cos(angle);
pl(j)=sqrt(e1)/sqrt(u1)*sqrt(1-e0*u0*sin(angle)^2/(e1*u1));
CC(j)=(X(1,1)+X(1,2)*pl(j))*p0(j)+(X(2,1)+X(2,2)*pl(j));
BB(j)=(X(1,1)+X(1,2)*pl(j))*p0(j)-(X(2,1)+X(2,2)*pl(j));
r(j) =BB(j)/CC(j);
t(j) =2*p0(j)/CC(j);
R(j)=abs(r(j)^2);
T(j)=pl(j)/p0(j)*abs(t(j)^2);
A(j)=1-T(j)-R(j);
end
plot(lambda(1:5000),T(1:5000),'k');
hold on
评论0