高速フーリエ変換

MATLABによる震度算出方法。算出方法そのものは国土交通省 気象庁で公開されているものとなる。FFT(高速フーリエ変換)で周波数領域に変換してバンドパスフィルタと通し、IFFT(逆高速フーリエ変換)で再度、時間領域に戻す。計算式があっていても、演算順序やFPUの違い、FPUモードの違いで演算誤差がでる。