samllpig1984 发表于 2009-11-3 16:47

关于hhspectrum的结果成图。

我用小波分析信号间断点的相关信息，让后想用hhspectrun分析间断点的时频信息；
新图像：（有四张图的）表示的小波分析的第1，2 层近似系数和1，2 层的细节系数，看出间断点的位置在》50左右的位置。
两张图的是对第2层细节系数（四张图的最后一个）的振幅和频率的分析。我想得到的图像是横坐标是频率，纵坐标是振幅，可图像却像第三张所示，不知道是什么东西；
s=leleccum(1575:1720);
=wavedec(s,5,'db3');
a1=wrcoef('a',c,l,'db3',1);
a2=wrcoef('a',c,l,'db3',2);
d1=wrcoef('d',c,l,'db3',1);
d2=wrcoef('d',c,l,'db3',2);
subplot(411);plot(a1);
subplot(412);plot(a2);
subplot(413);plot(d1);
subplot(414);plot(d2);a1,a2,d1,d2为第1，2 层近似系数和1，2层   =hhspectrum(d2);细节系数
然后就是画出A，f的图像
plot(f,A);（最后一张的命令）；为什么不正确，请大侠致电一下

真人 发表于 2009-11-4 14:42

plot(f,abs(A));试一下

samllpig1984 发表于 2009-11-5 15:18

看了hhspectrum的源代码，A本来就是abs函数求的，所以plot(f,A)和plot(f,abs(A))是一样的
页: [1]
查看完整版本: 关于hhspectrum的结果成图。