問題タブ [dft]

Java で小さな離散フーリエ変換を記述して、明確な 400 Hz 正弦波信号 (pcm 符号付きショートとして 1 秒) で振幅スペクトルを見つけようとしています。

最初に、複素数の DFT を計算します。

次に、マグニチュード スペクトルを計算します。

しかし、何らかの理由で、16000 サンプルすべてをチェックすることを選択した場合、「マーカー」で 400 Hz の結果しか得られません。しかし、800 では 800/2 = 400 Hz が最大周波数として表示されるため、800 サンプルのみを選択した場合でも 400 Hz でピークが表示されるべきではありませんか?

800 サンプルを選択すると 20 Hz になり、1600 サンプルを選択すると 40 Hz になり、これは常に 1/40 * サンプル レートになります。

私は一体何を見逃しているのですか、それとも間違っていたのですか? 結果は奇妙です..

複素数値で逆 DFT を実行すると、オーディオ信号を再構築できることに注意してください。

これは正しいです？p1はソース画像と宛先画像であり、グレースケール画像で作業していて「ith」チャンネルp2にアクセスしているため、3 つのチャンネルはすべて同じ値を持っています。pp1[i]

FFTW を使用するコードを KissFFT を使用するように移植しようとしています。
コードはfftwf_plan_r2r_2d()with をFFTW_REDFT01使用します。

KissFFT での同等の呼び出しは何ですか?

この呼び出し ( with ) が DCT と同等である場合、代わりにOpenCVFFTW_REDFT01などの直接 DCT 変換を使用できますか? 反射や対称化など、必要な入力データの変更はありますか?cv::dct

私は自分の2D DFTを書く必要があり、現在これを使用しています

小さい画像では問題なく動作しますが、ピクセルの寸法または配列が大きくなると、非常に遅くなります。誰でも助けることができますか？

これim3は、ピクセル値が格納されている配列でありim、アルゴリズムの実行後に値を格納する配列です

学位プロジェクトとして音符検出器を開発したいのですが、ゼロからやりたいと思っています。そのオーディオ音楽ファイルからすべての情報を抽出し、結果として振幅を与える「.wav」ファイルのコードを書きました。

次に、DFT のコードを作成しました。軸 (実数/虚数) の 1 つが振幅/大きさで、もう 1 つが位相である複素数として出力されます。

ここで、周波数 (ベクトルではなくヘルツ) で答えが欲しいので、DFT が適切な出力を提供するかどうかを確認できます。DFT 出力を周波数に変換するにはどうすればよいですか?

これを C 言語でコーディングする必要があり、組み込みライブラリを使用したくない

上記で利用できる適切な実装があるかどうかを知る必要があります。ライブラリ/リソースはほとんど見つかりませんでしたが、それらのほとんどは異なるcppファイル内に分散されているため、複雑すぎて使用できませんでした。

のように、この投稿では、いくつかのリソースについて静かに言及していますが、どれも私には当てはまりません。

わずか 1 cpp 以内のコンパクトな実装が必要です。

多項式の次数は最大 10^6 の範囲です。、係数はすべて非常に小さく、最大 10^1 です。