k-Wave

1. zhangjinao@hotmail.com
   

   Member
   Joined: Jul '18
   Posts: 2

   I computed beam pattern pressure P using the following two ways:

   i) P = extractAmpPhase(sensor_data.p, 1/kgrid.dt, transducer.freq);
   P = reshape(extracted_P_amp, kgrid.Nx, kgrid.Ny, kgrid.Nz);
   P = P(:,:,kgrid.Nz/2);

   ii) sensor_P = reshape(sensor_data.p, kgrid.Nx, kgrid.Ny, kgrid.Nz, time_points);
   sensor_P_xy = sensor_P(:,:,kgrid.Nz/2,:);
   sensor_P_xy = reshape(sensor_P_xy, kgrid.Nx, kgrid.Ny, time_points);
   [freq, amp_spect] = spect(sensor_P_xy, 1/kgrid.dt, 'Dim', 3);
   [f1_value, f1_index] = findClosest(freq, transducer.freq);
   P = amp_spect(:, :, f1_index);

   I wonder if the above can offer the same results?

   Posted 6 years ago #

2. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,643

   Hi zhangjinao,

   I haven't run your code, but if you open extractAmpPhase, that's more or less exactly how it works. However, by default, it also zero pads the FFT and windows the input signal.

   Brad.

   Posted 6 years ago #

RSS feed for this topic

Reply

You must log in to post.