用Python求取语音的基音频率

要求取语音的基音频率，可以通过自相关函数实现。下面是一个用Python计算语音基音频率的例子：

import numpy as np
import scipy.signal as signal

# 读取语音文件
sample_rate, signal_data = wavfile.read('speech.wav')

# 帧长和帧移
frame_length = int(sample_rate * 0.02)
frame_step = int(sample_rate * 0.01)

# 分帧
frames = np.array([signal_data[i:i+frame_length] for i in range(0, len(signal_data)-frame_length, frame_step)])

# 求取自相关函数
corr = np.array([np.correlate(frame, frame, mode='full') for frame in frames])
corr = corr[:, len(corr[0])//2:]

# 找到每一帧的基音周期
periods = []
for i in range(len(corr)):
    peak = signal.find_peaks(corr[i], height=0)[0]
    if len(peak) > 1:
        period = (peak[1] - peak[0]) / sample_rate
        periods.append(period)
    else:
        periods.append(0)

# 计算基音频率
fundamental_freq = np.array(periods)
fundamental_freq[fundamental_freq > 0] = 1 / fundamental_freq[fundamental_freq > 0]
fundamental_freq = fundamental_freq * sample_rate

# 输出基音频率
print(fundamental_freq)

需要注意的是，这个方法只适用于单音高语音的基音频率提取。对于多音高复杂语音，基音频率的提取需要使用更加复杂的算法。