views:

3792

answers:

5

I'd like to compute power spectral density of time series; do some bandpass, lowpass, and highpass filtering; maybe some other basic stuff.

Is there a nice open-source Java library to do this?

I've hunted a bit without success (e.g., Googling "power spectral density java" or "signal processing java" and clicking through links, looking in Apache Commons, Sourceforge, java.net, etc.).

There are lots of applets, books, tutorials, commercial products, etc., that don't meet my needs.

Update: I found org.apache.commons.math.transform for Fourier transforms. This doesn't implement power spectral density, bandpass, etc., but it is something.

A: 

It looks pretty sparse. Try Signalgo or jein or the Intel Signal Processing Library, although I think the last one is just a JNI wrapper.

I saw a lot of those applets you were talking about. I think you may be able to get the JARs for them and use the class APIs inside. May have to use eclipse and jad to decompile and figure out what they do, though, due to lack of documentation. Try the source on this page for example.

John Ellinwood
Thanks for the refs. Honestly, I appreciate them.signalgo is "alpha (i.e. not tested at all!!)"jein has no source code that I can find.Intel won't be open source.
dfrankow
+2  A: 

I found the book Java Digital Signal Processing and its example source code. You might look through the code to see if it fits your needs.

You can also check out DSP Laboratory.

As duffymo and basszero mentioned in the comments, there have been changes to Java since the publication of Java DSP that may impact some of the code examples. In particular, the (relatively) new Concurrency Utilties package might prove useful.

Bill the Lizard
Wow, this book is vintage 1997. I'm not sure I'd recommend it. The signal processing ideas have not changed, but Java sure has.
duffymo
What changes have been made that you think will have an effect on digital signal processing?
Bill the Lizard
@duffymo: The ONLY advance in Java since then is some of the concurrency primitives (since dsp BEGS for multithreading). Generally dsp is math on primitive arrays. I'm w/ Bill.
basszero
I hadn't thought of that, and I should have. Concurrency is definitely worth mentioning. Thanks, guys.
Bill the Lizard
Good points both - I stand corrected. Thanks.
duffymo
More important for us, a book is unlikely to have open source licensing, so we won't be able to use it for our endeavors.
dfrankow
Of all the possibilities anyone has listed, thus far "DSP Laboratory" (a dead open source project) looks like it might have the right spirit, but I can't tell from the code. I don't understand its FOB, FTB, FTJ, etc.
dfrankow
My apologies, the book has GPL licensing on its source!
dfrankow
A: 

I found another resource, although it's not a library: http://www.dickbaldwin.com/tocdsp.htm. It's just a basic discussion of signal processing and Fourier transform, with some Java examples. See for example tutorials 1478, 1482, 1486. Not sure what the license on the code is.

dfrankow
+2  A: 

My first suggestion is to not do your DSP implementation in Java. My second suggestion would to roll your own simple DSP implementations yourself in Java.


Why not to use Java:

I work in a company who's main selling point is DSP related products (about 300 Million revenue per year)... and none of our DSP is in Java... so forgive me when I am hesitant to read about someone who wants to implement DSP in java.

If you are going to be doing non-trivial DSP then you shouldn't be using Java. The reason that DSP is so painful to implement in Java is because all the good DSP implementations use low level memory management tricks, pointers (crazy amounts of pointers), large raw data arrays, etc.

Why to use Java:

If you are doing simple DSP stuff roll your own Java implementation. Simple DSP things like PSD and filtering are both relatively easy to implement (easy implementation but they won't be fast) because there is soo many implementation examples and well documented theory online.

In my case I implemented a PSD function in Java once because I was graphing the PSD in a Java GUI so it was easiest to just take the performance hit in Java and have the PSD computed in the java GUI and then plot it.


How to implement a PSD:

The PSD is usually just the magnitude of the FFT displayed in dB. Start by copying the C code from Numerical Recipes section talking about FFT. Convert the FFT code to Java. (If your data array is real use Numerical Recipes: FFT of Single Real Function. Look for function: void realft(float data[], unsigned long n, int isign) ). Take the 10 * log10(abs(data[])) of the output from realft(). Now you have your own PSD in Java and you have learned something about DSP/FFT.

Before you get upset about suggesting converting C code to Java let me tell you that I have already implemented this specific code to Java and it works decently (the code doesn't use any fancy C tricks that would make the Java implementation horrible).


How to implement lowpass, bandpass filtering:

The easiest implementation (not the most computationally efficient) would in my opinion be using an FIR filter and doing time domain convolution.

Convolution is very easy to implement it is two nested for loops and there are literally millions of example code on the net.

The FIR filter will be the tricky part if you don't know anything about filter design. The easiest method would be to use Matlab to generate your FIR filter and then copy the coefficents into java. I suggest using firpmord() and firpm() from Matlab. Shoot for -30 to -50 dB attenuation in the stopband and 3 dB ripple in the passband.

Trevor Boyd Smith
I'm not working on this anymore, but it would not be appropriate to start from Numerical Recipes for a commercial app. Apache Commons Math has an FFT: http://commons.apache.org/math/userguide/transform.html.
dfrankow
A: 

Check out MARF

Cary