0

信号のエンベロープを取得しようとしています。古いトピックで、古いコードを見つけましabs(hilbert(signal))た。

1/3 オクターブ バンドを使用して時間領域で加速度信号をフィルター処理し、上記のコードを使用すると、信号は次のようになります。

赤 --> ヒルベルト ブルー --> シグナル

結果
(出典: 1.ii.gl )

私はこのようなものをもっと期待していました:

例
(ソース: tugraz.at )

封筒がうまく入らないようです。

1/3 オクターブバンドのコードは次のとおりです。

Fs = 5000;
N= 4;     
Fc1= 1425;  
Fc2= 1795;  
h= fdesign.bandpass('N,F3dB1,F3dB2', N, Fc1, Fc2, Fs);
Hd= design(h, 'butter');

私は何を間違っていますか?

4

0 に答える 0