matlab对.wav文件进行短时傅里叶变换计算短时能量、过零率的代码

当使用MATLAB对.wav文件进行短时傅里叶变换计算短时能量和过零率时，可以按照以下步骤进行：

1. 读取.wav文件：

filename = 'your_file.wav';
[y, Fs] = audioread(filename);

其中，`y`是音频数据，`Fs`是采样率。

2. 分帧：

frameLength = 256; % 设置帧长
overlap = 0.5; % 设置重叠率
frameShift = frameLength * (1 - overlap); % 计算帧移
numFrames = floor((length(y) - frameLength) / frameShift) + 1; % 计算帧数

frames = zeros(frameLength, numFrames); % 存储分帧结果
for i = 1:numFrames
    index = (i-1) * frameShift + 1;
    frames(:, i) = y(index:index+frameLength-1);
end

3. 对每一帧进行短时傅里叶变换：

nfft = 512; % 设置FFT点数
spectrogram = zeros(nfft/2+1, numFrames); % 存储短时傅里叶变换结果

for i = 1:numFrames
    frame = frames(:, i);
    spectrum = abs(fft(frame, nfft));
    spectrogram(:, i) = spectrum(1:nfft/2+1);
end

4. 计算短时能量：

energy = sum(abs(frames).^2);

5. 计算过零率：

zeroCrossingRate = sum(abs(diff(sign(frames)))) / (2 * frameLength);

以上是MATLAB对.wav文件进行短时傅里叶变换计算短时能量和过零率的代码。