Hilbert+Huang 源代码

%样本点数 N=4096; %正常时心音信号 [x,fxs,bits] = wavread('normal',N); %肺动脉瓣狭窄时心音信号 [y,fys,bits] = wavread('pulmonarystenosis',N); %设置…

%产生非平稳信号 sig=real(amgauss(128).*fmlin(128)); %%时域波形 figure(1) plot(sig,'LineWidth',2); xlabel('时间 t'); ylabel('幅值 A'); %设置窗…

%产生非平稳信号 sig=real(amgauss(128).*fmlin(128)); %%时域波形 figure(1) plot(sig,'LineWidth',2); xlabel('时间 t'); ylabel('幅值 A'); %设置窗…

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 …

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 x = hilbert(xr,n) if nargin0 & 2*fix(n/2)==n h([1 n/2+1]) = 1; h(2:n/2) = 2; elseif n>0 h(1) = 1; h(2…

function [hf, ha] = fah(data,dt) % The function FAH calculates an improved Hilbert frequency and amplitude % of data(n,k), whe…

function [nhf, nha, hex, hE, hN] = nhfam5(data,dt) % The function NHFAM5 calculates the smoothed modified Hilbert frequency a…