%啁啾光纖光柵模擬%
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