Scilabで地震の震度を計算する

　地震の震度は単にその加速度の大きさで決めるのかと思っていたら、決められたフィルターを通した上に計算式を使って出す必要がある。会社で作ってと言われたので、Scilabで作ってみた。
　決められた周波数特性になるように窓を作成し、地震波形をフーリエ変換して周波数軸にしたうえで作った窓をかけて、逆フーリエをする。単にFFTフィルターだ。FFTフィルターは個人的にはOrijinというソフトで使っているのしか知らない。位相がずれないので何か使う時もあるかもと以前作ったのが役に立った。