私はいくつかの高調波でさまざまな周波数のトーンを再生するいくつかのサウンドファイルを生成しています。
最終的には、これらのサウンドは小さなスピーカーを備えたデバイスで再生されます。Python：サウンドファイルの周波数解析
私はスピーカーの周波数応答曲線を持っており、Pythonで次の操作を実行したい：
gnuplot誰もが知っています：
あなたはPylab/Matplotlibは言及していませんが、動作します。以下はその例です（シングルチャンネル信号を想定）。
x, fs, nbits = audiolab.wavread('schubert.wav')
audiolab.play(x, fs)
N = 4*fs # four seconds of audio
X = scipy.fft(x[:N])
Xdb = 20*scipy.log10(scipy.absolute(X))
f = scipy.linspace(0, fs, N, endpoint=False)
pylab.plot(f, Xdb)
pylab.xlim(0, 5000) # view up to 5 kHz
Y = X*H
y = scipy.real(scipy.ifft(Y))
Hとは何ですか？あなたはNを意味しましたか？ – waffleman
いいえ、私は 'H'を意味しました。申し訳ありませんが、私は怠け者なので、一歩一歩をスキップしました。 'h 'がインパルス応答であれば、' H = scipy.fft（h） 'となります。あなたの質問の項目2では、伝達関数の周波数応答が直接与えられているので、 'H'になります。 –
numpyとmatPlotLibを使用できます。以下のコードのようなもの：
spectrum = numpy.fft.fft(signal)
frequencies = numpy.fft.fftfreq(len(spectrum))
pylab.plot(frequencies,spectrum)
pylab.show()
fftスペクトルのグラフが表示されます。
SciPyの紹介には、FFTを行い、結果をプロットすることが含まれています。http://www.scipy.org/Getting_Started –
あなたのロゴはおいしく見えます。 –