3Dパワースペクトル、つまり周波数シェルの平均パワーを計算しようとしています。私はシェル密度の計算を正しく行っていると思いますが、私は各シェルの頻度をどのように決定するか分かりません。fftwで3Dパワースペクトルを計算する
サンプリングレートFs
が各次元で同じであり、各次元の元のサンプルの長さも同じ値N
であるとします。シェルの「インデックス」はidx = sqrt(i*i + j*j + k*k)
です。ここで、i、j、およびkは各方向の範囲です。このシェルの頻度をどのように計算するのですか?
「周波数シェル」とは何ですか?これはプログラミングの質問よりも物理学の質問のほうに聞こえます。一般に、実数値の空間領域ボリュームデータセットの離散フーリエ変換は、複素数値の周波数領域ボリュームデータセットを生成する。あなたは何か違うものがほしいと思うようです。 – MooseBoys
あなたが何を求めようとしているかを完全には確認していません。 FFTW計画の設定、計画の実行、または結果の解釈方法を尋ねようとしていますか? – RyanP
@MooseBoysつまり、各周波数の密度値(f x f *)を平均したいと思っています。 3Dのためシェルです。2Dのサークルになります。 – cgreen