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;
????????%直流自耦合系數(shù)
????????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