Matlab模拟布拉格光纤光栅反射谱.zip

function [f]

=transmit_fiber(L,kappa,num,n_eff,lambda,lambda_Brag,deltaneff,b)

%-------------±äÁ¿¶¨Òå-------------------------------

% n_eff ----------------¹âÕ¤ÓÐЧÕÛÉäÂÊ

% L ----------------¹âÕ¤³¤¶È

% lambda_Brag ----------¹âÕ¤ÖÐÐIJ¨³¤

% kappa ----------------¹â²¨Ä£Ê½»¥ñîºÏϵÊý

*(1./lambda(num)-1./lambda_Brag)+(2*pi*deltaneff)/lambda(num);

s(num)=sqrt(kappa.^2-delta(num).^2);

kappa=(pi*b*deltaneff)/lambda(num);

s11(num)=cosh(s(num)*L)-i*(delta(num)/s(num))*sinh(s(num)*L);

s12(num)=-i*(kappa/s(num))*sinh(s(num)*L);

s21(num)=i*(kappa/s(num))*sinh(s(num)*L);

s22(num)=cosh(s(num)*L)+i*(delta(num)/s(num))*sinh(s(num)*L);

f=[s11(num) s12(num);s21(num) s22(num)];

end

n_eff=1.456;

F=[1 0;0 1];

for num=1:500

f=transmit_fiber(L,kappa,num,n_eff,lambda,lambda_Brag,deltaneff,b);

f=f*F;

r3(num)=f(2,1)/f(1,1);

R3(num)=(abs(-r3(num)))^2;%·´ÉäÂÊ

end

plot(lambda*1e9,R3,'b:');

title('FBG·´ÉäÆ×');