--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/test/test_ao_noisegen2D_2.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,116 @@
+% A test script for ao/noisegen2D
+% test ao/noisegen2D with input the model and output the filter
+% 
+% DESCRIPTION: Run noisegen2D with test data and test procedure accuracy
+% 
+% L. Ferraioli 10-11-08
+% 
+% $Id: test_ao_noisegen2D_2.m,v 1.1 2010/05/03 19:04:45 luigi Exp $
+% 
+
+%% General use variables and vectors
+
+userdir = 'C:\Users\Luigi'; % You should set your own dir
+
+f = logspace(-6,log10(5),300);
+fs = 10;
+
+Nt = fs*10000;
+
+
+%% CSD Noise models
+
+fundir = fullfile(userdir,'ltp_data_analysis\MDCs\MDC1_UTN');
+cf = cd;
+cd(fundir);
+[TF,CSD] = mdc1_tf_models(plist('f',f,'fs',fs));
+cd(cf);
+
+
+%% Plot Model Tfs
+
+iplot(TF(1,1),TF(1,2),TF(2,1),TF(2,2))
+
+%% Plot Model CSD
+
+iplot(CSD(1,1),CSD(1,2),CSD(2,1),CSD(2,2))
+
+%% Noise generation - 3 inputs
+
+pl = plist(...
+    'fs', fs, ...
+    'Iunits', '', ...
+    'Ounits', 'm', ...
+    'MaxIter', 80, ...
+    'PoleType', 2, ...
+    'MinOrder', 15, ...
+    'MaxOrder', 45, ...
+    'Weights', 3, ...
+    'FITTOL', 1e-3,...
+    'MSEVARTOL', 1e-1,...
+    'UseSym', 0,...
+    'Plot', false,...
+    'Disp', false);
+
+fil = noisegen2D(CSD(1,1),CSD(1,2),CSD(2,2), pl);
+
+%% filtering data
+
+a1 = ao(plist('tsfcn', 'randn(size(t))', 'fs', 10, 'nsecs', 1e5));
+a2 = ao(plist('tsfcn', 'randn(size(t))', 'fs', 10, 'nsecs', 1e5));
+
+a = matrix(a1,a2,plist('shape',[2 1]));
+
+b = filter(a,fil);
+
+%% Checking results and starting data
+
+iplot(b.objs)
+
+%% Make psd
+
+plpsd = plist('navs',8,'order',1,'olap',50);
+bxx1 = b.objs(1).psd(plpsd);
+bxx2 = b.objs(2).psd(plpsd);
+
+iplot(bxx1,CSD(1,1),bxx2,CSD(2,2))
+
+%% Noise generation - 4 inputs
+
+pl = plist(...
+    'fs', fs, ...
+    'Iunits', '', ...
+    'Ounits', 'm', ...
+    'MaxIter', 80, ...
+    'PoleType', 2, ...
+    'MinOrder', 15, ...
+    'MaxOrder', 45, ...
+    'Weights', 3, ...
+    'FITTOL', 1e-3,...
+    'MSEVARTOL', 1e-1,...
+    'UseSym', 0,...
+    'Plot', false,...
+    'Disp', false);
+
+fil = noisegen2D(CSD(1,1),CSD(1,2),CSD(2,1),CSD(2,2), pl);
+
+%% filtering data
+
+a1 = ao(plist('tsfcn', 'randn(size(t))', 'fs', 10, 'nsecs', 1e5));
+a2 = ao(plist('tsfcn', 'randn(size(t))', 'fs', 10, 'nsecs', 1e5));
+
+a = matrix(a1,a2,plist('shape',[2 1]));
+
+b = filter(a,fil);
+
+%% Checking results and starting data
+
+iplot(b.objs)
+
+%% Make psd
+
+plpsd = plist('navs',8,'order',1,'olap',50);
+bxx1 = b.objs(1).psd(plpsd);
+bxx2 = b.objs(2).psd(plpsd);
+
+iplot(bxx1,CSD(1,1),bxx2,CSD(2,2))
\ No newline at end of file