See example in documentation <fft>. As noted if you arrange your data where each channel is a column, then FFT will operate and return the FFT for each channel directly; you can get the PSDs for each by continuing the same pattern on the output arrays as is shown in the example for one channel.
Just give it a go and let us know where you run into a particular issue...
