Hilbert+Huang 源代码

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

function A=myhilb(n, m) %MYHILB 本函数用来演示MATLAB语言的函数编写方法。 % A=MYHILB(N, M) 将产生一个N行M列的Hilbert矩阵A； % A=MYHILB(N) 将产生一个NxN的方Hil…

%样本点数 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=myhilb(n, m) %MYHILB 本函数用来演示MATLAB语言的函数编写方法。 % A=MYHILB(N, M) 将产生一个N行M列的Hilbert矩阵A； % A=MYHILB(N) 将产生一个NxN的方Hil…

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 seismic=s_hilbert(seismic,varargin) % Function computes Hilbert transform or instantaneous amplitude of seismic data % Wr…

% Chapter 7: Example 7.20 % Freq. Samp. Tech.: Hilbert Transformer % M = 51; alpha = (M-1)/2; Dw = 2*pi/M; l = 0:M-1; w…

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