aboutsummaryrefslogtreecommitdiff
path: root/gaussian_focus.m
diff options
context:
space:
mode:
authorMatt Argao <mcargao@email.wm.edu>2013-03-31 22:35:49 -0400
committerMatt Argao <mcargao@email.wm.edu>2013-03-31 22:35:49 -0400
commit4d08c162d5ce4ac8b7eaeb9d985bfef62c886e99 (patch)
tree45733367770ddc4585faa09976b403a858d5d8e5 /gaussian_focus.m
parent2f3d161985ae2e57ed966981e58eb117d3a49c71 (diff)
downloadmode_match-4d08c162d5ce4ac8b7eaeb9d985bfef62c886e99.tar.gz
mode_match-4d08c162d5ce4ac8b7eaeb9d985bfef62c886e99.zip
Added descriptions to each function.v3.0
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