대부분의 계수가 0 일 때 IIR 필터를 수행하는 방법

Question

I 파이썬에서 오디오 효과를 어떻게해야합니까? 예를 들어, 간단한 에코 효과 화학식 : 이것은 IIR 필터이며 scipy.signal에 lfilter()에 의해 계산 될 수

y[n] = x[n] + k*y[n-1000] 

:

import numpy as np 
import time 
import scipy.signal as signal 

pulse = np.zeros(10000) 
pulse[0] = 1.0 

a = np.zeros(1000) 
a[[0,999]] = 1, -0.7 

start = time.clock() 
out = signal.lfilter([1], a, pulse) 
print time.clock() - start 

import pylab as pl 
pl.plot(out) 
pl.show() 

문제는 다음 계수의 대부분 의 a가 0이고 필터가 매우 빠르게 계산 될 수 있지만, lfilter()는 이것을 인식 할 수 없으며 모든 제로 계수를 사용합니다.

나는이 간단한 예제에 대한 특정 계산을 코딩 할 수 있지만 일반적인 해결책을 찾고 있습니다.

Answer 1

이 시도 : 내 노트북에

import scipy 
import scipy.signal as sig 
import time 

# Input signal. 
x = scipy.randn(50000) 

# Filter coefficients. 
a = scipy.zeros(1001) 
a[[0,-1]] = [1, -0.7] 

# Method using lfilter. 
start = time.clock() 
y0 = sig.lfilter([1], a, x) 
end = time.clock() - start 
print end 

# Method using for loop. 
start = time.clock() 
y1 = x 
for i in range(1000, y1.size): 
    y1[i] += 0.7*y1[i-1000] 
end = time.clock() - start 
print end 

# Check that both outputs are equal. 
print scipy.square(y0-y1).sum() 

: 방법 1 0.38 초, 방법 0.13초을 2

참고 : N 샘플의 지연 들어, a[N]을하지 설정해야합니다 a[N-1].