Commit 6aeb7726 authored by henryjandrew's avatar henryjandrew
Browse files

Merge branch 'develop' of https://git.rwth-aachen.de/ita/VAMatlab into develop

parents f7becf20 a0f6dcbb
......@@ -5,44 +5,69 @@
% to hard drive.
free_field_only = false;
direct_field_only = false; % with shadowing of edge
diffracted_field_only = false;
%% Connect and set up simple scene
va = VA( 'localhost' );
va.add_search_path( 'D:/Users/andrew/dev/ITASuite/VA/VACore/data' );
va.add_search_path( 'D:/Users/andrew/dev/ITASuite/VA/VAMatlab/matlab' );
va.add_search_path( 'C:/dev/VA/VACore/data' );
va.add_search_path( 'D:/Users/stienen/dev/VA/VACore/data' );
va.add_search_path( 'C:/ITASoftware/Raven/RavenDatabase/SoundDatabase' );
assert( ~( free_field_only && diffracted_field_only ) )
assert( ~( free_field_only && direct_field_only ) )
assert( ~( diffracted_field_only && direct_field_only ) )
c = va.get_homogeneous_medium_sound_speed();
receiver_pos = [ 3 0 3 ]; % OpenGL coordinates
f = ita_ANSI_center_frequencies;
% House corner
w = itaInfiniteWedge( [ 1 0 0 ], [ 0 0 -1 ], [ 0 0 0 ] ); % OpenGL coordinates
L = va.create_sound_receiver( 'VA_Listener' );
receiver_pos = [ 1 1.7 2 ]; % OpenGL coordinates
va.set_sound_receiver_position( L, receiver_pos )
%H = va.create_directivity_from_file( '$(DefaultHRIR)' );
H = va.create_directivity_from_file( 'ITA_Artificial_Head_5x5_44kHz_128.v17.ir.daff' );
%% Connect and set up simple scene
va = VA;
va.set_sound_receiver_directivity( L, H );
try
va.connect;
dry_run = ~va.get_connected;
va.add_search_path( 'D:/Users/andrew/dev/ITASuite/VA/VACore/data' );
va.add_search_path( 'D:/Users/andrew/dev/ITASuite/VA/VAMatlab/matlab' );
va.add_search_path( 'C:/dev/VA/VACore/data' );
va.add_search_path( 'D:/Users/stienen/dev/VA/VACore/data' );
va.add_search_path( 'C:/ITASoftware/Raven/RavenDatabase/SoundDatabase' );
va.add_search_path( 'C:\Users\jonas\sciebo\ITA\Lehre\Masterarbeiten\2019 Henry Andrew\car_pass-by_corner' );
c = va.get_homogeneous_medium_sound_speed();
L = va.create_sound_receiver( 'VA_Listener' );
va.set_sound_receiver_position( L, receiver_pos )
%H = va.create_directivity_from_file( '$(DefaultHRIR)' );
H = va.create_directivity_from_file( 'ITA-Kunstkopf_HRIR_AP11_Pressure_Equalized_3x3_256.v17.ir.daff' );
S = va.create_sound_source( 'VA_Source' );
%X = va.create_signal_source_buffer_from_file( 'musting_idle.wav' );
X = va.create_signal_source_buffer_from_file( 'Full song - Vulfpeck - Dean Town.wav' );
va.set_sound_receiver_directivity( L, H );
S = va.create_sound_source( 'VA_Source' );
va.set_signal_source_buffer_playback_action( X, 'play' )
va.set_signal_source_buffer_looping( X, true );
va.set_sound_source_signal_source( S, X )
%X = va.create_signal_source_buffer_from_file( 'WelcomeToVA.wav' );
X = va.create_signal_source_buffer_from_file( 'LKW Leerlauf 2.wav' );
%X = va.create_signal_source_buffer_from_file( 'Full song - Vulfpeck - Dean Town.wav' );
va.set_signal_source_buffer_playback_action( X, 'play' )
va.set_signal_source_buffer_looping( X, true );
va.set_sound_source_signal_source( S, X )
catch
dry_run = true;
c = 343;
S = -1;
L = -1;
end
% House corner
w = itaFiniteWedge( [ 1 0 0 ], [ 0 0 -1 ], [ 0 -2 0 ], 4 ); % OpenGL coordinates
%% Example for a synchronized scene update & audio processing simulation/auralization
timestep = 128 / 44100; % here: depends on block size and sample rate
manual_clock = 0;
va.set_core_clock( 0 );
if ~dry_run
va.set_core_clock( 0 );
end
spatialstep = 0.1;
disp( [ 'Resulting sound source speed: ' num2str( spatialstep / timestep ) ' m/s' ] )
......@@ -54,27 +79,33 @@ x = linspace( -1, 1, numsteps ) * 50; % motion from x = -50m to x = 50m
n0_of_paths = 2;
H_log = []; % Diffracted field
IL_log = []; % Insertion loss
h = waitbar( 0, 'Hold on, running auralization' );
direct_deleted = false; %set to true when the direct path has been deleted
for n = 1:length( x )
% Modify scene as you please (position has no real effect for prototype generic path renderer)
source_pos = [ x( n ) 1.1 -3 ]; % OpenGL coordinates
%distance = sum( abs( source_pos - receiver_pos ) );
distance = norm(source_pos - receiver_pos);
va.set_sound_source_position( S, source_pos );
source_pos = [ x( n ) 0 -3 ]; % OpenGL coordinates
if ~dry_run
va.set_sound_source_position( S, source_pos );
end
%% Direct sound path
% Manually create direct sound path and diffracted sound path
shadow_zone = ita_diffraction_shadow_zone(w, source_pos, receiver_pos); %is receiver in shadow zone?
shadow_zone = ita_diffraction_shadow_zone( w, source_pos, receiver_pos ); %is receiver in shadow zone?
if n > 1 && last_shadow_zone ~= shadow_zone
fprintf( 'Shadow zone has been crossed on frame %i\n', n )
end
last_shadow_zone = shadow_zone;
% Direct sound path (id 0)
distance = norm( source_pos - receiver_pos );
prop_path_direct = struct();
prop_path_direct.ID = 0;
prop_path_direct.path = strcat( 'direct_', num2str(S), '_to_', num2str(L) ); %sound path ID
prop_path_direct.gain = 1/distance;
prop_path_direct.identifier = 'direct_path';
prop_path_direct.gain = 1 / distance;
prop_path_direct.delay = distance / c; %delay from sound emitting from source to being received at listener
prop_path_direct.frequencies = f; %Frequencies corresponding to the mags
values = ones(1,31);
......@@ -84,61 +115,73 @@ for n = 1:length( x )
prop_path_direct.delete = false; %set to true when an existing path should be deleted
prop_path_direct.position = source_pos;
if free_field_only
prop_path_direct.audible = true;
prop_path_direct.audible = true; % Free field: always audible
else
prop_path_direct.audible = ~shadow_zone;
prop_path_direct.audible = ~shadow_zone; % Direct field: may be occluded
end
% Diffracted sound path (id 1)
%prop_path_direct.audible = false;
apex = w.get_aperture_point( source_pos, receiver_pos );
detour = norm( source_pos - apex ) + norm(apex - receiver_pos );
H_direct = exp( -1i * 2 * pi * f / c / distance ) / distance;
prop_path_diffracted = struct();
%prop_path_diffracted.distance = detour;
prop_path_diffracted.delay = detour / c;
prop_path_diffracted.frequencies = f;
[H, D, A] = ita_diffraction_utd( w, source_pos, receiver_pos, f, c );
%% Diffracted sound path
rho = norm(apex - source_pos);
prop_path_diffracted.gain = A / rho;
prop_path_diffracted.frequency_magnitudes = abs( D );
%prop_path_diffracted.frequency_magnitudes = values;
apex = w.get_aperture_point( source_pos, receiver_pos );
detour = norm( source_pos - apex ) + norm( apex - receiver_pos );
[ H, D, A ] = ita_diffraction_utd( w, source_pos, receiver_pos, f, c );
rho = norm( apex - source_pos );
prop_path_diffracted = struct();
prop_path_diffracted.source = S;
prop_path_diffracted.receiver = L;
prop_path_diffracted.ID = 1;
prop_path_diffracted.path = strcat( 'diff1_', num2str(S), '_to_', num2str(L) );
prop_path_diffracted.identifier = 'diffracted_path';
prop_path_diffracted.delete = false;
prop_path_diffracted.position = apex;
prop_path_diffracted.audible = true;
prop_path_diffracted.delay = detour / c;
prop_path_diffracted.frequencies = f;
prop_path_diffracted.gain = A / rho;
prop_path_diffracted.frequency_magnitudes = abs( D );
% @todo at some point: load all patqhs from pre-calculated simulation -> ITAPropagationPathSim (C++) output
% D:/Users/andrew/dev/ITASuite/VA/VACore/tests/BinauralOutdoorNoiseRendererTest.VACore.ini
%disp(strcat('**diff dirn: ', num2str(apex),' Delay = ',num2str(distance/c),'** Direct sound: ',num2str(source_pos),' Delay = ',num2str(detour/c),'**'));
% Update wave fronts in renderer
%% Update wave fronts in renderer
paths_update = struct();
paths_update.prop_path_0 = prop_path_direct;
if ~free_field_only
paths_update.prop_path_1 = prop_path_diffracted;
if ~diffracted_field_only
paths_update.direct_path = prop_path_direct;
end
if ~free_field_only && ~direct_field_only
paths_update.diffracted_path = prop_path_diffracted;
end
va.set_rendering_module_parameters( 'MyBinauralOutdoorNoise', paths_update );
% Increment core clock
manual_clock = manual_clock + timestep;
va.call_module( 'manualclock', struct( 'time', manual_clock ) );
% Process audio chain by incrementing one block
va.call_module( 'virtualaudiodevice', struct( 'trigger', true ) );
if ~dry_run
rends = va.get_rendering_modules();
for rend = 1:numel( rends )
if strcmpi( rend.class, 'BinauralOutdoorNoise' )
va.set_rendering_module_parameters( rend.id, paths_update );
end
end
% Increment core clock
manual_clock = manual_clock + timestep;
va.call_module( 'manualclock', struct( 'time', manual_clock ) );
% Process audio chain by incrementing one block
va.call_module( 'virtualaudiodevice', struct( 'trigger', true ) );
end
waitbar( n / numsteps )
H_log = [ H_log; H ];
if shadow_zone
IL_log = [ IL_log; H ./ H_direct ];
else
IL_log = [ IL_log; ( H + H_direct ) ./ H_direct ];
end
end
close( h )
va.disconnect
disp( 'Stop VA to export simulation results from rendering module(s)' )
if ~dry_run
va.disconnect
disp( 'Stop VA to export simulation results from rendering module(s)' )
end
......@@ -6,10 +6,13 @@
% Requires ITA-Toolbox, obtain from http://www.ita-toolbox.org
%ppa_folder = '../../../ITAGeometricalAcoustics/ITAPropagationPathSim/tests/CombinedModel/UrbanTrajectory';
ppa_folder = '../../../dist/win32-x64.vc12/bin/UrbanTrajectory';
if isfolder( ppa_folder ) && ~exist( ppa_folder, 'dir' )
if ~exist( ppa_folder, 'dir' )
ppa_folder = '../../../ITAGeometricalAcoustics/ITAPropagationPathSim/tests/CombinedModel/UrbanTrajectory';
end
if ~exist( ppa_folder, 'dir' )
error( 'Could not find trajectory source folder "%s".', ppa_folder )
end
ppa_diffraction_model = 'utd';
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment