elif c == 'a':
[fs, dataI, dataQ] = IOs.readWaveFile(filename)
sp.fftPlot(dataI + 1j * dataQ, fs=fs)
elif c == 'c':
fs = 1000
t = np.arange(0, 1, 1 / fs)
x = 2 * np.cos(2 * np.pi * 100 * t) + 0.01 * np.cos(
2 * np.pi * 200 * t) + 0.005 * np.cos(2 * np.pi * 300 * t)
#plt.plot(x)
#plt.grid()
#plt.show()
#mLib.fftPlot(x, fs=1000)
sd.Thd(x, fs)
sd.ThdN(x, fs)
elif c == 'w':
sRate = 44100.0
t = np.arange(0, 2.0, 1 / sRate)
samples = np.sin(2 * np.pi * 440.0 * t)
IOs.writeWaveFile('test1.wav', samples, int(sRate))
sampleRate = 44100
samples = sd.KarplusStrong(1 * sampleRate, sampleRate, 391)
IOs.writeWaveFile('test2.wav', samples, int(sampleRate))
samples -= np.mean(samples)
windowed = samples * signal.blackmanharris(len(samples))
