--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/test/example_2.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,64 @@
+% A test script for the AO implementation.
+%
+% M Hewitson 08-02-07
+%
+% $Id: example_2.m,v 1.27 2009/02/02 15:20:38 hewitson Exp $
+%
+function example_2()
+  
+  
+  % Make test AOs
+  
+  nsecs = 10;
+  fs    = 100;
+  
+  pl = plist('nsecs', nsecs, 'fs', fs, 'tsfcn', 'randn(size(t))');
+  
+  a1 = ao(pl);
+  a2 = ao(pl);
+  a3 = ao(pl);
+  a4 = ao(pl);
+  a5 = ao(pl);
+  
+  % Subtract two data
+  a6 = a5 - a4;
+  a6.setName;
+  
+  % Make LPSD of each
+  
+  % Window function
+  w = specwin('Kaiser', 1000, 250);
+  
+  % parameter list for lpsd
+  pl = plist('Kdes', 100, 'Kmin', 10, 'Jdes', 1000, 'Win', w);
+  
+  % use lpsd
+  a7  = lpsd(a1, pl); a7.setName;
+  a8  = lpsd(a2, pl); a8.setName;
+  a9  = lpsd(a3, pl); a9.setName;
+  a10 = lpsd(a6, pl); a10.setName;
+  
+  % some manipulation and plot
+  a11 = a7+a8;
+  a12 = a9.*a10;
+  a13 = a11./a12;
+  a14 = sqrt(a13);
+  a15 = a14.^3 ./ a13;
+  
+  iplot(a15)
+  
+  % plot history
+  plot(a15.hist)
+  
+  save(a15, 'a15.xml');
+  
+  % Reproduce from history
+  
+  a_out = rebuild(a15);
+  
+  % the last object is always a1
+  iplot(a_out)
+  plot(a_out.hist)
+  
+  
+end
\ No newline at end of file