clear all

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%         TSDATA OBJECTS
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Make test AOS

p      = plist('waveform', 'noise', 'fs', 10, 'nsecs', 10);
tsao1  = ao(p);
p      = plist('waveform', 'sine wave', 'fs', 10, 'nsecs', 10, 'f', 1, 'phi', 0);
tsao2  = ao(p);

tsvec = [tsao1 tsao2];
tsmat = [tsao1 tsao2; tsao1 tsao2];

%% Default plot
miplot(tsao1, tsao2)
miplot(tsvec)
miplot(tsmat)

%% Change colors and line styles

pl = plist('Colors', {'g', 'k', 'm'}, 'LineStyles', {'--'}, 'LineWidths', [1 4]);
miplot(tsao1, tsao2, pl);

%% Change arrangement to single plots

pl = plist('Arrangement', 'single');
miplot(tsao1, tsao2, pl);

%% Change arrangement to subplots

% Also override the second legend text and the first line style
pl = plist('Arrangement', 'subplots', 'LineStyles', {'--'}, 'Legends', {'', 'My Sine Wave'});
miplot(tsao1, tsao2, pl);

%% Changel ylabel

% Overide the Y-labels
pl = plist('Arrangement', 'subplots', 'YLabels', {'signal 1', 'signal 2'}, 'Legends', {'', 'My Sine Wave'});
miplot(tsao1, tsao2, pl);

%% Changel xlabel

% Overide the X-labels
pl = plist('Arrangement', 'subplots', 'XLabels', {'', 'Time-stamps'}, 'Legends', {'', 'My Sine Wave'});
miplot(tsao1, tsao2, pl);

%% No legends

pl = plist('Arrangement', 'subplots', 'Legends', 'off');
miplot(tsao1, tsao2, pl);

%% Check math functions

pl = plist('Arrangement', 'subplots', 'YMaths', 'y.^2', 'XMaths', {'', 'log(x)'});
miplot(tsao1, tsao2, pl);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%         FSDATA OBJECTS
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%% Make test AOs

fsd1 = resp(miir(plist('fs', 100)), plist('f', logspace(-4, log10(50), 1000)));
fsd1 = set(fsd1, 'yunits', 'm/V');
fsd2 = ltpda_pwelch(tsao1+tsao2);
fsd3 = fsd2*10;
fsd4 = resp(miir(plist('type', 'highpass')));
fsd4 = set(fsd4, 'yunits', 'm/V');

%% Default plot

miplot(fsd1, fsd2, fsd3, fsd4)

%% Subplots

pl = plist('Arrangement', 'subplots');
miplot(fsd1, fsd2, fsd3, fsd4, pl)

%% Single plots

pl = plist('Arrangement', 'single');
miplot(fsd1, fsd2, fsd3, fsd4, pl)

%% Use math function on y-data

pl = plist('Ymaths', {'abs(y)'});
miplot(fsd1, pl);

%% Change colors and line styles

pl = plist('Colors', {'g', 'k', 'm'}, 'LineStyles', {'--'}, 'LineWidths', [1 4]);
miplot(fsd1, fsd2, pl);

%% Change arrangement to subplots

% Also override the second legend text and the first line style
pl = plist('Arrangement', 'subplots', 'LineStyles', {'--'}, 'Legends', {'', 'My Sine Wave'});
miplot(fsd1, fsd2, pl);

%% Changel ylabel

% Overide the Y-labels
pl = plist('Arrangement', 'subplots', 'YLabels', {'signal 1 mag', '', 'signal 2'}, 'Legends', {'', 'My Sine Wave'});
miplot(fsd1, fsd2, pl);

%% Changel xlabel

% Overide the X-labels
pl = plist('Arrangement', 'subplots', 'XLabels', {'', 'Freq'}, 'Legends', {'', 'My Sine Wave'});
miplot(fsd1, fsd2, pl);

%% No legends

pl = plist('Arrangement', 'subplots', 'Legends', 'off');
miplot(fsd1, fsd2, pl);

%% Set Y scales

pl = plist('Arrangement', 'subplots', 'YScales', {'', 'log', 'lin'});
miplot(fsd1, fsd2, pl);

%% Set X scales

pl = plist('Arrangement', 'subplots', 'XScales', {'', 'lin', 'lin'});
miplot(fsd1, fsd2, pl);

%% Set Y ranges

pl = plist('Arrangement', 'subplots', 'YRanges', {[0.1 10], [], [1e-3 100]});
miplot(fsd1, fsd2, pl);


%% Set X ranges

pl = plist('Arrangement', 'subplots', 'YRanges', [0.001, 10], 'XRanges', {[0.1 10], [0 10], [0.1 10]});
miplot(fsd1, fsd2, pl);


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%         XYDATA OBJECTS
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%% Make test AOs

xy1  = ao(xydata(1:100, cos(1:100)));
xy1  = set(xy1, 'xunits', 'NA', 'yunits', 'NA');
xy2  = ao(xydata(1:100, sin(1:100)));
xy2  = set(xy2, 'xunits', 'NA', 'yunits', 'NA');


%% Default plot
miplot(xy1, xy2)

%% Use Math function on y-data

pl = plist('Ymaths', {'abs(y)', 'y.^2'}, 'Xmaths', 'log(x)');
miplot(xy1, xy2, pl);


%% Change colors and line styles

pl = plist('Colors', {'g', 'k', 'm'}, 'LineStyles', {'--'}, 'LineWidths', [1 4]);
miplot(xy1, xy2, pl);

%% Change arrangement to single plots

pl = plist('Arrangement', 'single');
miplot(xy1, xy2, pl);

%% Change arrangement to subplots

% Also override the second legend text and the first line style
pl = plist('Arrangement', 'subplots', 'LineStyles', {'--'}, 'Legends', {'', 'My Sine Wave'});
miplot(xy1, xy2, pl);

%% Changel ylabel

% Overide the Y-labels
pl = plist('Arrangement', 'subplots', 'YLabels', {'signal 1', 'signal 2'}, 'Legends', {'', 'My Sine Wave'});
miplot(xy1, xy2, pl);

%% Changel xlabel

% Overide the X-labels
pl = plist('Arrangement', 'subplots', 'XLabels', {'', 'Time-stamps'}, 'Legends', {'', 'My Sine Wave'});
miplot(xy1, xy2, pl);

%% No legends

pl = plist('Arrangement', 'subplots', 'Legends', 'off');
miplot(xy1, xy2, pl);



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%         CDATA OBJECTS
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%% Make test AOS

c1  = ao(randn(1,10));
c2  = ao(1:10);


%% Default plot
miplot(c1, c2)

%% Change colors and line styles

pl = plist('Colors', {'g', 'k', 'm'}, 'LineStyles', {'--'}, 'LineWidths', [1 4]);
miplot(c1, c2, pl);

%% Change arrangement to single plots

pl = plist('Arrangement', 'single');
miplot(c1, c2, pl);

%% Change arrangement to subplots

% Also override the second legend text and the first line style
pl = plist('Arrangement', 'subplots', 'LineStyles', {'--'}, 'Legends', {'', 'My Sine Wave'});
miplot(c1, c2, pl);

%% Changel ylabel

% Overide the Y-labels
pl = plist('Arrangement', 'subplots', 'YLabels', {'signal 1', 'signal 2'}, 'Legends', {'', 'My Sine Wave'});
miplot(c1, c2, pl);

%% Changel xlabel

% Overide the X-labels
pl = plist('Arrangement', 'subplots', 'XLabels', {'', 'Time-stamps'}, 'Legends', {'', 'My Sine Wave'});
miplot(c1, c2, pl);

%% No legends

pl = plist('Arrangement', 'subplots', 'Legends', 'off');
miplot(c1, c2, pl);


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%         XYZ DATA OBJECTS
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%% Test data
nsecs = 100;
fs    = 100;
p     = plist('waveform', 'noise', 'fs', fs, 'nsecs', nsecs);
n     = ao(p);
p     = plist('waveform', 'chirp', 'fs', fs, 'nsecs', nsecs, 'f0', 0, 'f1', 50, 't1', nsecs);
s     = ao(p);

a = s+n;

% Make spectrogram

sxx = spectrogram(a, plist('Nfft', 4*fs));
rxx = spectrogram(s, plist('Nfft', 4*fs));

%% Plot

pl = plist('Arrangement', 'subplots', 'YMaths', {'log10(y)', 'log10(y)'}, ...
          'Zmaths', {'20*log10(z)', '20*log10(z)'}, ...
          'Xmaths', {'log10(x)'});

miplot(sxx, rxx, pl)


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%         MIXED OBJECTS
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%% XY and TSDATA
miplot(xy1, tsao1, xy2, tsao2, c1)