samer@0: % apowspec1 - Arrow to compute power spectra (automatic sizing version)
samer@0: %
samer@0: % apowspec1 ::
samer@0: %    (N:natural -> [[N]]) ~'analysis window function'
samer@0: % -> arrow( {[[N]]}, {[[dftbins(N)]]}, empty).
samer@0:  
samer@0: function o=apowspec1(windowfn)
samer@0: 	o=arrf(@mk,1,1);
samer@0: 	function f=mk(sizes)
samer@0: 		N=sizes{1}(1);	
samer@0: 		W=spdiag(windowfn(N));
samer@0: 		f=@(x)powspec(W*x,N);
samer@0: 	end
samer@0: end