signal.freqz
したがって，filtfiltのgain はこの複素振幅の４乗になる．2乗が振幅をgainにするためで，もう一回の２乗がfiltfiltでは２回フィルターをかけるためである．
周波数はradians/sampleの単位で0～piで表現されるので，通常のfrequency=1/period に変換するには，1/np.pi0*.5を乗じる．
デジタルフィルターの係数を与えてその応答を，複素振幅で求める．
b,a=signal.butter(4,frq_pb)
w,h=signal.freqz(b,a)
frq=w/np.pi*0.5
scipy.signal.freqz obtains filter response in complex amplitude.
#Example
b,a=signal.butter(4,frq_pb)
w,h=signal.freqz(b,a)
frq=w/np.pi*0.5
Thus, the gain of filtfilt is expressed by quadruple of this amplitude.
Frequency, w, is given by radians/sample from 0-pi, and can be converted by multiplying 1/np.pi0*.5.