sampledgrating.asv

matlab例程

lamda=1.50:0.0001:1.60;
i=0;
dim=length(lamda);
r=ones(1,dim);    
%k1=ones(1,dim);
%k2=ones(1,dim);
lamdaB=1.55;    
Kapa=0.01;
Neff=3.21;
dataN12=Kapa*lamdaB/2;
period=8;
gburstl=5;
ls=60;
gratingp=lamdaB/Neff/4;
t=gburstl/gratingp*2;

hom=ls-gburstl;
for lam=1.50:0.0001:1.60;
i=i+1;        


n1=Neff+dataN12/2;
n2=Neff-dataN12/2;
k1=(2*pi*(Neff+dataN12/2)/lam)+j*0.001;
k2=(2*pi*(Neff-dataN12/2)/lam)+j*0.001;
k=(2*pi*(Neff)/lam)
matrixhom=[exp(-j*k*hom),0;0,exp(j*k*hom)];

Tgrating1=1/n1*[Neff,dataN12;dataN12,Neff]*[exp(j*k1*gratingp),0;0,exp(-j*k1*gratingp)]*1/n2*[Neff,-dataN12;-dataN12,Neff]*[exp(j*k2*gratingp),0;0,exp(-j*k2*gratingp)];
matrixbur=Tgrating1^t;
matrixsam=matrixbur*matrixhom;
r(i)=Tgrating(2)/Tgrating(4);
end
plot(lamda,abs(r));