Hilbert+Huang 源代码

function [DH,DW] = hilbfilt(N, F, GF, W, delay) %HILBFILT Desired frequency response for Hilbert filters. % CREMEZ(N,F,'hilbfil…

function[A,B]=myhilb(n,m) % MYHILB 生成一个Hilbert矩阵 % [A,B]=myhilb(n,m) % where % n,m are size of the Hilbert matrix,if only one…

function [A,f,tt] = hhspectrum(imf,t,l,aff) % [A,f,tt] = HHSPECTRUM(imf,t,l,aff) computes the Hilbert-Huang spectrum % % inputs: …

% HSPEC: Hilbert Amplitude Spectrum % % [S,freq]=hspec(imf,N); % % S - Time-frequency-amplitude matrix % Columns are …

%%------和hspec.m不同的是，hspec1.m令其返回能量，频率和时间间隔 % HSPEC: Hilbert Amplitude Spectrum % % [S,freq]=hspec(imf,N); % % S - Time-fre…

function [A,f,tt] = hhspectrum(imf,t,l,aff) % [A,f,tt] = HHSPECTRUM(imf,t,l,aff) computes the Hilbert-Huang spectrum % % inputs: …

% HSPEC: Hilbert Amplitude Spectrum % % [S,freq]=hspec(imf,N); % % S - Time-frequency-amplitude matrix % Columns are …

%%------和hspec.m不同的是，hspec1.m令其返回能量，频率和时间间隔 % HSPEC: Hilbert Amplitude Spectrum % % [S,freq]=hspec(imf,N); % % S - Time-fre…

function xl=loweq(x,ts,f0) % xl=loweq(x,ts,f0) %LOWEQ returns the lowpass equivalent of the signal x % f0 is t…

function u=lssb_mod(m,ts,fc) % u=lssb_mod(m,ts,fc) %LSSB_MOD takes signal m sampled at ts and carrier % freq. fc as input a…