sesan10.m

calculating fiber eigen modes

%求解两层、3阶第1个模式的色散曲线   @@@(程序4）

%运用割线法

clear all;

l=3;
wl1=500e-9;
wlN=746.37492e-9;
deltawl=0.5e-9;
wl=wl1:deltawl:wlN;

N=1+fix((wlN-wl1)./deltawl);
neff=zeros(1,N);
neff1=1.45749154001038;
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,'m')
grid on;

xlabel('归一化频率')
ylabel('归一化有效折射率')


hold on;