forked from sfstoolbox/sfs-matlab
-
Notifications
You must be signed in to change notification settings - Fork 0
/
time_response_point_source.m
98 lines (88 loc) · 4.29 KB
/
time_response_point_source.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
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
function varargout = time_response_point_source(X,xs,conf)
%TIME_RESPONSE_POINT_SOURCE simulates the time response for a point source at
%the given listener position
%
% Usage: [s,t] = time_response_point_source(X,xs,conf)
%
% Input parameters:
% X - listener position / m
% xs - position of point source / m
% conf - configuration struct (see SFS_config)
%
% Output parameters:
% s - simulated time response
% t - corresponding time axis / s
%
% TIME_RESPONSE_POINT_SOURCE(X,xs,conf) simulates the time response of the
% sound field at the given position X. The sound field is simulated for a
% point source at the given source position xs.
%
% See also: sound_field_imp_wfs, freq_response_wfs, time_response_nfchoa
%*****************************************************************************
% The MIT License (MIT) *
% *
% Copyright (c) 2010-2017 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. *
% *
% http:https://sfstoolbox.org [email protected] *
%*****************************************************************************
%% ===== Checking of input parameters ==================================
nargmin = 3;
nargmax = 3;
narginchk(nargmin,nargmax);
isargposition(X);
isargxs(xs);
isargstruct(conf);
%% ===== Configuration ==================================================
% Plotting result
useplot = conf.plot.useplot;
fs = conf.fs;
showprogress = conf.showprogress;
% Disable progress bar for the sound field function
conf.showprogress = 0;
% Disable plotting, otherwise the sound_field_imp fails
conf.plot.useplot = 0;
%% ===== Computation ====================================================
% Get the position of the loudspeaker from point source position.
% NOTE: its directivity [0 -1 0] will be ignored
x0 = [xs 0 -1 0 1];
% Generate time axis (0-500 samples)
t = (0:500)';
s = zeros(1,length(t));
for ii = 1:length(t)
if showprogress, progress_bar(ii,length(t)); end
% Calculate sound field at the listener position
p = sound_field_imp(X(1),X(2),X(3),x0,'ps',dirac_imp(),t(ii),conf);
s(ii) = real(p);
end
% Return parameter
if nargout>0, varargout{1}=s; end
if nargout>1, varargout{2}=t; end
%% ===== Plotting ========================================================
if nargout==0 || useplot
figure;
figsize(conf.plot.size(1),conf.plot.size(2),conf.plot.size_unit);
plot(t/fs*1000,s);
ylabel('amplitude / dB');
xlabel('time / ms');
end