http://www.k-wave.org/forum/topic/extractampphase-and-spect Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 03:11: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/extractampphase-and-spect#post-6579 Thu, 13 Sep 2018 15:57:59 +0000 Bradley Treeby 6579@http://www.k-wave.org/forum/ <p>Hi zhangjinao,</p> <p>I haven't run your code, but if you open <code>extractAmpPhase</code>, that's more or less exactly how it works. However, by default, it also zero pads the FFT and windows the input signal.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/extractampphase-and-spect#post-6534 Thu, 26 Jul 2018 08:23:19 +0000 zhangjinao@hotmail.com 6534@http://www.k-wave.org/forum/ <p>I computed beam pattern pressure P using the following two ways:</p> <p>i) P = extractAmpPhase(sensor_data.p, 1/kgrid.dt, transducer.freq);<br /> P = reshape(extracted_P_amp, kgrid.Nx, kgrid.Ny, kgrid.Nz);<br /> P = P(:,:,kgrid.Nz/2);</p> <p>ii) sensor_P = reshape(sensor_data.p, kgrid.Nx, kgrid.Ny, kgrid.Nz, time_points);<br /> sensor_P_xy = sensor_P(:,:,kgrid.Nz/2,:);<br /> sensor_P_xy = reshape(sensor_P_xy, kgrid.Nx, kgrid.Ny, time_points);<br /> [freq, amp_spect] = spect(sensor_P_xy, 1/kgrid.dt, 'Dim', 3);<br /> [f1_value, f1_index] = findClosest(freq, transducer.freq);<br /> P = amp_spect(:, :, f1_index);</p> <p>I wonder if the above can offer the same results? </p>