ment">% mean rate for each channel</span>
0024 meansp=Msp(:,ones(1,K),ones(1,size(H,1)));  <span class="comment">% add taper and frequency indices to meansp</span>
0025 meansp=permute(meansp,[3,2,1]); <span class="comment">% permute to get meansp with the same dimensions as H</span>
0026 data=data(:,:,ones(1,K));<span class="comment">% add taper indices to the data</span>
0027 data=permute(data,[1 3 2]); <span class="comment">% permute data to be of the same dimensions as H</span>
0028 data_proj=data.*