Skip to content

Commit d983b5e

Browse files
committed
Added missing references for mono 2.5D WFS ps driving function
1 parent 544ea78 commit d983b5e

File tree

1 file changed

+9
-5
lines changed

1 file changed

+9
-5
lines changed

SFS_monochromatic/driving_functions_mono/driving_function_mono_wfs_ps.m

Lines changed: 9 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -25,6 +25,9 @@
2525
% Synthesis Revisited", AES124
2626
% E. Verheijen (1997) - "Sound Reproduction by Wave Field Synthesis", PhD
2727
% thesis, TU Delft
28+
% D. Opperschall (2002) - "Realisierung eines Demonstrators für
29+
% Punktquellen und ebene Wellen für ein Wellenfeldsynthese-System",
30+
% Master thesis, Universität Erlangen-Nürnberg
2831
% F. Völk (2010) - "Psychoakustische Experimente zur Distanz mittels
2932
% Wellenfeldsynthese erzeugter Hörereignisse", DAGA, p.1065-66
3033
% S. Spors, J. Ahrens (2010) - "Analysis and Improvement of
@@ -226,8 +229,7 @@
226229
%
227230
elseif strcmp('opperschall',driving_functions)
228231
% --- Opperschall -------------------------------------------------
229-
% Opperschall, Equation (3.14)
230-
% Note: Driving function with only one stationary phase
232+
% Driving function with only one stationary phase
231233
% approximation, reference to one point in field
232234
%
233235
% 2.5D correction factor
@@ -238,14 +240,16 @@
238240
%
239241
g0 = sqrt( vector_norm(x0-xref,2) ./ (vector_norm(xs-x0,2) + vector_norm(x0-xref,2)) );
240242
% ______
241-
% |2pi iw (x0-xs) nx0
243+
% | i w (x0-xs) nx0
242244
% D_2.5D(x0,w) = g0 _ |------ ------------- e^(-i w/c |x0-xs|)
243-
% \| c |x0-xs|^(3/2)
245+
% \|2pi c |x0-xs|^(3/2)
246+
%
247+
% see Opperschall (2002), p.14 (3.1), (3.14), (3.15)
244248
%
245249
% r = |x0-xs|
246250
r = vector_norm(x0-xs,2);
247251
% driving signal
248-
D = sqrt(2*pi*1i*omega/c) * g0 .* vector_product(x0-xs,nx0,2) ./ r.^(3/2) .* exp(-1i*omega/c .* r);
252+
D = sqrt(1i*omega/(2*pi*c)) * g0 .* vector_product(x0-xs,nx0,2) ./ r.^(3/2) .* exp(-1i*omega/c .* r);
249253
%
250254
elseif strcmp('volk2010',driving_functions)
251255
% --- Voelk 2010 --------------------------------------------------

0 commit comments

Comments
 (0)