func.m

来自「光纤中纤芯基模的有效折射率计算程序,二分法求解！」· M 代码 · 共 12 行

function y=func(lamda,Neffco)
n2=1.450;       % 包层折射率
n1=1.458;       % 纤芯折射率
a1=2.625e-6;    % 纤芯半径
a2=6.25e-5;     % 包层半径
V=2*pi*a1*sqrt(n1^2-n2^2)./lamda;
b=(Neffco.^2-n2^2)/(n1^2-n2^2);
m=V.*sqrt(1-b);
n=V.*sqrt(b);
y=m.*besselj(1,m)./besselj(0,m)-n.*besselk(1,n)./besselk(0,n);