view m-toolbox/test/test_ao_iplotyy.m @ 12:86aabb42dd84 database-connection-manager

Use utils.repository utilities
author Daniele Nicolodi <nicolodi@science.unitn.it>
date Mon, 05 Dec 2011 16:20:06 +0100
parents f0afece42f48
children
line wrap: on
line source

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