ltpda: testing/utp_1.1/utps/ao/utp_ao

annotate testing/utp_1.1/utps/ao/utp_ao_firwhiten.m @ 52:daf4eab1a51e database-connection-manager tip

Fix. Default password should be [] not an empty string

author	Daniele Nicolodi <nicolodi@science.unitn.it>
date	Wed, 07 Dec 2011 17:29:47 +0100
parents	409a22968d5e
children

rev	line source
44 409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	1 % UTP_AO_PSD a set of UTPs for the ao/firwhiten method
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	2 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	3 % M Hewitson 06-08-08
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	4 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	5 % $Id: utp_ao_firwhiten.m,v 1.21 2011/11/16 08:16:08 mauro Exp $
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	6 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	7
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	8 % <MethodDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	9 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	10 % The firwhiten method of the ao class computes the spectral density of time-series AOs.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	11 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	12 % </MethodDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	13
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	14 function results = utp_ao_firwhiten(varargin)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	15
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	16 % Check the inputs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	17 if nargin == 0
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	18
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	19 % Some keywords
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	20 class = 'ao';
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	21 mthd = 'firwhiten';
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	22
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	23 results = [];
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	24 disp('******************************************************');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	25 disp(['**** Running UTPs for ' class '/' mthd]);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	26 disp('******************************************************');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	27
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	28 % Test AOs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	29 [at1,at2,at3,at4,at5,at6,atvec,atmat] = eval(['get_test_objects_' class]);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	30
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	31 % Exception list for the UTPs:
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	32 [ple1,ple2,ple3,ple4,ple5,ple6] = get_test_ples();
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	33
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	34 % Run the tests
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	35 results = [results utp_01]; % getInfo call
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	36 results = [results utp_02]; % Vector input
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	37 results = [results utp_03]; % Matrix input
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	38 results = [results utp_04]; % List input
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	39 results = [results utp_05]; % Test with mixed input
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	40 results = [results utp_06]; % Test history is working
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	41 results = [results utp_07]; % Test the modify call works
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	42 results = [results utp_08]; % Test the data shape
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	43 results = [results utp_09]; % Test with complex plist
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	44 results = [results utp_10]; % Test the spectral falttening
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	45 results = [results utp_11(mthd, at1, ple1)]; % Test plotinfo doesn't disappear
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	46 results = [results utp_12(mthd, at1, ple1)]; % Test errors are cleared
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	47
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	48 disp('Done.');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	49 disp('******************************************************');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	50
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	51 elseif nargin == 1 % Check for UTP functions
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	52 if strcmp(varargin{1}, 'isutp')
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	53 results = 1;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	54 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	55 results = 0;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	56 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	57 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	58 error('### Incorrect inputs')
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	59 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	60
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	61 %% UTP_01
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	62
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	63 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	64 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	65 % Tests that the getInfo call works for this method.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	66 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	67 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	68 function result = utp_01
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	69
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	70
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	71 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	72 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	73 % Test that the getInfo call works for no sets, all sets, and each set
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	74 % individually.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	75 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	76 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	77
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	78 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	79 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	80 % Call for no sets
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	81 io(1) = eval([class '.getInfo(''' mthd ''', ''None'')']);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	82 % Call for all sets
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	83 io(2) = eval([class '.getInfo(''' mthd ''')']);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	84 % Call for each set
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	85 for kk=1:numel(io(2).sets)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	86 io(kk+2) = eval([class '.getInfo(''' mthd ''', ''' io(2).sets{kk} ''')']);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	87 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	88 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	89 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	90 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	91 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	92 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	93 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	94
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	95 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	96 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	97 % 1) Check that getInfo call returned an minfo object in all cases.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	98 % 2) Check that all plists have the correct parameters.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	99 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	100 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	101
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	102 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	103 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	104 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	105 % check we have minfo objects
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	106 if isa(io, 'minfo')
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	107 prefs = getappdata(0, 'LTPDApreferences');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	108 defaultWinType = char(prefs.getMiscPrefs.getDefaultWindow);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	109
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	110 %%% SET 'None'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	111 pn = 1;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	112 if ~isempty(io(pn).sets), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	113 if ~isempty(io(pn).plists), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	114 %%% Check all Sets
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	115 pn = 2;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	116 if ~any(strcmpi(io(pn).sets, 'Default')), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	117 if numel(io(pn).plists) ~= numel(io(pn).sets), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	118 %%%%%%%%%% SET 'Default'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	119 pn = 3;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	120 if io(pn).plists.nparams ~= 7, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	121 % Check key
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	122 if ~io(pn).plists.isparam('nfft'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	123 if ~io(pn).plists.isparam('bw'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	124 if ~io(pn).plists.isparam('hc'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	125 if ~io(pn).plists.isparam('win'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	126 if ~io(pn).plists.isparam('order'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	127 if ~io(pn).plists.isparam('firwin'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	128 if ~io(pn).plists.isparam('ntaps'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	129 % Check default value
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	130 if ~isequal(io(pn).plists.find('nfft'), -1), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	131 if ~isequal(io(pn).plists.find('bw'), 20), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	132 if ~isequal(io(pn).plists.find('hc'), .8), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	133 if ~strcmpi(io(pn).plists.find('win'), defaultWinType), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	134 if ~isequal(io(pn).plists.find('order'), -1), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	135 if ~strcmpi(io(pn).plists.find('win'), defaultWinType), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	136 if ~isequal(io(pn).plists.find('ntaps'), 256), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	137 % Check options
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	138 if ~isequal(io(pn).plists.getOptionsForParam('nfft'), {-1}), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	139 if ~isequal(io(pn).plists.getOptionsForParam('bw'), {20}), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	140 if ~isequal(io(pn).plists.getOptionsForParam('hc'), {.8}), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	141 if ~isequal(io(pn).plists.getOptionsForParam('win'), specwin.getTypes), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	142 if ~isequal(io(pn).plists.getOptionsForParam('order'), {-1 0 1 2 3 4 5 6 7 8 9}), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	143 if ~isequal(io(pn).plists.getOptionsForParam('firwin'), specwin.getTypes), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	144 if ~isequal(io(pn).plists.getOptionsForParam('ntaps'), {256}), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	145 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	146 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	147 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	148 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	149 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	150
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	151 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	152 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	153 end % END UTP_01
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	154
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	155 %% UTP_02
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	156
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	157 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	158 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	159 % Tests that the firwhiten method works with a vector of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	160 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	161 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	162 function result = utp_02
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	163
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	164 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	165 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	166 % Test that the firwhiten method works for a vector of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	167 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	168 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	169
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	170 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	171 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	172 avec = [at5 at5 at6];
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	173 N = 512;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	174 fwin = specwin('Hanning', N+1);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	175 [out, outf, outxx] = firwhiten(avec, plist('Ntaps', N, 'FIRwin', fwin));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	176 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	177 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	178 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	179 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	180 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	181 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	182
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	183 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	184 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	185 % 1) Check the number of elements in 'out'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	186 % 2) Check the number of filters (outf) and noise-floor estimates (outxx)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	187 % 3) Check that each output AO contains the correct data.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	188 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	189 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	190
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	191 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	192 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	193 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	194 % Check we have the correct number of outputs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	195 if numel(out) ~= numel(avec), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	196 if numel(outf) ~= numel(avec), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	197 if numel(outxx) ~= numel(avec), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	198 % Check the output data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	199 prefs = getappdata(0, 'LTPDApreferences');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	200 swin = char(prefs.getMiscPrefs.getDefaultWindow);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	201 for kk = 1:numel(avec)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	202 pl_psd = plist('Nfft', length(avec(kk).y), 'Win', specwin(swin), 'Order', -1, 'Scale', 'ASD');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	203 axx = psd(avec(kk), pl_psd);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	204 nxx = smoother(axx, plist('width', 20, 'hc', 0.8));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	205 w = 1./nxx;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	206 ff = mfir(w, plist('Win', fwin, 'N', N));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	207 res = filter(avec(kk), ff);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	208 if ~isequal(res.x, out(kk).x), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	209 if ~isequal(res.y, out(kk).y), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	210 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	211 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	212 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	213 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	214 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	215
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	216 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	217 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	218 end % END UTP_02
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	219
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	220 %% UTP_03
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	221
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	222 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	223 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	224 % Tests that the firwhiten method works with a matrix of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	225 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	226 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	227 function result = utp_03
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	228
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	229 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	230 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	231 % Test that the firwhiten method works for a matrix of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	232 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	233 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	234
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	235 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	236 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	237 amat = [at1 at5 at6; at5 at6 at1];
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	238 N = 512;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	239 fwin = specwin('Hanning', N+1);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	240 [out, outf, outxx] = firwhiten(amat, plist('Ntaps', N, 'FIRwin', fwin));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	241 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	242 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	243 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	244 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	245 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	246 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	247
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	248 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	249 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	250 % 1) Check the number of elements in 'out'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	251 % 2) Check the number of filters (outf) and noise-floor estimates (outxx)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	252 % 3) Check that each output AO contains the correct data.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	253 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	254 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	255
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	256 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	257 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	258 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	259 % Check we have the correct number of outputs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	260 if numel(out) ~= numel(amat), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	261 if numel(outf) ~= numel(amat), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	262 if numel(outxx) ~= numel(amat), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	263 % Check the output data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	264 prefs = getappdata(0, 'LTPDApreferences');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	265 swin = char(prefs.getMiscPrefs.getDefaultWindow);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	266 for kk = 1:numel(amat)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	267 pl_psd = plist('Nfft', length(amat(kk).y), 'Win', specwin(swin), 'Order', -1, 'Scale', 'ASD');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	268 axx = psd(amat(kk), pl_psd);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	269 nxx = smoother(axx, plist('width', 20, 'hc', 0.8));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	270 w = 1./nxx;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	271 ff = mfir(w, plist('Win', fwin, 'N', N));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	272 res = filter(amat(kk), ff);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	273 if ~isequal(res.x, out(kk).x), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	274 if ~isequal(res.y, out(kk).y), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	275 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	276 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	277 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	278 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	279 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	280
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	281 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	282 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	283 end % END UTP_03
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	284
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	285 %% UTP_04
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	286
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	287 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	288 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	289 % Tests that the firwhiten method works with a list of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	290 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	291 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	292 function result = utp_04
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	293
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	294 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	295 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	296 % Test that the firwhiten method works for a list of AOs as input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	297 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	298 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	299
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	300 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	301 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	302 N = 512;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	303 fwin = specwin('Hanning', N+1);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	304 [out, outf, outxx] = firwhiten(at1,at5,at6, plist('Ntaps', N, 'FIRwin', fwin));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	305 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	306 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	307 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	308 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	309 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	310 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	311
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	312 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	313 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	314 % 1) Check the number of elements in 'out'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	315 % 2) Check the number of filters (outf) and noise-floor estimates (outxx)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	316 % 3) Check that each output AO contains the correct data.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	317 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	318 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	319
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	320 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	321 aoin = [at1, at5, at6];
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	322 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	323 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	324 % Check we have the correct number of outputs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	325 if numel(out) ~= 3, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	326 if numel(outf) ~= 3, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	327 if numel(outxx) ~= 3, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	328 % Check the output data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	329 prefs = getappdata(0, 'LTPDApreferences');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	330 swin = char(prefs.getMiscPrefs.getDefaultWindow);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	331 for kk = 1:numel(aoin)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	332 pl_psd = plist('Nfft', length(aoin(kk).y), 'Win', specwin(swin), 'Order', -1, 'Scale', 'ASD');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	333 axx = psd(aoin(kk), pl_psd);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	334 nxx = smoother(axx, plist('width', 20, 'hc', 0.8));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	335 w = 1./nxx;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	336 ff = mfir(w, plist('Win', fwin, 'N', N));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	337 res = filter(aoin(kk), ff);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	338 if ~isequal(res.x, out(kk).x), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	339 if ~isequal(res.y, out(kk).y), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	340 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	341 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	342 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	343 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	344 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	345
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	346 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	347 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	348 end % END UTP_04
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	349
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	350 %% UTP_05
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	351
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	352 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	353 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	354 % Tests that the firwhiten method works with a mix of different shaped AOs as
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	355 % input.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	356 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	357 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	358 function result = utp_05
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	359
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	360 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	361 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	362 % Test that the firwhiten method works with an input of matrices and vectors
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	363 % and single AOs.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	364 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	365 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	366
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	367 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	368 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	369 N = 512;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	370 fwin = specwin('Hanning', N+1);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	371 pl = plist('Ntaps', N, 'FIRwin', fwin);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	372 [out, outf, outxx] = firwhiten(at1,[at5 at6],at5,[at5 at1; at6 at1],at6, pl);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	373 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	374 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	375 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	376 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	377 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	378 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	379
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	380 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	381 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	382 % 1) Check the number of elements in 'out'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	383 % 2) Check the number of filters (outf) and noise-floor estimates (outxx)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	384 % 3) Check that each output AO contains the correct data.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	385 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	386 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	387
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	388 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	389 aoin = [at1, reshape([at5 at6], 1, []), at5, reshape([at5 at1; at6 at1], 1, []), at6];
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	390 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	391 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	392 % Check we have the correct number of outputs
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	393 if numel(out) ~= 9, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	394 if numel(outf) ~= 9, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	395 if numel(outxx) ~= 9, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	396 % Check the output data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	397 prefs = getappdata(0, 'LTPDApreferences');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	398 swin = char(prefs.getMiscPrefs.getDefaultWindow);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	399 for kk = 1:numel(aoin)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	400 pl_psd = plist('Nfft', length(aoin(kk).y), 'Win', specwin(swin), 'Order', -1, 'Scale', 'ASD');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	401 axx = psd(aoin(kk), pl_psd);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	402 nxx = smoother(axx, plist('width', 20, 'hc', 0.8));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	403 w = 1./nxx;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	404 ff = mfir(w, plist('Win', fwin, 'N', N));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	405 res = filter(aoin(kk), ff);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	406 if ~isequal(res.x, out(kk).x), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	407 if ~isequal(res.y, out(kk).y), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	408 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	409 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	410 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	411 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	412 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	413
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	414 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	415 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	416 end % END UTP_05
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	417
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	418 %% UTP_06
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	419
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	420 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	421 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	422 % Tests that the firwhiten method properly applies history.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	423 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	424 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	425 function result = utp_06
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	426
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	427 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	428 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	429 % Test that the result of applying the firwhiten method can be processed back
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	430 % to an m-file.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	431 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	432 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	433
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	434 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	435 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	436 N = 512;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	437 fwin = specwin('Hanning', N+1);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	438 pl = plist('Ntaps', N, 'FIRwin', fwin);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	439 out = firwhiten(at5, pl);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	440 mout = rebuild(out);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	441 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	442 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	443 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	444 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	445 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	446 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	447
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	448 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	449 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	450 % 1) Check that the last entry in the history of 'out' corresponds to
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	451 % 'firwhiten'.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	452 % 2) Check that the re-built object is the same object as 'out'.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	453 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	454 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	455
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	456 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	457 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	458 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	459 % Check the last step in the history of 'out'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	460 if ~strcmp(out.hist.methodInfo.mname, 'firwhiten'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	461 % Check the re-built object
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	462 if ~eq(mout, out, ple4), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	463 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	464 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	465 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	466 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	467
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	468 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	469 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	470 end % END UTP_06
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	471
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	472 %% UTP_07
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	473
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	474 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	475 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	476 % Tests that the firwhiten method can modify the input AO.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	477 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	478 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	479 function result = utp_07
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	480
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	481 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	482 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	483 % Test that the firwhiten method can modify the input AO by calling
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	484 % with no output.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	485 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	486 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	487
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	488 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	489 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	490 N = 512;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	491 fwin = specwin('Hanning', N+1);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	492 pl = plist('Ntaps', N, 'FIRwin', fwin);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	493 % copy at1 to work with
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	494 ain = ao(at1);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	495 % modify ain
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	496 aout = ain.firwhiten(pl);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	497 ain.firwhiten(pl);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	498 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	499 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	500 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	501 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	502 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	503 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	504
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	505 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	506 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	507 % 1) Check that 'at1' and 'ain' are now different.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	508 % 2) Check that 'ain' is firwhiten(at1).
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	509 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	510 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	511
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	512 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	513 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	514 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	515 % Check that firwhiten modified the input by comparing to the copy
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	516 if eq(ao(at1), ain, ple1), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	517 % Check that firwhiten doesn't modified the input for the function notation
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	518 if ~eq(aout, ain, ple1), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	519 % Check that the modified input is the firwhiten of the copy
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	520 % Check the output data of ain
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	521 prefs = getappdata(0, 'LTPDApreferences');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	522 swin = char(prefs.getMiscPrefs.getDefaultWindow);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	523 pl_psd = plist('Nfft', length(at1.y), 'Win', specwin(swin), 'Order', -1, 'Scale', 'ASD');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	524 axx = psd(at1, pl_psd);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	525 nxx = smoother(axx, plist('width', 20, 'hc', 0.8));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	526 w = 1./nxx;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	527 ff = mfir(w, plist('Win', fwin, 'N', N));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	528 res = filter(at1, ff);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	529 if ~isequal(res.x, ain.x), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	530 if ~isequal(res.y, ain.y), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	531 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	532 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	533 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	534 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	535
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	536 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	537 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	538 end % END UTP_07
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	539
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	540 %% UTP_08
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	541
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	542 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	543 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	544 % Tests that the firwhiten method keeps the data shape of the input object.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	545 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	546 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	547 function result = utp_08
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	548
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	549 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	550 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	551 % Test that the firwhiten method keeps the data shape of the input object. The
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	552 % input AO must be an AO with row data and an AO with column data.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	553 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	554 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	555
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	556 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	557 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	558 N = 512;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	559 fwin = specwin('Hanning', N+1);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	560 pl = plist('Ntaps', N, 'FIRwin', fwin);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	561 out1 = firwhiten(at5, pl);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	562 out2 = firwhiten(at6, pl);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	563 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	564 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	565 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	566 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	567 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	568 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	569
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	570 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	571 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	572 % 1) Check that the shpe of the data doesn't change.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	573 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	574 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	575
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	576 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	577 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	578 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	579 % Check the shape of the output data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	580 if size(out1.data.y,2) ~= 1, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	581 if size(out2.data.y,1) ~= 1, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	582 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	583 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	584 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	585 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	586
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	587 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	588 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	589 end % END UTP_08
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	590
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	591 %% UTP_09
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	592
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	593 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	594 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	595 % Tests that the firwhiten method with a complex plist.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	596 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	597 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	598 function result = utp_09
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	599
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	600 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	601 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	602 % Test that the result of applying the firwhiten method with a complex plist
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	603 % can be processed back to a m-file.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	604 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	605 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	606
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	607 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	608 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	609 N = 512;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	610 Nfft = 100;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	611 order = 1;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	612 fwin = specwin('Hanning', N+1);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	613 swin = specwin('Kaiser', 10, 100);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	614 pl = plist('Ntaps', N, 'FIRwin', fwin, 'win', swin, 'Nfft', Nfft, 'order', order);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	615 [out, outf, outxx] = firwhiten(at5, pl);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	616 mout = rebuild(out);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	617 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	618 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	619 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	620 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	621 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	622 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	623
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	624 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	625 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	626 % 1) Check the output data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	627 % 2) Check the output filter
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	628 % 3) Check the noise-floor estimates
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	629 % 4) Check that the re-built object is the same object as 'out'.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	630 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	631 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	632
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	633 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	634 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	635 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	636 % Check the last step in the history of 'out'
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	637 if ~strcmp(out.hist.methodInfo.mname, 'firwhiten'), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	638 % Check the re-built object
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	639 if ~eq(mout, out, ple2), atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	640
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	641 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	642 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	643 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	644 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	645
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	646 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	647 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	648 end % END UTP_09
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	649
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	650 %% UTP_10
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	651
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	652 % <TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	653 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	654 % Test the spectral falttening capability of firwhiten method.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	655 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	656 % </TestDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	657 function result = utp_10
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	658
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	659 % <SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	660 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	661 % Test that the application of the firwhiten method enhances the
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	662 % spectral flatness of input data.
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	663 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	664 % </SyntaxDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	665
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	666 try
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	667 % <SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	668 % Making test data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	669 fs = at5.fs;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	670 pzm = pzmodel(2, [pz(0.1, 2) pz(0.5, 1)] , [pz(1) pz(2.55, 2)]);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	671 ft = miir(pzm,plist('fs',fs));
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	672 af = filter(at5, ft); % Colored noise
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	673
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	674 % Whitening
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	675 N = 512;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	676 Nfft = 100;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	677 order = 1;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	678 fwin = specwin('Hanning', N+1);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	679 swin = specwin('Kaiser', 10, 100);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	680 pl = plist('Ntaps', N, 'FIRwin', fwin, 'win', swin, 'Nfft', Nfft, 'order', order);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	681 out = firwhiten(af, pl);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	682
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	683 % </SyntaxCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	684 stest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	685 catch err
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	686 disp(err.message)
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	687 stest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	688 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	689
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	690 % <AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	691 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	692 % 1) Calculate PSD of input and whitened data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	693 % 2) Compare relative spectral flatness coefficients
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	694 %
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	695 % </AlgoDescription>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	696
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	697 atest = true;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	698 if stest
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	699 % <AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	700 % Check the output data
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	701 % pl_psd = plist('Nfft', Nfft, 'Win', swin, 'Order', order, 'Scale', 'PSD');
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	702 axx = af.psd;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	703 awxx = out.psd;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	704
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	705 % Claculating flatness
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	706 sf1 = utils.math.spflat(axx.data.y);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	707 sf2 = utils.math.spflat(awxx.data.y);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	708
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	709 % Checking flatness
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	710 if sf1>sf2, atest = false; end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	711
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	712 % </AlgoCode>
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	713 else
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	714 atest = false;
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	715 end
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	716
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	717 % Return a result structure
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	718 result = utp_prepare_result(atest, stest, dbstack, mfilename);
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	719 end % END UTP_10
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	720
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	721
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	722
409a22968d5e Add unit tests Daniele Nicolodi <nicolodi@science.unitn.it> parents: diff changeset	723 end

Mercurial > hg > ltpda

annotate testing/utp_1.1/utps/ao/utp_ao_firwhiten.m @ 52:daf4eab1a51e database-connection-manager tip