clear all;

%% Make some test data

amps = [1 2 3 4 5 6 7];
freqs = [0.005 0.05 0.1 0.2 0.3 0.4 0.5 0.7];
nsecs = [400 200 100 100 100 100 100 100];
gaps  = [50 50 150 50 50 50 50 50];
toffs(1) = gaps(1);
for kk=2:numel(nsecs)
  toffs(kk) = sum(gaps(1:kk)) + sum(nsecs(1:kk-1));
end

sigs = ao(plist('waveform', 'sine wave', ...
  'nsecs', nsecs, ...
  'fs', 10, ...
  'A', amps, ...
  'f', freqs, ...
  'toff', toffs, 'yunits', 'V'))

sigs.zeropad(plist('N', 1000))
sigs = sigs + ao(plist('tsfcn', '10*randn(size(t))', 'fs', 10, 'nsecs', sigs.nsecs));
sigs.setName;
% sigs.iplot

%% filter the data
pzm = pzmodel(1, {[0.05 2]}, {[0.5 2]});
pzm.setIunits('V');
pzm.setOunits('m');
f = miir(pzm, plist('fs', 10));
r = f.resp(plist('f', logspace(-3,log10(5), 1000)));
% iplot(r)

out = filter(sigs, f) + ...
  ao(plist('tsfcn', '1*randn(size(t))', 'fs', 10, 'nsecs', sigs.nsecs));
out.setName;
iplot(sigs, out)

%% Create the start times
startTimes = time;
t0 = sigs.t0;
for ii=1:numel(toffs) 
  startTimes(ii) = t0 + toffs(ii);
end
startTimes

%% Swept-sine measurement


pl = plist('start times', startTimes, ...
            'durations', nsecs, ...
            'input', sigs, 'frequencies', freqs);
          
t1 = quasiSweptSine(out, pl)          


t1.dy ./ abs(t1.y)

%% Plot abs/phase

ppl = plist('markers', {'', 'o'}, 'linestyles', {'', 'none'},...
  'yranges', {[1e-3,1e2], [-180,180]})
iplot(r, t1, ppl)

%% Plot real/imag

ppl = plist('markers', {'', 'o'}, 'linestyles', {'', 'none'},...
  'YERRL', {[], t1.y-t1.dy}, ...
  'COMPLEXPLOTTYPE', 'realimag', 'xscales', {'all', 'log'}, ...
  'yscales', {'lin', 'lin'})
iplot(r, t1, ppl)

%% Plot abs/rad

ppl = plist('markers', {'', 'o'}, 'linestyles', {'', 'none'},...
  'YERRL', {[], t1.y-t1.dy}, ...
  'COMPLEXPLOTTYPE', 'absrad', 'xscales', {'all', 'log'}, ...
  'yscales', {'lin', 'lin'})
iplot(r, t1, ppl)

%% Measure tf with ltfe

Te = ltfe(sigs, out);

iplot(Te)

Te = tfe(sigs, out, plist('navs', 100));

iplot(Te)