cw_angular_spectrum()

Description: This function uses the Angular Spectrum Approach to quickly calculate 3D continuous-wave pressures given an intial pressure field.

Usage

pressure = cw_angular_spectrum(p0, coord_grid, medium, f0, nfft);
pressure = cw_angular_spectrum(p0, coord_grid, medium, f0, nfft, type); 

Arguments

Output Parameters

Notes

To choose the number of FFT terms to use, the simplest method is to zero-pad the source pressure and increase the number of terms until an acceptable error is reached. Numbers of terms that are powers of two (e.g. 512, 1024) usually result in the fastest calculations.