2
ゼロパディングの奇数長FFTの場合の実数値のリンギング
だから私は、信号と逆変換の周波数応答をゼロで埋める周波数領域補間器を書き込もうとしています。対処しなければならない2つのケースがあります。 偶数長の応答ビンを分割する必要があるのは、あいまいだからです。そこで、スペクトルの負の部分をコピーし、その間にゼロを追加します。Fs/2Fs/2F_s/2n*(interp-1)-1 奇数長の応答ビンがないため、正/負の周波数を分割し、それらの間にゼロを挿入するだけです。Fs/2Fs/2F_s/2n*(interp-1) ゼロパディングを行うコードはこちらにあります // Copy negative frequency components to end of buffer and zero out middle // inp - input buffer of complex floats // n - transform size // interp - interpolation amount void zero_pad_freq(cfloat_t *inp, size_t n, size_t interp) { if ((n % 2) == 0) { …
13
fft
interpolation
c