<html>
<body>
Sorry - I'll join there.<br><br>
- Ray<br><br>
At 10:00 AM 12/8/2015, you wrote:<br><br>
<br>
<blockquote type=cite class=cite cite="">On Tue, Dec 8, 2015 at 9:30 AM,
R Schumacher
<<a href="mailto:rays@blue-cove.com">rays@blue-cove.com</a>>
wrote:<br>

<dl>
<dd>We have a function which describes a frequency response correction to
piezo devices we use. To flatten the FFT, it is similar to:<br>

<dd>Cdis_t = .5<br>

<dd>N = 8192<br>

<dd>for n in range(8192):<br>

<dd>  B3 = n * 2560 / N<br>

<dd>  Fc(n) = 1 /
((B3/((1/(Cdis_t*2*pi))**2+B3**2)**0.5)*(-0.01*log(B3) +
1.04145))<br><br>

<dd>In practice it really only matters for low frequencies.<br><br>

<dd>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.<br>

<dd>So, can one use<br>

<dd>
<a href="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</a>
<br>

<dd>and, how does one generate b,a from the given Fourrier domain
flattening function?<br><br>

</dl><br>
This should go to either
<a href="mailto:scipy-user@scipy.org">scipy-user@scipy.org</a> or
<a href="mailto:scipy-dev@scipy.org">scipy-dev@scipy.org</a><br><br>
Chuck <br>
_______________________________________________<br>
NumPy-Discussion mailing list<br>
NumPy-Discussion@scipy.org<br>
<a href="https://mail.scipy.org/mailman/listinfo/numpy-discussion" eudora="autourl">
https://mail.scipy.org/mailman/listinfo/numpy-discussion</a></blockquote>
</body>
</html>