drdz.m

光纤光栅的反射谱 包含两段代码 在MATLAB中仿真

function drdz=drdz(z,r) 

global delta Neff delta_neff_average tlamda Pg Pe epsrate

drdz=-i*2.0*r*(delta+2*pi*delta_neff_average./tlamda-0.5*(-2.0*pi/Pg)*(1.0-1.0/(1+epsrate*(1-Pe)*z)^2))-i*(1+r*r)*pi*delta_neff_average./tlamda;