FarField javítások

RTX class/FarField.m

@@ -7,17 +7,20 @@ classdef FarField
     methods (Static)
         function field = CalculateFarField(ntheta, nfi, apertureSize, wavelength, histogram)
             field = zeros(1, ntheta);
-            theta = linspace(-pi/6, pi/6, 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
                 er = [cos(theta(i)) sin(theta(i)) 0];
                 for j = 1:nr
+                    area = r(j) * c1;
+                    ind = 2*nr-j+1;
                     for k = 1:nfi
                         rv = [0 r(j)*sin(fi(k)) r(j)*cos(fi(k))];
-                        area = (apertureSize/(2*nr)) * r(j)*2*pi/nfi;
-                        field(i) = field(i) + histogram(2*nr-j+1)*exp(1i*2*pi*rv*er'/wavelength)*area;
+                        field(i) = field(i) + histogram(ind)*exp(c2*rv*er')*area;
                     end
                 end
             end