diff options
| -rw-r--r-- | Gaussian2D.m | 13 |
1 files changed, 7 insertions, 6 deletions
diff --git a/Gaussian2D.m b/Gaussian2D.m index 586ea36..7063b06 100644 --- a/Gaussian2D.m +++ b/Gaussian2D.m @@ -1,15 +1,16 @@ function [z] = Gaussian2D(p,X,Y); +% beam position (center of Gaussian) cx = p(1); -cy = p(2); +cy = p(2); +% waiste size 1/e in amplitude, 1/e^2 in intensity wx = p(3); wy = p(4); -amp = p(5); -theta = p(6); -background=p(7); +amp = p(5); % beam amplitude or peak intensity +theta = p(6); % if beam is ellipticle then it is major axis angle +background=p(7); + Xn = (X-cx)*cos(theta) - (Y-cy)*sin(theta); Yn = (X-cx)*sin(theta) + (Y-cy)*cos(theta); -%z = amp*(exp(-2*(Xn).^2./(wx^2)-2*(Yn).^2./(wy^2))) + background; - z = amp*(exp(-2*((Xn).^2./(wx^2)+(Yn).^2./(wy^2)))) + background; |