窗口做短时傅里叶变换

可以使用Python中的scipy库来进行短时傅里叶变换。具体步骤如下：

1. 读取音频文件，将其转换为一维数组。
2. 定义窗口大小和步长，将音频信号分成多个窗口。
3. 对每个窗口进行加窗操作，常用的窗口函数有汉宁窗、汉明窗等。
4. 对每个窗口进行傅里叶变换，得到频谱。
5. 将所有窗口的频谱拼接起来，得到整个音频信号的短时傅里叶变换结果。

下面是一个简单的示例代码：

import numpy as np
from scipy import signal
from scipy.fft import fft, fftshift
import matplotlib.pyplot as plt
from scipy.io import wavfile

# 读取音频文件
fs, data = wavfile.read('audio.wav')

# 定义窗口大小和步长
win_size = 1024
hop_size = 512

# 分帧
frames = signal.frame(data, win_size, hop_size)

# 加窗
window = signal.hann(win_size)
frames *= window

# 傅里叶变换
spectrums = np.abs(fft(frames, axis=1))

# 拼接频谱
stft = np.concatenate(spectrums, axis=0)

# 绘制频谱图
plt.imshow(stft.T, aspect='auto', origin='lower')
plt.xlabel('Time')
plt.ylabel('Frequency')
plt.show()