%% Make test data for whiten1D
fs = 100;
a = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'nsecs', 1000,'yunits','V'));

% filter
pzm = pzmodel(1e-1, {[0.1 3]}, {[10 2]});
ft = miir(pzm,plist('fs',fs));

af = filter(a, ft);
save(af,'whiten.mat');