% changes data stored as a row vector to a column vector</span>
0018 [N,C]=size(data); <span class="comment">% size of data</span>
0019 K=size(tapers,2); <span class="comment">% size of tapers</span>
0020 tapers=tapers(:,:,ones(1,C)); <span class="comment">% add channel indices to tapers</span>
0021 H=fft(tapers,nfft,1); <span class="comment">% fourier transform of the tapers</span>
0022 Nsp=sum(data,1); <span class="comment">% number of spikes in each channel</span>
0023 Msp=Nsp'./N; <span class="com