Commit e0c0f059 authored by Dipl.-Ing. Jonas Stienen's avatar Dipl.-Ing. Jonas Stienen
Browse files

Using get_aperture_point2 because old impl is faulty

parent afa0bd83
function aperture_point = get_aperture_point( obj, source_pos, receiver_pos )
% get_aperture_point Returns aperture point on wedge (closest point on wedge
% between source and receiver)
aperture_point = obj.get_aperture_point2( source_pos, receiver_pos );
return
% ... the rest is faulty!
assert( numel( source_pos ) == 3 )
assert( numel( receiver_pos ) == 3 )
......
......@@ -8,7 +8,7 @@ function diffr_field = ita_diffraction_utd_approx( wedge, source_pos, receiver_p
% [1] Tsingos, Funkhouser et al. - Modeling Acoustics in Virtual Environments using the Uniform Theory of Diffraction
%
% Example:
% att = ita_diffraction_utd_approximated( wedge, source_pos, receiver_pos, frequency_vec )
% att = ita_diffraction_utd_approximated( wedge, source_pos, receiver_pos, frequency_vec, speed_of_sound )
%
%% Assertions
if nargin < 6
......
......@@ -14,7 +14,7 @@ freq = [100, 200, 400, 800, 1600, 3200, 6400, 12800, 24000]';
alpha_d = linspace( 0, inf_wdg.opening_angle, 100 );
% Set different receiver positions rotated around the aperture
rcv_positions = ita_diffraction_align_points_around_aperture( inf_wdg, rcv_start_pos, alpha_d, apex_point, ref_face );
rcv_positions = ita_diffraction_align_points_around_aperture( inf_wdg, rcv_start_pos, alpha_d, apex_point, 1 );
rcv_positions = rcv_positions(2 : end-1, :);
% wavenumber
......
......@@ -18,4 +18,4 @@ end
omega = 2 * pi * f;
T = 1i .* omega .* m_dash ./ ( 2 .* Z + 1i .* omega .* m_dash ) - 1;
% R = 10 log20( 1 + ( omega m' / ( 2 rh0_0 c ) ^2 )
% R = 10 log10( 1 + ( omega m' / ( 2 rh0_0 c ) ^2 )
Supports Markdown
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