lecal.asv

双光纤耦合器的建模仿真

%数值积分函数lecal：le 计算式
function y=lecal(z)
h0=5.0e-6;%z=0处两光纤纤芯距，工艺参数
R=2;%光纤曲率半径
Y=3.0e-6;%Y为2根光纤垂直纸面的位移,工艺参数
W=0.5851;
a=5.0e-6; %光纤线芯半径
hz=((h0+(z.^2)./R).^2+Y.^2).^(0.5);
y=(h0./hz).^(0.5).*exp(-1*W.*(hz-h0)./a);