diff options
Diffstat (limited to 'fitness.m')
| -rw-r--r-- | fitness.m | 60 |
1 files changed, 40 insertions, 20 deletions
@@ -3,25 +3,26 @@ function [Energy, Waist, Penalty] = fitness( q_0, q_final, x_final, optics_posit x0 = 0;
Np=20; % # of pts between start position and second lens
N_collimated = 10; % # of pts between second lens and third lens to be collimated region
+ Energy = 0;
x1=optics_positions(1);
x2=optics_positions(2);
x3=optics_positions(3);
- x_array1=linspace(x0,x2,Np);
- x_array2=linspace(x2, x3, N_collimated);
- x = cat(2,x_array1,x_array2);
-
- q_f_trial_forward = gbeam_propagation(x,q_0,x0,optics_placer(optics_positions, optics_focal_length));
- [Waist_trial_forward, Radius_trial_forward] = q2wr(q_f_trial_forward, lambda);
- q_f_trial_backward = gbeam_propagation(x,q_final,x_final,optics_placer(optics_positions, optics_focal_length));
- [Waist_trial_backward, Radius_trial_backward] = q2wr(q_f_trial_backward, lambda);
-
- Energy = 0;
- Penalty_waist_mismatch = sum(abs((Waist_trial_forward-Waist_trial_backward)./min(Waist_trial_forward, Waist_trial_backward)))/Np;
-
- Energy = 1e-2*Penalty_waist_mismatch;
-
+% x_array1=linspace(x0,x2,Np);
+% x_array2=linspace(x2, x3, N_collimated);
+% x = cat(2,x_array1,x_array2);
+%
+% q_f_trial_forward = gbeam_propagation(x,q_0,x0,optics_placer(optics_positions, optics_focal_length));
+% [Waist_trial_forward, Radius_trial_forward] = q2wr(q_f_trial_forward, lambda);
+% q_f_trial_backward = gbeam_propagation(x,q_final,x_final,optics_placer(optics_positions, optics_focal_length));
+% [Waist_trial_backward, Radius_trial_backward] = q2wr(q_f_trial_backward, lambda);
+%
+% Energy = 0;
+% Penalty_waist_mismatch = sum(abs((Waist_trial_forward-Waist_trial_backward)./min(Waist_trial_forward, Waist_trial_backward)))/Np;
+%
+% Energy = 1e-2*Penalty_waist_mismatch;
+%
% penalty calculation
% do not put lenses too close to each other and end positions
lens_size=0.03;
@@ -58,14 +59,33 @@ function [Energy, Waist, Penalty] = fitness( q_0, q_final, x_final, optics_posit Energy = Energy + penalty_lenses_outside_optical_path;
% make collimated region between 2nd and 3rd lens
- %intialize intermediate points between lenses
- q_intermediate = q_f_trial_forward((x2<x) & (x<x3));
- lambda_over_waist_sq = (-imag(1./q_intermediate)); %with numerical factor
- coef = 1e-3;
- penalty_not_collimated_beam = coef * sum(exp((std(lambda_over_waist_sq)/mean(lambda_over_waist_sq)/abs(optics_positions(2) - optics_positions(3)))));
+ [w0, r0] = q2wr(q_0, lambda);
+ [ w, w_pos ] = self_gbeam_propagation( w0, optics_positions, optics_focal_length, x0, lambda );
+
+ coef = 1;
+ d_object = ((optics_positions(end) - w_pos(end - 1))^2)^.5;
+ d_lens = optics_positions(end) - optics_positions(end - 1);
+ penalty_not_collimated_beam = coef * exp(-(d_object/d_lens)^2);
+
Energy = Energy + penalty_not_collimated_beam;
- Penalty = [ penalty_lenses_too_closely_spaced; penalty_lenses_too_closely_spaced; penalty_not_collimated_beam];
+% % waist at end matches desired waist
+ coef = 10;
+ waist_desired = 3.709E-5;
+
+ penalty_waist = coef *((waist_desired-w(end))^2)^.5;
+
+ Energy = Energy + penalty_waist;
+
+ %intialize intermediate points between lenses
+% q_intermediate = q_f_trial_forward((x2<x) & (x<x3));
+% lambda_over_waist_sq = (-imag(1./q_intermediate)); %with numerical factor
+%
+% coef = 1e-3;
+% penalty_not_collimated_beam = coef * sum(exp((std(lambda_over_waist_sq)/mean(lambda_over_waist_sq)/abs(optics_positions(2) - optics_positions(3)))));
+% Energy = Energy + penalty_not_collimated_beam;
+%
+ Penalty = [ penalty_lenses_too_closely_spaced; penalty_lenses_outside_optical_path; penalty_not_collimated_beam; penalty_waist];
end
|