Я намагаюся додати документацію для всіх віконних функцій в scipy.signal, і я застряг на Slepian (такий же, як DPSS?) Та узагальнених гауссових вікнах, про які я раніше ніколи не чув.
Є дві змінні, які є параметрами форми певного типу, p
в узагальненому гауссовому та width
сліпійському. ( sig
Здається, сигма, стандартне відхилення.)
2 питання:
Замість мене зворотної інженерії та здогадувань, хтось може пояснити, як називаються ці змінні і що вони роблять?
Чи можете ви пояснити, для чого ці вікна корисні чи де вони використовуються?
def general_gaussian(M, p, sig, sym=True):
"""Return a window with a generalized Gaussian shape.
The Gaussian shape is defined as ``exp(-0.5*(x/sig)**(2*p))``, the
half-power point is at ``(2*log(2)))**(1/(2*p)) * sig``.
"""
if M < 1:
return np.array([])
if M == 1:
return np.ones(1, 'd')
odd = M % 2
if not sym and not odd:
M = M + 1
n = np.arange(0, M) - (M - 1.0) / 2.0
w = np.exp(-0.5 * (n / sig) ** (2 * p))
if not sym and not odd:
w = w[:-1]
return w
def slepian(M, width, sym=True):
"""Return the M-point slepian window.
"""
if (M * width > 27.38):
raise ValueError("Cannot reliably obtain slepian sequences for"
" M*width > 27.38.")
if M < 1:
return np.array([])
if M == 1:
return np.ones(1, 'd')
odd = M % 2
if not sym and not odd:
M = M + 1
twoF = width / 2.0
alpha = (M - 1) / 2.0
m = np.arange(0, M) - alpha
n = m[:, np.newaxis]
k = m[np.newaxis, :]
AF = twoF * special.sinc(twoF * (n - k))
[lam, vec] = linalg.eig(AF)
ind = np.argmax(abs(lam), axis=-1)
w = np.abs(vec[:, ind])
w = w / max(w)
if not sym and not odd:
w = w[:-1]
return w
Можливі матчі:
Функція dpss_windows nipy використовує NW
"стандартизовану половину пропускної здатності, що відповідає 2NW = BW * f0 = BW * N / dt, але з dt, прийнятим як 1"
Dpss Matlab використовує time_halfbandwidth
Це те саме вікно? Це time_halfbandwidth
те саме, що і width
?
Це визначення DPSS має "бажану частоту відсічення основної долі в радіанах на секунду".
Узагальнений нормальний розподіл має β (рівний вдвічі p
?), Який просто називають параметром форми, з нормальним розподілом для β = 1 і розподілом Лапласа для β = 2.