glass1

Contents

load class catalog
match glass name
check dispersion coefficients
plot dispersion

load class catalog

clear all
load glass
disp(glass)

        name: {103x1 cell}
          nd: [103x1 double]
    disperse: [103x6 double]
          ne: [103x1 double]
          vd: [103x1 double]
          ve: [103x1 double]

match glass name

name = 'N-BK7';
idx = strmatch('N-BK7',glass.name);
fprintf('glass %s nd %.8g\n',name,glass.nd(idx));

glass N-BK7 nd 1.5168

check dispersion coefficients

wvl = [ 0.48613 0.58756 0.65627];
a = glass.disperse(idx,:);
rn = sellmeier(a,wvl);
vd = (rn(2)-1)/(rn(1)-rn(3));
fprintf('vd (from table) %g\n',glass.vd(idx));
fprintf('nd %.8g vd %g\n',rn(2),vd);

vd (from table) 64.17
nd 1.5168001 vd 64.1664

plot dispersion

w = linspace(0.4,0.7,100);
plot(w,sellmeier(a,w),'k','LineWidth',2);
set(gca,'FontSize',14);
grid;
xlabel('wavelength   \lambda (\mu m)')
ylabel('index','FontSize',14)
title('Refractive index vs \lambda for N-BK7 glass');