From 4d08c162d5ce4ac8b7eaeb9d985bfef62c886e99 Mon Sep 17 00:00:00 2001 From: Matt Argao Date: Sun, 31 Mar 2013 22:35:49 -0400 Subject: Added descriptions to each function. --- self_gbeam_propagation.m | 38 +++++--------------------------------- 1 file changed, 5 insertions(+), 33 deletions(-) (limited to 'self_gbeam_propagation.m') diff --git a/self_gbeam_propagation.m b/self_gbeam_propagation.m index 10bcf05..824be5a 100644 --- a/self_gbeam_propagation.m +++ b/self_gbeam_propagation.m @@ -1,11 +1,11 @@ function [ w, w_pos ] = self_gbeam_propagation( w0, x_lens, f, x0, lambda ) -%SELF_GBEAM_PROPAGATION Summary of this function goes here -% Detailed explanation goes here +%Beam propagation based on Sidney A. Self's "Focusing of spherical Gaussian +%beams". Applied Optics, Vol. 22, Issue 5, pp. 658-661 (1983) n_lens = 3; -w = w0; -w_pos = x0; +w = w0; %Initial waist +w_pos = x0; %Initial waist position for i = 1:n_lens x_from_next_lens = x_lens(i) - w_pos(i); @@ -14,32 +14,4 @@ for i = 1:n_lens w_pos(i+1) = x_lens(i) + s_new; end -% x_pos= x_lens; -% x(1)=w_pos(1); -% waist_prop(1)=w(1); -% -% x_temp=linspace(x0 ,x_pos(1)); -% waist_prop=[waist_prop, w(1)*sqrt(1+((x_temp-w_pos(1))/(pi*w(1)^2/lambda)).^2)]; -% x=[x x_temp]; -% -% for i=1:n_lens -% x_temp=linspace(x_pos(i),x_pos(i+1)); -% waist_prop=[waist_prop, w(i)*sqrt(1+((x_temp-w_pos(i))/(pi*w(i)^2/lambda)).^2)]; -% x=[x x_temp]; -% end -% -% x_temp=linspace(x_pos(end),x_pos(end)+abs(f(end))); -% waist_prop=[waist_prop, w(end)*sqrt(1+((x_temp-w_pos(end))/(pi*w(end)^2/lambda)).^2)]; -% x=[x x_temp]; -% -% -% figure(1) -% plot(x,waist_prop,'r',x,-waist_prop,'r'); -% - - - - - -end - +end \ No newline at end of file -- cgit v1.2.3