aboutsummaryrefslogtreecommitdiff
path: root/gaussian_focus.m
diff options
context:
space:
mode:
Diffstat (limited to 'gaussian_focus.m')
-rw-r--r--gaussian_focus.m17
1 files changed, 11 insertions, 6 deletions
diff --git a/gaussian_focus.m b/gaussian_focus.m
index ff444ee..cceeda0 100644
--- a/gaussian_focus.m
+++ b/gaussian_focus.m
@@ -1,9 +1,14 @@
function [ w, s ] = gaussian_focus( w0, s0, f, lambda )
-%GAUSSIAN_FOCUS Summary of this function goes here
-% Detailed explanation goes here
+%Focusing of Gaussian beams based on Sidney A. Self's "Focusing of
+%spherical Gaussian beams" Applied Optics, Vol. 22, Issue 5, pp. 658-661 (1983)
+% Input: lambda = wavelength;
+% w0 = waist before the lens;
+% s0 = distance of waist before lens (s>0 before lens)
+% f = focal length of lens (f>0 for converging lens)
+% Output: w = new waist;
+% s = position of waist from lens
-zR = pi*w0^2/lambda;
-s = f*(1+(s0/f-1)/((s0/f-1)^2+(zR/f)^2));
+zR = pi*w0^2/lambda; %Rayleigh Range
+s = f*(1+(s0/f-1)/((s0/f-1)^2+(zR/f)^2)); %Eq. (9b)
w = w0/sqrt((1-s0/f)^2+(zR/f)^2);
-end
-
+end \ No newline at end of file