matlab对各种光纤光栅的仿真

matlab对各种光纤光栅的特性参数的仿真代码。

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）+（i*dc/sm）*sinh（sm*l）]*exp（-i*dc*l）;
    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）;
    end
    if t==1
    B（::t）=[T1 T2;T3 T4];
    else
        B（::t）=B（::t-1）*[T1T2;T3T4];
    end
end
A=B（::500）;
%C（x）=（abs（-A（21）*A（12）/A（22）+A（11）））^2;%投射谱
C（x）=（abs（A（21）/A（22）））^2;%反射谱
end

x=1:2000;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1108  2010-08-30 20:50  matlab 仿真\相移布拉格光栅反射谱\xiangyi.m

     文件        822  2010-08-30 20:31  matlab 仿真\均匀布拉格光栅反射普、透射谱\main.m

     文件       1133  2010-08-31 09:52  matlab 仿真\切趾光栅反射谱\qiezhi.m

     文件      18297  2010-08-31 10:10  matlab 仿真\切趾光栅反射谱\qiezhi1.jpg

     文件      18284  2010-08-31 09:48  matlab 仿真\切趾光栅反射谱\weiqiezhi.jpg

     文件        909  2010-08-31 10:52  matlab 仿真\啁啾光栅反射谱\zhoujiuguangshan.m

     文件      19361  2010-08-31 10:58  matlab 仿真\啁啾光栅反射谱\zhoujiuguangshanp.jpg

     文件        896  2010-08-31 16:50  matlab 仿真\取样光栅反射谱\quyangguangshan.m

     文件      37464  2010-08-31 16:51  matlab 仿真\取样光栅反射谱\2.jpg

     文件      36149  2010-08-31 17:09  matlab 仿真\取样光栅反射谱\2queji.jpg

     文件      17771  2010-09-01 12:13  matlab 仿真\取样切趾光栅反射谱\36.jpg

     文件       1189  2010-09-01 12:12  matlab 仿真\取样切趾光栅反射谱\new.m

     目录          0  2010-08-30 20:55  matlab 仿真\相移布拉格光栅反射谱

     目录          0  2010-08-30 20:55  matlab 仿真\均匀布拉格光栅反射普、透射谱

     目录          0  2010-08-31 11:01  matlab 仿真\切趾光栅反射谱

     目录          0  2010-08-31 11:01  matlab 仿真\啁啾光栅反射谱

     目录          0  2010-08-31 15:47  matlab 仿真\取样光栅反射谱

     目录          0  2010-09-01 12:14  matlab 仿真\取样切趾光栅反射谱

     目录          0  2010-08-30 20:55  matlab 仿真

     文件       1139  2011-06-13 20:30  matlab+仿真\matlab 仿真\切趾光栅反射谱\qiezhi.m

     文件      18297  2010-08-31 10:10  matlab+仿真\matlab 仿真\切趾光栅反射谱\qiezhi1.jpg

     文件      18284  2010-08-31 09:48  matlab+仿真\matlab 仿真\切趾光栅反射谱\weiqiezhi.jpg

     文件      37464  2010-08-31 16:51  matlab+仿真\matlab 仿真\取样光栅反射谱\2.jpg

     文件      36149  2010-08-31 17:09  matlab+仿真\matlab 仿真\取样光栅反射谱\2queji.jpg

     文件        896  2010-08-31 16:50  matlab+仿真\matlab 仿真\取样光栅反射谱\quyangguangshan.m

     文件      17771  2010-09-01 12:13  matlab+仿真\matlab 仿真\取样切趾光栅反射谱\36.jpg

     文件       1189  2010-09-01 12:12  matlab+仿真\matlab 仿真\取样切趾光栅反射谱\new.m

     文件        909  2010-08-31 10:52  matlab+仿真\matlab 仿真\啁啾光栅反射谱\zhoujiuguangshan.m

     文件      19361  2010-08-31 10:58  matlab+仿真\matlab 仿真\啁啾光栅反射谱\zhoujiuguangshanp.jpg

     文件        822  2010-08-30 20:31  matlab+仿真\matlab 仿真\均匀布拉格光栅反射普、透射谱\main.m

............此处省略12个文件信息