Mercurial > hg > ltpda

diff m-toolbox/test/test_ao_zeropad.m @ 0:f0afece42f48

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Import.
author Daniele Nicolodi <nicolodi@science.unitn.it>
date Wed, 23 Nov 2011 19:22:13 +0100
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/test/test_ao_zeropad.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,25 @@
+
+
+% Compute time-series
+Nsecs = 1e2; % Number of seconds desired
+fs    = 10;
+
+s = ao(plist('tsfcn', 'sin(2*pi*1.137*t)', 'fs', fs, 'Nsecs', Nsecs));
+n = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'Nsecs', Nsecs));
+b = n + s;
+
+% Zero pad x3
+bz = zeropad(b, plist('Factor', 10));
+
+% Zero pad with N samples
+cz = zeropad(b, plist('N', 100));
+
+iplot(b,bz,cz)
+
+% Spectra
+bxx = abs(fft(b)).^2;
+bzxx = abs(fft(bz)).^2;
+czxx = abs(fft(cz)).^2;
+
+iplot(bzxx, bxx, czxx)
+