http://www.k-wave.org/forum/topic/how-to-define-pressure-waves-for-multiple-sources Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 23:08:27 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/how-to-define-pressure-waves-for-multiple-sources#post-148 Thu, 12 May 2011 15:27:35 +0000 maren 148@http://www.k-wave.org/forum/ <p>Hi Brad,<br /> many thanks for your response and also for making this work public,<br /> it s pretty clear to me now,<br /> best,<br /> georgios </p> http://www.k-wave.org/forum/topic/how-to-define-pressure-waves-for-multiple-sources#post-146 Wed, 11 May 2011 23:29:10 +0000 Bradley Treeby 146@http://www.k-wave.org/forum/ <p>Hi Georgios,</p> <p>Thanks for your question. Matrices in MATLAB are actually stored in memory as a one-dimensional sequence of numbers. The numbers are stored column by column moving top to bottom and left to right. Instead of using the row and column index to access a particular position in a matrix, we can also use its position in this sequence (this is called the linear index). From the <a href="http://www.mathworks.com/help/techdoc/math/f1-85462.html">MATLAB documentation</a>:</p> <p>"Matrix A</p> <pre><code>A = [2 6 9; 4 2 8; 3 5 1] A = 2 6 9 4 2 8 3 5 1</code></pre> <p>is actually stored in memory as the sequence</p> <p><code>2, 4, 3, 6, 2, 5, 9, 8, 1</code></p> <p>The element at row 3, column 2 of matrix A (value = 5) can also be identified as element 6 in the actual storage sequence. To access this element, you have a choice of using the standard A(3,2) syntax, or you can use A(6), which is referred to as linear indexing."</p> <p>This ordering is used to order the source points (and the sensor points) in k-Wave. For example, if your binary mask for <code>source.p_mask</code> was equal to a circle:</p> <pre><code>source.p_mask = 0 1 0 1 0 1 1 0 1 0 1 0</code></pre> <p>The actual ordering of the source points would follow:<br /> <pre><code>source.p_mask = 0 3 0 1 0 5 2 0 6 0 4 0</code></pre> <p>This is the ordering used in assigning the source terms to <code>source.p</code> etc.</p> <p>I hope that helps - let me know if you still have any questions. We will work on making the documentation clearer for the next release.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/how-to-define-pressure-waves-for-multiple-sources#post-144 Wed, 11 May 2011 14:55:16 +0000 maren 144@http://www.k-wave.org/forum/ <p>Hello,<br /> I am having some difficulty understanding how k-Wave handles different time varying pressures for different sources. Assuming I have 8 sources distributed circularly in the grid, I need to define 8 different pressure signals, one for each source. This I have done by making source.p a (number_of_sources)x(length(kgrid.t_array)) array.<br /> Given that sources are defined just by 1 s on the grid array (there is no ordering), I am having some problem understanding how can I match the rows of the source.p array to the actual sources on the grid. To put it more simply, how does the software know which source signal corresponds to each source?<br /> In the documentation it is mentioned that source.p can be " a matrix of time series following the source elements using MATLAB's standard column-wise linear matrix index ordering" which I am finding hard to decipher.<br /> Any help would be greatly appreciated as I am stuck on this for a couple of days.<br /> best regards,<br /> georgios </p>