6 May
2009
6 May
'09
7:59 p.m.
Le mercredi 06 mai 2009 à 10:21 -0400, josef.pktd@gmail.com a écrit :
On Wed, May 6, 2009 at 10:00 AM, Talbot, Gerry
wrote: Sorry, I guess I wasn't clear, I meant:
for n in xrange(1,N): y[n] = A*x[n] + B*y[n-1]
So y[n-1] is the result from the previous loop iteration.
I was using scipy.signal for this but I have to look up what I did exactly. I think either signal.correlate or using signal.lti.
Josef
Isn't it what scipy.signal.lfilter does ?
y=scipy.signal.lfilter([A],[1,-B],x)
You may be careful with initial conditions...
--
Fabrice Silva