+
−mc
+
−
+
−
+
−%% Make test AOs
+
−
+
−nsecs = 10000;
+
−fs    = 10;
+
−
+
−pl = plist('nsecs', nsecs, 'fs', fs, 'tsfcn', 'sin(2*pi*7.433*t) + randn(size(t))');
+
−
+
−a1 = ao(pl);
+
−a1.setName('a1');
+
−
+
−%% Make two filters
+
−
+
−lp = miir(plist('type', 'lowpass', 'fs', fs, 'fc', 0.1, 'order', 5));
+
−hp = miir(plist('type', 'highpass', 'fs', fs, 'fc', 1, 'order', 2));
+
−
+
−%% Serial
+
−
+
−a1lp = filter(a1, lp);
+
−a1lphp = filter(a1lp, hp);
+
−
+
−a1f = filter(a1, [lp hp], plist('bank', 'serial'));
+
−
+
−iplot(a1lphp.psd, a1f.psd, plist('LineStyles', {'', '--'}))
+
−
+
−%% Parallel
+
−
+
−a1lp = filter(a1, lp);
+
−a1lphp = a1lp + filter(a1, hp);
+
−
+
−a1f = filter(a1, [lp hp]);
+
−
+
−iplot(a1lphp.psd, a1f.psd, plist('LineStyles', {'', '--'}))
+
−
+
−%% Test filter history serial
+
−
+
−nsecs = 100;
+
−fs    = 10;
+
−
+
−pl = plist('nsecs', nsecs, 'fs', fs, 'tsfcn', 'sin(2*pi*0.433*t)');
+
−
+
−a1 = ao(pl);
+
−a1.setName('a1');
+
−
+
−segs = split(a1, plist('split_type', 'chunks', 'n', 2));
+
−
+
−filts = [lp hp];
+
−
+
−% full
+
−[a1lps, filts] = filter(a1, filts, plist('bank', 'serial'));
+
−
+
−% segement 1
+
−[a1lps1, filts] = filter(segs(1), filts, plist('bank', 'serial'));
+
−
+
−% segement 2
+
−a1lps2 = filter(segs(2), filts, plist('bank', 'serial'));
+
−
+
−iplot(a1lps1, a1lps2, a1lps, plist('Xranges', [49 51], 'markers', {'o', 'o', ''}))
+
−
+
−%% Test filter history parallel
+
−
+
−nsecs = 100;
+
−fs    = 10;
+
−
+
−pl = plist('nsecs', nsecs, 'fs', fs, 'tsfcn', 'sin(2*pi*0.433*t)');
+
−
+
−a1 = ao(pl);
+
−a1.setName('a1');
+
−
+
−segs = split(a1, plist('split_type', 'chunks', 'n', 2));
+
−
+
−filts = [lp hp];
+
−
+
−% full
+
−[a1lps, filts] = filter(a1, filts);
+
−
+
−% segement 1
+
−[a1lps1, filts] = filter(segs(1), filts);
+
−
+
−% segement 2
+
−a1lps2 = filter(segs(2), filts);
+
−
+
−iplot(a1lps1, a1lps2, a1lps, plist('Xranges', [49 51], 'markers', {'o', 'o', ''}))
+
−
+
−%% Test Serial with different units
+
−
+
−a1.setYunits('m');
+
−
+
−% This should output N
+
−iunits1 = unit('m');
+
−ounits1 = unit('V');
+
−iunits2 = unit('V');
+
−ounits2 = unit('N');
+
−
+
−lp = miir(plist('type', 'lowpass', 'fs', fs, 'fc', 0.1, 'order', 5,'iunits', iunits1, 'ounits', ounits1));
+
−hp = miir(plist('type', 'highpass', 'fs', fs, 'fc', 1, 'order', 2,'iunits', iunits2, 'ounits', ounits2));
+
−
+
−a1f = filter(a1, [lp hp], plist('bank', 'serial'));
+
−
+
−%% Test Parallel with different units
+
−
+
−% This must not to be allowed
+
−iunits1 = unit('m');
+
−ounits1 = unit('V');
+
−iunits2 = unit('m');
+
−ounits2 = unit('V');
+
−
+
−lp = miir(plist('type', 'lowpass', 'fs', fs, 'fc', 0.1, 'order', 5,'iunits', iunits1, 'ounits', ounits1));
+
−hp = miir(plist('type', 'highpass', 'fs', fs, 'fc', 1, 'order', 2,'iunits', iunits2, 'ounits', ounits2));
+
−
+
−a1f = filter(a1, [lp hp], plist('bank', 'parallel'));
+
−
+
−