FarField.m

classdef FarField
    
    properties

    end
    
    methods (Static)
        function field = CalculateFarField(ntheta, nfi, apertureSize, wavelength, histogram)
            field = zeros(1, ntheta);
            theta = linspace(-pi/18, pi/18, ntheta);
            fi = linspace(0, 2*pi, nfi);
            nr = floor(length(histogram)/2);
            r = linspace(0, apertureSize/2, nr);
            c1 = (apertureSize/(2*nr))*2*pi/nfi;
            c2 = 1i*2*pi/wavelength;
            for i = 1:ntheta
                c3 = c2*sin(theta(i));
                for j = 1:nr
                    area = r(j) * c1;
                    c4 = c3*r(j);
                    ind = 2*nr-j+1;
                    for k = 1:nfi
                        field(i) = field(i) + histogram(ind)*exp(c4*sin(fi(k)))*area;
                    end
                end
            end
        end
    end
end