90 Optical Fiber Communication Systems with MATLAB
®
and Simulink
®
Models
G1=0.7028554; % Sellmeier’s coefficients for germanium
G2=0.4146307;
% doped silica (concentration B in table)
G3=0.8974540;
lambda1=0.0727723e-6;
%Wavelengths for germanium doped silica
lambda2=0.1143085e-6;
lambda3=9.896161e-6;
c = 299792458;
%Speed of light
π
= 3.1415926;
%Greek letter
π
a = 2.4e-6;
%Core radius
delta = 0.0043;
%Greek letter delta (ref. index difference
%between core and cladding)
% Calculating the refractive index
% -----------------------------------
npow2oflambda = 1 + (G1.*lambda.^2./(lambda.^2.-lambda1*lambda1)) ...
+ (G2.*lambda.^2./(lambda.^2.-lambda2*lambda2)) ...
+ (G3.*lambda.^2./(lambda.^2.-lambda3*lambda3));
noflambda=sqrt(npow2oflambda);
pointer=find(lambda==1.550e-6); ...