啁啾光纤光栅的反射谱和时延仿真

clear all
clc
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;
%      （2）高斯变迹
%       deltaneff = 0.00005*exp（（-64*（-L/2+i*L/N）^4）/L^4）;
        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