回复 #121 mofei 的帖子
已经对了啊，还有什么问题吗？ 出现的频率不对啊，应该是在60Hz和90Hz处出现，但我那副图中两个频率都小于50Hz啊。回复 #123 mofei 的帖子
y1=2*sin(60*pi*t);y2=5*sin(90*pi*t);
sin(2*pi*f*t)中的f才是频率，这样你就知道上面两个信号的频率到底是多少了。 不好意思，是我看错了。谢谢您。 原帖由 破凰 于 2007-4-27 21:51 发表 http://www.chinavib.com/forum/images/common/back.gif
form求边际谱时所用程序是没有问题的，用的是矩形积分公式。
他所得结果不正确的原因是：输入的应是调用了toimage后的结果，而不是调用了hhspectrum后的结果。
下面给一段程序，大家可以去试下。边际谱的分析结果是 ...
下面是我仿真的图
第一张是用你的程序作出来的
第二张是用nspab函数和mspc函数作出来的
好像有区别呀
为什么呢？
Hilbert边际谱和FT变换后的幅频谱
这得出的Hilbert边际谱和FT变换后的幅频谱为什么会有这么大的区别呢，到底哪个幅值才是真正的实际幅值呢？有参考价值吗？ 程序如下：
load shuju
fs=5120;
N=4096;
a1=a(1:N,1);
a2=abs(fft(a1))*2/N;
f=fs*(0:N/2-1)/N;
n=length(f);
subplot(211)
plot(f,a2(1:n))
xlabel('频率 / Hz');
ylabel('幅值');
title('FT后的幅频图')
imf=emd(a1);
=hhspectrum(imf);
=toimage(A,fa,tt,length(tt));
for k=1:size(E,1)
bjp(k)=sum(E(k,:))*1/fs;
end
f=(0:N-3)/N*(fs/2);
subplot(212)
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
title('Hilbert边际谱') 没人回答吗 EMD分解的IMF能量和原信号其能量是不相等的。所以边际谱能量不能和FFT的能量相比。边际谱能量只能说明某个信号存在，能量相对于其他的大小。
求助：EMD分解后，单个IMF的能量谱怎么做？
EMD分解后，单个IMF的能量谱怎么做？回复 113楼 chiu 的帖子
你的那个为什么乘以fs，有没有找到答案呢，为什么啊，我很不明白，为什么是对toimage的每一行进行求和？麻烦你一定要回一下，谢谢 这得出的Hilbert边际谱和FT变换后的幅频谱为什么会有这么大的区别呢，到底哪个幅值才是真正的实际幅值呢？有参考价值吗？ 程序如下：
load shuju
fs=5120;
N=4096;
a1=a(1:N,1);
a2=abs(fft(a1))*2/N;
f=fs*(0:N/2-1)/N;
n=length(f);
subplot(211)
plot(f,a2(1:n))
xlabel('频率 / Hz');
ylabel('幅值');
title('FT后的幅频图')
imf=emd(a1);
=hhspectrum(imf);
=toimage(A,fa,tt,length(tt));
for k=1:size(E,1)
bjp(k)=sum(E(k,:))*1/fs;
end
f=(0:N-3)/N*(fs/2);
subplot(212)
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
title('Hilbert边际谱')
附件
shuju.mat (14.75 KB)
2008-12-25 13:18, 下载次数: 46
所用数据
0.JPG (32.74 KB)
2008-12-25 13:18
图形 原因是什么啊？？？？怎末没人回答？？