Mercurial > hg > ltpda
view m-toolbox/test/test_ao_fngen.m @ 29:54f14716c721 database-connection-manager
Update Java code
author | Daniele Nicolodi <nicolodi@science.unitn.it> |
---|---|
date | Mon, 05 Dec 2011 16:20:06 +0100 |
parents | f0afece42f48 |
children |
line wrap: on
line source
% Tests the fngen method of the AO class. % % M Hewitson 18-04-08 % % $Id: test_ao_fngen.m,v 1.3 2008/09/19 14:03:50 ingo Exp $ % clear fs = 12; nsecs = 1e3; % Make a sine-wave s = ao(plist('tsfcn', 'randn(size(t)) + sin(2*pi*1e-1*t)', 'fs', fs, 'nsecs', nsecs)); % Add noise and signal pl = plist('type', 'highpass', ... 'order', 1, ... 'gain', 2.0, ... 'fs', 10, ... 'fc', 0.2); a = filter(s, miir(pl)); % Compute target PSD Nfft = fs*100; axx = pwelch(a, plist('Nfft', Nfft, 'Scale', 'PSD', 'win', specwin('Hanning', 10))); % Compute new time-series Nsecs = 1e4; % Number of seconds to produce swin = specwin('Kaiser', 10, 150); % Window for blending b = fngen(axx, plist('Win', swin, 'Nsecs', Nsecs)); %% Output from noise generator iplot(b) sxx = pwelch(a,b, plist('Nfft', 100*fs, 'Scale', 'AS')); % iplot(sxx(1)./sxx(2), plist('YScales', {'All', 'lin'})) iplot(sxx)