马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我利用HHT处理一个分段正弦信号,信号为y1 = sin(2*pi*10*t1),y2 = sin(2*pi*30*t2);采样频率为512,共1024个点。下面是用HHT做的分析。(当然这个信号本身就是一个IMF不需要在EMD分解),我做的三维谱图和时频图都对,但是那个边际谱图幅值误差很大,绘制边际谱的程序是利用论坛中提到的那个。有没有高手帮我找找毛病,改进一下。谢谢了!程序如下:
function plot_bjp(imf,fs,N)
if size(imf,1) > 1
[A,fa,tt] = hhspectrum(imf(1:end-1, :));
else
[A,fa,tt] = hhspectrum(imf);
end
[E, tt1] = toimage(A,fa,tt,length(tt));
for k = 1:size(E,1)
bjp(k) = sum(E(k,:))/size(E,2);
end
f = (0:N-3)/N*(fs/2);
figure
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
title('边际谱','color','r');
set(gca,'Xlim',[0,fs/2]);
分析结果如下:
|