spiketrain_powerspectrum

  SPIKETRAIN_POWERSPECTRUM - Compute power spectrum of a spiketrain using pmtm

   [Pxx,Pxxc,F] = SPIKETRAIN_POWERSPECTRUM(SPIKETIMES, TIME_RESOLUTION)

   Use the Tomson Multi-Taper method (see help pmtm) to calculate the 
   power spectrum of a list of SPIKETIMES, sampled with time resolution
   TIME_RESOLUTION (e.g., 0.001 for 1/1000 of the units of SPIKETIMES).

   One may modify the behavior of the function by passing name/value
   pairs:
   Parameter (default value) :  Description
   ------------------------------------------------------------------------
   NFFT ([])                 :  Number of points to include in the FFT
                             :     (see HELP PMTM)


   Example:
     dt = 0.001;
     t = 0:dt:100;  % 100 seconds long
     rate = 10*rectify(cos(2*pi*0.5*t)*dt); % 
     spiketimes = t(find(rand(size(rate))<rate)); 
     [Pxx,Pxxc,F] = spiketrain_powerspectrum(spiketimes,0.01);
     figure;
     loglog(F,Pxx); 
     hold on;
     loglog(F,Pxxc(:,1),'r-');
     loglog(F,Pxxc(:,2),'r-');
     title('Tomson Multitaper Power Spectral Density Estimate');
     xlabel('Frequency (Hz)');
     ylabel('Power/frequency (dB/Hz)');