view m-toolbox/test/test_fft.m @ 13:e05504b18072 database-connection-manager

Move more functions to utils.repository
author Daniele Nicolodi <nicolodi@science.unitn.it>
date Mon, 05 Dec 2011 16:20:06 +0100
parents f0afece42f48
children
line wrap: on
line source

% Test fft() operator for AOs.
%
% M Hewitson 19-04-07
%
% $Id: test_fft.m,v 1.14 2009/07/24 13:38:19 luigi Exp $
%
function test_fft()
  
  
  % Make test AOs
  nsecs = 10;
  fs    = 1000;  
  pl = plist('nsecs', nsecs, 'fs', fs, 'tsfcn', 'sin(2*pi*7.433*t) + randn(size(t))');
  a1 = ao(pl);
  a1 = a1.setName;
  
  % Take plain one sided
  a2 = fft(a1,plist('type','plain'));
  iplot(abs(a2), plist('XScales', {'All', 'lin'}));
  
  % Take abs of one-sided
  a3 = abs(fft(a1,plist('type','one'))).^2;
  a4 = psd(a1, plist('Nfft', 10000, 'Win', specwin('Rectangular', 10000)));
  
  % Plot ratio
  rat = a3./a4;
  
  iplot(rat, plist('YScales', 'lin'))
  iplot(a3,a4)
  
  % Two-sided fft
  a5 = fft(a1, plist('type', 'two'));
  iplot(abs(a5), plist('XScales', {'All', 'lin'}));
  
end
% END