sellmeier.m

来自「利用电磁场的源激发方法来计算光子晶体波导例如光子晶体光纤」· M 代码 · 共 35 行

function epsilon = sellmeier(lambda)
global C
lambda = lambda*1e6;
% Corning
% B1 = 0.68374049400;
% B2 = 0.42032361300;
% B3 = 0.58502748000;
% C1 = 0.00460352869;
% C2 = 0.01339688560;
% C3 = 64.49327320000;
% 
% epsilon = 1+lambda.^2.*(B1./(lambda.^2-C1) + B2./(lambda.^2-C2) + B3./(lambda.^2-C3));

% Agrawal
% b1 = 0.6961663;
% b2 = 0.4079426;
% b3 = 0.897494;
% w1 = 2*pi*C/0.0684043;
% w2 = 2*pi*C/0.1162414;
% w3 = 2*pi*C/9.896161;
% w = 2*pi*C./lambda;
% epsilon = 1+B1*w1^2./(w1.^2-w.^2) + B2*w2^2./(w2.^2-w.^2) + B3*w3^2./(w1.^2-w.^2);

% Dyott
a1 = 0.2045154578;
a2 = 0.06451676258;
a3 = 0.1311583151;
z1 = 0.06130807320;
z2 = 0.1108859848;
z3 = 8.964441861;
alpha = a1*lambda.^2./(lambda.^2-z1^2) + a2*lambda.^2./(lambda.^2-z2^2) + a3*lambda.^2./(lambda.^2-z3^2);
epsilon = (1+2*alpha)./(1-alpha);