切趾光纤光栅仿真

neff=3.5;
k0=500;
L=1e-2;
pe=221e-9; 
wl=1546e-9;
l=L/500; 
z=zeros（1501）;
B=zeros（22500）;
C=zeros（1100）;
z=0:2e-5:1e-2; 
for x=1:2000
    wl=wl+0.001e-9;
    cc=2*pi*neff/wl;
    dc=cc-pi/pe;
    sm=sqrt（（k^2-dc^2））;
for t=1:500
    k=k0*exp（-8*（z（t）/L）^2）;
    if t==150
    T1=[cosh（sm*l）+（i*dc/sm）*sinh（sm*l）]*exp（-i*dc*l）*exp（i*pi）;
    T2=i*（k/sm）*sinh（sm*l）*exp（-i*dc*（z（t）+z（t+1）））;
    T3=-（i*k/sm）*sinh（sm*l）*exp（i*dc*（z（t）+z（t+1）））;
    T4=[cosh（sm*l）-（i*dc/sm）*sinh（sm*l）]*exp（i*dc*l）*exp（-i*pi）;
    else
    T1=[cosh（sm*l）