diff options
Diffstat (limited to 'gbeam_propagation_froward_only.m')
| -rw-r--r-- | gbeam_propagation_froward_only.m | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/gbeam_propagation_froward_only.m b/gbeam_propagation_froward_only.m index ba8c1cf..612d30f 100644 --- a/gbeam_propagation_froward_only.m +++ b/gbeam_propagation_froward_only.m @@ -2,6 +2,9 @@ function q = gbeam_propagation_froward_only(x_pos, q_in, x_in, optics_elements) % calculate the 'q' parameter of the Gaussian beam propagating through optical % 'optics_elements' only in the positive direction along 'x' axis at points 'x_pos' % takes the gaussian beam with initial q_in parameter at x_in +% Note!: if optics element position coinside with a x_pos at some place +% the reported value will be the one calculated after the element +% so the resulting value will coinside with submitted q_in % % all x_pos must be to the right of x_in % x_pos must be monotonic! @@ -38,3 +41,24 @@ function q = gbeam_propagation_froward_only(x_pos, q_in, x_in, optics_elements) q(i)=q_last_calc; endfor end + +%!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.8098e-96 + 1.3453e+02i; +%! x0=0; +%! optics={lns1,lns2,lns3}; +%! x=[0, 0.1, 0.2, 0.3, 0,4, 0.5, 0.6, 0.7]; +%! q_test = [1.8098e-96 + 1.3453e+02i, 1.0000e-01 + 1.3453e+02i, -1.5260e-01 + 1.7310e-04i, -5.2600e-02 + 1.7310e-04i, -5.2600e-02 + 1.7310e-04i, 3.4371e+00 + 1.8961e-03i, 3.4371e+00 + 1.8961e-03i, 3.4371e+00 + 1.8961e-03i, 3.4371e+00 + 1.8961e-03i]; +%! +%! q=gbeam_propagation(x,q0,x0,optics); +%! assert(q,q_test, -1e-4) + + |