sesan9.m

来自「calculating fiber eigen modes」· M 代码 · 共 40 行

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

%运用割线法

clear all;

l=1;
wl1=500e-9;
wlN=694.453188e-9;
deltawl=5e-9;
wl=wl1:deltawl:wlN;

N=1+fix((wlN-wl1)./deltawl);
neff=zeros(1,N);
neff1=1.45294294183074;
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('归一化有效折射率')


hold on;