diff options
| author | Matt Argao <mcargao@email.wm.edu> | 2013-02-06 15:04:23 -0500 |
|---|---|---|
| committer | Matt Argao <mcargao@email.wm.edu> | 2013-02-06 15:04:23 -0500 |
| commit | 83b8b9b72f2d9fe27e98e7ad2ad27d545dea9873 (patch) | |
| tree | ba8e0f84c9591c68a3f73502a4a26596b5070277 /solution_stability.m | |
| parent | 8c4278397279db6644624f79179c05bc300a78da (diff) | |
| download | mode_match-83b8b9b72f2d9fe27e98e7ad2ad27d545dea9873.tar.gz mode_match-83b8b9b72f2d9fe27e98e7ad2ad27d545dea9873.zip | |
Added subplot of solution's stability under propagation plot.
Diffstat (limited to 'solution_stability.m')
| -rw-r--r-- | solution_stability.m | 7 |
1 files changed, 3 insertions, 4 deletions
diff --git a/solution_stability.m b/solution_stability.m index f3363b1..324621f 100644 --- a/solution_stability.m +++ b/solution_stability.m @@ -1,11 +1,10 @@ -function [ hist_h, hist_x ] = solution_stability( q_0, q_final, x_final, optics_positions, optics_focal_length, lambda ) +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 -n = 10000; lens_displacement = 10^-3; -x = lens_displacement*randn(n,3); +x = lens_displacement*randn(n_hist,3); -for i=1:n +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 ); end |