mc

%% Time-series data

p      = plist('waveform', 'noise', 'fs', 10, 'nsecs', 10);
tsao1  = ao(p);
tsao1.setT0('2008-10-10 12:00:00')

p      = plist('waveform', 'sine wave', 'fs', 10, 'nsecs', 10, 'A', 10, 'f', 1, 'phi', 0);
tsao2  = tsao1+ao(p);
tsao2.setT0('2008-10-10 12:00:02')

iplotyy(tsao1, tsao2, plist('RightLineStyle', '--'))

%% Spectra



p      = plist('waveform', 'noise', 'fs', 10, 'nsecs', 10);
tsao1  = ao(p);
tsao1.setT0('2008-10-10 12:00:00')

p      = plist('waveform', 'sine wave', 'fs', 10, 'nsecs', 10, 'A', 10, 'f', 1, 'phi', 0);
tsao2  = ao(p);
tsao2.setT0('2008-10-10 12:00:02')


fs1 = tsao1.psd;
fs2 = tsao2.psd;

iplotyy(fs1, fs2);

iplotyy(fs1, fs2, plist('LeftYscale', 'lin'));

iplotyy(fs1, fs2, plist('LeftYscale', 'lin', 'RightYscale', 'lin'));

iplotyy(fs1, fs2, plist('LeftYscale', 'lin', 'RightYscale', 'lin', 'LeftXscale', 'lin', 'RightXscale', 'lin'));

iplotyy(fs1, fs2, plist('LeftXscale', 'lin', 'RightXscale', 'lin'));

%% Complex spectra

pzm = pzmodel(1,1,100);
pr1 = pzm.resp;

pzm = pzmodel(1,100,1);
pr2 = pzm.resp;


iplotyy(abs(pr1), phase(pr1))
iplotyy(abs(pr2), phase(pr2))

