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

R Schumacher rays at blue-cove.com
Tue Dec 8 13:12:03 EST 2015

```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
><<mailto:rays at blue-cove.com>rays at 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 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 really only matters for low frequencies.
>
>I suggested that we might be able to do a time
>domain correction as a forward-reverse FFT
>filter using the function, but another said it
>can also be applied in the time domain using a bilinear transform.
>So, can one use
><http://docs.scipy.org/doc/scipy-0.16.0/reference/generated/scipy.signal.bilinear.html>http://docs.scipy.org/doc/scipy-0.16.0/reference/generated/scipy.signal.bilinear.html
>and, how does one generate b,a from the given
>Fourrier domain flattening function?
>
>
>This should go to either
><mailto:scipy-user at scipy.org>scipy-user at scipy.org
>  or <mailto:scipy-dev at scipy.org>scipy-dev at scipy.org
>
>Chuck
>_______________________________________________
>NumPy-Discussion mailing list
>NumPy-Discussion at scipy.org
>https://mail.scipy.org/mailman/listinfo/numpy-discussion
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/numpy-discussion/attachments/20151208/f1aa7813/attachment.html>
```