trig2 demo
Contents
define optical system
cv = [ 0 0.034654 -0.005408 0];
th = [ 0 1 49.42301 0];
rn = [ 1 1.5 1 1 ];
trace oblique ray through edge of pupil
y = linspace(-1,1,11)';
Uchief = -4*pi/180;
Uchief
yo = trig_ray(0,Uchief,cv,th,rn);
yo
ye = trig_ray(y*6.25*cos(Uchief), Uchief, cv, th, rn) - yo;
format short g
[y ye]
Uchief =
-0.069813
yo =
-3.4963
ans =
-1 0.15813
-0.8 0.094939
-0.6 0.053566
-0.4 0.027695
-0.2 0.011616
0 0
0.2 -0.012293
0.4 -0.030394
0.6 -0.059594
0.8 -0.10553
1 -0.1744
trace full fan
y = linspace(-1,1,101)';
ye = trig_ray(y*6.25*cos(Uchief), Uchief, cv, th, rn) - yo;
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 3*y.^2 4*y.^3 5*y.^4 6*y.^5];
a = c\d
a =
6.3113
0.64446
6.0187
-0.016692
0.18143
rms deviation from fit
diff = d - c*a;
rms = sqrt(sum(diff.*diff)/length(diff))
rms =
0.0006876