关于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);(最后一张的命令);为什么不正确,请大侠致电一下 plot(f,abs(A));试一下 看了hhspectrum的源代码,A本来就是abs函数求的,所以plot(f,A)和plot(f,abs(A))是一样的
页:
[1]