function PhaseFiber_by_TransmissionMatrix(Fai)
lamda=1e-9*linspace(1540,1560,1000);
R=Transmission_FBG2(Fai);
subplot(2,1,1)
plot(lamda*1e9,R);
title('相移光栅的反射谱线');
xlabel('波长 /nm');
ylabel('反射率');
grid on
subplot(2,1,2)
plot(lamda*1e9,1-R)
title('相移光栅的透射谱线');
xlabel('波长 /nm');
ylabel('透射率');
end
function [F1]=Transmission_FBG1(n,lamda,lamda_B,v,dn,n_eff,j)
delta=2*n_eff*pi*(1./lamda-1./lamda_B);
k=pi*dn/lamda_B;
q=sqrt(k.^2-delta.^2);
L(1)=220e-6;
f11(j,1)=(cosh(q(j)*L(1))-i*delta(j)/q(j)*sinh(q(j)*L(1)));
f12(j,1)=-(i*k/q(j)*sinh(q(j)*L(1)));
f21(j,1)=(i*k/q(j)*sinh(q(j)*L(1)));
f22(j,1)=(cosh(q(j)*L(1))+i*delta(j)/q(j)*sinh(q(j)*L(1)));
F1=