% MATLAB version '9.7.0.1190202 (R2019b)' load('../data/time_series.mat') nw = 4; fs = 1000; [pxx, f] = pmtm(time_series, nw, [], fs); save('../data/matlab_multitaper', 'pxx', 'f', 'fs', 'nw');