【光学】基于matlab单相移点光纤光栅仿真【含Matlab源码 2652期】.zip

clear all; clc %单相移点光纤光栅 lamda=1549.6*1e-9:0.00020*1e-9:1550.6*1e-9; % 波长范围 avNeff=0.00018; % 平均有效折射率变化 Neff=1.45; lamda_B=1550*1e-9; %反射波长 v=1; L1=0.0052; L2=0.0072; M=5001; % M is sampled FBG number for o=1:M p1=pi; %相位大小 P111=exp(-1i*p1/2); % phase delay matrix P122=exp(1i*p1/2); P112=0; P121=0; matrixP1=[P111,P112;P121,P122]; k=pi*v*avNeff/lamda(o); sigma=2*pi*Neff*(1/lamda(o)-1/lamda_B)+2*pi*avNeff/lamda(o); s=sqrt(k.^2-sigma.^2); F111=cosh(s*L1)-1i*sigma./s.*sinh(s*L1); % one section FBG matrix F122=cosh(s*L1)+1i*sigma./s.*sinh(s*L1); F112=-1i*k./s.*sinh(s*L1); F121=-F112; matrixF1=[F111,F112;F121,F122]; F11=cosh(s*L2)-1i*sigma./s.*sinh(s*L2); % one section FBG matrix F22=cosh(s*L2)+1i*sigma./s.*sinh(s*L2); F12=-1i*k./s.*sinh(s*L2); F21=-F12; matrixF=[F11,F12;F21,F22]; matrixT1=matrixF; matrixT2=matrixF1;%*matrixP1*matrixF;%*matrixP2*matrixF1; %R(o)=(abs(matrixF(2)/matrixF(1))).^2; matrixT=matrixF1*matrixP1*matrixF; r21(o)=matrixT1(3)/matrixT1(1); r22(o)=matrixT2(2)/matrixT2(1); r2(o)=matrixT(2)/matrixT(1); R21(o)=abs(r21(o)).^2; R22(o)=abs(r21(o)).^2; R2(o)=abs(r2(o)).^2; end q=matrixT; plot(lamda/1e-9,1-R2); xlabel('Wavelength （nm）') ylabel('透射率'); figure(3); plot(lamda/1e-9,10*log10(1-R2)); xlabel('Wavelength （nm）') ylabel('透射率/DB'); figure(2); plot(lamda/1e-9,R2) xlabel('Wavelength （nm）');ylabel('反射率'); figure(4); plot(lamda/1e-9,10*log10(R2)); xlabel('Wavelength （nm）') ylabel('反射率/DB');