-
Notifications
You must be signed in to change notification settings - Fork 38
/
Copy pathssr_generic_wfs.m
78 lines (71 loc) · 3.94 KB
/
ssr_generic_wfs.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
function ir = ssr_generic_wfs(xs,src,conf)
%SSR_GENRIC_WFS impulse response for the generic renderer of the SoundScape
%Renderer
%
% Usage: ir = ssr_generic_wfs(xs,src,conf)
%
% Input parameters:
% xs - virtual source position / m
% src - source type: 'pw' -plane wave
% 'ps' - point source
% 'fs' - focused source
% conf - configuration struct (see SFS_config)
%
% Output parameters:
% ir - impulse response for the desired loudspeaker array
%
% SSR_GENERIC_WFS(xs,src,conf) calculates an impulse response for a virtual
% source at xs for the loudspeakers of a WFS array. Every loudspeaker of
% the array is represented by one column in the impulse response.
% For the generic renderer it is of importance to know what position the first
% loudspeaker of your array will have. For example, the included circular
% array has its first loudspeaker at phi=0deg which is on the x-axis. If you
% have another setup you have to provide it with conf.secondary_sources.x0.
%
% See also: generic_nfchoa, brs_wfs, driving_function_imp_wfs
%*****************************************************************************
% The MIT License (MIT) *
% *
% Copyright (c) 2010-2019 SFS Toolbox Developers *
% *
% Permission is hereby granted, free of charge, to any person obtaining a *
% copy of this software and associated documentation files (the "Software"), *
% to deal in the Software without restriction, including without limitation *
% the rights to use, copy, modify, merge, publish, distribute, sublicense, *
% and/or sell copies of the Software, and to permit persons to whom the *
% Software is furnished to do so, subject to the following conditions: *
% *
% The above copyright notice and this permission notice shall be included in *
% all copies or substantial portions of the Software. *
% *
% THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR *
% IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, *
% FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL *
% THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER *
% LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING *
% FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER *
% DEALINGS IN THE SOFTWARE. *
% *
% The SFS Toolbox allows to simulate and investigate sound field synthesis *
% methods like wave field synthesis or higher order ambisonics. *
% *
% https://sfs.readthedocs.io [email protected] *
%*****************************************************************************
%% ===== Checking of input parameters ==================================
nargmin = 3;
nargmax = 3;
narginchk(nargmin,nargmax);
isargxs(xs);
isargchar(src);
isargstruct(conf);
%% ===== Configuration ==================================================
N = conf.N;
%% ===== Main ============================================================
% Secondary sources
x0 = secondary_source_positions(conf);
% Create empty impulse response for all secondary sources
ir = zeros(N,size(x0,1));
[x0,idx] = secondary_source_selection(x0,xs,src);
x0 = secondary_source_tapering(x0,conf);
% Driving signals for the active speakers
ir(:,idx) = driving_function_imp_wfs(x0,xs,src,conf);