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