在光波导理论中,求解波导色散曲线的常用数值方法之一是效折射率法,本文给出了有效折射率求解思路,并给出了具体的matlab程序,可供学习集成光学的学生参考使用。

global V b;

% a:d = 1:1

figure(1);

% x方向偏振

NTEx = linspace(1.5100, 1.5370, 2000);

for n = 0:1

dTE = yTE_DispersionFun(NTEx, n);

for m = 0:1

k = 1;

for i = 1:2000

if(NTEx(i) <= 1.5360)

NTMe = linspace(1.5100, NTEx(i), 4000);

aTM = xTM_DispersionFun(NTMe, NTEx(i), m);

for j = 1:4000

if(abs(aTM(j) - dTE(i)) < 2e-2) V(k) = 2*dTE(i)*sqrt(1.5360^2 -1.5100^2);

b(k) = (NTMe(j)^2 - 1.5100^2)/(1.5360^2 -1.5100^2);

k = k+1;

end;

end;

end;

end;

plot(V, b,'r');

hold on;

pause;

clear V b;

end;

end;

% y方向偏振

NTMx = linspace(1.5095, 1.5360, 2000);

for n = 0:1

dTM = yTM_DispersionFun(NTMx, n);

for m = 0:1

k=1;