啁啾光纤光栅的光学特征matlab

L=0.04;neff=1.45;C=10*1e-9;N=50;M=1501;
deltaneff = 0.0001;
lamda1=1549;lamda2=1551;
lamda=linspace（lamda1lamda2M） * 1e-9;
deltalamda=（lamda2-lamda1）/M * 1e-9 ;
tic
for k=1:M
F=[10;01];
for i =1:N
% （1）均匀
deltaneff = 0.00005;

lamda_D =（1550 -C*L/2+C*i * L/N） * 1e-9;
sigma= 2 * pi * neff * （1/lamda（k） -1/lamda_D） + 2 * pi * deltaneff/lamda （k）+ （4 * pi * neff） * C * （ -L/2+i *L/N）/lamda_D^ 2;
kac = pi * deltaneff/lamda （k）;
RB= sqrt（kac^2 -sigma^2）;
F=F*[cosh（RB* L/N）-j*（sigma/RB）*sinh（RB*L/N）-j*（kac/RB）*sinh（RB*L/N）;j*（kac/RB）*sinh（RB*L/N）cosh（RB*L/N）+j*（sigma/RB）*sinh（RB*L/N）];
end
R（k） = （abs（ -F（3）/F（1）））^2;
Q（k）=phase（（ -F（3）/F（1）））;
end
toc
tao（1）=Q（1） ;tao（2） =Q（2） ;tao（3） =Q（3）;
for i=4:M
if （abs（Q（i -1）-Q（i））<=1）
tao（i） = （（lamda1+i * 0.001）^2 * 1e-18/（2 *pi * 3e-4） * （Q（i-1）-Q（i））/deltalamda）;
else
tao（i） = （（lamda1+i * 0.001）^ 2 * 1e-18/（2 * pi * 3e-4） * （Q（i-3） -Q（i-2））/deltalamda）;
end
end