啁啾光纤光栅

%啁啾光纤光栅模拟%
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-4）*（Q（i-1）-Q（i））/deltala