Mercurial > hg > ltpda
view src/ltpda_dft/test_ltpda_lpsd_new.m @ 10:75007001cbfe database-connection-manager
Check for binary only objects
author | Daniele Nicolodi <nicolodi@science.unitn.it> |
---|---|
date | Mon, 05 Dec 2011 16:20:06 +0100 |
parents | f0afece42f48 |
children |
line wrap: on
line source
% 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'))