http://www.k-wave.org/forum/topic/2d-photoacoustic-simulation-about-a-homogeneous-phantom Support for the k-Wave MATLAB toolbox en-US Thu, 14 May 2026 20:16:13 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/2d-photoacoustic-simulation-about-a-homogeneous-phantom#post-8858 Mon, 10 Jul 2023 02:40:39 +0000 Wenson 8858@http://www.k-wave.org/forum/ <p>Dear all,<br /> I am going to do the 2D photoacoustic simulation about a rectangle homogeneous phantom.<br /> But I found the channel data of the simulation is weird, like two lines, so as the beamforming image.<br /> Here is my MATLAB code:<br /> <pre><code>% Transducer Setting f0 = 11e6; % Center frequency [Hz] fs = 4*f0; % Sampling frequency [Hz] soundv = 1540; lambda = soundv/f0; BW = 60; % 2D Simulation Setting Nx = 128; dx = 0.12e-3; %[m] F_number = 1; depth = F_number*Nx*dx; dz = lambda/3; Nz = 2*round(depth/dz); kgrid = kWaveGrid(Nx,dx,Nz,dz); % Medium medium.sound_speed = soundv; % Source source.p0 = zeros(Nx,Nz); hypo = zeros(Nx,Nz); hypo(:, (Nz/2-50 : Nz/2+50)) = 1e6; % Rectangle phantom source.p0 = hypo; % Sensor sensor.mask = zeros(Nx,Nz); sensor.mask(:,1) = 1; senor.frequency_response = [f0,BW]; % 2D Simulation sensor_data_2D = kspaceFirstOrder2D(kgrid,medium,source,sensor);</code></pre> <p>Is there any stupid mistake I miss?<br /> Thanks. </p>