关于MATLAB的问题

首先要定义函数，然后用函数名调用。这个说起来有点复杂，建议你百度matlab函数的使用或者看书上的函数定义，你就会明白的。

网页链接 网页链接这两个经验挺好的

2014版的自带函数名称换了。由defaultStream换成getGlobalStream。亲测好用。

你打开eemd函数，对一下部分，进行如下修改，即可。

function savedState = set_seed(seedNo)

defaultStream = RandStream.getGlobalStream %￥￥￥￥￥原程序是，RandStream.defaultStream。新版matlab中没有这个函数了，换成了getGlobalStream。因此，我也给换成了getGlobalStream￥￥￥￥￥￥%%%%%

savedState = defaultStream.State% 自

rand('seed',seedNo) % 己

randn('seed',seedNo)% 注

% 释

end % 处

% 已经保存运行了

function return_seed(savedState)

RandStream.getGlobalStream.State = savedState%￥￥￥￥￥原程序是，RandStream.defaultStream。新版matlab中没有这个函数了，换成了getGlobalStream。因此，我也给换成了getGlobalStream￥￥￥￥￥￥5%%%%%

end

这个分解是基于希尔伯特-黄变换和希尔伯特变换而来，通过黄变换滤除信号里局域的直流成分，短时内是纯交流成分。而之所以不直接使用希尔伯特变换后利用复信号的d(phi)/dt=w (phi是角度）的方式来求信号的瞬时频率，是因为信号可能存在非频率波动的成分，也就是说可能信号的波动是因幅值引起，这样求得的频率可能是负值，而希尔伯特变换的信号频率不存在负的，所以才需要以上的黄变换这个分解imf的过程。但这个分解精度纯粹因信号本身而异，存在一定的风险。不过据我认为，信号在离散的情况下，黄变换并不是必须的，只需要在希尔伯特变换之前加一些简单处理即可，这种方式得到的瞬时频率与用了imf后得到的再matlab下得到的图像基本一样。楼主可自己编写程序，这个程序包不是必须

---