S变换 matlab源码

 
function [sttf] = st（timeseries）

TRUE = 1;
FALSE = 0;
%%% DEFAULT PARAMETERS  [change these for your particular application]
verbose = TRUE;          
removeedge= FALSE;
analytic_signal =  FALSE;
factor = 1;

% use defaults for input variables
% calculate the sampled time and frequency values from the two sampling rates
t = （0:500-1）*1;
spe_nelements =ceil（（250 - 0+1）/1）   ;
f = （0 + [0:spe_nelements-1]*1）/（1*500）;

% The actual S Transform function is here:
st = strans（timeseries025011verboseremoveedgeanalytic_signalfactor）; 
% this function is below thus nicely encapsulated

%WRITE switch statement on nargout
% if 0 then plot amplitude spectrum
return


%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^


function st = strans（timeseriesminfreqmaxfreqsamplingratefreqsamplingrateverboseremoveedgeanalytic_signalfactor）; 
% Returns the Stockwell Transform STOutput of the time-series
% Code by R.G. Stockwell.
% Reference is “Localization of the Complex Spectrum: The S Transform“
% from IEEE Transactions on Signal Processing vol. 44. number 4
% April 1996 pages 998-1001.
%
%-------Inputs Returned------------------------------------------------
%         - are all taken care