http://www.k-wave.org/forum/topic/creating-circular-source Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 01:21:56 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/creating-circular-source#post-4894 Fri, 12 Dec 2014 12:29:48 +0000 Bradley Treeby 4894@http://www.k-wave.org/forum/ <p>Hi gondolio,</p> <p>The input signal is defined by assigning the time series input/s to <code>source.p</code>. This can be a single time series (in which case it is applied to all the grid points that form part of the source), or a matrix of time series following the source elements using MATLAB's standard column-wise linear matrix index ordering. If you don't want to use the <code>focus</code> function, then just assign <code>input_signal</code> directly to <code>source.p</code>.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/creating-circular-source#post-4882 Tue, 09 Dec 2014 23:34:48 +0000 gondolio 4882@http://www.k-wave.org/forum/ <p>Hello Bradley,</p> <p>I am trying to model a circular transducer without the makeTransducer() function and I'm unsure how to do that without using the focus() function as shown in examples. Do you have any suggestions? Here is what I have so far:</p> <pre><code>source_strength = 0.076e6; % [Pa] tone_burst_freq = 1.5e6; % [Hz] tone_burst_cycles = 60; source_radius = 15; % [grid points] source_shape = makeDisc(Ny, Nz, Ny/2, Nz/2, source_radius); source_z_pos = 1; % [grid points] source.p_mask = zeros(Nx, Ny, Nz); source.p_mask(source_z_pos,:,:) = source_shape; input_signal = toneBurst(1/kgrid.dt, tone_burst_freq, tone_burst_cycles,&#39;Plot&#39;,true,&#39;Envelope&#39;,&#39;Rectangular&#39;, &#39;SignalLength&#39;,(1/kgrid.dt)*200e-6); input_signal = (source_strength./(medium.sound_speed*medium.density)).*input_signal;</code></pre> <p>So what do I do with the input_signal from here. I believe I need to use it in a function to obtain source.p? I don't want the beam pattern to be focused. Thanks! </p>