我正试图获取信号的包络。在一个旧主题中,我发现了一些旧代码abs(hilbert(signal))
。
使用1/3倍频程频带对时域中的加速度信号进行滤波后,使用上面的代码,我的信号如下所示:
红色 - >希尔伯特蓝 - >信号
result http://1.ii.gl/gaNtcE9To.png
我期待更像这样的事情:
example http://www.emt.tugraz.at/publications/diplomarbeiten/da_hoebenreich/img56.gif
好像我的信封不太合适。
这是我的1/3倍频程代码:
Fs = 5000;
N= 4;
Fc1= 1425;
Fc2= 1795;
h= fdesign.bandpass('N,F3dB1,F3dB2', N, Fc1, Fc2, Fs);
Hd= design(h, 'butter');
我做错了什么?