diff options
Diffstat (limited to 'solution_stability.m')
| -rw-r--r-- | solution_stability.m | 9 |
1 files changed, 5 insertions, 4 deletions
diff --git a/solution_stability.m b/solution_stability.m index 324621f..359335c 100644 --- a/solution_stability.m +++ b/solution_stability.m @@ -1,12 +1,13 @@ -function [ hist_h, hist_x ] = solution_stability( q_0, q_final, x_final, optics_positions, optics_focal_length, lambda, n_hist ) -%SOLUTION_STABILITY Summary of this function goes here -% Detailed explanation goes here +function [ hist_h, hist_x ] = solution_stability( q_0, q_final, x_final, optics_positions, optics_focal_length, lambda, n_hist, self_flag ) +%Plots histogram of n_hist number of sample points +%Area under plot is equal to the stability of the solution. + lens_displacement = 10^-3; x = lens_displacement*randn(n_hist,3); for i=1:n_hist xd= x(i,:); - fitness_array(i) = fitness( q_0, q_final, x_final, optics_positions + xd, optics_focal_length, lambda ); + fitness_array(i) = fitness( q_0, q_final, x_final, optics_positions + xd, optics_focal_length, lambda, self_flag ); end [hist_h, hist_x] = hist(fitness_array,100); |