🚀 go-pugleaf

In comp.dsp <1992Aug5.132704.1@zodiac.dsto.oz.au>
Witold Waldman <st_waldman@zodiac.dsto.oz.au> writes

> I would like to obtain source code in Fortran, C or
> Pascal for computing the Hilbert transform of a band-pass
> filter response function. Note that only magnitude response
> data is available for the band-pass filter, and the Hilbert
> transform would be used to compute the corresponding
> phase response of the band-pass filter.

It seems to me that you are trying to generate complex data from
real only measurements. You have two options: you can implement the
Hilbert Transform with analog components ,(nothing more than a 90
degree phase shifter) or digitally perform the transform. Each has its
own advantages but since you asked for source code I assume you've already
decided on the digital approach. You can just use a generic FFT to do this
as follows: perform FFT on your data and then multiply the array by the
"signum" function ( -I for 0<i<N/2-1 , +I for N/2<i<N-1 ) and then do an
inverse FFT to get the "quadrature" or imaginary part of your signal.
This may not be the most efficient way of doing the transform and in fact
I would be interseted in better ways if anyone knows of any. Hope this helps.

Jacob Hirbawi
JcbHrb@CERF.net