http://www.k-wave.org/forum/topic/comparison-of-modelling-functions-example Support for the k-Wave MATLAB toolbox en-US Tue, 02 Jun 2026 11:23:08 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/comparison-of-modelling-functions-example#post-5148 Thu, 16 Jul 2015 12:45:37 +0000 bencox 5148@http://www.k-wave.org/forum/ <p>Hi Suryakant,</p> <p>Check your sensor.mask, which seems to be defined twice. I suspect you need to remove the line <code>sensor.mask = cart2grid(kgrid, sensor.mask);</code></p> <p>Ben </p> http://www.k-wave.org/forum/topic/comparison-of-modelling-functions-example#post-5146 Thu, 16 Jul 2015 11:58:47 +0000 Suryakant 5146@http://www.k-wave.org/forum/ <p>Can anybody please help us? </p> http://www.k-wave.org/forum/topic/comparison-of-modelling-functions-example#post-5141 Sat, 11 Jul 2015 13:26:22 +0000 Suryakant 5141@http://www.k-wave.org/forum/ <p>Hey,</p> <p>We are trying to use this code to produce a graph between pressure vs. time using kspaceSecondOrder function and kspaceFirstOrder2D function but we are not getting the result. Can anyone please help us. Following is the code :</p> <p>Nx = 128;<br /> Ny = 128;<br /> Nz = 128;<br /> dx = .1*10^-3;<br /> dy = .1*10^-3;<br /> t_end = 6*10^-6;<br /> medium.sound_speed = 1500;<br /> medium.alpha_coeff = .75;<br /> medium.alpha_power = 1.5;<br /> disc_magnitude = 5;<br /> disc_x_pos = 50;<br /> disc_y_pos = 50;<br /> disc_radius = 8;<br /> disc_1 = disc_magnitude*makeDisc(Nx, Ny, disc_x_pos, disc_y_pos, disc_radius);<br /> disc_magnitude = 3;<br /> disc_x_pos = 80;<br /> disc_y_pos = 60;<br /> disc_radius = 5;<br /> disc_2 = disc_magnitude*makeDisc(Nx, Ny, disc_x_pos, disc_y_pos, disc_radius);<br /> source.p0 = disc_1 + disc_2;<br /> sensor_x_pos = Nx/2;<br /> sensor_y_pos = Ny/2;<br /> sensor_radius = Nx/2-22;<br /> sensor_arc_angle = 3*pi/2;<br /> sensor.mask =makeCircle(Nx,Ny,sensor_x_pos,sensor_y_pos,sensor_radius,sensor_arc_angle);<br /> kgrid = makeGrid(Nx, dx, Ny,dy);<br /> kgrid.t_array = makeTime(kgrid,medium.sound_speed,[],t_end);<br /> sensor.mask = cart2grid(kgrid, sensor.mask);<br /> sensor_data = kspaceSecondOrder(kgrid, medium, source, sensor, 'ExpandGrid', true);</p> <p>sensor_data_reordered =reorderSensorData(kgrid, sensor , sensor_data);</p> <p>figure;<br /> imagesc(sensor_data_reordered, [-1,1]);<br /> colormap(getColorMap);<br /> ylabel('Sensor Position');<br /> xlabel('Time Step');<br /> colorbar;</p> <p>[~,scale,prefix] = scaleSI(max(kgrid.t_array));<br /> plot(kgrid.t_array*scale,sensor_data.p); </p>