crclayton Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 22:31:02 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/curved-medium-boundaries#post-5761 Tue, 22 Nov 2016 16:44:37 +0000 Bradley Treeby 5761@http://www.k-wave.org/forum/ <p>Hi crclayton,</p> <p>There are various functions included in k-Wave for creating other geometries, for example, <code>makeLine</code>, <code>makeCircle</code>, etc. There's a full list in the documentation under "Functions - By Category". If you want to do something more complex, you'll need to code it yourself, I'm afraid.</p> <p>Hope that helps,</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/curved-medium-boundaries#post-5743 Tue, 15 Nov 2016 21:23:22 +0000 crclayton 5743@http://www.k-wave.org/forum/ <p>Problem elaborated on in a more mathematical sense on stackoverflow here: <a href="https://stackoverflow.com/questions/40618988/how-to-generate-a-matrix-with-a-function-pattern" rel="nofollow">https://stackoverflow.com/questions/40618988/how-to-generate-a-matrix-with-a-function-pattern</a> </p> http://www.k-wave.org/forum/topic/curved-medium-boundaries#post-5740 Tue, 15 Nov 2016 02:06:38 +0000 crclayton 5740@http://www.k-wave.org/forum/ <p>I'm trying to expand upon the Snell's Law example: <a href="http://www.k-wave.org/documentation/example_tvsp_snells_law.php" rel="nofollow">http://www.k-wave.org/documentation/example_tvsp_snells_law.php</a></p> <p>But instead of a single straight boundary, I would like the boundaries to be diagonal to model a waveguide that tapers into a narrower point like so: <a href="https://i.imgur.com/KzxoVBc.png" rel="nofollow">https://i.imgur.com/KzxoVBc.png</a></p> <p>Is there a way to assign a <code>medium.sound_speed</code> such that it has some curvature? Currently I'm using the following code to generate two parallel boundaries, but ultimately I would like to tilt the boundaries so they act as a funnel.</p> <p><code>medium.sound_speed(:, 64+width/2 : 64+width/2+1) = c02;</code><br /> <code>medium.sound_speed(:, 64-width/2-1 : 64-width/2) = c02;</code></p> <p>Is there a way to do this?</p> <p>Thank you very much in advance.</p> <p>Edit:</p> <p>I managed to create a rudimentary funnel that tapers at a 45 degree then spreads out again, but it's far from elegant:</p> <pre><code>width = 18; medium.sound_speed(1:60, 64+width/2 : 64+width/2+1) = c02; medium.sound_speed(1:60, 64-width/2-1 : 64-width/2) = c02; for i=1:1:60 medium.sound_speed(i, i) = c02; end for i=1:1:60 medium.sound_speed(i, 128-i) = c02; end width=5; medium.sound_speed(60:80, 64+width/2 : 64+width/2+1) = c02; medium.sound_speed(60:80, 64-width/2-1 : 64-width/2) = c02;</code></pre>