Mercurial > hg > ishara

diff arrows/dsp/apowspec.m @ 0:672052bd81f8

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Initial partial import.
author samer
date Wed, 19 Dec 2012 22:38:28 +0000
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/arrows/dsp/apowspec.m	Wed Dec 19 22:38:28 2012 +0000
@@ -0,0 +1,19 @@
+% apowspec - Arrow to compute power spectra
+%
+% apowspec ::
+%    [[N]] ~'analysis window'
+% -> arrow( {[[N]]}, {[[dftbins(N)]]}, empty).
+% 
+% apowspec ::
+%    [[M]] ~'analysis window',
+%    N     ~'size for FFT'
+% -> arrow( {[[M]]}, {[[dftbins(N)]]}, empty).
+function [o,fx]=apowspec(window,N)
+	W=spdiag(window);
+	if nargin<2, N=length(window); end
+	fx=(0:dftbins(N)-1)/N;
+	o=arr(@ps);
+	function y=ps(x)
+		y=powspec(W*x,N);
+	end
+end