140 likes | 263 Views
STORM WP7 Software library for nonlinear analysis of fluctuations in space plasma time series. P. Kovács, A. Koppán C . Munteanu , M. Echim , STORM Annual Meeting, Graz, 25-26 November, 2013. Outline. STORM related Matlab functions ASCII data import tool
E N D
STORM WP7Software libraryfornonlinearanalysis of fluctuationsinspaceplasmatime series P. Kovács, A. Koppán C. Munteanu, M. Echim, STORM Annual Meeting, Graz, 25-26 November, 2013
Outline • STORM related Matlab functions • ASCII data import tool • CDF file import tool • Database of thedeliverables (??) • STORM Annual Meeting, Graz, 25-26 November, 2013
STORMMatlab functions PSD Psd_Pwelch.m Psd_Plot.m Kurt_vs_Time.m Pdf_plot.m STR_Func.m Pmodel_Fit.m STR_plot.m PowerLFit.m FigParam_Save.m PLFit_vs_Time.m LogMean.m Increment.m GapFill.m • function [Pxx, f] = Psd_Pwelch(indata, segments, overlap, fs , varargin) • % Psd_Pwelch ; Computes the Power Spectral Density (PSD) function of • % vector "indata" using the Welch algorithm (see Welch, 1967 or Matlab • % Documentation). By default, the window segments are multiplied by Hamming • % window of length given as floor(length(indata)/segments). However, • % Gaussian, Hann and rectangular alternative windows can also be selected. • % • % [Pxx, f] = Psd_Pwelch(indata, segments, overlap, fs) • % • % INPUTS: • % indata = input time series • % segments = number of overlapping segments in which separate PSDs are • % computed (Final PSD is given by the mean of the segment PSDs) • % overlap = overlap of the segments in per cent ([0 100]) • % fs = sapling frequency in Hz • % varargin = Type of window with which the segments of input time-series • % are multiplied before the PSD computation. Its value can be 'rectwin' • % for rectangular, 'gausswin' for Gaussian or 'hann' foHann windows. If • % varargin is not given, the Hamming window is used, by default. • % • % OUTPUTS: • % [Pxx f] = Frequency, f, and the corresponding PSD, Pxx, of the input • % time-series. (The first element of the spectrum corresponding to f = • % 0 Hz is discarded) • % • % NOTE: • % The conventional PSD function can be obtained by setting segments = 1 • % and overlap = 0 . • % • % REQUIRED MATLAB TOOLBOX: • % - Signal Processing Toolbox • window = floor(length(indata)/segments); • L_overlap = floor(window*overlap/100); • if(numel(varargin)) • switch(varargin{1}) • case'rectwin' • window = rectwin(window); • case'gausswin' • window = gausswin(window); • case'hann' • window = hann(window); • end • end • [Pxx, f] = pwelch(indata(:,2), window, L_overlap, [], fs); • Pxx = Pxx(2:length(Pxx)); • f = f(2:length(f)); PDF STR General
ASCII Data Import Possible Input timeformats: yy mm dd HH MM SS yyDoY HH MM SS yyyy-mm-ddTHH:MM:SS
ASCII Data Import data14354×5 double filename ‚u01010minsh.dat’ header5×4 char unit 4×2 char
CDF File Import TIME_RESOLUTION ??? Time_resolution ??? time_resolution
CDF File Import – Gaphandling FillMethods LinearInterpolation Fillwithvalue Resample STORM Annual Meeting, Graz, 25-26 November, 2013
Unevendatasampling STORM Annual Meeting, Graz, 25-26 November, 2013
Database of deliverables STORM Annual Meeting, Graz, 25-26 November, 2013
Browsinginthedatabase SearchOptions: Mission Instrument Data Type Period