http://www.k-wave.org/forum/topic/large-scale-modelling-for-doa-problem Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 01:13:44 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/large-scale-modelling-for-doa-problem#post-5142 Tue, 14 Jul 2015 16:48:30 +0000 zia.roghani 5142@http://www.k-wave.org/forum/ <p>Hello AVB,<br /> Run the code below. It has a source. A linear sensors array. perform some DOA algorithm like CAPON or Delay-and-sum on the recorded data you will get the direction of source. It has a linear array of there sensors placed at the possition (-0.4667,-99.8667) (0,-99.8667) and ( 0.4667,-99.8667). </p> <p>clear all;<br /> clc;<br /> close all;<br /> PML_size = 5;</p> <p>C = 1400;<br /> f_max = 1500;<br /> dx = C/(2*f_max);<br /> dy=dx;<br /> grid_length = 200;<br /> Ny = round(grid_length/dy);<br /> Nx = Ny;<br /> kgrid = makeGrid (Nx,dx,Nx,dy);<br /> co_ordinate = sqrt(kgrid.x.^2 + kgrid.y.^2);<br /> mid_pt = find(kgrid.x==0);<br /> mid_pt = mid_pt(1);</p> <p>% define the properties of the propagation medium<br /> medium.sound_speed = C; % [m/s]<br /> medium.alpha_coeff = 0.75; % [dB/(MHz^y cm)]<br /> medium.alpha_power = 1.5; </p> <p>% create the time array<br /> cfl = 0.2; %Courant-Friedrichs-Lewy number<br /> f_sig=1000;<br /> FS=10*f_sig;<br /> dt=1/FS;<br /> kgrid.t_array= 0:dt:0.5;</p> <p>% define a time varying sinusoidal source<br /> A=5; % [Pa]<br /> source.p = A*sin(2*pi*f_sig*kgrid.t_array);</p> <p>%define a tone burst<br /> %sampling_freq = 1/dt; % [Hz]<br /> %tone_burst_freq = 1.5e3; % [Hz]<br /> %tone_burst_cycles = 3;<br /> %source.p = source_mag*toneBurst(sampling_freq, tone_burst_freq, tone_burst_cycles);</p> <p>% define a single source point<br /> source.p_mask = zeros(Nx, Ny);<br /> xs=mid_pt;<br /> ys=mid_pt;<br /> source.p_mask(xs, ys) = 1;</p> <p>%define sensor point<br /> incr=1;% 0.4667m spacing<br /> sensor.mask = zeros(Nx, Ny);<br /> y_sens = 10;<br /> x_sens=(mid_pt-incr):incr:(mid_pt+incr);<br /> sensor.mask(x_sens, y_sens) = 1;</p> <p>%define the acoustic parameters to record<br /> sensor.record = {'p'};<br /> sensor_data = kspaceFirstOrder2D(kgrid, medium, source, sensor, 'PMLInside',true, 'PlotPML', true );</p> <p>%plot the recorded data</p> <p>figure;<br /> subplot(4, 1, 1);<br /> [t_sc scale prefix] = scaleSI(max(kgrid.t_array(:)));<br /> plot(source.p, 'b-');<br /> xlabel(['Time [' prefix 's]']);<br /> ylabel('Signal Amplitude');<br /> axis tight;<br /> title('Source Pressure Signal');</p> <p>subplot(4, 1, 2), plot(kgrid.t_array*scale, sensor_data.p(1, :), 'r-');<br /> xlabel(['Time [' prefix 's]']);<br /> ylabel('Signal Amplitude');<br /> axis tight;<br /> title('Sensor-1 Pressure Signal');</p> <p>subplot(4, 1, 3), plot(kgrid.t_array*scale, sensor_data.p(2, :), 'r-');<br /> xlabel(['Time [' prefix 's]']);<br /> ylabel('Signal Amplitude');<br /> axis tight;<br /> title('Sensor Pressure Signal');</p> <p>subplot(4, 1, 4), plot(kgrid.t_array*scale, sensor_data.p(3, :), 'r-');<br /> xlabel(['Time [' prefix 's]']);<br /> ylabel('Signal Amplitude');<br /> axis tight;<br /> title('Sensor Pressure Signal'); </p> http://www.k-wave.org/forum/topic/large-scale-modelling-for-doa-problem#post-5036 Wed, 08 Apr 2015 16:12:14 +0000 bencox 5036@http://www.k-wave.org/forum/ <p>Hi AVB, </p> <p>What exactly are you stuck on?</p> <p>Ben </p> http://www.k-wave.org/forum/topic/large-scale-modelling-for-doa-problem#post-5024 Fri, 27 Mar 2015 14:56:30 +0000 AVB 5024@http://www.k-wave.org/forum/ <p>Hello!<br /> I am trying to use k-wave toolkit for modelling direction of arrival data for multiple acoustic sensors.</p> <p>Could you tell what is the best way to model 30 kHz acoustic waves from omnidirectional transciver 100m far away from several acoustic rescievers (located 20 to 40 cm far from each other)? </p>