k-Wave Toolbox

• Getting Started
• Run the First Example
• Build Your Own Simulation
• Troubleshooting
• Reference List
• Examples
• Initial Value Problems
• Example: Homogenous Propagation Medium
• Example: Using A Binary Sensor Mask
• Example: Defining A Sensor Mask By Opposing Corners
• Example: Loading External Image Maps
• Example: Heterogeneous Propagation Medium
• Example: Saving Movie Files
• Example: Recording The Particle Velocity
• Example: Defining A Gaussian Sensor Frequency Response
• Example: Comparison Of Modelling Functions
• Example: Setting An Initial Pressure Gradient
• Example: Simulations In One Dimension
• Example: Simulations In Three Dimensions
• Example: Photoacoustic Waveforms in 1D, 2D and 3D
• Time Varying Source Problems
• Example: Monopole Point Source In A Homogeneous Propagation Medium
• Example: Dipole Point Source In A Homogeneous Propagation Medium
• Example: Simulating Transducer Field Patterns
• Example: Steering A Linear Array
• Example: Snell's Law And Critical Angle Reflection
• Example: The Doppler Effect
• Example: Diffraction Through A Slit
• Example: Simulations In Three-Dimensions
• Sensor Directivity
• Example: Focussed Detector in 2D
• Example: Focussed Detector in 3D
• Example: Modelling Sensor Directivity in 2D
• Example: Modelling Sensor Directivity in 3D
• Example: Sensor Element Directivity in 2D
• Example: Focussed 2D Array with Directional Elements
• Photoacoustic Image Reconstruction
• Example: 2D FFT Reconstruction For A Line Sensor
• Example: 3D FFT Reconstruction For A Planar Sensor
• Example: 2D Time Reversal For A Line Sensor
• Example: 2D Time Reversal For A Circular Sensor
• Example: 3D Time Reversal For A Planar Sensor
• Example: 3D Time Reversal For A Spherical Sensor
• Example: Image Reconstruction With Directional Sensors
• Example: Image Reconstruction With Bandlimited Sensors
• Example: Iterative Image Improvement Using Time Reversal
• Example: Attenuation Compensation Using Time Reversal
• Example: Attenuation Compensation Using Time Variant Filtering
• Example: Automatic Sound Speed Selection
• Diagnostic Ultrasound Simulation
• Example: Defining An Ultrasound Transducer
• Example: Simulating Ultrasound Beam Patterns
• Example: Using An Ultrasound Transducer As A Sensor
• Example: Simulating B-mode Ultrasound Images
• Example: Simulating B-mode Images Using A Phased Array
• Numerical Analysis
• Example: Controlling The Absorbing Boundary Layer
• Example: Source Smoothing
• Example: Filtering A Delta Function Input Signal
• Example: Modelling Power Law Absorption
• Example: Modelling Nonlinear Wave Propagation
• Example: Optimising k-Wave Performance
• Using The C++ Code
• Example: Running C++ Simulations
• Example: Saving Input Files in Parts
• Elastic Wave Propagation
• Example: Explosive Source In A Layered Medium
• Example: Plane Wave Absorption
• Example: Shear Waves And Critical Angle Reflection
• Example: Simulations In Three Dimensions
• Functions - By Category
• Functions - Alphabetical List
• Release Notes
• License

k-Wave Toolbox

fitPowerLawParams

Fit power law absorption parameters for highly absorbing media

Syntax

[a0_fit, y_fit] = fitAbsorptionValues(a0, y, c0, f_min, f_max)
[a0_fit, y_fit] = fitAbsorptionValues(a0, y, c0, f_min, f_max, plot_fit)

Description

fitPowerLawParams calculates the absorption parameters that should be defined in the simulation functions given the desired power law absorption behaviour defined by a0 and y. This takes into account the actual absorption behaviour exhibited by the fractional Laplacian wave equation.

This fitting is required when using large absorption values or high frequencies, as the fractional Laplacian wave equation solved in kspaceFirstOrder1D, kspaceFirstOrder2D, kspaceFirstOrder3D, and kspaceSecondOrder no longer encapsulates absorption of the form a = a0*f^y.

The returned values should be used to define the medium.alpha_coeff and medium.alpha_power within the simulation functions. The absorption behaviour over the frequency range f_min:f_max will then follow the power law defined by a0 and y.

Inputs

a0	desired power law absorption prefactor [dB/(MHz^y cm)]
y	desired power law exponent
c0	medium sound speed [m/s]
f_min	minimum frequency of interest [Hz]
f_max	maximum frequency of interest [Hz]

Optional Inputs

plot_fit

boolean controlling whether the final fit is displayed (default = false)

Outputs

a0_fit	power law absorption prefactor that should be used to define medium.alpha_coeff in the simulation functions
y_fit	power law exponent that should be used to define medium.alpha_power in the simulation functions

See Also

kspaceFirstOrder1D, kspaceFirstOrder2D, kspaceFirstOrder3D

findClosest

flyThrough

© 2009-2014 Bradley Treeby and Ben Cox.