可以运行的EMD工具箱针对hhspectrum错误-instfreq.m

可以运行的EMD工具箱针对hhspectrum错误-instfreq.m
在论坛上载了那个最经典的EMD，但是运行后发现出现错误。但是好像很多人也只是随便发，不管是否可以运行，可能仅仅是为了挣麦片。
与同学研究了下，发现是缺少了一个instfreq的函数，于是就从网站上载了一个函数给有需要的人参考下。

EMD中文件夹的install不用管它。可以用一个例子测试下
t=1:500;
t=t*1/2000;
x=sin;
y=0.4*sin;
z=x y;
plot
imf=emd;
emd_visu,imf)
[A,f,tt]=hhspectrum;
[im,tt]=toimage;
disp_hhs;
colormap)


如果单纯用EMD包，运行到[A,f,tt]=hhspectrum;就会出错，hhspectrum函数需要有instfreq函数包才行。



希望对大家有帮助，为了方便，也提供了一个emd instfreq函数的可以下载

function [fnormhatt]=instfreq（xtLtrace）;
%INSTFREQ Instantaneous frequency estimation.
%	[FNORMHATT]=INSTFREQ（XTLTRACE） computes the instantaneous 
%	frequency of the analytic signal X at time instant（s） T using the
%	trapezoidal integration rule.
%	The result FNORMHAT lies between 0.0 and 0.5.
% 
%	X : Analytic signal to be analyzed.
%	T : Time instants	        （default : 2:length（X）-1）.
%	L : If L=1 computes the （normalized） instantaneous frequency 
%	    of the signal X defined as angle（X（T+1）*conj（X（T-1）） ;
%	    if L>1 computes a Maximum Likelihood estimation of the
%	    instantaneous frequency of the deterministic part of the signal
%	    blurried in a white gaussian noise.
%	    L must be an integer       	（default : 1）.
%	TRACE : if nonzero the progression of the algorithm is shown
%	                                （default : 0）.
%	FNORMHAT : Output （normalized） instantaneous frequency.
%	T : Time instants.
%
%	Examples : 
%	 x=fmsin（700.050.3525）; [instft]=instfreq（x）; plot（tinstf）
%	 N=64; SNR=10.0; L=4; t=L+1:N-L; x=fmsin（N0.050.3540）;
%	 sig=sigmerge（xhilbert（randn（N1））SNR）;
%	 plotifl（t[instfreq（sigtL）instfreq（xt）]）; grid;
%	 title （‘theoretical and estimated instantaneous frequencies‘）;
%
%	See also  KAYTTH SGRPDLAY.

%	F. Auger March 1994 July 1995.
%	Copyright （c） 1996 by CNRS （France）.
%
%	------------------- CONFIDENTIAL PROGRAM -------------------- 
%	This program can not be used without the authorization