L = 10*1e-3;%光栅长度
neff = 1.447;%光栅有效折射率
v = 1;%光栅可见度
delta_neff = 2*10^-4;%光栅有效折射率变化
lambda_B = 1550e-9;%FBG中心波长
lambda = 1e-9*linspace(1548,1552,1000);%波长范围
delta = 2*pi*neff*(1./lambda-1/lambda_B); %布拉格波长谐量
sigma = 2*pi./lambda*delta_neff+delta; %直流自耦合系数
kappa = pi./lambda*v*delta_neff; %交流耦合系数
s = sqrt(kappa.^2-sigma.^2);
R = (kappa.*kappa.*(sinh(s.*L)).^2)./(s.^2.*(cosh(s.*L)).^2+sigma.^2.*(sinh(s.*L).^2)); %反射率
subplot(1,2,1) ; plot(lambda, R) ; title(' FBG反射谱' ) ;
xlabel(' 波长nm' ) ; ylabel(' 反射率' ) ; grid on;
L1 = 10*1e-3;%光栅长度
x1 = linspace(-L1, L1, 1000);
f1 = exp(-8*(x1/L1) .^2);
neff1 = 1.447;%光栅有效折射率
v1 = 1;%光栅可见度
delta_neff1 = 2*10^-4;%光栅有效折射率变化
delta_neff1 = 2*10^-4*f1;
lambda_B1 = 1550e-9;%FBG中心波长
lambda1 = 1e-9*linspace(1549,1551,1000);%波长范围
delta1 = 2*pi*neff1*(1./lambda1-1/lambda_B1); %布拉格波长谐量
sigma1 = 2*pi./lambda1.*delta_neff1+delta1; %直流自耦合系数
kappa1 = pi./lambda1.*v1.*delta_neff1; %交流耦合系数
kappa2 = kappa1.*f1;
s2 = sqrt(kappa2.^2-sigma1.^2);
R1 = (kappa2.*kappa2.*(sinh(s2.*L1)).^2)./(s2.^2.*(cosh(s2.*L1)).^2+sigma1.^2.*(sinh(s2.*L1).^2)); %反射率
subplot(1,2,2) ; plot(lambda1, R1, ' r' ) ;
title(' 高斯函数切趾(G=8) ' ) ;
