numpy.hanning

numpy.hanning(M) [source]

Возвращает окно Ханна.

Окно Ханна представляет собой обрезной профиль, сформированный с помощью взвешенного косинуса.

Параметры:	M : int Количество точек в выходном окне. Если значение равно или меньше нуля, возвращается пустой массив.
Возвращает:	out : ndarray, форма(M,) Окно, с максимальным значением, нормализованным к единице (значение единица появляется только если M нечетное).

См. также

bartlett, blackman, hamming, kaiser

Примечания

Окно Ханна определяется как

Окно Ханна названо в честь Юлиуса фон Ханна, австрийского метеоролога. Также известен как окно Косинус. Некоторые авторы предпочитают называть его окном Ханны, чтобы избежать путаницы с очень похожим окном Хэмминга.

Большинство ссылок на окно Ханна взяты из литературы по обработке сигналов, где оно используется в качестве одной из многих функций окон для сглаживания значений. Оно также известно как функция аподизации (что означает «удаление стопы», т.е. сглаживание разрывов в начале и конце выборочного сигнала) или функция окон.

Ссылки

[R25]

Blackman, R.B. и Tukey, J.W., (1958) Измерение спектра мощностей, Dover Publications, Нью-Йорк.

[R26]

E.R. Kanasewich, «Анализ временных последовательностей в геофизике», The University of Alberta Press, 1975, стр. 106-108.

[R27]

Википедия, «Функция окна», http://en.wikipedia.org/wiki/Window_function

[R28]

W.H. Press, B.P. Flannery, S.A. Teukolsky и W.T. Vetterling, «Численные рецепты», Cambridge University Press, 1986, страница 425.

Примеры

>>> np.hanning(12)
array([ 0.        ,  0.07937323,  0.29229249,  0.57115742,  0.82743037,
        0.97974649,  0.97974649,  0.82743037,  0.57115742,  0.29229249,
        0.07937323,  0.        ])

Постройте окно и его частотную характеристику:

>>> from numpy.fft import fft, fftshift
>>> window = np.hanning(51)
>>> plt.plot(window)
[<matplotlib.lines.Line2D object at 0x...>]
>>> plt.title("Hann window")
<matplotlib.text.Text object at 0x...>
>>> plt.ylabel("Amplitude")
<matplotlib.text.Text object at 0x...>
>>> plt.xlabel("Sample")
<matplotlib.text.Text object at 0x...>
>>> plt.show()

>>> plt.figure()
<matplotlib.figure.Figure object at 0x...>
>>> A = fft(window, 2048) / 25.5
>>> mag = np.abs(fftshift(A))
>>> freq = np.linspace(-0.5, 0.5, len(A))
>>> response = 20 * np.log10(mag)
>>> response = np.clip(response, -100, 100)
>>> plt.plot(freq, response)
[<matplotlib.lines.Line2D object at 0x...>]
>>> plt.title("Frequency response of the Hann window")
<matplotlib.text.Text object at 0x...>
>>> plt.ylabel("Magnitude [dB]")
<matplotlib.text.Text object at 0x...>
>>> plt.xlabel("Normalized frequency [cycles per sample]")
<matplotlib.text.Text object at 0x...>
>>> plt.axis('tight')
(-0.5, 0.5, -100.0, ...)
>>> plt.show()