def amstim(msec=1000, modfreq=20, carrierfreq=200, samplingfreq=1000, amplitude=1.0, carriershape=numpy.sin, modshape=numpy.sin): """create an amplitude-modulated periodic wave""" ### carrier = Basic.wavegen(freq_hz=carrierfreq, samplingfreq_hz=samplingfreq, duration_msec=msec, container=wav(), waveform=carriershape) stim = Basic.ampmod(carrier, freq_hz=modfreq, waveform=modshape) stim.autoscale(amplitude) return stim
def zap(msec=10, samplingfreq=44100, amplitude=0.95, shape='sine'): """create a bipolar pulse""" ### f = 1000.0 / float(msec) if shape in ['sine', 'sin', 'square']: carrier = Basic.wavegen(freq_hz=f, samplingfreq_hz=samplingfreq, duration_msec=msec, container=wav(), waveform=numpy.sin) stim = Basic.ampmod(carrier, freq_hz=f) else: raise ValueError, 'unrecognized shape "%s"' % shape if shape == 'square': y = stim.y y[numpy.where(y < 0)] = -1 y[numpy.where(y > 0)] = +1 stim.autoscale(amplitude) return stim
def zap(msec=10, samplingfreq=44100, amplitude=0.95, shape='sine'): """create a bipolar pulse"""### f = 1000.0 / float(msec) if shape in ['sine', 'sin', 'square']: carrier = Basic.wavegen(freq_hz=f, samplingfreq_hz=samplingfreq, duration_msec=msec, container=wav(), waveform=numpy.sin) stim = Basic.ampmod(carrier, freq_hz=f) else: raise ValueError, 'unrecognized shape "%s"' % shape if shape == 'square': y = stim.y y[numpy.where(y < 0)] = -1 y[numpy.where(y > 0)] = +1 stim.autoscale(amplitude) return stim
def amstim(msec=1000, modfreq=20, carrierfreq=200, samplingfreq=1000, amplitude=1.0, carriershape=numpy.sin, modshape=numpy.sin): """create an amplitude-modulated periodic wave"""### carrier = Basic.wavegen(freq_hz=carrierfreq, samplingfreq_hz=samplingfreq, duration_msec=msec, container=wav(), waveform=carriershape) stim = Basic.ampmod(carrier, freq_hz=modfreq, waveform=modshape) stim.autoscale(amplitude) return stim