sesantong.m

calculating fiber eigen modes

function ss=sesantong(wl1,wlN,neff1,l)
%求解两层色散曲线的通用程序 函数 
%运用割线法
%wl1为最小波长
%wlN为此模式的截止波长
%neff1为wl1时的有效折射率
%l为模式的阶数

deltawl=1e-9;
wl=wl1:deltawl:wlN;

N=1+fix((wlN-wl1)./deltawl);
neff=zeros(1,N);
neff(1)=neff1;



for i=2:1:N
   u=secant(neff(i-1),neff(i-1)-0.001,wl(i),l);
    neff(i)=u;
end
f=3.*10.^8./wl;

r1=4e-6;
nc=1.458;
n=1.45;
k=2.*pi./wl;
beta=k.*neff;
b=(beta./k-n)./(nc-n);%归一化
v=r1.*k.*sqrt(nc.^2-n.^2);
plot(v,b,'k')
grid on;

xlabel('归一化频率')
ylabel('归一化有效折射率')
axis([0,12,0,1]);