小波与hilbert变换结合进行信号包络谱分析
这是我自己搞的一个仿真信号用小波与hilbert变换进行信号包络谱分析的程序大家看一下,给点意见啊。我是结合书上的编的,有几个地方一直不懂,
1,许多论文上说hilbert变换是这样进行的:原始信号X(t)带通滤波,进行hilbert变换得到x^(t),作为虚部,
然后用w(t)= x(t)+jx^(t)作为解析信号。包络谱的求法是:a(t)=X(t)和x^(t)的平方和再开放.但是这个程序里面怎么没有这个功能??
2,plot((0:nfft/2-1)/nfft*fs,p(1:nfft/2));这句程序有点难懂,麻烦哪位编程高手解释一下??
程序如下:
t=0:0.005:1*pi;
fs=10000
s=4*sin(2*200*pi*t).*(sin(2*4500*pi*t))+25*(sin(2*4500*pi*t));
subplot(411);plot(t,s)
=wavedec(s,1,'db10');
d1=wrcoef('d',c,l,'db10');
a1=0;
subplot(412);plot(d1);title('重构高频信号');
y=hilbert(d1);
y1=abs(y);
ydata=y-mean(y);
nfft=1024;
p=abs(fft(y1,nfft));
figure(2);
plot((0:nfft/2-1)/nfft*fs,p(1:nfft/2));
xlabel('频率');
ylabel('功率谱'); 1 叙述的没错,我机子上没装matlab不好查,再仔细看看hilbert.m,估计是问题不大
2 这句是画图,肯定都明白;我想不明白的是为什么/2,这是因为频谱是对称的,我们只要画正半轴,所以/2 ydata=y-mean(y);
这句是做什么用的?? plot((0:nfft/2-1)/nfft*fs,p(1:nfft/2));这句话好像是把你的采样信号转换为频率,并且保证横纵坐标的点数是一样的 不错,你这个程序根本就不能实现包络分析的功能,求包络可以用三次样条插值来拟合,
至于包络点你已经知道怎么求解了
那就很容易得到它的包络图了 要正半轴呀 “ 这句是画图,肯定都明白;我想不明白的是为什么/2,这是因为频谱是对称的,我们只要画正半轴,所以/2”
因为数据的采集都是一个转换过程,采样都需要遵守一个采样定理,因为有这个采样定理,所以你最后要是不将频谱/2,你得到的频谱图就是一个对称的。
页:
[1]