我有两个信号a和b作为两个1d numpy数组。我正在尝试将两个信号的光谱偏差计算为两个单独的值，即平均光谱偏差和最大光谱偏差。

到目前为止，我已经尝试过：

取各个功率谱的对数以计算最大和平均频谱偏差

   a = np.asarray(signal_a)
   b = np.asarray(signal_b)

   power_spectrum_a = np.abs(np.fft.fft(a))**2 
   power_spectrum_b - np.abs(np.fft.fft(b))**2

   spectral_dev = np.abs(20*np.log(power_spectrum_a/power_spectrum_b)))

   max_spectral_dev = np.max(spectral_dev_in_db)
   avg_spectral_dev = np.mean(spectral_dev_in_db)

我做对了吗？