trig1a

Contents

define optical system
trace ray through edge of pupil
ray fan
fit wavefront
rms deviation from fit

define optical system

cv = [ 0 0.034654 -0.005408 0];
th = [ 0 1 49.42301 0];
rn = [ 1 1.5 1 1 ];

trace ray through edge of pupil

y =  linspace(0,1,6)';
ye = trig_ray(y*6.25, 0, cv, th, rn);
[y ye]

ans =

            0            0
          0.2  -0.00083418
          0.4   -0.0067075
          0.6     -0.02283
          0.8    -0.054764
            1     -0.10864

ray fan

y = linspace(-1,1,101)';
ye = trig_ray(y*6.25, 0, cv, th, rn);
plot(y,ye);
xlabel('relative aperture');
ylabel('transverse ray error');

fit wavefront

lambda = 0.00055;
ua = -0.125;
epsilon = abs(lambda/ua);

d = - ye/epsilon;
c = [ 2*y 4*y.^3 6*y.^5 ];

a = c\d

rms deviation from fit

diff = d - c*a;
rms = sqrt(sum(diff.*diff)/length(diff))

rms =

   0.00051646