當前位置: 首頁>>代碼示例 >>用法及示例精選 >>正文


Python SciPy windows.flattop用法及代碼示例


本文簡要介紹 python 語言中 scipy.signal.windows.flattop 的用法。

用法:

scipy.signal.windows.flattop(M, sym=True)#

返回一個平頂窗口。

參數

M int

輸出窗口中的點數。如果為零,則返回空數組。當它為負數時會拋出異常。

sym 布爾型,可選

當為 True(默認)時,生成一個對稱窗口,用於濾波器設計。當為 False 時,生成一個周期窗口,用於頻譜分析。

返回

w ndarray

最大值歸一化為 1 的窗口(如果 M 為偶數且 sym 為 True,則不會出現值 1)。

注意

平頂窗用於對頻域中的信號幅度進行準確測量,與其他窗相比,從頻率倉中心到其邊的扇形誤差最小。這是一個 5th-order 餘弦窗口,優化了 5 個項以使主瓣最大程度地平坦。 [1]

參考

[1]

D'Antona、Gabriele 和 A. Ferrero,“測量係統的數字信號處理”,Springer Media,2006 年,第70 DOI:10.1007/0-387-28666-7

例子

繪製窗口及其頻率響應:

>>> import numpy as np
>>> from scipy import signal
>>> from scipy.fft import fft, fftshift
>>> import matplotlib.pyplot as plt
>>> window = signal.windows.flattop(51)
>>> plt.plot(window)
>>> plt.title("Flat top window")
>>> plt.ylabel("Amplitude")
>>> plt.xlabel("Sample")
>>> plt.figure()
>>> A = fft(window, 2048) / (len(window)/2.0)
>>> freq = np.linspace(-0.5, 0.5, len(A))
>>> response = 20 * np.log10(np.abs(fftshift(A / abs(A).max())))
>>> plt.plot(freq, response)
>>> plt.axis([-0.5, 0.5, -120, 0])
>>> plt.title("Frequency response of the flat top window")
>>> plt.ylabel("Normalized magnitude [dB]")
>>> plt.xlabel("Normalized frequency [cycles per sample]")
scipy-signal-windows-flattop-1_00.png
scipy-signal-windows-flattop-1_01.png

相關用法


注:本文由純淨天空篩選整理自scipy.org大神的英文原創作品 scipy.signal.windows.flattop。非經特殊聲明,原始代碼版權歸原作者所有,本譯文未經允許或授權,請勿轉載或複製。