diff --git a/RTX class/FarField.m b/RTX class/FarField.m index 2ae571a9284e8b6d5f4639f635bdcd91caf27b6c..c6453f28c330d80d066fa12ab8ef0001396aca21 100644 --- a/RTX class/FarField.m +++ b/RTX class/FarField.m @@ -14,13 +14,12 @@ classdef FarField c1 = (apertureSize/(2*nr))*2*pi/nfi; c2 = 1i*2*pi/wavelength; for i = 1:ntheta - er = [cos(theta(i)) sin(theta(i)) 0]; + c3 = sin(theta(i)); 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))]; - field(i) = field(i) + histogram(ind)*exp(c2*rv*er')*area; + field(i) = field(i) + histogram(ind)*exp(c2*c3*r(j)*sin(fi(k)))*area; end end end