RIAA補正カーブを精密に計算する | 創造の館音楽苦楽部

レコードの特性評価に必要なRIAAカーブの計算方法がネットにあまりない。あっても間違っていて結果がおかしい。そこで、エクセルで一発で出せる計算式を作った。

RIAA規格の時定数から、s領域の伝達関数は次のようである。

ここから振幅特性を計算できる。次は、dBの振幅特性を直接得るエクセルの計算式。

=20*LOG10(SQRT((1+(0.000318*2*PI()*ｆ)^2)/((1-0.000075*0.00318*(2*PI()*ｆ)^2)^2+((0.000075+0.00318)*2*PI()*ｆ)^2)))+19.911

ｆは周波数(Hz)、最後の19.911は１kHzでゼロにするためのオフセット補正。

pythonで離散領域のフィルター係数を求める場合は、次のようにする。

from scipy import signal
 b,a=signal.bilinear([0,3.18e-4,1],[2.385e-7,0.003255,1],samprate)

samprateはサンプリング周波数(Hz)。