numpy.fft.ifft()の戻り値の配列に意図せずにnanが含まれてしまい、困ってます。

Question

前提・実現したいこと

numpy.fft.ifft()を使いwavファイルからケプストラムを得るために、フーリエ逆変換をしようとしているのですが、numpy.nan,numpy.inf,numpy.NINFを含まない配列を渡したつもりだったのにnanを含む配列が返ってきてしまいます。

該当のソースコード

import numpy as np

#wavはあるwavファイルから0.04秒分を切り出してint16の一次元配列に変換したものです。各要素の範囲を-1から1に
#正規化してあります。

hanningWindow = np.hanning(len(wav))
wav = wav * hanningWindow[:,None]

n=2048

Plog=10 * np.log10(np.abs(np.fft.fft(wav,n))**2) #対数スペクトル化


nan1=np.where(Plog==np.inf,(sys.float_info.max),Plog) #正の無限大を有限に

nan2=np.where(nan1==np.NINF,-(sys.float_info.max),nan1) #負の無限大を有限に

print("isfinite ",np.all(np.isfinite(nan2))) #nan2の要素全てがnan、inf、NINFでない、有限の数のみなのかを
                                             #チェックしているつもりです

nan3=np.fft.ifft(nan2) #フーリエ逆変換

print("isnan ",np.any(np.isnan(nan3))) #nan3の要素のうち、一つにでもnanが含まれているかを
                                       #チェックしているつもりです

発生している問題(出力結果)

C:...\teratailyou.py:32: RuntimeWarning: divide by zero encountered in log10
  Plog=10 * np.log10(np.abs(np.fft.fft(wav,n))**2)
isfinite  True
isnan  True

試したこと

• numpy.fft.ifft()にnan,inf,NINFを含む配列を与えるとnanを含む配列が返ってくることを知り、

それらを有限の数に置き換えてみました。その結果前述の結果となってしまい、以降何も思いつけませんでした。
どうかご助力ください。

Comments

[t_obara]

置き換える前に、エラーになっているのでは？

[redflute_gumi]

log10の中での0除算エラーはlog10(0)の場合に発生し、これはNINFを返すため問題ないと考えていました。この考え方が誤っているということでしょうか。
　それとも、非有限数の置き換え時や非有限数の判定がおかしい、ということでしょうか。
　私は未熟でいただいたアドバイスの意味がよくわからず…。もう少し詳しくご教示くださいませんか？

[t_obara]

10 * np.log10(np.abs(np.fft.fft(wav,n))**2)
この計算を分解して、一つずつ確認してみては？

Answer 1

無限大を入れてエラーになるなら、float_info.maxすなわち無限大よりわずかに小さい(※)数を入れてもFFT内部の計算で無限大が生じて同様にNaNになるでしょう。
n=2048だとfloat_info.max/2048なら大丈夫かな? ちょっと自信がありませんが、まあとにかく、そのような無限大よりだいぶ小さい数を入れればNaNは出なくなると思います。
答えが正しいかどうかは分かりませんが。

※正確には、「浮動小数点数で無限大と扱われる最小の数」よりわずかに小さい数