Playing time-locked square wave using sound() function

2 views (last 30 days)
YCR on 8 Mar 2021
Answered: Gabriele Bunkheila on 26 Mar 2021
I'm a graduate student interested in playing a 1 KHz square wave, at a 50% duty cycle, as a sound mask for a type of stimulation that human participants will receive. This mask will need to last for 0.7 seconds.
I have found examples of how to use the square() function to generate the wave and even graph it, but I was wondering how it may be used with sound() to generate the noise directly from the computer? My other option would be to create an external trigger for a waveform signal generator, but I feel like generating the noise strictly from MATLAB is possible.
Any suggestions?

Accepted Answer

Star Strider
Star Strider on 8 Mar 2021
Try this:
Fs = 1E4;
t = linspace(0, 0.7, 0.7*Fs);
s = ((1+sign(sin(2*pi*t*1E3)))/2).'*[1 1]; % Create Binaural Signal
sound(s, Fs)
plot(t, s)
plot(t, s)
xlim([0 0.07])

More Answers (2)

Gabriele Bunkheila
Gabriele Bunkheila on 26 Mar 2021
Apologies for the late response. If you haven't done so yet, I recommend taking a look at audioOscillator and possibly its joint use with audioDeviceWriter. You can find some code to get started in this seemingly unrelated example.
One of the advantages of audioOscillator over other (still valid) legacy signal-generation functions pointed out in previous answers is that it enables generating the signal samples live, while allowing to tune parameters like duty cycle or tone frequency during execution (for example, via user interfaces). I don't understand that is a requirement for you here, but I thought I'd point that out for future visitors.

Walter Roberson
Walter Roberson on 8 Mar 2021
You can use sound() as Star Strider shows, but the timing on it is terrible for the purpose of stimulation studies.
AudioPlayer offers more control over timing, but it is still weak for the purpose of stimulation studies.
Audio System Toolbox offers finer control than either of those.
However... for stimulation studies, I recommend that you look at Psychtoolbox, a third-party product specifically designed for psychophysical experiments. It has visual and audio facilities, and tries hard to be precise, such as having a number of different internal strategies to work with different hardware setups to provide frame-accurate timings.
Walter Roberson
Walter Roberson on 8 Mar 2021
sound() is one of the worst choices for accurate timings. It constructs internal audio mechanisms and runs them whenever it is ready. It returns to the command line after it has started playing, before it has finished, but you cannot query the state.

Sign in to comment.


Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!