未知の信号のフーリエ変換からの奇妙な結果、何が悪いのですか?


7

工学を勉強して久しぶりですが、間違いは許してください。

デジタル加速度計からデータをサンプリングしました。20 Hzでサンプリングされ、1秒の期間に縮小されます。次に、データをフィルタリングしてノイズを除去します。システムは安静時に0を出力し、アクティビティ中はさらに多く出力します。

1秒に1回サンプリングされた1024サンプルの信号のフーリエ変換を調べています。私はこれを見ています:ここに画像の説明を入力してください

これは、フーリエ変換の実数と虚数の結果の実数と虚数のプロットです。 ここに画像の説明を入力してください

信号は周波数領域でらせんのように見えます。システムは周波数領域内の制限サイクルの周りで振動しているようです。

これは可能ですか?フーリエ変換はデジタル信号に対して完全に間違っていますか?この種の分析にはDFTのみを使用する必要がありますか?私のあるウィンドウサイズ間違っては? ウィンドウサイズが正しくないためにエイリアシングが発生していますか?

前回のエンジニアリングのクラスから数年が経ちましたが、この件に関する助力をいただければ幸いです。工学、信号解析、変換に関する本を3冊もらったが、この種の知識が戻ってくるまでには時間がかかる。


私はこれを注意深く見ていませんでしたが、プロットを見たときの私のひざの反応は、エイリアシングが起こっているということです。

いいえ、フーリエ変換はデジタル信号に問題はありません。はい、DFTは離散信号でのみ使用できます。
内部石

時間の経過に伴う加速度のグラフは、期待どおりのものに見えますか?
endolith

加速度計のグラフは良さそうです。私は再確認しましたが、プロットしているのは、特定の移動平均を超える加速度の変化率です。期待どおりに動作します。センサーに触れると、加速度が一瞬現れ、その後消えます。
Alex Stone、

回答:


7

私の推測では、ここには何も表示されません。一緒に移動。

「信号は[周波数領域でらせん状に見える]とのことですが、このらせんが本当に何であるかは、複雑な指数関数であると私は断言します。そして、信号に時間遅延を与えると複雑な指数が得られるので、これには何の問題もありません。f(xa)f^(ξ)e2πiaξ。らせんパターンは、取得している信号のパワーのほとんどが、T=0。(とにかくそこにあるべきではありません!)

目を見張ると、周波数プロットの振動の周期は〜28サンプルのようです。20hzのサンプリングレートで、私が言っていることが真実であれば、時間領域信号のほとんどのエネルギーは、おおよそt2820=1.4秒。近かったですか?

ここでの本当の問題はあなたの陰謀だろう。生のリアル/ imag DFT出力をプロットしているようです。しないでください。代わりに、実際の値またはimagの値からマグニチュード/位相を計算し、それをプロットします。


1
ありがとうございました !らせんをクリーンアップしましたが、毎秒25回振動しています。25 Hzで信号のサンプリングを開始した場合、信号の大部分は1秒のエポック内に集中しますか?
Alex Stone

プロットが奇妙に見えると聞いたことがあります。正弦波または正弦波の合計を使用して単純なFTを実行すると、予想どおり基本周波数にピークが発生します。取得したデータで同じことをしようとすると、このような奇妙な結果が得られます。ビンの大きさを詳しく見ていきます。
Alex Stone、

1
そのFFTプロットはかなり正常です。ある瞬間に1つの大きなスパイクがあり、他の場所ではゼロに近い時間領域データがある場合、そのデータに対してFFTを実行すると、周波数プロットでそのらせんパターン(複素指数)が得られるはずです。
David Cary

3

考えられるスターター:

エイリアシングを回避するために、ADCへの入力がサンプリングレートの半分よりもかなり低くなるようにローパスフィルター処理されていることを確認してください。20 Hzのサンプルでは、​​10 Hzの「納屋のドア」の無限カットオフローパスフィルターか、やや低い位置にある実物が必要です。シート21にはコンポーネントのエイリアシングの感触がありますが、そうではないかもしれません。

カットオフ周波数が10 Hz未満の1秒のサンプルでは、​​数サンプルしか得られません。私はあなたが本当にしていることを完全に逃しているかもしれません。

サンプリングされたウィンドウに含まれる完全な波形サイクル以外を処理するには、ウィンドウ処理が必要になる場合があります。FFTパスに含まれるサイクルの任意の部分を持ついくつかのサンプルと周波数成分について、存在しない強力な成分を生成できます。

FFTは帯域内ノイズを処理する必要があります。「ノイズフィルタリング」は、データフィルタリングの場合もあります。上記のようにナイキストレートフィルターが必要ですが、残りの通過帯域のすべてが正当な信号である可能性があります。


2
彼は、「デジタル加速度計」と言った。これは、ADCを内蔵し、セルフクロックサンプリング(たとえば、ADXL345)を備えた加速度計を意味する。このような加速度計には、おそらくアンチエイリアシングフィルターが組み込まれています。
11

はい、できれば加速度計にアンチエイリアシングフィルターが付いているといいのですが、OPが「1秒の期間に折りたたまれた」前に、OPが再びフィルターをかけ、意味のある信号を取得する必要があります。
David Cary

2

表示されているのは、1秒の「エポック」ウィンドウがデータと同期していないという事実です。FFT結果の位相はウィンドウのエッジを基準にしているため、ウィンドウのエッジが信号とのさまざまな位相関係に移動すると回転します。

実際に位相または相対位相を気にする場合は、ウィンドウオフセットを信号の周期にロックします。フェーズを気にしない場合は、複雑な結果の大きさを計算し、それを使用してください。


意味あり。現在のウィンドウがデータとともにスライドします。ウィンドウを増やすか減らすかについて何か提案はありますか?
Alex Stone
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.