Fbg1 n=1.45; T=1.55e-6/(2*n); L=2e-3; DeltaN=5e-4; wavelength=1.545e-6:0.01e-9:1.555e-6; % 起始波长 间隔 % 周期,中心波长 1550nm %一般在 1e-5~1e-3 % 有效折射率 % 光栅长度 k=pi*DeltaN./wavelength; % 交流耦合系数，2*pi/wavelength 波数即频率% 所有向量的乘 除要用点乘（.*）、点除（./） Beta=2*pi*n./wavelength; delta=Beta-pi/T; g=sqrt(k.^2-delta.^2); R=sinh(g*L).^2./(cosh(g*L).^2-(delta./k).^2); % 模式间的失谐量 % 传输常数 %Beta=2*n*pi./wavelength-pi/T g=abs(sqrt(k.^2-Beta.^2)); R=k.^2.*(sinh(g.*L)).^2./((Beta).^2.*(sinh(g.*L)).^2+g.^2.*(cosh(g.*L)).^2); plot(wavelength); ylabel('R'); xlabel('wavelength') hold on fbg2 n = 1.5; T = ((1*1.55)/(2*n))*10^5261(-6); L = 20*10^(-3); n1 = 5*10^(-5); wave0 = 1.5498*10^(-6); wave = wave0:0.001*10^(-8):1.5502*10^(-6); % 间隔太大，取小点 w = (pi*n1)./wave; beita3 = 2*n*pi./wave-pi/T; 点除（./） S = abs(sqrt(w.^2-(beita3).^2)); r = w.^2.*(sinh(S.*L)).^2./((beita3).^2.*(sinh(S.*L)).^2+S.^2.*(cosh(S.*L)).^2); plot(wave,r,'b-'); axis([1.5498*10^(-6) 1.5502*10^(-6) 0 1]); % 所有向量的开方（.^） % 所有向量的乘 4102 除要用点乘（1653.*）、 fbg3 function sampling_opticfiber(delta_n,T,L_all,Period) %---------------- sampling_opticfiber(3e-5,0.5,5e-2,1e-3)------------------------------- % M-----------------取样个数 % delta_n-----------调制折射率变化 % T-----------------取样光栅占空比 

% L_all-------------取样光栅总长度 % L_grating---------取样光栅节点长度 % Period------------取样光栅的取样周期%------------------------------------------------ disp('调制折射率变化为：');disp(delta_n); disp('取样光栅占空比：');disp(T); disp('取样光栅总长度:');disp(L_all); disp('取样光栅周期:');disp(Period); lambda=1e-9*linspace(1548,1552,1024); n_eff=1.45; L_grating=Period*T; M=L_all/Period; disp('取样个数：');disp(M); lambda_Brag=1550e-9; L_normal=Period-L_grating; for j=1:1024 F1=Optic_Fiber(lambda,lambda_Brag,delta_n,n_eff,j,L_grating); phai=Normal_OpticFiber(L_normal,n_eff,j,lambda); F1=(phai*F1*phai)^M; R(j)=(abs(-F1(2,1)/F1(1,1)))^2; end plot(lambda*1e9,R); axis([1548 1552 0 1]); title('The Reflective of Sampling Grating Fiber'); xlabel('Wavelength /nm'); ylabel('Reflective'); axis on; hold on end %-----------------取样光栅正常段--------------------------- function [phai]=Normal_OpticFiber(L_normal,n_eff,j,lambda) beita(j)=2*pi/lambda(j)*n_eff; phai=[exp(-i*beita(j)*L_normal),0;0,exp(i*beita(j)*L_normal)]; end %-----------------------取样光栅均匀段------------------------ function [F1]=Optic_Fiber(lambda,lambda_Brag,delta_n,n_eff,j,L_grating) delta=1/2 *2*pi*n_eff*(1./lambda-1/lambda_Brag); k=pi./lambda*delta_n; sigma1=2*pi./lambda*delta_n; sigma2=delta+sigma1; s=sqrt(k.^2-delta.^2); f11(j)=(cosh(s(j)*L_grating)-i*sigma2(j)/s(j)*sinh(s(j)*L_grating)); f12(j)=-(i*k(j)/s(j)*sinh(s(j)*L_grating)); f21(j)=(i*k(j)/s(j)*sinh(s(j)*L_grating)); f22(j)=(cosh(s(j)*L_grating)+i*sigma2(j)/s(j)*sinh(s(j)*L_grating)); 

F1=[f11(j) f12(j);f21(j) f22(j)]; end