电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

Hilbert边际谱

Hilbert边际谱_第1页
1/11
Hilbert边际谱_第2页
2/11
Hilbert边际谱_第3页
3/11
1、Hilbert边际谱 我觉得既然已经做出 EMD 了,也就是得到了 IMF。这个时候就是做 hilbert 幅值谱,然后对它积分就可以了。程序不是很难搞到吧! 我是用 hspec 画谱图的,自己又在后面添加了求边际谱的代码 for k=1:size(E) bjp(k)=sum(E(k,:))*1/fs; %fs 为采样频率; end figure plot(bjp); xlabel('频率 / Hz'); ylabel('幅值'); 比如我用两个正弦信号作仿真 fs=1000; t=1/fs:1/fs:1; y1=2*sin(40*pi*t); y2=5*sin(80*pi*t); y=[y1,y2]; % 信号 画出来的图很粗糙,更不用说对实际信号分析了,所以大家看看如何来修正?黄文章中边际谱对实际信号分析是很好的一条曲线 我用 hhspectrum 算了一下谱图,同时求了一下边际谱,边际谱程序基本想法同 form。 结果也不太好,20HZ 处还行,40HZ 就有些问题了,见附图 你自己再用这个试试 我没有用 rilling 的 hhspectru m nspab: function h1= nspab(data,nyy,minw,maxw,dt) % The function NSPAB generates a smoothed HHT spectrum of data(n,k) % in time-frequency space, where % n specifies the length of time series, and % k is the number of IMF components. % The frequency-axis range is prefixed. % Negative frequency sign is reversed. % % MATLAB Library function HILBERT is used to calculate the Hilbert transform. % % Example, [h,xs,w] = nspab(lod78_p',200,0,0.12,1,3224). % % Functions CONTOUR or IMG can be used to view the spectrum, % for example contour(xs,w,h) or img(xs,w,h). % % Calling sequence- % [h,xs,w] = nspab(data,nyy,minw,maxw,t0,t1) % % Input- % data - 2-D matrix data(n,k) of IMF components % nyy - the frequency resolution % minw - the minimum frequency % maxw - the maximum frequency % t0 - the start time % t1 - the end time % Output- % h - 2-D matrix of the HHT spectrum, where % the 1st dimension specifies the number of frequencies, % the 2nd dimension specifies the number of time values % xs - vector that specifies the time-axis values % w - vector that spec...

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

Hilbert边际谱

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部