Mercurial > hg > ltpda

view m-toolbox/test/test_ao_filtSubtract.m @ 52:daf4eab1a51e database-connection-manager tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix. Default password should be [] not an empty string
author Daniele Nicolodi <nicolodi@science.unitn.it>
date Wed, 07 Dec 2011 17:29:47 +0100
parents f0afece42f48
children
line wrap: on
line source

% 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)