FBG反射谱透射率的MATlab仿真，基于耦合模理论

L=0.01;
neff=1.447;
c=1*1e-7/（2*neff）;
deltaneff=0.0004;
N=100;
M=2000;
lamda1=1548;
lamda2=1552;
lamda=linspace（lamda1*1e-9lamda2*1e-9M）;
deltalamda=（lamda2*1e-9-lamda1*1e-9）/（M-1）;
for k=1:M
    F=[10;01];
    for n=1:N

deltaneff_z=deltaneff;
  %g=4;
   % deltaneff_z=deltaneff*exp（-g*g*（（2*n-1-N）/（2*N））^2）;
  %  H=0.5;
  %deltaneff_z=deltaneff*（1+H*cos（pi*（2*n-1-N）/N））/（1+H）;
    lamdaD=1550*1e-9-L*neff*c+2*neff*c*（L*n/N-0.5*L/N）;
    sigma=2*pi*neff*（1/lamda（k）-1/lamdaD）+2*pi*deltaneff_z/lamda（k）+（4*pi*neff）*c*（2*（neff+deltaneff））*（-L/2+n*L/N）/lamdaD^2;
    kac=pi*deltaneff_z/lamda（k）;
    RB=sqrt（kac^2-sigma^2）+eps;
    F=F*[cosh（RB*L/N）-（sigma/RB）*sinh（RB*