+
−% Test script for ao/filtSubtract
+
−%
+
−% M Nofrarias 14-05-09
+
−%
+
−% $Id: test_ao_filtSubtract.m,v 1.1 2009/08/06 09:12:12 miquel Exp $
+
−%
+
−
+
−%% Make test AOs
+
−
+
−% clear 
+
−
+
−nsecs = 1000;
+
−fs    = 10;
+
−
+
−
+
−pl1 = plist('nsecs', nsecs, 'fs', fs, 'tsfcn', 'randn(size(t))');
+
−a1 = ao(pl1);
+
−
+
−
+
−pl2 = plist('nsecs', nsecs, 'fs', fs, 'tsfcn', 'sin(2*pi*0.01*t)');
+
−a2 = ao(pl2);
+
−
+
−% linear coupling
+
−a3 = a1 + a2
+
−a3.setName('input')
+
−
+
−%%  Filter Subtraction
+
−
+
−w = specwin('Kaiser', 100, 90);
+
−% plist lcohere
+
−pll = plist(...
+
−  'type','ms',...
+
−  'Kdes', 100,...
+
−  'win',w,...
+
−  'Jdes', 200,...
+
−  'Order', 1,...
+
−  'variance','no');
+
−
+
−% plist zDomainFit fit
+
−plf = plist('fs',fs,...
+
−    'AutoSearch','on',...
+
−    'StartPoles',[],...
+
−    'StartPolesOpt','real',...
+
−    'maxiter',500,...
+
−    'minorder',3,...
+
−    'maxorder',4,...
+
−    'FITTOL',1e-8,...
+
−    'MSEVARTOL',1e-8,...
+
−    'Plot','off',...
+
−    'ForceStability','on',...
+
−    'direct term','off',...
+
−    'CheckProgress','off');
+
−  
+
−% plist split
+
−plb = plist('times',[10 1000],...
+
−  'filtspl',[0.1  5]);
+
−
+
−a4 = filtSubtract(a3,a1,pll,plf,plb);
+
−
+
−
+
−%% Plot
+
−
+
−iplot(a1,a3)
+
−
+
−iplot(a4)
+
−
+
−%% Test Subtraction
+
−
+
−b = lpsd(a3,a4)
+
−b.iplot
+
−
+
−%% History
+
−
+
−plot(a4.hist)
+
−
+
−
+
−
+
−
+
−
+
−