+
−mc;
+
−
+
−% Compute time-series
+
−Nsecs = 1e3; % Number of seconds desired
+
−fs    = 10;
+
−
+
−s = ao(plist('tsfcn', 'sin(2*pi*1.137*t)', 'fs', fs, 'Nsecs', Nsecs));
+
−n = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'Nsecs', Nsecs));
+
−b = n + s;
+
−
+
−% Make FFT
+
−tic
+
−bxx = fft(b);
+
−toc
+
−
+
−% Make DFT
+
−tic
+
−dxx = dft(b);
+
−toc
+
−
+
−% Make DFT at a particular frequency
+
−dsx = dft(b, plist('f', 1.137));
+
−
+
−% Plot
+
−iplot(bxx, dxx, dsx, plist('Markers', {'none', 'none', 'o'}));
+
−
+
−iplot(bxx./dxx)