汉明窗频谱分割程序matlab

汉明窗是一种常用的窗函数，可以用于谱分析。以下是使用Matlab实现汉明窗频谱分割的示例程序：

% 读入音频文件
[y,fs]=audioread('test.wav');

% 设置分帧的参数
frameSize=256; % 帧长为256个采样点
frameShift=128; % 帧移为128个采样点

% 分帧
nFrames=floor((length(y)-frameSize)/frameShift)+1; % 计算帧数
frames=zeros(frameSize,nFrames); % 存储所有帧
for i=1:nFrames
    startIdx=(i-1)*frameShift+1; % 当前帧的起始位置
    frames(:,i)=y(startIdx:startIdx+frameSize-1).*hamming(frameSize); % 加上汉明窗
end

% 计算每一帧的能量谱
nFFT=512; % FFT长度为512个点
spectrums=zeros(nFFT/2+1,nFrames); % 存储所有帧的能量谱
for i=1:nFrames
    spectrum=abs(fft(frames(:,i),nFFT)).^2; % 以2为底的FFT
    spectrums(:,i)=spectrum(1:nFFT/2+1); % 取FFT结果的前一半
end

% 进行频谱分割
threshold=mean(spectrums(:)); % 计算能量谱的平均值
segments=zeros(nFrames,1); % 存储分割结果
for i=1:nFrames
    if mean(spectrums(:,i))>threshold % 如果当前帧的能量谱大于平均值
        segments(i)=1; % 则将该帧标记为语音帧
    end
end

% 在原始音频信号上画出分割结果
t=(0:length(y)-1)/fs;
subplot(2,1,1);
plot(t,y);
axis([0 max(t) -1 1]);
xlabel('Time (s)');
title('Original Signal');
subplot(2,1,2);
plot(t,segments);
axis([0 max(t) -0.5 1.5]);
xlabel('Time (s)');
title('Speech Segments');

在程序中，首先读入音频文件，并设置了分帧的参数。然后，使用`hamming`函数加上汉明窗对每一帧进行加窗，再计算每一帧的能量谱。接着，根据能量谱的平均值进行频谱分割，将能量谱高于平均值的帧标记为语音帧。最后，在原始音频信号上画出分割结果。