matlab对光栅的仿真代码

matlab对三种光栅：长周期光栅、均匀光栅等的仿真代码

clc
clear;
np=0.5308e-6;       %光栅周期%
L=80000*np;        %光栅长度2mm%
n1=1.46;          %纤芯原始折射率%
v=1;              %边缘可见度%
DeltaN=1e-4;
Neff=n1;
wl_D=2*Neff*np

h=50;
SectionN=2*h+1;

Ch=1e-9;
z=L/（2*h+1）;

a=pi/4;
T=[exp（-i*a）0;0exp（i*a）];       %相移参数

X=[1;0];

wavelength=[1545:0.01:1554.99]*1e-9;      % 初始化输入波长区间%


for j=1:1000
    F=zeros（22）;
    M=[10;01];                           %初始化矩阵
    k=pi*v*DeltaN/wavelength（j）;    
    for n=-h:h
        Lamda=（wl_D+Ch*n/SectionN）;%每一section是均匀的，对应一个lamda.
        delta=2*pi*DeltaN/wavelength（j）+2*pi*Neff*（1/wavelength（j）-1/Lamda）;
        g=sqrt（k^2-delta^2）;
        A=cosh（g*z）-i*delta*sinh（g*z）/g;
        B=-i*k*sinh（g*z）/g;
        C=cosh（g*z）+i*delta*sinh（g*z）/g;
        F=[AB;-BC];
        M=F*M;
    end
    Y=M*X;
    R（j）=abs（Y（21）/Y（11））^2;
end

plot（wavelengthR）;
axis（[1.54800e-6 1.5520e-6 0 1]）;
title（‘反射率‘）;
xlabel（‘wavelength‘）;
ylabel（‘Reflectivity‘）;
grid on;   

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

     文件       1049  2009-11-01 15:39  光栅MATLAB程序\FBG\FBG\CFBG\CFBG.asv

     文件       1049  2009-11-01 15:39  光栅MATLAB程序\FBG\FBG\CFBG\CFBG.m

     文件       1138  2009-10-10 11:02  光栅MATLAB程序\FBG\FBG\CFBG\CFBG_delay.asv

     文件       1165  2009-10-10 11:05  光栅MATLAB程序\FBG\FBG\CFBG\CFBG_delay.m

     文件       1087  2009-10-29 12:29  光栅MATLAB程序\FBG\FBG\CFBG\CFBG_EquivalentMirror.asv

     文件       1102  2009-10-29 20:39  光栅MATLAB程序\FBG\FBG\CFBG\CFBG_EquivalentMirror.m

     文件        836  2009-10-27 15:47  光栅MATLAB程序\FBG\FBG\CFBG\CFBG_VI.m

     文件        197  2009-09-26 16:04  光栅MATLAB程序\FBG\FBG\CFBG\func_ODE3.m

     文件        586  2009-07-09 14:23  光栅MATLAB程序\FBG\FBG\CFBG\piece.m

     文件        744  2009-09-03 19:53  光栅MATLAB程序\FBG\FBG\CFBG\VI_matrix0.m

     文件        661  2009-07-15 15:35  光栅MATLAB程序\FBG\FBG\CFBG\VI_matrix1.m

     文件        665  2009-04-07 13:07  光栅MATLAB程序\FBG\FBG\LFBG\Bar.m

     文件        862  2009-05-22 11:36  光栅MATLAB程序\FBG\FBG\LFBG\Cross.m

     文件        946  2009-06-05 18:41  光栅MATLAB程序\FBG\FBG\LFBG\Cross_ODE.m

     文件        112  2009-06-05 18:35  光栅MATLAB程序\FBG\FBG\LFBG\func2.m

     文件        879  2009-09-09 15:40  光栅MATLAB程序\FBG\FBG\UFBG\FBGODE.m

     文件        117  2009-05-27 16:42  光栅MATLAB程序\FBG\FBG\UFBG\func.m

     文件        971  2009-10-31 16:13  光栅MATLAB程序\FBG\FBG\UFBG\UFBG.asv

     文件        917  2009-10-31 16:18  光栅MATLAB程序\FBG\FBG\UFBG\UFBG.m

     文件        944  2009-10-31 16:47  光栅MATLAB程序\FBG\FBG\UFBG\UFBG_Delay.m

     文件        998  2009-10-31 17:17  光栅MATLAB程序\FBG\FBG\UFBG\UFBG_Delay1.asv

     文件        997  2009-10-31 17:17  光栅MATLAB程序\FBG\FBG\UFBG\UFBG_Delay1.m

     文件       1611  2009-06-16 13:48  光栅MATLAB程序\FBG\FBG\UFBG\UFBG_Delay_matrix.m

     文件     168327  2009-10-31 16:17  光栅MATLAB程序\FBG\FBG\UFBG\UFBG_L.fig

     文件       1091  2010-09-03 23:32  光栅MATLAB程序\FBG\FBG\UFBG\UFBG_matrix.m

     文件        743  2009-10-31 15:38  光栅MATLAB程序\FBG\FBG\UFBG\UFBG_ODE.asv

     文件        743  2010-07-04 22:21  光栅MATLAB程序\FBG\FBG\UFBG\UFBG_ODE.m

     文件        801  2009-10-30 17:06  光栅MATLAB程序\FBG\FBG\UFBG\UFBG_Symbolic.m

     文件        942  2009-10-27 16:00  光栅MATLAB程序\FBG\FBG\UFBG\UFBG_transmission_Delay.m

     文件       1136  2009-09-03 20:02  光栅MATLAB程序\FBG\FBG\UFBG\UFBG_VI.m

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