--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/src/ltpda_dft/test_ltpda_lpsd_new.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,69 @@
+% function test_ltpda_lpsd_new()
+% A test script for the AO implementation of lpsd.
+% 
+% M Hewitson 02-02-07
+% 
+% $Id: test_ltpda_lpsd_new.m,v 1.6 2008/02/05 13:33:11 hewitson Exp $
+% 
+
+clear all;
+
+%% Make test AOs
+
+%% poles and zeros
+p1 = pole(0.1, 2);
+p2 = pole(0.5, 3);
+p3 = pole(1);
+
+z1 = zero(0.5, 30);
+z2 = zero(1, 2);
+
+gain =1; 
+%% convert to get a and b coefficients
+%pole zero model
+pzm = pzmodel(gain, [p1 p2 p3], [z1 z2]);
+% parameter list for ltpda_noisegen
+pl = plist();
+pl = append(pl, param('nsecs', 100000));
+pl = append(pl, param('fs', 10));
+pl = append(pl, param('pzmodel', pzm));
+%% calling the noisegenerator
+[a1, pl1, pl2]  = ltpda_noisegen(pl);
+a1 = set(a1, 'yunits', 'V');
+
+%% Make LPSD of each
+
+% Window function
+w = specwin('Kaiser', 10, 150);
+
+% parameter list for lpsd
+pl = plist();
+pl = append(pl, param('Kdes', 100));
+pl = append(pl, param('Kmin', 2));
+pl = append(pl, param('Jdes', 1000));
+pl = append(pl, param('Win', w));
+pl = append(pl, param('Order', 1));
+
+plm = plist('M-FILE ONLY', 1);
+
+%% use old lpsd
+tic
+a4 = ltpda_lpsd(a1, pl, plm);
+told = toc;
+
+
+%% use new lpsd
+tic
+a5 = ltpda_lpsd(a1, pl);
+tnew = toc;
+
+sprintf('#### Old method: %f s', told)
+sprintf('#### New method: %f s', tnew)
+
+told/tnew
+
+%% add and plot
+iplot(a4, a5, plist('Legends', {'old','new'}))
+
+iplot(a4./ a5, plist('YScales', 'lin'))
+