http://www.k-wave.org/forum/topic/correct-post-processing-for-applying-transducer-bandwidth Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 00:43:31 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/correct-post-processing-for-applying-transducer-bandwidth#post-4713 Thu, 11 Sep 2014 14:24:11 +0000 Bradley Treeby 4713@http://www.k-wave.org/forum/ <p>Hi Sebastion,</p> <p>The <code>gaussianFilter</code> function applies a simple non-causal or zero-phase filter using the FFT. If you navigate to the "Filtering A Delta Function Input Signal Example" under the "Numerical Analysis" heading in the html help, then open the third m-file (under the "Filtering the input time series" title), and add the following code around line 65:</p> <pre><code>case 4 % filter the input signal freq = 10e6; bandwidth = 50; source_func_filtered = gaussianFilter(source_func, 1/dt, freq, bandwidth);</code></pre> <p>You can then compare the output of <code>gaussianFilter</code> with other types of filtering, and hopefully get a handle on what's happening. In particular, you'll notice that when using <code>gaussianFilter</code> or <code>applyFilter</code>/<code>filterTimeSeries</code> with the option <code>&#39;ZeroPhase&#39;, true</code>, the main features of the signal will still be in roughly the same place, but there might be oscillations <em>ahead</em> of time, i.e., the filter is not causal. </p> <p>If you compare this to <code>applyFilter</code>/<code>filterTimeSeries</code> with the option <code>&#39;ZeroPhase&#39;, false</code>, you'll notice there is a translation of the main features of the signal to ensure nothing appears ahead of time, i.e., the causal nature of the filter causes phase distortion. </p> <p>As you've noticed, the exact nature of the filtered output will depend on the filter shape. The best filter to use depends largely on what you're trying to achieve.</p> <p>Hope that helps,</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/correct-post-processing-for-applying-transducer-bandwidth#post-4694 Mon, 01 Sep 2014 10:21:23 +0000 bastlwastl 4694@http://www.k-wave.org/forum/ <p>Hi,</p> <p>I'm working on an ultrasonic simulation and i would like to compare the result to actual data from our experiments. The general impression seems to be correct, but in detail, the frequency bandwidth of the used transducer plays an important role.<br /> As far as i know, the manual addresses this issue only quickly and in this forum, i found these 2 topics:<br /> <a href="http://www.k-wave.org/forum/topic/simple-3d-simulation-received-signal-question" rel="nofollow">http://www.k-wave.org/forum/topic/simple-3d-simulation-received-signal-question</a><br /> <a href="http://www.k-wave.org/forum/topic/sensor-specifications" rel="nofollow">http://www.k-wave.org/forum/topic/sensor-specifications</a></p> <p>The second suggest to use <code>gaussianfilter</code> which is applying a frequency filtering (as far as i see). What seems to be correct compared to the spectra of the transducer provided by the manufacturer, results in clearly incorrect signals. Have a look here:<br /> <a href="http://homepages.physik.uni-muenchen.de/~Sebastian.Lehrack/gaussian_filter_spectra_3.5MHz.png" rel="nofollow">http://homepages.physik.uni-muenchen.de/~Sebastian.Lehrack/gaussian_filter_spectra_3.5MHz.png</a><br /> <a href="http://homepages.physik.uni-muenchen.de/~Sebastian.Lehrack/compared_signals_filtered.png" rel="nofollow">http://homepages.physik.uni-muenchen.de/~Sebastian.Lehrack/compared_signals_filtered.png</a><br /> The gaussian filter results in a signal which differs from the data we measured. It also produces a signal appearing before the actual signal. In a second step, we applied Matlab's standard butterworth low and high pass filter and the result agrees much better with the measured data but yet not perfectly. However, it seems to me that the <code>gaussianfilter</code> is to simple and to strict regarding higher frequencies. And although it is not indicated, a transducer can have higher resonances, at least theoretically.</p> <p>Arising from this, what is the intended or recommended way to simulate the detector bandwidth?</p> <p>Greetings,<br /> Sebastian </p>