k-Wave

1. gondolio
   

   Member
   Joined: Dec '14
   Posts: 3

   Hello Bradley,

   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:

   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,'Plot',true,'Envelope','Rectangular',
   'SignalLength',(1/kgrid.dt)*200e-6);
   
   input_signal = (source_strength./(medium.sound_speed*medium.density)).*input_signal;

   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!

   Posted 9 years ago #

2. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,643

   Hi gondolio,

   The input signal is defined by assigning the time series input/s to source.p. 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 focus function, then just assign input_signal directly to source.p.

   Brad.

   Posted 9 years ago #

RSS feed for this topic

Reply

You must log in to post.