Re: [Numpy-discussion] Q: Use of scipy.signal.bilinear

2015-12-08 Thread R Schumacher
Sorry - I'll join there. - Ray At 10:00 AM 12/8/2015, you wrote: On Tue, Dec 8, 2015 at 9:30 AM, R Schumacher <r...@blue-cove.com> wrote: We have a function which describes a frequency response correction to piezo devices we use. To flatten the FFT, it is similar

Re: [Numpy-discussion] Q: Use of scipy.signal.bilinear

2015-12-08 Thread Charles R Harris
On Tue, Dec 8, 2015 at 9:30 AM, R Schumacher wrote: > We have a function which describes a frequency response correction to > piezo devices we use. To flatten the FFT, it is similar to: > Cdis_t = .5 > N = 8192 > for n in range(8192): > B3 = n * 2560 / N > Fc(n) = 1 / ((B3/((1/(Cdis_t*2*pi))*

[Numpy-discussion] Q: Use of scipy.signal.bilinear

2015-12-08 Thread R Schumacher
We have a function which describes a frequency response correction to piezo devices we use. To flatten the FFT, it is similar to: Cdis_t = .5 N = 8192 for n in range(8192): B3 = n * 2560 / N Fc(n) = 1 / ((B3/((1/(Cdis_t*2*pi))**2+B3**2)**0.5)*(-0.01*log(B3) + 1.04145)) In practice it reall