http://www.k-wave.org/forum/topic/%e6%b0%94%e6%b3%a1%e7%bc%ba%e9%99%b7%e6%a8%a1%e6%8b%9f Support for the k-Wave MATLAB toolbox en-US Fri, 15 May 2026 03:48:33 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/%e6%b0%94%e6%b3%a1%e7%bc%ba%e9%99%b7%e6%a8%a1%e6%8b%9f#post-8491 Fri, 01 Apr 2022 02:36:47 +0000 sudelong 8491@http://www.k-wave.org/forum/ <p>I would like to ask about simulating bubble defects when using sinusoidal filtering to display values that exceed the maximum variable allowed by the program. Please ask where the error occurred</p> <p>% create the computational grid<br /> Nx = 110; % number of grid points in the x (row) direction<br /> Ny = 110; % number of grid points in the y (column) direction<br /> dx = 10e-3/Nx; % grid point spacing in the x direction [m]<br /> dy = dx; % grid point spacing in the y direction [m]<br /> kgrid = kWaveGrid(Nx, dx, Ny, dy);<br /> % create initial pressure distribution using makeDisc<br /> disc_x_pos = 55; % [grid points]<br /> disc_y_pos = 55; % [grid points]<br /> disc_radius = 10; % [grid points]<br /> disc = makeDisc(Nx, Ny, disc_x_pos, disc_y_pos, disc_radius);<br /> % define the properties of the propagation medium<br /> medium.sound_speed(disc == 1) = 340 ; % [m/s]<br /> medium.density(disc == 1) = 1290 ; % [kg/m^3]%<br /> % create the time array<br /> kgrid.makeTime(medium.sound_speed);<br /> % define a single source point<br /> source.p_mask = zeros(Nx, Ny);<br /> source.p_mask(29 , 40) = 1;<br /> % define a time varying sinusoidal source<br /> source_freq =2.5e6; % [Hz]<br /> source_mag = 10; % [Pa]<br /> source.p = source_mag * sin(2 * pi * source_freq * kgrid.t_array);<br /> % filter the source to remove high frequencies not supported by the grid<br /> source.p = filterTimeSeries(kgrid, medium, source.p); </p>