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

if(abs(aTM(j) - 2*dTE(i)) < 1e-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;

for i = 1:2000

NTEe = linspace(1.5100, NTMx(i), 4000);

aTE = xTE_DispersionFun(NTEe, NTMx(i), m);

for j = 1:4000

if(abs(aTE(j) - 2*dTM(i)) < 1e-2)

V(k) = 2*dTM(i)*sqrt(1.5360^2 - 1.5100^2);

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

k = k+1;

end;

end;

end;

plot(V,b,'b');

hold on;

pause;