FBG.rar_FBG反射谱_光纤FBG_布拉格光栅_布拉格反射_矩阵传输法FBG

clear; clc; close all; L=0.015; Neff=1.460; C=0; v=1; Period=530.08*1e-9; deltaN=0.00015; SectionN=100; D=125e-6;%光栅直径125um E=74.52;%杨思模量 V=0.17;%泊松比 P11=0.121;%光弹系数 P12=0.270;%光弹系数 Lamda=[1547:0.01:1549]*1e-9;%计算波长范围 ii=0; r=0; for Lamda=[1547:0.01:1549]*1e-9 F=[1,0;0,1];%初始化传输矩阵 LamdaD=2*Neff*Period; for i=1:SectionN;%计算相应波长下的传输矩阵 delta=2*pi*Neff*(1/Lamda-1/LamdaD);%归一化频率失谐量 sigma=2*pi*deltaN/Lamda;%直流耦合系数 sigma1=delta+sigma;%平均自耦合系数 Kac=pi*V*deltaN/Lamda;%交流耦合系数 RB=(Kac^2-sigma1^2)^0.5; F=F*[cosh(RB*L/SectionN)-j*(sigma1/RB)*sinh(RB*L/SectionN),... -j*(Kac/RB)*sinh(RB*L/SectionN);j*(Kac/RB)*sinh(RB*L/SectionN),... cosh(RB*L/SectionN)+j*(sigma1/RB)*sinh(RB*L/SectionN)];% 传输矩阵输矩阵 end r=F(3)/F(1); ii=ii+1; R(ii)=(abs(r))^2; r=0; end Lamda=[1547:0.01:1549]*1e-9; plot(Lamda,R,'b');