qq_33670784 2020-04-17 11:17 采纳率: 50%
浏览 1328
已结题

语音处理快速傅里叶变换xf = np.fft.rfft(xs,fft_size)/fft_size之后xf长度不变

语音处理快速傅里叶变换xf = np.fft.rfft(xs,fft_size)/fft_size之后信号长度没有变为fft_size/2 +1 长度仍然是xs的长度,按理来说频率个数应为:np.linspace(0, sample_rate/2, fft_size/2+1) 但是这样就会报错:
ValueError: x and y must have same first dimension, but have shapes (257,) and (512, 257)

如果用人工生成的周期信号这样处理不会有上述问题,但是使用wav的语音信号输入,就会出来经过rfft函数信号长度不变的情况,有没有哪位大佬知道这是为啥吗?

def plot_freq(signal, sample_rate, fft_size=512):
xs= signal[:fft_size]
xf = np.fft.rfft(xs,fft_size)/fft_size
freqs = np.linspace(0, sample_rate//2, fft_size/2+1)
xfp = 20 * np.log10(np.clip(np.abs(xf), 1e-20, 1e100))
print(len(xf))

print(len(freqs))
plt.figure(figsize=(20, 5))
plt.plot(freqs, xfp)
plt.xlabel('Freq(hz)')
plt.ylabel('dB')
plt.grid()
plt.show()


  • 写回答

1条回答 默认 最新

  • dabocaiqq 2020-04-17 15:56
    关注
    评论

报告相同问题?

悬赏问题

  • ¥50 导入文件到网吧的电脑并且在重启之后不会被恢复
  • ¥15 (希望可以解决问题)ma和mb文件无法正常打开,打开后是空白,但是有正常内存占用,但可以在打开Maya应用程序后打开场景ma和mb格式。
  • ¥15 绘制多分类任务的roc曲线时只画出了一类的roc,其它的auc显示为nan
  • ¥20 ML307A在使用AT命令连接EMQX平台的MQTT时被拒绝
  • ¥20 腾讯企业邮箱邮件可以恢复么
  • ¥15 有人知道怎么将自己的迁移策略布到edgecloudsim上使用吗?
  • ¥15 错误 LNK2001 无法解析的外部符号
  • ¥50 安装pyaudiokits失败
  • ¥15 计组这些题应该咋做呀
  • ¥60 更换迈创SOL6M4AE卡的时候,驱动要重新装才能使用,怎么解决?