scipy.signal.windows.

flattop

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

返回平頂窗。

參數:

M整數

輸出視窗中的點數。若為零，則返回空陣列。當為負數時，會拋出例外。

sym布林值，選用

當為 True（預設值）時，產生對稱視窗，用於濾波器設計。當為 False 時，產生週期性視窗，用於頻譜分析。

返回:

wndarray

視窗，最大值正規化為 1（但如果 M 為偶數且 sym 為 True，則不會出現值 1）。

註解

平頂窗用於在頻域中精確測量訊號振幅，與其他視窗相比，從頻率倉位的中心到邊緣的扇形誤差最小。這是一個 5 階餘弦視窗，其 5 個項經過最佳化，使主瓣最大程度地平坦。[1]

參考文獻

[1]

D’Antona, Gabriele, and A. Ferrero, “Digital Signal Processing for Measurement Systems”, Springer Media, 2006, p. 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]")

返回在新分頁中開啟