%---------------------------------------------------------------------------------------------
% PROGRAM: LGModes
% AUTHOR: Shaowen Xuan
% DATE: 6/7/16
%
% INPUT ARGUMENTS:
% ----------------
% z = filed that need to be decomposed
% dz = distance between neighbour elements in z
% w = the beam effective radius
% r = radial position vector (or matrix generated by meshgrid)
% theta = azimuthal position vector (or matrix generated by meshgrid)
% mmat = m meshgrid
% nmat = n meshgrid
%
% OUTPUT ARGUMENTS:
% ----------------
% coeffs = the coefficients of LG modes
%---------------------------------------------------------------------------------------------
function [ coeffs ] = LGDecompose( z, dz, w, r, theta, mmat, nmat)
%LGDecompose Summary of this function goes here
% Detailed explanation goes here
coeffs = zeros(size(mmat));
[rows, columns] = size(mmat);
u = 1;
for row=1:rows
for column=1:columns
z_LGMode = LGModes(mmat(row,column),nmat(row,column),u,w,r,theta);
coeff = sum(sum(z.*conj(z_LGMode)))*dz*dz;
coeffs(row,column) = coeff;
% coeffs(row,column) = real(coeff);
end
end
end
