From 5e73c8b9e65d95f0b9da403eed686ff0e8127d0e Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Frey=20Bal=C3=A1zs?= <frey.balazs96@gmail.com>
Date: Mon, 8 Jun 2020 15:23:52 +0200
Subject: [PATCH] =?UTF-8?q?FarField=20jav=C3=ADt=C3=A1sok?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
---
RTX class/FarField.m | 9 ++++++---
1 file changed, 6 insertions(+), 3 deletions(-)
diff --git a/RTX class/FarField.m b/RTX class/FarField.m
index a1ce28d..2ae571a 100644
--- a/RTX class/FarField.m
+++ b/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
--
GitLab