問題タブ [fftpack]
For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.
fortran - fftpack5.1を使用したフーリエ変換の不具合
離散フーリエ変換を計算するサブルーチンを含む Fortran 90 で FFTPACK5.1 を使用すると問題が発生します。私はなんとかそれをインストールしてルーチンを使用しましたが、周波数Aの単純な正弦波ですべてが問題ないかどうかを確認すると、 A (周波数空間、スペクトル) ではなく2Aでゼロ以外の係数が得られます。スペクトルに変化があり、その理由がわかりません。周波数軸のステップを正しく計算することはほぼ確実です（ただし、疑問があります）。
N は元の正弦波のポイント数、Fech はサンプル周波数で、周波数軸のステップを df(i)=Fech(i-1)/N として計算します。
私はrfft1fルーチンを使用しているので、経験があり、私の問題を知っている人がいれば、ここで何が問題なのかを理解するのは本当に素晴らしいことです。
これが私のコードです：
単純な正弦波では、周波数が 0.5 (cf コード) のディラックを持つことを期待していますが、代わりに、周波数領域で 1. のディラックが得られます。その上、スペクトルは奇妙に見えます...これが私が得たものです:
python - Python で同等の MATLAB fftfilt
MATLAB で作成した次の関数を Python に変換しようとしています。
そこで、MATLAB で定義されたこの関数を Python に実装しようとしています。
関数fftfiltは MATLAB の fftfilt のクローンであると考えていましたが、実行すると次のエラーが発生しました
だから、私の質問は: Python で MATLAB fftfilt に相当するものはありますか? 私の関数の目的はoutput_phase、位相信号の速い変動を修正し、n*90 度のシフトを修正することです。
python - このスクリプトを実行するとコンピュータがフリーズするのはなぜですか?
SciPy を使用して Python でスクリプトを作成し、信号に対して短時間のフーリエ変換を実行しました。1000 のタイムポイントを持つシグナルで実行したところ、問題なく実行されました。100 万のタイムポイントを持つ信号で実行すると、コンピューターがフリーズしました (コンピューターが応答せず、オーディオが再生されている場合、コンピューターはスキップとループのバズを出力します)。これは、私が試みた3回すべてで一貫して発生しました。何時間もかかるスクリプトを書いたことがありますが、実際にコンピューターがフリーズするスクリプトに遭遇したことはありません。理由はありますか？スクリプトは以下に掲載されています。