diff options
Diffstat (limited to 'gbeam_propagation.m')
| -rw-r--r-- | gbeam_propagation.m | 32 |
1 files changed, 32 insertions, 0 deletions
diff --git a/gbeam_propagation.m b/gbeam_propagation.m index 524d4ae..cd7c47c 100644 --- a/gbeam_propagation.m +++ b/gbeam_propagation.m @@ -26,8 +26,16 @@ function q = gbeam_propagation(x_pos, q_in, x_in, optics_elements) % we need to flip all optics elements around x_in as well for i=1:length(optics_elements_backw) optics_elements_backw{i}.x=x_in-optics_elements_backw{i}.x; + % there is a tricky case: if position of any optics element coincides + % with x_in i.e it is 0 in backwards x coordinates (or at index 1) then + % we need to apply abcd matrix of that element in advance. + if ( optics_elements_backw{i}.x == 0) + q_in_backw(1) = q_afteer_element( q_in_backw(1), optics_elements_backw{i}.abcd); + end + end + q_backw = gbeam_propagation_froward_only(x_backw, q_in_backw, 0, optics_elements_backw); % now we need to flip the radius of curvature again q_backw = -real(q_backw) + 1i*imag(q_backw); @@ -39,3 +47,27 @@ function q = gbeam_propagation(x_pos, q_in, x_in, optics_elements) endfunction +%!test +%! lambda= 1.064E-6 ; +%! f1=0.1526; +%! lns1.abcd=abcd_lens( f1) ; +%! lns1.x= 0.2; +%! f2=0.019; +%! lns2.abcd=abcd_lens( f2) ; +%! lns2.x= lns1.x+ f1 + f2; +%! lns3.abcd=abcd_lens( 0.0629) ; +%! lns3.x= 0.500; +%! q0=-1.5260e-01 + 1.7310e-04i; +%! x0=0.2; +%! optics={lns1,lns2,lns3}; +%! x=[0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7]; +%! q0=1.8098e-96 + 1.3453e+02i; +%! x0=0; +%! q_test = gbeam_propagation(x,q0,x0,optics); +%! +%! % test beam is the same but starting point is taken midway +%! % and coincides with optics element{1}.x since (x(3) = 0.2) +%! q=gbeam_propagation(x,q_test(3),x(3),optics); +%! assert(q,q_test, -1e-6) + + |