k-Wave

1. kai.rieger
   

   Member
   Joined: Nov '18
   Posts: 3

   Hello,

   I try to simulate a pressure source where all source points should have a certain amplitude. I calculated the pressure profile of my ultrasonic transducer and stored it in a 2D array with values between 0 and 1.
   Unfortunately I received the following error message:
   "If using a labelled source.p_mask, the source labels must be monotonically increasing and start from 1."

   If I understand the error message correctly, then it doesn't work the way I try. Is there any other way to do this?

   I use the following commands:

   %transducer is a 2D array with values between 0-1.
   source.p_mask = transducer;
   source.p = toneBurst(source_envelope, tone_burst_freq, tone_burst_cycles, 'Plot', true);

   Thank you very much for your help!

   Best Regards,
   Kai

   Posted 4 years ago #

2. kai.rieger
   

   Member
   Joined: Nov '18
   Posts: 3

   Hello,

   a small addition to my contribution, when I use source.p0, I can use the presure profile but than it is just a puls and not a toneBurst.

   Best Regards,
   Kai

   Posted 4 years ago #

3. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,640

   Hi Kai,

   The input to source.p_mask must be binary, i.e., a matrix of 1s and 0s. What is the value of transducer in your script?

   Brad.

   Posted 4 years ago #

4. kai.rieger
   

   Member
   Joined: Nov '18
   Posts: 3

   Hello Brad,
   the transducer is a matrix of 1s and 0s created with the command makeLine. That means there is a hard transition at the edge of the source (1 -> 0) and the produces side lobes. The real source has a smooth transition at the edges and therefore no side lobes. For me, it would be nice to chose values between 0 and 1 to build a realistic source.

   Best Regards,
   Kai

   Posted 4 years ago #

5. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,640

   Hi Kai,

   You can do this by creating source.p to have the same number of time series as elements in your source mask. For example, you can pass a vector or matrix of amplitudes to createCWSignals and then reshape this to assign to source.p.

   Hope that helps,

   Brad.

   Posted 4 years ago #

RSS feed for this topic

Reply

You must log in to post.