diff options
| -rw-r--r-- | fitter_check.m | 51 |
1 files changed, 19 insertions, 32 deletions
diff --git a/fitter_check.m b/fitter_check.m index 51e1480..5b0f8e3 100644 --- a/fitter_check.m +++ b/fitter_check.m @@ -1,9 +1,9 @@ +%Permute all possible lens combinations out of set of lenses lens_set = [.075, .203, .05, .03]; lens_set = [.075,.203]; lens_permutations = pick(lens_set,3,'or'); -% ########################################## -% Sample Solution +%Pre-defined Constants lambda= 1.064E-6 ; Ltot= 1.010675025828971 ; r0= 1.0E+100 ; @@ -12,42 +12,29 @@ x0= 0 ; wf= 3.709E-5 ; rf= 1.0E+100 ; xf= Ltot; - q0=wr2q(w0,r0,lambda); -x0=0; qf=wr2q(wf,rf,lambda); -xf=Ltot; -% -% optics={lns1,lns2,lns3}; -% figure(1) -% w_final_handmade = solution_visualization(q0,x0, qf, xf, optics, lambda); -% title('Hand made'); -% ########################################## +%End list + +%Mode match +[ possible_lens_placement, possible_lens_set, possible_sample_energy] = mode_match( q0, qf, Ltot, lambda, lens_permutations ); -[ possible_lens_placement, possible_lens_set, possible_sample_energy, n_possible_lens_placement, index ] = mode_match( q0, qf, Ltot, lambda, lens_permutations ); +%Remove similar solutions +n_truncate = 4; +[ possible_lens_placement_uniq, possible_lens_placement, possible_lens_set, index ] = remove_similar_soln( possible_sample_energy, possible_lens_placement, possible_lens_set, index, n_truncate ); %Visualize five best solutions -for n_graph = 1:n_possible_lens_placement - figure(n_graph) - w_final_trial = solution_visualization(q0,x0, qf, xf, optics_placer(possible_lens_placement(index(n_graph),:), possible_lens_set(index(n_graph),:)), lambda); - - str1=sprintf('\n f_1 = %0.4f, x_1 = %0.4f\n',possible_lens_set(index(n_graph),1),possible_lens_placement(index(n_graph),1)); - str2=sprintf('f_2 = %0.4f, x_2 = %0.4f\n',possible_lens_set(index(n_graph),2),possible_lens_placement(index(n_graph),2)); - str3=sprintf('f_3 = %0.4f, x_3 = %0.4f\n',possible_lens_set(index(n_graph),3),possible_lens_placement(index(n_graph),3)); - tstr='Other Solutions'; - title([tstr, str1, str2, str3]); -end +n_visualizations = 5; +pick_visualization( possible_lens_placement_uniq, possible_lens_placement, possible_lens_set, index, n_visualizations, q0, qf, Ltot, lambda ); -possible_lens_placement(index(1:n_graph),:) -possible_lens_set(index(1:n_graph),:) possible_sample_energy(index(1:n_graph),:) -%Visualize fitness function for fixed f2 and f3 -lens_set = [.075, .075, .203]; -f2= 0.40361319425309 ; -f3= 0.80361319425309 ; - -fitness_simplified=@(x) fitness(q0, qf, Ltot, [x, f2, f3], lens_set, lambda ); -figure(6) -ezplot(fitness_simplified, [0,Ltot])
\ No newline at end of file +% %Visualize fitness function for fixed f2 and f3 +% lens_set = [.075, .075, .203]; +% f2= 0.40361319425309 ; +% f3= 0.80361319425309 ; +% +% fitness_simplified=@(x) fitness(q0, qf, Ltot, [x, f2, f3], lens_set, lambda ); +% figure(6) +% ezplot(fitness_simplified, [0,Ltot])
\ No newline at end of file |