+
−% Model and generation
+
−
+
−
+
−
+
−fs = 1;
+
−
+
−p1 = pzmodel(1e-8,[pz(1e-3)  pz(0.1, 4)],...
+
−  [ pz(2e-1) pz(3e-1)]);
+
−
+
−pl3 = plist();
+
−pl3 = append(pl3, param('pzmodel', p1));
+
−pl3 = append(pl3, param('fs',fs));
+
−pl3 = append(pl3, param('nsecs', 3e3));
+
−
+
−a = ao(pl3);
+
−%% psd
+
−
+
−pl = plist(...
+
−  'Scale','ASD',...
+
− ... 'Win',  getappdata(0, 'ltpda_default_spectral_window'), ...
+
−  'Nfft', 1000,...
+
−  'Order', 1);
+
−
+
−% pl = plist('Kdes', 100,...
+
−%   'Jdes', 100000, ...
+
−%   'Order', 1, ...
+
−%   'Scale', 'ASD');
+
−
+
−ca = psd(a,pl);
+
−
+
−%% plot
+
−
+
−rp1 = resp(p1,plist('f1',1e-3,'f2',100,'nf',1000));
+
−
+
−iplot(rp1,ca)
+
−
+
−%% gaps
+
−by = mod((1:length(a.data.y))/10, 20)>19;
+
−b = ao(tsdata(double(~by), 10));
+
−c = a.*b.';
+
−
+
−
+
−%%
+
−d = c.gapfillingoptim(c, p1);
+
−%%
+
−iplot(psd(a, pl), psd(c, pl), psd(d, pl), resp(p1))
+
−iplot(a,c,d)
+
−
+
−