logogram.m

来自「老外写的小波变换的工具箱」· M 代码 · 共 29 行

function [t,g] = LogoGram(ser)
% LogoGram -- Time Series Spectrum Estimation
%  Usage
%    [t,g] = LogoGram(ser)
%  Inputs
%    ser      zero-mean time series, dyadic length
%  Outputs
%    t        ordinate of logogram
%    g        abscissa of logogram
%
n = length(ser);
Ij = abs(fft(ser)).^2 ./ ( n);
Yj = log(Ij) + .57721;
Yj(1) = Yj(1) - .57721 + (log(2) + .57721)/pi ;
Yj(n/2+1) = Yj(n/2+1) - .57721 + (log(2) + .57721)/pi ;
g = Yj .* sqrt(6)/pi;
t = 2 * pi * (0: (n-1)) ./ n;

    
    
 
 
%
%  Part of Wavelab Version 850
%  Built Tue Jan  3 13:20:42 EST 2006
%  This is Copyrighted Material
%  For Copying permissions see COPYING.m
%  Comments? e-mail [e-mail]